An official website of the United States government
The development of novel drug delivery systems, which are revolutionizing modern medicine, is benefiting from studies on microorganisms’ swimming. In this paper we consider a model microorganism (a squirmer) enclosed in a viscous droplet to investigate the effects of medium heterogeneity or geometry on the propulsion speed of the caged squirmer. We first consider the squirmer and droplet to be spherical (no shape effects) and derive exact solutions for the equations governing the problem. For a squirmer with purely tangential surface velocity, the squirmer is always able to move inside the droplet (even when the latter ceases to move as a result of large fluid resistance of the heterogeneous medium). Adding radial modes to the surface velocity, we establish a new condition for the existence of a co-swimming speed (where squirmer and droplet move at the same speed). Next, to probe the effects of geometry on propulsion, we consider the squirmer and droplet to be in Newtonian fluids. For a squirmer with purely tangential surface velocity, numerical simulations reveal a strong dependence of the squirmer's speed on shapes, the size of the droplet and the viscosity contrast. We found that the squirmer speed is largest when the droplet size and squirmer's eccentricity are small, and the viscosity contrast is large. For co-swimming, our results reveal a complex, non-trivial interplay between the various factors that combine to yield the squirmer's propulsion speed. Taken together, our study provides several considerations for the efficient design of future drug delivery systems.
more »
« less
This content will become publicly available on May 9, 2026
Hydrodynamic interactions between a sedimenting squirmer and a planar wall
The hydrodynamic interactions between a sedimenting microswimmer and a solid wall have ubiquitous biological and technological applications. A plethora of gravity-induced swimming dynamics near a planar no-slip wall provide a platform for designing artificial microswimmers that can generate directed propulsion through their translation–rotation coupling near a wall. In this work, we provide exact solutions for a squirmer (a model swimmer of spherical shape with a prescribed slip velocity) facing either towards or away from a planar wall perpendicular to gravity. These exact solutions are used to validate a numerical code based on the boundary integral method with an adaptive mesh for distances from the wall down to 0.1 % of the squirmer radius. This boundary integral code is then used to investigate the rich gravity-induced dynamics near a wall, mapping out the detailed bifurcation structures of the swimming dynamics in terms of orientation and distance to the wall. Simulation results show that a squirmer may traverse the wall, move to a fixed point at a given height with a fixed orientation in a monotonic way or in an oscillatory fashion, or oscillate in a limit cycle in the presence of wall repulsion.
more »
« less
- Award ID(s):
- 1951600
- PAR ID:
- 10625106
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 1010
- ISSN:
- 0022-1120
- Subject(s) / Keyword(s):
- active matter swimming low-Reynolds-number flows
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We analyze a nonlinear partial differential equation system describing the motion of a microswimmer in a nematic liquid crystal environment. For the microswimmer’s motility, the squirmer model is used in which self-propulsion enters the model through the slip velocity on the microswimmer’s surface. The liquid crystal is described using the well-established Beris–Edwards formulation. In previous computational studies, it was shown that the squirmer, regardless of its initial configuration, eventually orients itself either parallel or perpendicular to the preferred orientation dictated by the liquid crystal. Furthermore, the corresponding solution of the coupled nonlinear system converges to a steady state. In this work, we rigorously establish the existence of the steady state and also the finite-time existence for the time-dependent problem in a periodic domain. Finally, we will use a two-scale asymptotic expansion to derive a homogenized model for the collective swimming of squirmers as they reach their steady-state orientation and speed.more » « less
-
The effect of slip surfaces on the laminar–turbulent separatrix of plane Poiseuille flow is studied by direct numerical simulation. In laminar flows, the inclusion of the slip surfaces results in a drag reduction of over 10 %, which is in good agreement with previous studies and the theory of laminar slip flows. Turbulence lifetimes, the likelihood that turbulence is sustained, is investigated for transitional flows with various slip lengths. We show that slip surfaces decrease the likelihood of sustained turbulence compared to the no-slip case, and the likelihood is further decreased as slip length is increased. A more deterministic analysis of the effects of slip surfaces on a transition to turbulence is performed by using nonlinear travelling-wave solutions to the Navier–Stokes equations, also known as exact coherent solutions. Two solution families, dubbed P3 and P4, are used since their lower-branch solutions are embedded on the boundary of the basin of attraction of laminar and turbulent flows (Park & Graham, J. Fluid Mech. , vol. 782, 2015, pp. 430–454). Additionally, they exhibit distinct flow structures – the P3 and P4 are denoted as core mode and critical-layer mode, respectively. Distinct effects of slip surfaces on the solutions are observed by the skin-friction evolution, linear growth rate and phase-space projection of transitional trajectories. The slip surface appears to modify the transition dynamics very little for the core mode, but quite considerably for the critical-layer mode. Most importantly, the slip surface promotes different transition dynamics – an early and bypass-like transition for the core mode and a delayed and H- or K-type-like transition for the critical-layer mode. We explain these distinct transition dynamics based on spatio-temporal and quadrant analyses. It is found that slip surfaces promote the prevalence of strong wall-toward motions (sweep-like events) near vortex cores close to the channel centre, inducing an early transition, while long sustained ejection events are present in the region of the $$\unicode[STIX]{x1D6EC}$$ -shaped vortex cores close to the critical layer, resulting in a delayed transition. This should motivate flow control strategies to fully exploit these distinct transition dynamics for transition to turbulence.more » « less
-
Abstract When swimming near a solid planar boundary, bio-inspired propulsors can naturally equilibrate to certain distances from that boundary. How these equilibria are affected by asymmetric swimming kinematics is unknown. We present here a study of near-boundary pitching hydrofoils based on water channel experiments and potential flow simulations. We found that asymmetric pitch kinematics do affect near-boundary equilibria, resulting in the equilibria shifting either closer to or away from the planar boundary. The magnitude of the shift depends on whether the pitch kinematics have spatial asymmetry (e.g. a bias angle, θ 0 ) or temporal asymmetry (e.g. a stroke-speed ratio, τ ). Swimming at stable equilibrium requires less active control, while shifting the equilibrium closer to the boundary can result in higher thrust with no measurable change in propulsive efficiency. Our work reveals how asymmetric kinematics could be used to fine-tune a hydrofoil’s interaction with a nearby boundary, and it offers a starting point for understanding how fish and birds use asymmetries to swim near substrates, water surfaces, and sidewalls.more » « less
-
The general system of images for regularized Stokeslets (GSIRS) developed by Cortez and Varela [] is used extensively to model Stokes flow phenomena such as microorganisms swimming near a boundary. Our collaborative team uses dynamically similar scaled macroscopic experiments to test theories for forces and torques on spheres moving near a boundary and uses these data and the method of regularized Stokeslets (MRS) created by Cortez [] to calibrate the GSIRS. We find excellent agreement between theory and experiments, which provides experimental validation of exact series solutions for spheres moving near an infinite plane boundary. We test two surface discretization methods commonly used in the literature: the 6-patch method and the spherical centroidal Voronoi tessellation (SCVT) method. Our data show that a discretization method, such as SCVT, that uniformly distributes points provides the most accurate results when the motional symmetry is broken by the presence of a boundary. We use theory and the MRS to find optimal values for the regularization parameter in free space for a given surface discretization and show that the optimal regularization parameter values can be fit with simple formulas when using the SCVT method. We also present a regularization function with higher-order accuracy when compared with the regularization function previously introduced by Cortez []. The simulated force and torque values compare very well with experiments and theory for a wide range of boundary distances. However, we find that for a fixed discretization of the sphere, the simulations lose accuracy when the gap between the edge of the sphere and the wall is smaller than the average distance between discretization points in the SCVT method. We also show an alternative method to calibrate the GSIRS to simulate sphere motion arbitrarily close to the boundary. Our computational parameters and methods along with our matlab and python implementations of the series solution of Lee and Leal [], MRS, and GSIRS provide researchers with important resources to optimize the GSIRS and other numerical methods, so that they can efficiently and accurately simulate spheres moving near a boundary.more » « less
