4,523 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
A fluctuating boundary integral method for Brownian suspensions
We present a fluctuating boundary integral method (FBIM) for overdamped
Brownian Dynamics (BD) of two-dimensional periodic suspensions of rigid
particles of complex shape immersed in a Stokes fluid. We develop a novel
approach for generating Brownian displacements that arise in response to the
thermal fluctuations in the fluid. Our approach relies on a first-kind boundary
integral formulation of a mobility problem in which a random surface velocity
is prescribed on the particle surface, with zero mean and covariance
proportional to the Green's function for Stokes flow (Stokeslet). This approach
yields an algorithm that scales linearly in the number of particles for both
deterministic and stochastic dynamics, handles particles of complex shape,
achieves high order of accuracy, and can be generalized to three dimensions and
other boundary conditions. We show that Brownian displacements generated by our
method obey the discrete fluctuation-dissipation balance relation (DFDB). Based
on a recently-developed Positively Split Ewald method [A. M. Fiore, F. Balboa
Usabiaga, A. Donev and J. W. Swan, J. Chem. Phys., 146, 124116, 2017],
near-field contributions to the Brownian displacements are efficiently
approximated by iterative methods in real space, while far-field contributions
are rapidly generated by fast Fourier-space methods based on fluctuating
hydrodynamics. FBIM provides the key ingredient for time integration of the
overdamped Langevin equations for Brownian suspensions of rigid particles. We
demonstrate that FBIM obeys DFDB by performing equilibrium BD simulations of
suspensions of starfish-shaped bodies using a random finite difference temporal
integrator.Comment: Submitted to J. Comp. Phy
Inertial Coupling Method for particles in an incompressible fluctuating fluid
We develop an inertial coupling method for modeling the dynamics of
point-like 'blob' particles immersed in an incompressible fluid, generalizing
previous work for compressible fluids. The coupling consistently includes
excess (positive or negative) inertia of the particles relative to the
displaced fluid, and accounts for thermal fluctuations in the fluid momentum
equation. The coupling between the fluid and the blob is based on a no-slip
constraint equating the particle velocity with the local average of the fluid
velocity, and conserves momentum and energy. We demonstrate that the
formulation obeys a fluctuation-dissipation balance, owing to the
non-dissipative nature of the no-slip coupling. We develop a spatio-temporal
discretization that preserves, as best as possible, these properties of the
continuum formulation. In the spatial discretization, the local averaging and
spreading operations are accomplished using compact kernels commonly used in
immersed boundary methods. We find that the special properties of these kernels
make the discrete blob a particle with surprisingly physically-consistent
volume, mass, and hydrodynamic properties. We develop a second-order
semi-implicit temporal integrator that maintains discrete
fluctuation-dissipation balance, and is not limited in stability by viscosity.
Furthermore, the temporal scheme requires only constant-coefficient Poisson and
Helmholtz linear solvers, enabling a very efficient and simple FFT-based
implementation on GPUs. We numerically investigate the performance of the
method on several standard test problems...Comment: Contains a number of corrections and an additional Figure 7 (and
associated discussion) relative to published versio
Stationary shapes of deformable particles moving at low Reynolds numbers
Lecture Notes of the Summer School ``Microswimmers -- From Single Particle
Motion to Collective Behaviour'', organised by the DFG Priority Programme SPP
1726 (Forschungszentrum J{\"{u}}lich, 2015).Comment: Pages C7.1-16 of G. Gompper et al. (ed.), Microswimmers - From Single
Particle Motion to Collective Behaviour, Lecture Notes of the DFG SPP 1726
Summer School 2015, Forschungszentrum J\"ulich GmbH, Schriften des
Forschungszentrums J\"ulich, Reihe Key Technologies, Vol 110, ISBN
978-3-95806-083-
Numerical electrokinetics
A new lattice method is presented in order to efficiently solve the
electrokinetic equations, which describe the structure and dynamics of the
charge cloud and the flow field surrounding a single charged colloidal sphere,
or a fixed array of such objects. We focus on calculating the electrophoretic
mobility in the limit of small driving field, and systematically linearise the
equations with respect to the latter. This gives rise to several subproblems,
each of which is solved by a specialised numerical algorithm. For the total
problem we combine these solvers in an iterative procedure. Applying this
method, we study the effect of the screening mechanism (salt screening vs.
counterion screening) on the electrophoretic mobility, and find a weak
non-trivial dependence, as expected from scaling theory. Furthermore, we find
that the orientation of the charge cloud (i. e. its dipole moment) depends on
the value of the colloid charge, as a result of a competition between
electrostatic and hydrodynamic effects.Comment: accepted for publication in Journal of Physics Condensed Matter
(proceedings of the 2012 CODEF conference
- …