18,454 research outputs found
Hydrodynamics of Suspensions of Passive and Active Rigid Particles: A Rigid Multiblob Approach
We develop a rigid multiblob method for numerically solving the mobility
problem for suspensions of passive and active rigid particles of complex shape
in Stokes flow in unconfined, partially confined, and fully confined
geometries. As in a number of existing methods, we discretize rigid bodies
using a collection of minimally-resolved spherical blobs constrained to move as
a rigid body, to arrive at a potentially large linear system of equations for
the unknown Lagrange multipliers and rigid-body motions. Here we develop a
block-diagonal preconditioner for this linear system and show that a standard
Krylov solver converges in a modest number of iterations that is essentially
independent of the number of particles. For unbounded suspensions and
suspensions sedimented against a single no-slip boundary, we rely on existing
analytical expressions for the Rotne-Prager tensor combined with a fast
multipole method or a direct summation on a Graphical Processing Unit to obtain
an simple yet efficient and scalable implementation. For fully confined
domains, such as periodic suspensions or suspensions confined in slit and
square channels, we extend a recently-developed rigid-body immersed boundary
method to suspensions of freely-moving passive or active rigid particles at
zero Reynolds number. We demonstrate that the iterative solver for the coupled
fluid and rigid body equations converges in a bounded number of iterations
regardless of the system size. We optimize a number of parameters in the
iterative solvers and apply our method to a variety of benchmark problems to
carefully assess the accuracy of the rigid multiblob approach as a function of
the resolution. We also model the dynamics of colloidal particles studied in
recent experiments, such as passive boomerangs in a slit channel, as well as a
pair of non-Brownian active nanorods sedimented against a wall.Comment: Under revision in CAMCOS, Nov 201
Large-scale simulation of steady and time-dependent active suspensions with the force-coupling method
We present a new development of the force-coupling method (FCM) to address
the accurate simulation of a large number of interacting micro-swimmers. Our
approach is based on the squirmer model, which we adapt to the FCM framework,
resulting in a method that is suitable for simulating semi-dilute squirmer
suspensions. Other effects, such as steric interactions, are considered with
our model. We test our method by comparing the velocity field around a single
squirmer and the pairwise interactions between two squirmers with exact
solutions to the Stokes equations and results given by other numerical methods.
We also illustrate our method's ability to describe spheroidal swimmer shapes
and biologically-relevant time-dependent swimming gaits. We detail the
numerical algorithm used to compute the hydrodynamic coupling between a large
collection () of micro-swimmers. Using this methodology, we
investigate the emergence of polar order in a suspension of squirmers and show
that for large domains, both the steady-state polar order parameter and the
growth rate of instability are independent of system size. These results
demonstrate the effectiveness of our approach to achieve near continuum-level
results, allowing for better comparison with experimental measurements while
complementing and informing continuum models.Comment: 37 pages, 21 figure
On the transition to turbulence of wall-bounded flows in general, and plane Couette flow in particular
The main part of this contribution to the special issue of EJM-B/Fluids
dedicated to Patrick Huerre outlines the problem of the subcritical transition
to turbulence in wall-bounded flows in its historical perspective with emphasis
on plane Couette flow, the flow generated between counter-translating parallel
planes. Subcritical here means discontinuous and direct, with strong
hysteresis. This is due to the existence of nontrivial flow regimes between the
global stability threshold Re_g, the upper bound for unconditional return to
the base flow, and the linear instability threshold Re_c characterized by
unconditional departure from the base flow. The transitional range around Re_g
is first discussed from an empirical viewpoint ({\S}1). The recent
determination of Re_g for pipe flow by Avila et al. (2011) is recalled. Plane
Couette flow is next examined. In laboratory conditions, its transitional range
displays an oblique pattern made of alternately laminar and turbulent bands, up
to a third threshold Re_t beyond which turbulence is uniform. Our current
theoretical understanding of the problem is next reviewed ({\S}2): linear
theory and non-normal amplification of perturbations; nonlinear approaches and
dynamical systems, basin boundaries and chaotic transients in minimal flow
units; spatiotemporal chaos in extended systems and the use of concepts from
statistical physics, spatiotemporal intermittency and directed percolation,
large deviations and extreme values. Two appendices present some recent
personal results obtained in plane Couette flow about patterning from numerical
simulations and modeling attempts.Comment: 35 pages, 7 figures, to appear in Eur. J. Mech B/Fluid
Modeling hydrodynamic self-propulsion with Stokesian Dynamics. Or teaching Stokesian Dynamics to swim
We develop a general framework for modeling the hydrodynamic self-propulsion (i.e., swimming) of bodies (e.g., microorganisms) at low Reynolds number via Stokesian Dynamics simulations. The swimming body is composed of many spherical particles constrained to form an assembly that
deforms via relative motion of its constituent particles. The resistance tensor describing the hydrodynamic interactions among the individual particles maps directly onto that for the assembly. Specifying a particular swimming gait and imposing the condition that the swimming body is force- and torque-free determine the propulsive speed. The body’s translational and rotational
velocities computed via this methodology are identical in form to that from the classical theory for the swimming of arbitrary bodies at low Reynolds number. We illustrate the generality of the method through simulations of a wide array of swimming bodies: pushers and pullers, spinners, the
Taylor=Purcell swimming toroid, Taylor’s helical swimmer, Purcell’s three-link swimmer, and an amoeba-like body undergoing large-scale deformation. An open source code is a part of the supplementary material and can be used to simulate the swimming of a body with arbitrary geometry and swimming gait
- …