-
Active particles disturb the fluid around them as force dipoles, or
stresslets, which govern their collective dynamics. Unlike swimming speeds, the
stresslets of active particles are rarely determined due to the lack of a
suitable theoretical framework for arbitrary geometry. We propose a general
method, based on the reciprocal theorem of Stokes flows, to compute stresslets
as integrals of the velocities on the particle's surface, which we illustrate
for spheroidal chemically-active particles. Our method will allow tuning the
stresslet of artificial swimmers and tailoring their collective motion in
complex environments.
-
A popular method to induce synthetic propulsion at the microscale is to use
the forces created by surface-produced gas bubbles inside the asymmetric body
of a catalytic swimmer (referred to in the literature as microrocket). Gas
bubbles nucleate and grow within the catalytic swimmer and migrate toward one
of its opening under the effect of asymmetric geometric confinement, thus
generating a net hydrodynamic force which propels the device. In this paper we
use numerical simulations to develop a joint chemical (diffusive) and
hydrodynamic (Stokes) analysis of the bubble growth within a conical catalytic
microrocket and of the associated bubble and microrocket motion. Our
computational model allows us to solve for the bubble dynamics over one full
bubble cycle ranging from its nucleation to its exiting the conical rocket and
therefore to identify the propulsion characteristics as function of all design
parameters, including geometry and chemical activity of the motor, surface
tension phenomena, and all physicochemical constants. Our results suggest that
hydrodynamics and chemistry partially decouple in the motion of the bubbles,
with hydrodynamics determining the distance travelled by the microrocket over
each cycle while chemistry setting the bubble ejection frequency. Our numerical
model finally allows us to identify an optimal microrocket shape and size for
which the swimming velocity (distance travelled per cycle duration) is
maximized.
-
A microscopic bubble of soluble gas always dissolves in finite time in an
under-saturated fluid. This diffusive process is driven by the difference
between the gas concentration near the bubble, whose value is governed by the
internal pressure through Henry's law, and the concentration in the far field.
The presence of neighbouring bubbles can significantly slow down this process
by increasing the effective background concentration and reducing the diffusing
flux of dissolved gas experienced by each bubble. We develop theoretical
modelling of such diffusive shielding process in the case of small microbubbles
whose internal pressure is dominated by Laplace pressure. We first use an exact
semi-analytical solution to capture the case of two bubbles, analyse in detail
the shielding effect and show that hydrodynamic effects are mostly negligible
except in the case of almost-touching bubbles. In order to tackle the case of
multiple bubbles, we then derive and validate two analytical approximate yet
generic frameworks, first using the method of reflections and then by proposing
a self-consistent continuum description. Using both modelling frameworks, we
examine the dissolution of regular 1D, 2D and 3D bubble lattices. Bubbles
located at the edge of the lattices dissolve first, while innermost bubbles
benefit from the diffusive shielding effect, leading to the inward propagation
of a dissolution front within the lattice. We show that diffusive shielding
leads to severalfold increases in the dissolution time which grows
logarithmically with the number of bubbles in 1D-lattices and algebraically in
two and three dimensions, scaling respectively as its square root and
2/3-power. We further illustrate the sensitivity of the dissolution patterns to
initial fluctuations in bubble size or arrangement in the case of large and
dense lattices, as well as non-intuitive oscillatory effects.
-
Janus particles with the ability to move phoretically in self-generated
chemical concentration gradients are model systems for active matter. Their
motion typically consists of straight paths with rotational diffusion being the
dominant reorientation mechanism. In this paper, we show theoretically that by
a suitable surface coverage of both activity and mobility, translational and
rotational motion can be induced arbitrarily in three dimensions. The resulting
trajectories are in general helical, and their pitch and radius can be
controlled by adjusting the angle between the translational and angular
velocity. Building on the classical mathematical framework for axisymmetric
self-phoretic motion under fixed-flux chemical boundary condition, we first
show how to calculate the most general three-dimensional motion for an
arbitrary surface coverage of a spherical particle. After illustrating our
results on surface distributions, we next introduce a simple intuitive patch
model to serve as a guide for designing arbitrary phoretic spheres.
-
Singular solutions of the Stokes equations play important roles in a variety
of fluid dynamics problems. They allow the calculation of exact flows, are the
basis of the boundary integral methods used in numerical computations, and can
be exploited to derive asymptotic flows in a wide range of physical problems.
The most fundamental singular solution is the flow's Green function due to a
point force, termed the Stokeslet. Its expression is classical both in free
space and near a flat surface. Motivated by problems in biological physics
occurring near corners, we derive in this paper the asymptotic behaviour for
the Stokeslet both near and far from a corner geometry by using complex
analysis on a known double integral solution for corner flows. We investigate
all possible orientations of the point force relative to the corner and all
corner geometries from acute to obtuse. The case of salient corners is also
addressed for point forces aligned with both walls. We use experiments on beads
sedimenting in corn syrup to qualitatively test the applicability of our
results. The final results and scaling laws will allow to address the role of
hydrodynamic interactions in problems from colloidal science to microfluidics
and biological physics.
-
Recent experiments proposed to use confined bacteria in order to generate
flows near surfaces. We develop a mathematical and a computational model of
this fluid transport using a linear superposition of fundamental flow
singularities. The rotation of a helical bacterial flagellum induces both a
force and a torque on the surrounding fluid, both of which lead to a net flow
along the surface. The combined flow is in general directed at an angle to the
axis of the flagellar filament. The optimal pumping is thus achieved when
bacteria are tilted with respect to the direction in which one wants to move
the fluid, in good agreement with experimental results. We further investigate
the optimal helical shapes to be used as micropumps near surfaces and show that
bacterial flagella are nearly optimal, a result which could be relevant to the
expansion of bacterial swarms.
-
The bacterium Helicobacter pylori causes ulcers in the stomach of humans by
invading mucus layers protecting epithelial cells. It does so by chemically
changing the rheological properties of the mucus from a high-viscosity gel to a
low-viscosity solution in which it may self-propel. We develop a two-fluid
model for this process of swimming under self-generated confinement. We solve
exactly for the flow and the locomotion speed of a spherical swimmer located in
a spherically symmetric system of two Newtonian fluids whose boundary moves
with the swimmer. We also treat separately the special case of an immobile
outer fluid. In all cases, we characterise the flow fields, their spatial
decay, and the impact of both the viscosity ratio and the degree of confinement
on the locomotion speed of the model swimmer. The spatial decay of the flow
retains the same power-law decay as for locomotion in a single fluid but with a
decreased magnitude. Independently of the assumption chosen to characterise the
impact of confinement on the actuation applied by the swimmer, its locomotion
speed always decreases with an increase in the degree of confinement. Our
modelling results suggest that a low-viscosity region of at least six times the
effective swimmer size is required to lead to swimming with speeds similar to
locomotion in an infinite fluid, corresponding to a region of size above
$\approx 25~\mu$m for Helicobacter pylori.
-
The low-Reynolds number hydrodynamics of slender ribbons is accurately
captured by slender-ribbon theory, an asymptotic solution to the Stokes
equation which assumes that the three length scales characterising the ribbons
are well separated. We show in this paper that the force distribution across
the width of an isolated ribbon located in a infinite fluid can be determined
analytically, irrespective of the ribbon's shape. This, in turn, reduces the
surface integrals in the slender-ribbon theory equations to a line integral
analogous to the one arising in slender-body theory to determine the dynamics
of filaments. This result is then used to derive analytical solutions to the
motion of a rigid plate ellipsoid and a ribbon torus and to propose a ribbon
resistive-force theory, thereby extending the resistive-force theory for
slender filaments.
-
It is known from the wave-like motion of microtubules in motility assays that
the piconewton forces that motors produce can be sufficient to bend the
filaments. In cellular phenomena such as cytosplasmic streaming, molecular
motors translocate along cytoskeletal filaments, carrying cargo which entrains
fluid. When large numbers of such forced filaments interact through the
surrounding fluid, as in particular stages of oocyte development in
$Drosophila~melanogaster$, complex dynamics are observed, but the detailed
mechanics underlying them has remained unclear. Motivated by these
observations, we study here perhaps the simplest model for these phenomena: an
elastic filament, pinned at one end, acted on by a molecular motor treated as a
point force. Because the force acts tangential to the filament, no matter what
its shape, this "follower-force" problem is intrinsically non-variational, and
thereby differs fundamentally from Euler buckling, where the force has a fixed
direction, and which, in the low Reynolds number regime, ultimately leads to a
stationary, energy-minimizing shape. Through a combination of linear stability
theory, analytical study of a solvable simplified "two-link" model, and
numerical studies of the full elastohydrodynamic equations of motion we
elucidate the Hopf bifurcation that occurs with increasing forcing of a
filament, leading to flapping motion analogous to the high Reynolds number
oscillations of a garden hose with a free end.
-
Many cells exploit the bending or rotation of flagellar filaments in order to
self-propel in viscous fluids. While appropriate theoretical modelling is
available to capture flagella locomotion in simple, Newtonian fluids,
formidable computations are required to address theoretically their locomotion
in complex, nonlinear fluids, e.g.~mucus. Based on experimental measurements
for the motion of rigid rods in non-Newtonian fluids and on the classical
Carreau fluid model, we propose empirical extensions of the classical Newtonian
resistive-force theory to model the waving of slender filaments in
non-Newtonian fluids. By assuming the flow near the flagellum to be locally
Newtonian, we propose a self-consistent way to estimate the typical shear-rate
in the fluid, which we then use to construct correction factors to the
Newtonian local drag coefficients. The resulting non-Newtonian resistive-force
theory, while empirical, is consistent with the Newtonian limit, and with the
experiments. We then use our models to address waving locomotion in
non-Newtonian fluids, and show that the resulting swimming speeds are
systematically lowered -- a result which we are able to capture asymptotically
and to interpret physically. An application of the models to recent
experimental results on the locomotion of Caenorhabditis elegans in polymeric
solutions shows reasonable agreement and thus captures the main physics of
swimming in shear-thinning fluids.
-
Inspired by recent experiments using synthetic microswimmers to manipulate
droplets, we investigate the low-Reynolds-number locomotion of a model swimmer
(a spherical squirmer) encapsulated inside a droplet of comparable size in
another viscous fluid. Meditated solely by hydrodynamic interactions, the
encaged swimmer is seen to be able to propel the droplet, and in some
situations both remain in a stable co-swimming state. The problem is tackled
using both an exact analytical theory and a numerical implementation based on
boundary element method, with a particular focus on the kinematics of the
co-moving swimmer and droplet in a concentric configuration, and we obtain
excellent quantitative agreement between the two. The droplet always moves
slower than a swimmer which uses purely tangential surface actuation but when
it uses a particular combination of tangential and normal actuations, the
squirmer and droplet are able to attain a same velocity and stay concentric for
all times. We next employ numerical simulations to examine the stability of
their concentric co-movement, and highlight several stability scenarios
depending on the particular gait adopted by the swimmer. Furthermore, we show
that the droplet reverses the nature of the far-field flow induced by the
swimmer: a droplet cage turns a pusher swimmer into a puller, and vice versa.
Our work sheds light on the potential development of droplets as self-contained
carriers of both chemical content and self-propelled devices for controllable
and precise drug deliveries.
-
Experiments involving the two-dimensional passive diffusion of colloidal
boomerangs tracked off their centre of mobility have shown striking
non-Gaussian tails in their probability distribution function [Chakrabarty et
al., Soft Matter 12, 4318 (2016)]. This in turn can lead to anomalous diffusion
characteristics, including mean drift. In this paper, we develop a general
theoretical explanation for these measurements. The idea relies on calculating
the two-dimensional probability densities at the centre of mobility of the
particle, where all distributions are Gaussian, and then transforming them to a
different reference point. Our model clearly captures the experimental results,
without any fitting parameters, and demonstrates that the one-dimensional
probability distributions may also exhibit strongly non-Gaussian tops. These
results indicate that the choice of tracking point can cause a considerable
departure from Gaussian statistics, potentially causing some common modelling
techniques to fail.
-
Phoretic swimmers provide new avenues to study non-equilibrium statistical
physics and are also hailed as a promising technology for bioengineering at the
cellular scale. Exact solutions for the locomotion of such swimmers have been
restricted so far to spheroidal shapes. In this paper we solve for the flow
induced by the canonical non-simply connected shape, namely an axisymmetric
phoretic torus. The analytical solution takes the form of an infinite series
solution, which we validate against boundary element computations. For a torus
of uniform chemical activity, confinement effects in the hole allow the torus
to act as a pump, which we optimize subject to fixed particle surface area.
Under the same constraint, we next characterize the fastest swimming Janus
torus for a variety of assumptions on the surface chemistry. Perhaps
surprisingly, none of the optimal tori occur in the limit where the central
hole vanishes.
-
One of the hallmarks of active matter is its rich nonlinear dynamics and
instabilities. Recent numerical simulations of phototactic algae showed that a
thin jet of swimmers, obtained from hydrodynamic focusing inside a Poiseuille
flow, was unstable to longitudinal perturbations with swimmers dynamically
clustering (Jibuti et al., Phys. Rev. E, 90, 2014). As a simple starting point
to understand these instabilities, we consider in this paper an initially
homogeneous one-dimensional line of aligned swimmers moving along the same
direction, and characterise its instability using both a continuum framework
and a discrete approach. In both cases, we show that hydrodynamic interactions
between the swimmers lead to instabilities in density for which we compute the
growth rate analytically. Lines of pusher-type swimmers are predicted to remain
stable while lines of pullers (such as flagellated algae) are predicted to
always be unstable.
-
Catalytic swimmers have attracted much attention as alternatives to
biological systems for examining collective microscopic dynamics and the
response to physico-chemical signals. Yet, understanding and predicting even
the most fundamental characteristics of their individual propulsion still
raises important challenges. While chemical asymmetry is widely recognized as
the cornerstone of catalytic propulsion, different experimental studies have
reported that particles with identical chemical properties may propel in
opposite directions. Here, we show that, beyond its chemical properties, the
detailed shape of a catalytic swimmer plays an essential role in determining
its direction of motion, demonstrating the compatibility of the classical
theoretical framework with experimental observations.
-
Artificial phoretic particles swim using self-generated gradients in chemical
species (self-diffusiophoresis) or charges and currents (self-electrophoresis).
These particles can be used to study the physics of collective motion in active
matter and might have promising applications in bioengineering. In the case of
self-diffusiophoresis, the classical physical model relies on a steady solution
of the diffusion equation, from which chemical gradients, phoretic flows and
ultimately the swimming velocity, may be derived. Motivated by disk-shaped
particles in thin films and under confinement, we examine the extension to two
dimensions. Because the two-dimensional diffusion equation lacks a steady state
with the correct boundary conditions, Laplace transforms must be used to study
the long-time behavior of the problem and determine the swimming velocity. For
fixed chemical fluxes on the particle surface, we find that the swimming
velocity ultimately always decays logarithmically in time. In the case of
finite Peclet numbers, we solve the full advection-diffusion equation
numerically and show that this decay can be avoided by the particle moving to
regions of unconsumed reactant. Finite advection thus regularizes the
two-dimensional phoretic problem.
-
Swimming microorganisms often have to propel in complex, non-Newtonian
fluids. We carry out experiments with self-propelling helical swimmers driven
by an externally rotating magnetic field in shear-thinning, inelastic fluids.
Similarly to swimming in a Newtonian fluid, we obtain for each fluid a
locomotion speed which scales linearly with the rotation frequency of the
swimmer, but with a prefactor which depends on the power index of the fluid.
The fluid is seen to always increase the swimming speed of the helix, up to 50%
faster and thus the strongest of such type reported to date. The maximum
relative increase for a fluid power index of around 0.6. Using simple scalings,
we argue that the speed increase is not due to the local decrease of the flow
viscosity around the helical filament but hypothesise instead that it
originates from confinement-like effect due to viscosity stratification around
the swimmer.
-
Cellular biology abound with filaments interacting through fluids, from
intracellular microtubules, to rotating flagella and beating cilia. While
previous work has demonstrated the complexity of capturing nonlocal
hydrodynamic interactions between moving filaments, the problem remains
difficult theoretically. We show here that when filaments are closer to each
other than their relevant length scale, the integration of hydrodynamic
interactions can be approximately carried out analytically. This leads to a set
of simplified local equations, illustrated on a simple model of two interacting
filaments, which can be used to tackle theoretically a range of problems in
biology and physics.
-
Numerical study of the hydrodynamics of thin sheets and ribbons presents
difficulties associated with resolving multiple length scales. To circumvent
these difficulties, asymptotic methods have been developed to describe the
dynamics of slender fibres and ribbons. However, such theories entail
restrictions on the shapes that can be studied, and often break down in regions
where standard boundary element methods are still impractical. In this paper we
develop a regularised stokeslet method for ribbons and sheets in order to
bridge the gap between asymptotic and boundary element methods. The method is
validated against the analytical solution for plate ellipsoids, as well as the
dynamics of ribbon helices and an experimental microswimmer. We then
demonstrate the versatility of this method by calculating the flow around a
double helix, and the swimming dynamics of a microscale "magic carpet".
-
It is now well established that nearby beating pairs of eukaryotic flagella
or cilia typically synchronize in phase. A substantial body of evidence
supports the hypothesis that hydrodynamic coupling between the active
filaments, combined with waveform compliance, provides a robust mechanism for
synchrony. This elastohydrodynamic mechanism has been incorporated into
`bead-spring' models in which flagella are represented by microspheres tethered
by radial springs as internal forces drive them about orbits. While these
low-dimensional models reproduce the phenomenon of synchrony, their parameters
are not readily relatable to those of flagella. More realistic models which
reflect the elasticity of the axonemes and active force generation take the
form of fourth-order nonlinear PDEs. While computational studies have shown
synchrony, the effects of hydrodynamic coupling between nearby filaments
governed by such models have been theoretically examined only in the regime of
interflagellar distances $d$ large compared to flagellar length $L$. Yet, in
many biological situations $d/L \ll 1$. Here, we first present an asymptotic
analysis of the hydrodynamic coupling between two filaments in the regime $d/L
\ll 1$, and find that the form of the coupling is independent of the details of
the internal forces that govern the motion of the filaments. The analysis is
like the localized induction approximation for vortex filament motion, extended
to the case of mutual induction. To understand how the coupling mechanism leads
to synchrony of extended objects, we introduce a heuristic model of flagellar
beating, a single fourth-order nonlinear PDE whose form is derived from
symmetry considerations, the physics of elasticity, and the overdamped nature
of the dynamics. Analytical and numerical studies of this model illustrate how
synchrony between two filaments is achieved through the asymptotic coupling.
-
Flagellated bacteria on nutrient-rich substrates can differentiate into a
swarming state and move in dense swarms across surfaces. A recent experiment
measured the flow in the fluid around an Escherichia coli swarm (Wu, Hosu and
Berg, 2011 Proc. Natl. Acad. Sci. USA 108 4147). A systematic chiral flow was
observed in the clockwise direction (when viewed from above) ahead of the swarm
with flow speeds of about $10~\mu$m/s, about 3 times greater than the radial
velocity at the edge of the swarm. The working hypothesis is that this flow is
due to the action of cells stalled at the edge of a colony that extend their
flagellar filaments outwards, moving fluid over the virgin agar. In this work
we quantitatively test his hypothesis. We first build an analytical model of
the flow induced by a single flagellum in a thin film and then use the model,
and its extension to multiple flagella, to compare with experimental
measurements. The results we obtain are in agreement with the flagellar
hypothesis. The model provides further quantitative insight into the flagella
orientations and their spatial distributions as well as the tangential speed
profile. In particular, the model suggests that flagella are on average
pointing radially out of the swarm and are not wrapped tangentially.
-
Swimming cells and microorganisms are a critical component of many biological
processes. In order to better interpret experimental studies of low Reynolds
number swimming, we combine experimental and numerical methods to perform an
analysis of the flow-field around the swimming nematode Caenorhabditis elegans.
We first use image processing and particle tracking velocimetry to extract the
body shape, kinematics, and flow-fields around the nematode. We then construct
a three-dimensional model using the experimental swimming kinematics and employ
a boundary element method to simulate flow-fields, obtaining very good
quantitative agreement with experiment. We use this numerical model to show
that calculation of flow shear rates using purely planar data results in
significant underestimates of the true three-dimensional value. Applying
symmetry arguments, validated against numerics, we demonstrate that the
out-of-plane contribution can be accounted for via incompressibility and
therefore simply calculated from particle tracking velocimetry. Our results
show how fundamental fluid mechanics considerations may be used to improve the
accuracy of measurements in biofluiddynamics.
-
Internal phoretic flows due to the interactions of solid boundaries with
local chemical gradients may be created using chemical patterning.
Alternatively, we demonstrate here that internal flows might also be induced by
geometric asymmetries of chemically-homogeneous surfaces. We characterise the
circulatory flow created in a cavity enclosed between two eccentric cylindrical
walls of uniform chemical activity. Local gradients of the diffusing solute
induce a slip flow along the surface of the cylinders, leading to a circulatory
bulk flow pattern which can be solved analytically in the diffusive limit. The
flow strength can be controlled by adjusting the relative positions of the
cylinders and an optimal configuration is identified. These results provide a
model system for tunable phoretic pumps.
-
The drag anisotropy of slender filaments is a critical physical property
allowing swimming in low-Reynolds number flows, and without it linear
translation is impossible. Here we show that, in contrast, net rotation can
occur under isotropic drag. We first demonstrate this result formally by
considering the consequences of the force- and torque-free conditions on
swimming bodies and we then illustrate it with two examples (a simple swimmers
made of three rods and a model bacterium with two helical flagellar filaments).
Our results highlight the different role of hydrodynamic forces in generating
translational vs.~rotational propulsion.
-
Ribbons are long narrow strips possessing three distinct material length
scales (thickness, width, and length) which allow them to produce unique shapes
unobtainable by wires or filaments. For example when a ribbon has half a twist
and is bent into a circle it produces a M\"obius strip. Significant effort has
gone into determining the structural shapes of ribbons but less is know about
their behavior in viscous fluids. In this paper we determine, asymptotically,
the leading-order hydrodynamic behavior of a slender ribbon in Stokes flows.
The derivation, reminiscent of slender-body theory for filaments, assumes that
the length of the ribbon is much larger than its width, which itself is much
larger than its thickness. The final result is an integral equation for the
force density on a mathematical ruled surface, termed the ribbon plane, located
inside the ribbon. A numerical implementation of our derivation shows good
agreement with the known hydrodynamics of long flat ellipsoids, and
successfully captures the swimming behavior of artificial microscopic swimmers
recently explored experimentally. We also study the asymptotic behavior of a
ribbon bent into a helix, that of a twisted ellipsoid, and we investigate how
accurately the hydrodynamics of a ribbon can be effectively captured by that of
a slender filament. Our asymptotic results provide the fundamental framework
necessary to predict the behavior of slender ribbons at low Reynolds numbers in
a variety of biological and engineering problems.