Computational mean-field modeling of confined active fluids

Abstract We present a new framework for the efficient simulation of the dynamics of active fluids in complex two- and three-dimensional microfluidic geometries. Focusing on the case of a suspension of microswimmers such as motile bacteria, we adopt a continuum mean-field model based on partial differential equations for the evolution of the concentration, polarization and nematic tensor fields, which are nonlinearly coupled to the Navier-Stokes equations for the fluid flow driven by internal active stresses. A level set method combined with an adaptive mesh refinement scheme on Quad-/Octree grids is used to capture complex domain shapes while refining the solution near boundaries or in the neighborhood of sharp gradients. A hybrid finite volumes/finite differences method is implemented in which the concentration field is treated using finite volumes to ensure mass conservation, while the polarization and nematic alignment fields are treated using a combination of finite differences and finite volumes for enhanced accuracy. The governing equations for these fields are solved along with the Navier-Stokes equations, which are evolved using an unconditionally stable projection solver. We illustrate the versatility and robustness of our method by analyzing spontaneous active flows in various two- and three-dimensional systems. Our results show excellent agreement with previous models and experiments and pave the way for further developments in active microfluidics.

[1]  Frédéric Gibou,et al.  A Multigrid Method on Non-Graded Adaptive Octree and Quadtree Cartesian Grids , 2013, J. Sci. Comput..

[2]  F. Bretherton.  The motion of rigid particles in a shear flow at low Reynolds number , 1962, Journal of Fluid Mechanics.

[3]  D. Saintillan.  Rheology of Active Fluids , 2018 .

[4]  T. Lubensky,et al.  Statistical mechanics and hydrodynamics of bacterial suspensions , 2009, Proceedings of the National Academy of Sciences.

[5]  A. Maude,et al.  Non-random Distribution of Bull Spermatozoa in a Drop of Sperm Suspension , 1963, Nature.

[6]  Anke Lindner,et al.  Living on the edge: transfer and traffic of E. coli in a confined flow. , 2015, Soft matter.

[7]  Michael J Shelley,et al.  Stability of active suspensions. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  ANALYTICAL STRUCTURE , 2020, Research Handbook on Unjust Enrichment and Restitution.

[9]  G. B. Jeffery.  The motion of ellipsoidal particles immersed in a viscous fluid , 1922 .

[10]  T. Aslam.  A partial differential equation approach to multidimensional extrapolation , 2004 .

[11]  Frédéric Gibou,et al.  A Sharp Computational Method for the Simulation of the Solidification of Binary Alloys , 2015, J. Sci. Comput..

[12]  Sriram Ramaswamy,et al.  Rheology of active-particle suspensions. , 2003, Physical review letters.

[13]  David Saintillan,et al.  Microfluidic flow actuation using magnetoactive suspensions , 2018 .

[14]  Frédéric Gibou,et al.  A stable projection method for the incompressible Navier-Stokes equations on arbitrary geometries and adaptive Quad/Octrees , 2015, J. Comput. Phys..

[15]  Donald L. Koch,et al.  Critical bacterial concentration for the onset of collective swimming , 2009, Journal of Fluid Mechanics.

[16]  E. Virga.  Instabilities , 2018, Variational Theories for Liquid Crystals.

[17]  K. Drescher,et al.  Direct measurement of the flow field around swimming microorganisms. , 2010, Physical review letters.

[18]  Michael Shelley,et al.  Instabilities and nonlinear dynamics of concentrated active suspensions , 2013 .

[19]  Chiu Fan Lee.  Active particles under confinement: aggregation at the wall and gradient formation inside a channel , 2013, 1304.5556.

[20]  Mohammad Mirzadeh,et al.  A second-order discretization of the nonlinear Poisson-Boltzmann equation over irregular geometries using non-graded adaptive Cartesian grids , 2011, J. Comput. Phys..

[21]  Donald L. Koch,et al.  Collective Hydrodynamics of Swimming Microorganisms: Living Fluids , 2011 .

[22]  Frank Losasso,et al.  Simulating water and smoke with an octree data structure , 2004, SIGGRAPH 2004.

[23]  Jörn Dunkel,et al.  Confinement stabilizes a bacterial suspension into a spiral vortex. , 2013, Physical review letters.

[24]  Jay X. Tang,et al.  Accumulation of microswimmers near a surface mediated by collision and rotational Brownian motion. , 2009, Physical review letters.

[25]  Raymond E. Goldstein,et al.  Directed collective motion of bacteria under channel confinement , 2016, 1603.01143.

[26]  Frédéric Gibou,et al.  Geometric integration over irregular domains with application to level-set methods , 2007, J. Comput. Phys..

[27]  P. Gao,et al.  Dynamics and structure of an apolar active suspension in an annulus , 2017, Journal of Fluid Mechanics.

[28]  M. Shelley,et al.  Instabilities and pattern formation in active particle suspensions: kinetic theory and continuum simulations. , 2008, Physical review letters.

[29]  Takuji Ishikawa,et al.  Orientational order in concentrated suspensions of spherical microswimmers , 2011, 1302.2288.

[30]  David Saintillan,et al.  Transport of a dilute active suspension in pressure-driven channel flow , 2015, Journal of Fluid Mechanics.

[31]  Rothschild,et al.  Non-random Distribution of Bull Spermatozoa in a Drop of Sperm Suspension , 1963, Nature.

[32]  Frédéric Gibou,et al.  A second order accurate projection method for the incompressible Navier-Stokes equations on non-graded adaptive grids , 2006, J. Comput. Phys..

[33]  Gerhard Gompper,et al.  Wall accumulation of self-propelled spheres , 2013 .

[34]  Frédéric Gibou,et al.  A second-order sharp numerical method for solving the linear elasticity equations on irregular domains and adaptive grids - Application to shape optimization , 2013, J. Comput. Phys..

[35]  M. Shelley.  The Dynamics of Microtubule/Motor-Protein Assemblies in Biology and Physics , 2016 .

[36]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[37]  Andrey Sokolov,et al.  Engineering bacterial vortex lattice via direct laser lithography , 2018, Nature Communications.

[38]  David Saintillan,et al.  Instability regimes in flowing suspensions of swimming micro-organisms , 2011 .

[39]  David Saintillan,et al.  The Dilute Rheology of Swimming Suspensions: A Simple Kinetic Model , 2010 .

[40]  Raymond E. Goldstein,et al.  Ferromagnetic and antiferromagnetic order in bacterial vortex lattices , 2013, Nature Physics.

[41]  J. Posner,et al.  Phoretic Self-Propulsion , 2017 .

[42]  Jay X. Tang,et al.  Accumulation of swimming bacteria near a solid surface. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  J. Joanny,et al.  Spontaneous flow transition in active polar gels , 2005, q-bio/0503022.

[44]  J. Dunkel,et al.  Fluid dynamics and noise in bacterial cell–cell and cell–surface scattering , 2011, Proceedings of the National Academy of Sciences.

[45]  M. Shelley,et al.  Emergence of coherent structures and large-scale flows in motile suspensions , 2012, Journal of The Royal Society Interface.

[46]  M. Graham,et al.  Transport and collective dynamics in suspensions of confined swimming particles. , 2005, Physical review letters.

[47]  Frédéric Gibou,et al.  Level-set simulations of soluble surfactant driven flows , 2017, J. Comput. Phys..

[48]  Maxime Theillard,et al.  Geometric control of active collective motion. , 2016, Soft matter.

[49]  M. Shelley,et al.  Instabilities, pattern formation and mixing in active suspensions , 2008 .

[50]  Seth Fraden,et al.  Transition from turbulent to coherent flows in confined three-dimensional active fluids , 2017, Science.

[51]  Raymond E. Goldstein,et al.  Nonlinear concentration patterns and bands in autochemotactic suspensions , 2018, Physical Review E.

[52]  Hanan Samet,et al.  Applications of spatial data structures - computer graphics, image processing, and GIS , 1990 .

[53]  H. H. Wensink,et al.  Meso-scale turbulence in living fluids , 2012, Proceedings of the National Academy of Sciences.

[54]  Frédéric Gibou,et al.  A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids , 2006, J. Comput. Phys..

[55]  S. Ramaswamy,et al.  Hydrodynamics of soft active matter , 2013 .

[56]  T. Powers,et al.  The hydrodynamics of swimming microorganisms , 2008, 0812.2887.

[57]  David A Boy,et al.  An Adaptive, Finite Difference Solver for the Nonlinear Poisson-Boltzmann Equation with Applications to Biomolecular Computations , 2013 .

[58]  Francis G. Woodhouse,et al.  Spontaneous circulation of confined active suspensions. , 2012, Physical review letters.

[59]  David Saintillan,et al.  Chaotic dynamics and oxygen transport in thin films of aerotactic bacteria , 2012 .

[60]  Daniel T. N. Chen,et al.  Spontaneous motion in hierarchically assembled active matter , 2012, Nature.

[61]  Michael J Shelley,et al.  Orientational order and instabilities in suspensions of self-locomoting rods. , 2007, Physical review letters.

[62]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[63]  Michael Shelley,et al.  Active suspensions and their nonlinear models , 2013 .

[64]  Gerhard Gompper,et al.  Run-and-tumble dynamics of self-propelled particles in confinement , 2015, 1503.06454.

[65]  M. Shelley,et al.  Analytical structure, dynamics, and coarse-graining of a kinetic model of an active fluid , 2017, 1703.00969.

[66]  Michael J Shelley,et al.  Multiscale modeling and simulation of microtubule-motor-protein assemblies. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[67]  Enkeleida Lushi,et al.  Fluid flows created by swimming bacteria drive self-organization in confined suspensions , 2014, Proceedings of the National Academy of Sciences.

[68]  Harold Auradou,et al.  Turning Bacteria Suspensions into Superfluids. , 2015, Physical review letters.