Microscopic Modelling of Active Bacterial Suspensions

We present two-dimensional simulations of chemotactic self-propelled bacteria swimming in a viscous fluid. Self-propulsion is modelled by a couple of forces of same intensity and opposite direction applied on the rigid bacterial body and on an associated region in the fluid representing the flagellar bundle. The method for solving the fluid flow and the motion of the bacteria is based on a variational formulation written on the whole domain, strongly coupling the fluid and the rigid particle problems: rigid motion is enforced by penalizing the strain rate tensor on the rigid domain, while incompressibility is treated by duality. This model allows to achieve an accurate description of fluid motion and hydrodynamic interactions in moderate to concentrated active suspensions. A mesoscopic model is also used, in which the size of the bacteria is supposed to be much smaller than the elements of fluid: the perturbation of the fluid due to propulsion and motion of the swimmers is neglected, and the fluid is only subjected to the buoyant forcing induced by the presence of the bacteria, which are denser than the fluid. Although this model does not accurately take into account hydrodynamic interactions, it is able to reproduce complex collective dynamics observed in concentrated bacterial suspensions, such as bioconvection. From a mathematical point of view, both models lead to a minimization problem which is solved with a standard Finite Element Method. In order to ensure robustness, a projection algorithm is used to deal with contacts between particles. We also reproduce chemotactic behaviour driven by oxygen: an advection-diffusion equation on the oxygen concentration is solved in the fluid domain, with a source term accounting for oxygen consumption by the bacteria. The orientations of the individual bacteria are subjected to random changes, with a frequency that depends on the surrounding oxygen concentration, in order to favor the direction of the concentration gradient.

[1]  H. Berg,et al.  Temporal comparisons in bacterial chemotaxis. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[2]  R. Glowinski.  Finite element methods for incompressible viscous flow , 2003 .

[3]  William S. Ryu,et al.  Real-Time Imaging of Fluorescent Flagellar Filaments , 2000, Journal of bacteriology.

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

[5]  T. Pedley,et al.  Hydrodynamic Phenomena in Suspensions of Swimming Microorganisms , 1992 .

[6]  Bertrand Maury,et al.  Apparent viscosity of a mixture of a Newtonian fluid and interacting particles , 2005 .

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

[8]  P. G. Ciarlet,et al.  Introduction a l'analyse numerique matricielle et a l'optimisation , 1984 .

[9]  Pierre Lubin,et al.  An adaptative augmented Lagrangian method for three-dimensional multimaterial flows , 2004 .

[10]  R. Temam,et al.  Mathematical Modeling in Continuum Mechanics: Index , 2000 .

[11]  Leonid Berlyand,et al.  Three-dimensional model for the effective viscosity of bacterial suspensions. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  R. Goldstein,et al.  Self-concentration and large-scale coherence in bacterial dynamics. , 2004, Physical review letters.

[13]  Bertrand Maury.  Numerical Analysis of a Finite Element/Volume Penalty Method , 2009, SIAM J. Numer. Anal..

[14]  H. Berg,et al.  Dynamics of bacterial swarming. , 2010, Biophysical journal.

[15]  Andrey Sokolov,et al.  Reduction of viscosity in suspension of swimming bacteria. , 2009, Physical review letters.

[16]  I. Tuval,et al.  Bacterial swimming and oxygen transport near contact lines. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[17]  Bertrand Maury,et al.  A time-stepping scheme for inelastic collisions , 2006, Numerische Mathematik.

[18]  Andrey Sokolov,et al.  Enhanced mixing and spatial instability in concentrated bacterial suspensions. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  A. Lefebvre.  Fluid-Particle simulations with FreeFem++ , 2007 .

[20]  A. Libchaber,et al.  Particle diffusion in a quasi-two-dimensional bacterial bath. , 2000, Physical review letters.

[21]  H. Berg.  Random Walks in Biology , 2018 .

[22]  J. Happel,et al.  Low Reynolds number hydrodynamics , 1965 .

[23]  Dale Kaiser,et al.  Bacterial Swarming: A Re-examination of Cell-Movement Patterns , 2007, Current Biology.

[24]  Howard C. Berg,et al.  E. coli in Motion , 2003 .

[25]  Sangtae Kim,et al.  Microhydrodynamics: Principles and Selected Applications , 1991 .

[26]  T. Pedley,et al.  Run and Tumble Chemotaxis in a Shear Flow: The Effect of Temporal Comparisons, Persistence, Rotational Diffusion, and Cell Shape , 2009, Bulletin of mathematical biology.

[27]  M. Hillairet.  Lack of Collision Between Solid Bodies in a 2D Incompressible Viscous Flow , 2007 .

[28]  R. Cortez,et al.  Fluid dynamics of self-propelled microorganisms, from individuals to concentrated populations , 2007 .

[29]  B. Maury,et al.  SIMULATION OF SELF-PROPELLED CHEMOTACTIC BACTERIA IN A STOKES FLOW , 2010 .

[30]  Fernando Peruani,et al.  Self-propelled particles with fluctuating speed and direction of motion in two dimensions. , 2007, Physical review letters.

[31]  Salima Rafaï,et al.  Effective viscosity of microswimmer suspensions. , 2009, Physical review letters.

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

[33]  B. Maury,et al.  A PENALTY METHOD FOR THE SIMULATION OF FLUID - RIGID BODY INTERACTION , 2005 .

[34]  M. Alber,et al.  Periodic reversal of direction allows Myxobacteria to swarm , 2009, Proceedings of the National Academy of Sciences.

[35]  D. Gérard-Varet,et al.  Regularity Issues in the Problem of Fluid Structure Interaction , 2008, 0805.2654.

[36]  I. Aranson,et al.  Effective shear viscosity and dynamics of suspensions of micro-swimmers from small to moderate concentrations , 2011, Journal of mathematical biology.

[37]  Leonid Berlyand,et al.  Effective viscosity of dilute bacterial suspensions: a two-dimensional model , 2008, Physical biology.

[38]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .