skip to main content


Title: Three-dimensional simulations of undulatory and amoeboid swimmers in viscoelastic fluids
Microorganisms often move through viscoelastic environments, as biological fluids frequently have a rich microstructure owing to the presence of large polymeric molecules. Research on the effect of fluid elasticity on the swimming kinematics of these organisms has usually been focused on those that move via cilia or flagellum. Experimentally, Shen (X. N. Shen et al. , Phys. Rev. Lett. , 2011, 106 , 208101) reported that the nematode C. elegans , a model organism used to study undulatory motion, swims more slowly as the Deborah number describing the fluid's elasticity is increased. This phenomenon has not been thoroughly studied via a fully resolved three-dimensional simulation; moreover, the effect of fluid elasticity on the swimming speed of organisms moving via euglenoid movement, such as E. gracilis , is completely unknown. In this study, we discuss the simulation of the arbitrary motion of an undulating or pulsating swimmer that occupies finite volume in three dimensions, with the ability to specify any differential viscoelastic rheological model for the surrounding fluid. To accomplish this task, we use a modified version of the Immersed Finite Element Method presented in a previous paper by Guido and Saadat in 2018 (A. Saadat et al. , Phys. Rev. E , 2018, 98 , 063316). In particular, this version allows for the simulation of deformable swimmers such that they evolve through an arbitrary set of specified shapes via a conformation-driven force. From our analysis, we observe several key trends not found in previous two-dimensional simulations or theoretical analyses for C. elegans , as well as novel results for the amoeboid motion. In particular, we find that regions of high polymer stress concentrated at the head and tail of the swimming C. elegans are created by strong extensional flow fields and are associated with a decrease in swimming speed for a given swimming stroke. In contrast, in two dimensions these regions of stress are commonly found distributed along the entire body, likely owing to the lack of a third dimension for polymer relaxation. A comparison of swim speeds shows that the calculations in two-dimensional simulations result in an over-prediction of the speed reduction. We believe that our simulation tool accurately captures the swimming motion of the two aforementioned model swimmers and furthermore, allows for the simulation of multiple deformable swimmers, as well as more complex swimming geometries. This methodology opens many new possibilities for future studies of swimmers in viscoelastic fluids.  more » « less
Award ID(s):
1803765
NSF-PAR ID:
10308900
Author(s) / Creator(s):
 ;  ;
Date Published:
Journal Name:
Soft Matter
Volume:
15
Issue:
24
ISSN:
1744-683X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Here, we use magnetically driven self-assembled achiral swimmers made of two to four superparamagnetic micro-particles to provide insight into how swimming kinematics develop in complex, shear-thinning fluids. Two model shear-thinning polymer fluids are explored, where measurements of swimming dynamics reveal contrasting propulsion kinematics in shear-thinning fluids vs a Newtonian fluid. When comparing the velocity of achiral swimmers in polymer fluids to their dynamics in water, we observe kinematics dependent on (1) no shear-thinning, (2) shear-thinning with negligible elasticity, and (3) shear-thinning with elasticity. At the step-out frequency, the fluidic environment's viscoelastic properties allow swimmers to propel faster than their Newtonian swimming speed, although their swimming gait remains similar. Micro-particle image velocimetry is also implemented to provide insight into how shear-thinning viscosity fluids with elasticity can modify the flow fields of the self-assembled magnetic swimmers. Our findings reveal that flow asymmetry can be created for symmetric swimmers through either the confinement effect or the Weissenberg effect. For pseudo-chiral swimmers in shear-thinning fluids, only three bead swimmers show swimming enhancement, while four bead swimmers always have a decreased step-out frequency velocity compared to their dynamics in water.

     
    more » « less
  2. Microorganisms are commonly found swimming in complex biological fluids such as mucus and these fluids respond elastically to deformation. These viscoelastic fluids have been previously shown to affect the swimming kinematics of these microorganisms in non-trivial ways depending on the rheology of the fluid, the particular swimming gait and the structural properties of the immersed body. In this report we put forth a previously unmentioned mechanism by which swimming organisms can experience a speed increase in a viscoelastic fluid. Using numerical simulations and asymptotic theory we find that significant swirling flow around a microscopic swimmer couples with the elasticity of the fluid to generate a marked increase in the swimming speed. We show that the speed enhancement is related to the introduction of mixed flow behind the swimmer and the presence of hoop stresses along its body. Furthermore, this effect persists when varying the fluid rheology and when considering different swimming gaits. This, combined with the generality of the phenomenon (i.e. the coupling of vortical flow with fluid elasticity near a microscopic swimmer), leads us to believe that this method of speed enhancement could be present for a wide range of microorganisms moving through complex fluids. 
    more » « less
  3. We conduct experiments with flexible swimmers to address the impact of fluid viscoelasticity on their locomotion. The swimmers are composed of a magnetic head actuated in rotation by a frequency-controlled magnetic field and a flexible tail whose deformation leads to forward propulsion. We consider both viscous Newtonian and glucose-based Boger fluids with similar viscosities. We find that the elasticity of the fluid systematically enhances the locomotion speed of the swimmer and that this enhancement increases with Deborah number. Using particle image velocimetry to visualize the flow field, we find a significant difference in the amount of shear between the rear and leading parts of the swimmer head. We conjecture that viscoelastic normal stresses lead to a net elastic forces in the swimming direction and thus a faster swimming speed.

     
    more » « less
  4. null (Ed.)
    We study bubble motion in a vertical capillary tube under an external flow. Bretherton ( J. Fluid Mech. , vol. 10, issue 2, 1961, pp. 166–188) has shown that, without external flow, a bubble can spontaneously rise when the Bond number ( ${Bo} \equiv \rho g R^2 / \gamma$ ) is above the critical value ${Bo}_{cr}=0.842$ , where $\rho$ is the liquid density, $g$ the gravitational acceleration, $R$ the tube radius and $\gamma$ the surface tension. It was then shown by Magnini et al. ( Phys. Rev. Fluids , vol. 4, issue 2, 2019, 023601) that the presence of an imposed liquid flow, in the same (upward) direction as buoyancy, accelerates the bubble and thickens the liquid film around it. In this work we carry out a systematic study of the bubble motion under a wide range of upward and downward external flows, focusing on the inertialess regime with Bond numbers above the critical value. We show that a rich variety of bubble dynamics occurs when an external downward flow is applied, opposing the buoyancy-driven rise of the bubble. We reveal the existence of a critical capillary number of the external downward flow ( ${Ca}_l \equiv \mu U_l/\gamma$ , where $\mu$ is the fluid viscosity and $U_l$ is the mean liquid speed) at which the bubble arrests and changes its translational direction. Depending on the relative direction of gravity and the external flow, the thickness of the film separating the bubble surface and the tube inner wall follows two distinct solution branches. The results from theory, experiments and numerical simulations confirm the existence of the two solution branches and reveal that the two branches overlap over a finite range of ${Ca}_l$ , thus suggesting non-unique, history-dependent solutions for the steady-state film thickness under the same external flow conditions. Furthermore, inertialess symmetry-breaking shape profiles at steady state are found as the bubble transits near the tipping points of the solution branches, which are shown in both experiments and three-dimensional numerical simulations. 
    more » « less
  5. We analyse the pressure-driven flow of the Oldroyd-B fluid in slowly varying arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $\Delta p$ . We first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ( $De$ ), we provide analytical expressions for the velocity, stress and the $q$ – $\Delta p$ relation in the weakly viscoelastic limit up to $O(De^2)$ . Furthermore, we exploit the reciprocal theorem derived by Boyko $\&$ Stone ( Phys. Rev. Fluids , vol. 6, 2021, L081301) to obtain the $q$ – $\Delta p$ relation at the next order, $O(De^3)$ , using only the velocity and stress fields at the previous orders. We validate our analytical results with two-dimensional numerical simulations in the case of a hyperbolic, symmetric contracting channel and find excellent agreement. While the velocity remains approximately Newtonian in the weakly viscoelastic limit (i.e. the theorem of Tanner and Pipkin), we reveal that the pressure drop strongly depends on the viscoelastic effects and decreases with $De$ . We elucidate the relative importance of different terms in the momentum equation contributing to the pressure drop along the symmetry line and identify that a pressure drop reduction for narrow contracting geometries is primarily due to gradients in the viscoelastic shear stresses. We further show that, although for narrow geometries the viscoelastic axial stresses are negligible along the symmetry line, they are comparable or larger than shear stresses in the rest of the domain. 
    more » « less