DeepAI AI Chat
Log In Sign Up

A fast Chebyshev method for the Bingham closure with application to active nematic suspensions

by   Scott Weady, et al.

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 2d-1 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.


page 11

page 16

page 19


Efficient and Accurate Adaptive Resolution for Weakly-Compressible SPH

In this paper we propose an accurate, and computationally efficient meth...

Learning invariance preserving moment closure model for Boltzmann-BGK equation

As one of the main governing equations in kinetic theory, the Boltzmann ...

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

We demonstrate the utility of deep learning for modeling the clustering ...

A neural network closure for the Euler-Poisson system based on kinetic simulations

This work deals with the modeling of plasmas, which are charged-particle...

Deep Neural Networks for Data-Based Turbulence Models

In this work, we present a novel data-based approach to turbulence model...