A fast Chebyshev method for the Bingham closure with application to active nematic suspensions
arXiv:2106.14817 · doi:10.1016/j.jcp.2021.110937
Abstract
Continuum kinetic theories provide an important tool for the analysis and simulation of particle suspensions. When those particles are anisotropic, the addition of a particle orientation vector to the kinetic description yields a dimensional theory which becomes intractable to simulate, especially in three dimensions or near states where the particles are highly aligned. Coarse-grained theories that track only moments of the particle distribution functions provide a more efficient simulation framework, but require closure assumptions. For the particular case where the particles are apolar, the Bingham closure has been found to agree well with the underlying kinetic theory; yet the closure is non-trivial to compute, requiring the solution of an often nearly-singular nonlinear equation at every spatial discretization point at every timestep. In this paper, we present a robust, accurate, and efficient numerical scheme for evaluating the Bingham closure, with a controllable error/efficiency tradeoff. To demonstrate the utility of the method, we carry out high-resolution simulations of a coarse-grained continuum model for a suspension of active particles in parameter regimes inaccessible to kinetic theories. Analysis of these simulations reveals that inaccurately computing the closure can act to effectively limit spatial resolution in the coarse-grained fields. Pushing these simulations to the high spatial resolutions enabled by our method reveals a coupling between vorticity and topological defects in the suspension director field, as well as signatures of energy transfer between scales in this active fluid model.
References in corpus (16)
- Defect annihilation and proliferation in active nematics
- Topological structure and dynamics of three dimensional active nematics
- Efficient numerical diagonalization of hermitian 3x3 matrices
- Structure formation in active networks
- Spontaneous Circulation of Confined Active Suspensions
- Multiscale polar theory of microtubule and motor-protein assemblies
- Vorticity, Defects and Correlations in Active Turbulence
- Universal scaling of active nematic turbulence
- Active turbulence in active nematics
- Dynamic self-assembly of microscale rotors and swimmers
- From microscopic theory to macroscopic theory: a systematic study on static modeling for liquid crystals
- Analytical structure, dynamics, and coarse-graining of a kinetic model of an active fluid
- Cascade or not cascade? Energy transfer and elastic effects in active nematics
- On the stupendous beauty of closure
- A Fast Algorithm for the Moments of Bingham Distribution
- Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere