Instabilities and nonlinear dynamics of concentrated active suspensions

Barath Ezhilan, Michael Shelley, David Saintillan

Research output: Contribution to journalArticle

Abstract

Suspensions of active particles, such as motile microorganisms and artificial microswimmers, are known to undergo a transition to complex large-scale dynamics at high enough concentrations. While a number of models have demonstrated that hydrodynamic interactions can in some cases explain these dynamics, collective motion in experiments is typically observed at such high volume fractions that steric interactions between nearby swimmers are significant and cannot be neglected. This raises the question of the respective roles of steric vs hydrodynamic interactions in these dense systems, which we address in this paper using a continuum theory and numerical simulations. The model we propose is based on our previous kinetic theory for dilute suspensions, in which a conservation equation for the distribution function of particle configurations is coupled to the Stokes equations for the fluid motion [D. Saintillan and M. J. Shelley, "Instabilities, pattern formation, and mixing in active suspensions," Phys. Fluids20, 123304 (2008)]10.1063/1.3041776. At high volume fractions, steric interactions are captured by extending classic models for concentrated suspensions of rodlike polymers, in which contacts between nearby particles cause them to align locally. In the absence of hydrodynamic interactions, this local alignment results in a transition from an isotropic base state to a nematic base state when volume fraction is increased. Using a linear stability analysis, we first investigate the hydrodynamic stability of both states. Our analysis shows that suspensions of pushers, or rear-actuated swimmers, typically become unstable in the isotropic state before the transition occurs; suspensions of pullers, or head-actuated swimmers, can also become unstable, though the emergence of unsteady flows in this case occurs at a higher concentration, above the nematic transition. These results are also confirmed using fully nonlinear numerical simulations in a periodic cubic domain, where pusher and puller suspensions are indeed both found to exhibit instabilities at sufficiently high volume fractions; these instabilities lead to unsteady chaotic states characterized by large-scale correlated motions and strong density fluctuations. While the dynamics in suspensions of pushers are similar to those previously reported in the dilute regime, the instability of pullers is novel and typically characterized by slower dynamics and weaker hydrodynamic velocities and active input power than in pusher suspensions at the same volume fraction.

Original languageEnglish (US)
Article number070607
JournalPhysics of Fluids
Volume25
Issue number7
DOIs
StatePublished - Jul 18 2013

Fingerprint

Volume fraction
Hydrodynamics
hydrodynamics
Linear stability analysis
Kinetic theory
Computer simulation
Unsteady flow
interactions
Microorganisms
Distribution functions
Conservation
flow stability
conservation equations
unsteady flow
Fluids
microorganisms
kinetic theory
Polymers
simulation
distribution functions

ASJC Scopus subject areas

  • Condensed Matter Physics

Cite this

Instabilities and nonlinear dynamics of concentrated active suspensions. / Ezhilan, Barath; Shelley, Michael; Saintillan, David.

In: Physics of Fluids, Vol. 25, No. 7, 070607, 18.07.2013.

Research output: Contribution to journalArticle

Ezhilan, Barath ; Shelley, Michael ; Saintillan, David. / Instabilities and nonlinear dynamics of concentrated active suspensions. In: Physics of Fluids. 2013 ; Vol. 25, No. 7.
@article{62fcafd3f8f74182a1ef4ac838e11d22,
title = "Instabilities and nonlinear dynamics of concentrated active suspensions",
abstract = "Suspensions of active particles, such as motile microorganisms and artificial microswimmers, are known to undergo a transition to complex large-scale dynamics at high enough concentrations. While a number of models have demonstrated that hydrodynamic interactions can in some cases explain these dynamics, collective motion in experiments is typically observed at such high volume fractions that steric interactions between nearby swimmers are significant and cannot be neglected. This raises the question of the respective roles of steric vs hydrodynamic interactions in these dense systems, which we address in this paper using a continuum theory and numerical simulations. The model we propose is based on our previous kinetic theory for dilute suspensions, in which a conservation equation for the distribution function of particle configurations is coupled to the Stokes equations for the fluid motion [D. Saintillan and M. J. Shelley, {"}Instabilities, pattern formation, and mixing in active suspensions,{"} Phys. Fluids20, 123304 (2008)]10.1063/1.3041776. At high volume fractions, steric interactions are captured by extending classic models for concentrated suspensions of rodlike polymers, in which contacts between nearby particles cause them to align locally. In the absence of hydrodynamic interactions, this local alignment results in a transition from an isotropic base state to a nematic base state when volume fraction is increased. Using a linear stability analysis, we first investigate the hydrodynamic stability of both states. Our analysis shows that suspensions of pushers, or rear-actuated swimmers, typically become unstable in the isotropic state before the transition occurs; suspensions of pullers, or head-actuated swimmers, can also become unstable, though the emergence of unsteady flows in this case occurs at a higher concentration, above the nematic transition. These results are also confirmed using fully nonlinear numerical simulations in a periodic cubic domain, where pusher and puller suspensions are indeed both found to exhibit instabilities at sufficiently high volume fractions; these instabilities lead to unsteady chaotic states characterized by large-scale correlated motions and strong density fluctuations. While the dynamics in suspensions of pushers are similar to those previously reported in the dilute regime, the instability of pullers is novel and typically characterized by slower dynamics and weaker hydrodynamic velocities and active input power than in pusher suspensions at the same volume fraction.",
author = "Barath Ezhilan and Michael Shelley and David Saintillan",
year = "2013",
month = "7",
day = "18",
doi = "10.1063/1.4812822",
language = "English (US)",
volume = "25",
journal = "Physics of Fluids",
issn = "1070-6631",
publisher = "American Institute of Physics Publising LLC",
number = "7",

}

TY - JOUR

T1 - Instabilities and nonlinear dynamics of concentrated active suspensions

AU - Ezhilan, Barath

AU - Shelley, Michael

AU - Saintillan, David

PY - 2013/7/18

Y1 - 2013/7/18

N2 - Suspensions of active particles, such as motile microorganisms and artificial microswimmers, are known to undergo a transition to complex large-scale dynamics at high enough concentrations. While a number of models have demonstrated that hydrodynamic interactions can in some cases explain these dynamics, collective motion in experiments is typically observed at such high volume fractions that steric interactions between nearby swimmers are significant and cannot be neglected. This raises the question of the respective roles of steric vs hydrodynamic interactions in these dense systems, which we address in this paper using a continuum theory and numerical simulations. The model we propose is based on our previous kinetic theory for dilute suspensions, in which a conservation equation for the distribution function of particle configurations is coupled to the Stokes equations for the fluid motion [D. Saintillan and M. J. Shelley, "Instabilities, pattern formation, and mixing in active suspensions," Phys. Fluids20, 123304 (2008)]10.1063/1.3041776. At high volume fractions, steric interactions are captured by extending classic models for concentrated suspensions of rodlike polymers, in which contacts between nearby particles cause them to align locally. In the absence of hydrodynamic interactions, this local alignment results in a transition from an isotropic base state to a nematic base state when volume fraction is increased. Using a linear stability analysis, we first investigate the hydrodynamic stability of both states. Our analysis shows that suspensions of pushers, or rear-actuated swimmers, typically become unstable in the isotropic state before the transition occurs; suspensions of pullers, or head-actuated swimmers, can also become unstable, though the emergence of unsteady flows in this case occurs at a higher concentration, above the nematic transition. These results are also confirmed using fully nonlinear numerical simulations in a periodic cubic domain, where pusher and puller suspensions are indeed both found to exhibit instabilities at sufficiently high volume fractions; these instabilities lead to unsteady chaotic states characterized by large-scale correlated motions and strong density fluctuations. While the dynamics in suspensions of pushers are similar to those previously reported in the dilute regime, the instability of pullers is novel and typically characterized by slower dynamics and weaker hydrodynamic velocities and active input power than in pusher suspensions at the same volume fraction.

AB - Suspensions of active particles, such as motile microorganisms and artificial microswimmers, are known to undergo a transition to complex large-scale dynamics at high enough concentrations. While a number of models have demonstrated that hydrodynamic interactions can in some cases explain these dynamics, collective motion in experiments is typically observed at such high volume fractions that steric interactions between nearby swimmers are significant and cannot be neglected. This raises the question of the respective roles of steric vs hydrodynamic interactions in these dense systems, which we address in this paper using a continuum theory and numerical simulations. The model we propose is based on our previous kinetic theory for dilute suspensions, in which a conservation equation for the distribution function of particle configurations is coupled to the Stokes equations for the fluid motion [D. Saintillan and M. J. Shelley, "Instabilities, pattern formation, and mixing in active suspensions," Phys. Fluids20, 123304 (2008)]10.1063/1.3041776. At high volume fractions, steric interactions are captured by extending classic models for concentrated suspensions of rodlike polymers, in which contacts between nearby particles cause them to align locally. In the absence of hydrodynamic interactions, this local alignment results in a transition from an isotropic base state to a nematic base state when volume fraction is increased. Using a linear stability analysis, we first investigate the hydrodynamic stability of both states. Our analysis shows that suspensions of pushers, or rear-actuated swimmers, typically become unstable in the isotropic state before the transition occurs; suspensions of pullers, or head-actuated swimmers, can also become unstable, though the emergence of unsteady flows in this case occurs at a higher concentration, above the nematic transition. These results are also confirmed using fully nonlinear numerical simulations in a periodic cubic domain, where pusher and puller suspensions are indeed both found to exhibit instabilities at sufficiently high volume fractions; these instabilities lead to unsteady chaotic states characterized by large-scale correlated motions and strong density fluctuations. While the dynamics in suspensions of pushers are similar to those previously reported in the dilute regime, the instability of pullers is novel and typically characterized by slower dynamics and weaker hydrodynamic velocities and active input power than in pusher suspensions at the same volume fraction.

UR - http://www.scopus.com/inward/record.url?scp=84881523458&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84881523458&partnerID=8YFLogxK

U2 - 10.1063/1.4812822

DO - 10.1063/1.4812822

M3 - Article

VL - 25

JO - Physics of Fluids

JF - Physics of Fluids

SN - 1070-6631

IS - 7

M1 - 070607

ER -