• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Your access to this publication is provided through the subscription of Carl von Ossietzky Universitat. Library Welcome Message Help

Library Welcome Message Help

close

Attention, institutional online journal administrators: if the institution name printed in this screen message is incorrect, you can request a change by sending an email to librarynames@aip.org.

Please include your institutional account number, the account administrator's email address and the preferred online display name.

FULL-TEXT OPTIONS:

Phys. Fluids 24, 091902 (2012); http://dx.doi.org/10.1063/1.4754873 (12 pages)

Interactions between active particles and dynamical structures in chaotic flow

Nidhi Khurana and Nicholas T. Ouellette

Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, USA

View MapView Map

(Received 10 April 2012; accepted 11 September 2012; published online 26 September 2012)

Using a simple model, we study the transport dynamics of active, swimming particles advected in a two-dimensional chaotic flow field. We work with self-propelled, point-like particles that are either spherical or ellipsoidal. Swimming is modeled as a combination of a fixed intrinsic speed and stochastic terms in both the translational and rotational equations of motion. We show that the addition of motility to the particles causes them to feel the dynamical structure of the flow field in a different way from fluid particles, with macroscopic effects on the particle transport. At low swimming speeds, transport is suppressed due to trapping on transport barriers in the flow; we show that this effect is enhanced when stochastic terms are added to the swimming model or when the particles are elongated. At higher speeds, we find that elongated swimmers tend be attracted to the stable manifolds of hyperbolic fixed points, leading to increased transport relative to swimming spheres. Our results may have significant implications for models of real swimming organisms in finite-Reynolds-number flows.

© 2012 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. MODEL
    1. Flow field
    2. Swimmers
  3. SPHERICAL PARTICLES
    1. Deterministic swimming
    2. Translational stochasticity
    3. Rotational stochasticity
  4. ELLIPSOIDAL PARTICLES
  5. SUMMARY AND CONCLUSIONS

KEYWORDS and PACS

PACS

ARTICLE DATA

PUBLICATION DATA

ISSN

1070-6631 (print)  
1089-7666 (online)

    References

  1. M. R. Maxey and J. J. Riley, “Equation of motion for a small rigid sphere in a nonuniform flow,” Phys. Fluids 26, 883 (1983).
  2. A. Babiano, J. H. E. Cartwright, O. Piro, and A. Provenzale, “Dynamics of a small neutrally buoyant sphere in a fluid and targeting in Hamiltonian systems,” Phys. Rev. Lett. 84, 5764 (2000). [MEDLINE]
  3. G. B. Jeffery, “The motion of ellipsoidal particles immersed in a viscous fluid,” Proc. R. Soc. London, Ser. A 102, 161 (1922).
  4. S. Parsa, J. S. Guasto, M. Kishore, N. T. Ouellette, and J. P. Gollub, “Rotation and alignment of rods in two-dimensional chaotic flow,” Phys. Fluids 23, 043302 (2011)PHFLE6000023000004043302000001.
  5. Z. Toroczkai and T. Tél, “Introduction: Active chaotic flow,” Chaos 12, 372 (2002)CHAOEH000012000002000372000001. [ISI] [MEDLINE]
  6. J. R. Howse, R. A. L. Jones, A. J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, “Self-motile colloidal particles: From directed propulsion to random walk,” Phys. Rev. Lett. 99, 048102 (2007). [ISI] [MEDLINE]
  7. R. Dreyfus, J. Baudry, M. L. Roper, M. Fermingier, H. A. Stone, and J. Bibette, “Microscopic artificial swimmers,” Nature (London) 437, 862 (2005).
  8. T. J. Pedley and J. O. Kessler, “Hydrodynamic phenomena in suspensions of swimming microorganisms,” Annu. Rev. Fluid Mech. 24, 313 (1992). [ISI]
  9. M. J. Lighthill, “On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers,” Commun. Pure Appl. Math. 5, 109 (1952).
  10. E. M. Purcell, “Life at low Reynolds number,” Am. J. Phys. 45, 3 (1977).
  11. E. Lauga and T. R. Powers, “The hydrodynamics of swimming microorganisms,” Rep. Prog. Phys. 72, 096601 (2009).
  12. X.-L. Wu and A. Libchaber, “Particle diffusion in a quasi-two-dimensional bacterial bath,” Phys. Rev. Lett. 84, 3017 (2000). [MEDLINE]
  13. P. T. Underhill, J. P. Hernandez-Ortiz, and M. D. Graham, “Diffusion and spatial correlations in suspensions of swimming particles,” Phys. Rev. Lett. 100, 248101 (2008). [MEDLINE]
  14. D. Saintillan and M. Shelley, “Instabilities, pattern formation, and mixing in active suspensions,” Phys. Fluids 20, 123304 (2008)PHFLE6000020000012123304000001.
  15. J. E. Colgate and K. M. Lynch, “Mechanics and control of swimming: a review,” IEEE J. Ocean. Eng. 29, 660 (2004).
  16. J. O. Kessler, “Hydrodynamic focusing of motile algal cells,” Nature (London) 313, 218 (1985). [Inspec] [ISI]
  17. W. M. Durham, J. O. Kessler, and R. Stocker, “Disruption of vertical motility by shear triggers formation of thin phytoplankton layers,” Science 323, 1067 (2009). [MEDLINE]
  18. W. M. Durham, E. Climent, and R. Stocker, “Gyrotaxis in a steady vortical flow,” Phys. Rev. Lett. 106, 238102 (2011).
  19. C. Torney and Z. Neufeld, “Transport and aggregation of self-propelled particles in fluid flows,” Phys. Rev. Lett. 99, 078101 (2007). [ISI]
  20. N. Khurana, J. Blawzdziewicz, and N. T. Ouellette, “Reduced transport of swimming particles in chaotic flow due to hydrodynamic trapping,” Phys. Rev. Lett. 106, 198104 (2011).
  21. Marcos, J. R. Seymour, M. Luhar, W. M. Durham, J. G. Mitchell, A. Macke, and R. Stocker, “Microbial alignment in flow changes ocean light climate,” Proc. Natl. Acad. Sci. U.S.A. 108, 3860 (2011).
  22. F. G. Schmitt and L. Seuront, “Intermittent turbulence and copepod dynamics: Increase in encounter rates through preferential concentration,” J. Mar. Syst. 70, 263 (2008).
  23. T. H. Solomon and J. P. Gollub, “Chaotic particle transport in time-dependent Rayleigh-Bénard convection,” Phys. Rev. A 38, 6280 (1988). [MEDLINE]
  24. T. Nishikawa, Z. Toroczkai, and C. Grebogi, “Advective coalescence in chaotic flows,” Phys. Rev. Lett. 87, 038301 (2001). [ISI] [MEDLINE]
  25. R. Festa, A. Mazzino, and M. Todini, “Dynamics of light particles in oscillating cellular flows,” Phys. Rev. E 80, 035301 (2009).
  26. H. Aref, “Stirring by chaotic advection,” J. Fluid Mech. 143, 1 (1984).
  27. J. M. Ottino, The Kinematics of Mixing: Stretching, Chaos, and Transport (Cambridge University Press, Cambridge, 1989).
  28. G. A. Voth, G. Haller, and J. P. Gollub, “Experimental measurements of stretching fields in fluid mixing,” Phys. Rev. Lett. 88, 254501 (2002).
  29. H. C. Berg and D. A. Brown, “Chemotaxis in Escherichia coli analysed by three-dimensional tracking,” Nature (London) 239, 500 (1972). [MEDLINE]
  30. M. Polin, I. Tuval, K. Drescher, J. P. Gollub, and R. E. Goldstein, “Chlamydomonas swims with two “gears” in a Eukaryotic version of run-and-tumble locomotion,” Science 325, 487 (2009). [MEDLINE]
  31. E. G. Altmann and A. Endler, “Noise-enhanced trapping in chaotic scattering,” Phys. Rev. Lett. 105, 244102 (2010).
  32. G. Haller, “Distinguished material surfaces and coherent structures in three-dimensional fluid flows,” Physica D 149, 248 (2001).

Figures (click on thumbnails to view enlargements)

FIG.1
(a) Snapshot of one unit cell of the flow field at zero phase (that is, for t such that sin Ωt = 0). Arrows show the local fluid velocity, and the shading shows the vorticity. As time progresses, the flow field oscillates sinusoidally in the horizontal direction. (b) Poincare section (at zero phase) for a fluid element initially in the chaotic sea, for B = 0.12 and Ω = 6.28. Only one quarter of the unit cell is shown (corresponding to the lower-left vortex in (a)); the rest of the unit cell is related by symmetry. The central empty region is a period-1 elliptic island; the surrounding empty regions are a period-3 island chain. (c) Finite-time Lyapunov exponent (FTLE) field at zero phase. Again, only one quarter of the unit cell is shown.

FIG.1 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.2
Chaotic diffusion coefficient D normalized by D0, the diffusion coefficient for fluid elements, as a function of the swimming speed vs, for spherical swimmers with σs = σr = 0. Error bars are computed from the statistical fluctuations between many sets of simulations. Two distinct regions of suppressed long-time transport are seen, corresponding to trapping by the period-3 islands and by the period-1 islands.20

FIG.2 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.3
Chaotic diffusion coefficient D normalized by D0 as a function of σs, for vs = 0, σr = 0, and α = 0. When compared with Fig. 2, transport is more strongly suppressed, but the distinct signatures of each type of elliptic island are no longer present.

FIG.3 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.4
Average time for a swimmer to cross a cell boundary as a function of its spatial location for (a) vs = 0.004 and σs = σr = 0 and (b) σs = 0.15 and vs = σr = 0. Only one quarter of the flow domain is shown. The shade/color bar gives the cell-crossing time in flow cycles. For the deterministic swimmer in (a), trapping is strong in the period-3 islands and on a small ring just inside the period-1 island. The purely stochastic swimmer can wander into the core of the period-1 island, where trapping is strongest.

FIG.4 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.5
Chaotic diffusion coefficients D in the two-dimensional parameter space spanned by vs and σs for σr = 0. The shade/color bar shows D relative to D0. The black line separates the regions of suppressed transport (D/D0 < 1) from those of enhanced transport (D/D0 > 1).

FIG.5 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.6
Chaotic diffusion coefficients D in the two-dimensional parameter space spanned by vs and σr for σs = 0. The shade/color bar shows D relative to D0. The black line separates the regions of suppressed transport (D/D0 < 1) from those of enhanced transport (D/D0 > 1). Transport is suppressed in larger region of parameter space than it was for purely translational stochasticity case in Fig. 5.

FIG.6 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.7
Spatially resolved maps of the average time to cross a cell boundary for vs = 0.004, and σr = (a) 0.1, (b) 1.0, (c) 3.0, and (d) 6.0. σs = 0 for all panels. The shade/color bar gives the cell-crossing times in flow cycles. The times are much longer than the comparable vs = 0.004, σr = 0 case shown in Fig. 4a.

FIG.7 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.8
(a) Chaotic diffusion coefficient D normalized by D0 as a function of eccentricity α for deterministic swimmers with vs = 0.002. Transport is much more strongly suppressed for ellipsoids of intermediate eccentricities than it is for spheres. (b)–(d) Probability density functions of swimmer position for (b) α = 0, (c) α = 0.5, and (d) α = 1. The strong suppression of transport for ellipsoidal particles is due to the formation of attractors.

FIG.8 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.9
(a) Chaotic diffusion coefficient D normalized by D0 as a function of eccentricity α for deterministic swimmers with vs = 0.08. For this speed, transport of ellipsoidal particles is strongly enhanced relative to spheres. (b)–(d) Probability density functions of swimmer position for (b) α = 0, (c) α = 0.5, and (d) α = 1. The strong enhancement of transport is due to the clustering of ellipsoids on the stable manifolds of the hyperbolic fixed points.

FIG.9 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.10
Mean nearest neighbor distance δNN scaled by the value for passive particles δ0 as a function of eccentricity for (a) vs = 0.002 and (b) vs = 0.08. High aspect ratio swimmers tend to be closer to each than spherical swimmers are, leading to enhanced encounter rates.

FIG.10 Download High Resolution Image (.zip file) | Export Figure to PowerPoint



Close
Google Calendar
ADVERTISEMENT

close