skip to main content


Title: Noise-Induced Aggregation of Swimmers in the Kolmogorov Flow
We investigate a model for the dynamics of ellipsoidal microswimmers in an externally imposed, laminar Kolmogorov flow. Through a phase-space analysis of the dynamics without noise, we find that swimmers favor either cross-stream or rotational drift, depending on their swimming speed and aspect ratio. When including noise, i.e., rotational diffusion, we find that swimmers are driven into certain parts of phase space, leading to a nonuniform steady-state distribution. This distribution exhibits a transition from swimmer aggregation in low-shear regions of the flow to aggregation in high-shear regions as the swimmer’s speed, aspect ratio, and rotational diffusivity are varied. To explain the nonuniform phase-space distribution of swimmers, we apply a weak-noise averaging principle that produces a reduced description of the stochastic swimmer dynamics. Using this technique, we find that certain swimmer trajectories are more favorable than others in the presence of weak rotational diffusion. By combining this information with the phase-space speed of swimmers along each trajectory, we predict the regions of phase space where swimmers tend to accumulate. The results of the averaging technique are in good agreement with direct calculations of the steady-state distributions of swimmers. In particular, our analysis explains the transition from low-shear to high-shear aggregation.  more » « less
Award ID(s):
1806355 1825379
NSF-PAR ID:
10337598
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Frontiers in Physics
Volume:
9
ISSN:
2296-424X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ . 
    more » « less
  2. Magnetically-actuated swimming microrobots are an emerging tool for navigating and manipulating materials in confined spaces. Recent work has demonstrated that it is possible to build such systems at the micro and nanoscales using polymer microspheres, magnetic particles and DNA nanotechnology. However, while these materials enable an unprecedented ability to build at small scales, such systems often demonstrate significant polydispersity resulting from both the material variations and the assembly process itself. This variability makes it difficult to predict, let alone optimize, the direction or magnitude of microswimmer velocity from design parameters such as link shape or aspect ratio. To isolate questions of a swimmer's design from variations in its physical dimensions, we present a novel experimental platform using two-photon polymerization to build a two-link, buoyant milliswimmer with a fully customizable shape and integrated flexible linker (the swimmer is underactuated, enabling asymmetric cyclic motion and net translation). Our approach enables us to control both swimming direction and repeatability of swimmer performance. These studies provide ground truth data revealing that neither the first order nor second order models currently capture the key features of milliswimmer performance. We therefore use our experimental platform to develop design guidelines for tuning the swimming speeds, and we identify the following three approaches for increasing speed: (1) tuning the actuation frequency for a fixed aspect ratio, (2) adjusting the aspect ratio given a desired range of operating frequencies, and (3) using the weaker value of linker stiffness from among the values that we tested, while still maintaining a robust connection between the links. We also find experimentally that spherical two-link swimmers with dissimilar link diameters achieve net velocities comparable to swimmers with cylindrical links, but that two-link spherical swimmers of equal diameter do not. 
    more » « less
  3. The rheology of suspensions of rings (tori) rotating in an unbounded low Reynolds number simple shear flow is calculated using numerical simulations at dilute particle number densities ( n ≪ 1 ). Suspensions of non-Brownian rings are studied by computing pair interactions that include hydrodynamic interactions modeled using slender body theory and particle collisions modeled using a short-range repulsive force. Particle contact and hydrodynamic interactions were found to have comparable influences on the steady-state Jeffery orbit distribution. The average tilt of the ring away from the flow-vorticity plane increased during pairwise interactions compared to the tilt associated with Jeffery rotation and the steady-state orbit distribution. Particle stresses associated with the increased tilt during the interaction were found to be comparable to the stresses induced directly by particle contact forces and the hydrodynamic velocity disturbances of other particles. The hydrodynamic diffusivity coefficients in the gradient and vorticity directions were also obtained and were found to be two orders of magnitude larger than the corresponding values in fiber suspensions at the same particle concentrations. Rotary Brownian dynamics simulations of isolated Brownian rings were used to understand the shear rate dependence of suspension rheology. The orbit distribution observed in the regime of weak Brownian motion, P e ≫ ϕ T − 3, was surprisingly similar to that obtained from pairwise interaction calculations of non-Brownian rings. Here, the Peclet number P e is the ratio of the shear rate and the rotary diffusivity of the particle and ϕ T is the effective inverse-aspect ratio of the particle (approximately equal to 2 π times the inverse of its non-dimensional Jeffery time period). Thus, the rheology results obtained from pairwise interactions should retain accuracy even for weakly Brownian rings ( n ≪ 1 and ϕ T − 3 ≪ P e ). 
    more » « less
  4. 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
  5. 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