1,051 research outputs found
A minimal model for acoustic forces on Brownian particles
We present a generalization of the inertial coupling (IC) [Usabiaga et al. J.
Comp. Phys. 2013] which permits the resolution of radiation forces on small
particles with arbitrary acoustic contrast factor. The IC method is based on a
Eulerian-Lagrangian approach: particles move in continuum space while the fluid
equations are solved in a regular mesh (here we use the finite volume method).
Thermal fluctuations in the fluid stress, important below the micron scale, are
also taken into account following the Landau-Lifshitz fluid description. Each
particle is described by a minimal cost resolution which consists on a single
small kernel (bell-shaped function) concomitant to the particle. The main role
of the particle kernel is to interpolate fluid properties and spread particle
forces. Here, we extend the kernel functionality to allow for an arbitrary
particle compressibility. The particle-fluid force is obtained from an imposed
no-slip constraint which enforces similar particle and kernel fluid velocities.
This coupling is instantaneous and permits to capture the fast, non-linear
effects underlying the radiation forces on particles. Acoustic forces arise
either because an excess in particle compressibility (monopolar term) or in
mass (dipolar contribution) over the fluid values. Comparison with theoretical
expressions show that the present generalization of the IC method correctly
reproduces both contributions. Due to its low computational cost, the present
method allows for simulations with many particles using a standard Graphical
Processor Unit (GPU)
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
Particle hydrodynamics: from molecular to colloidal fluids
A method for particle hydrodynamics based on an hybrid Eulerian-Lagrangian approach is presented. Particles are solved in the continuum space while the fluid is solved in an Eulerian mesh, and described by finite volume fluctuating hydrodynamics. This set-up is particulary suited for micron-size devices where the Reynolds number is small but thermal fluctuations are important. The fluid-particle coupling force is obtained by imposing zero relative (particle-fluid) velocity at discrete points representing the particle sites. In this work particles are described by an only site which neglect rotation. The momentum exchanged between fluid and particle is transfered instantaneously and this brings about several benefits such as a correct treatment of inertia and proper particle velocity fluctuations uniquely driven by the fluid thermal forces. The present scheme is designed for incompressible and compressible fluids at low Mach number. This is theoretically shown by analyzing the consistency between the Eulerian and Lagrangian momentum balance.A series of tests up to moderate Reynolds number and acoustic forces under ultrasound waves are also presented.
1 INTRODUCTIO
Swimming Efficiently by Wrapping
Single flagellated bacteria are ubiquitous in nature. They exhibit various swimming
modes using their flagella to explore complex surroundings such as soil and porous
polymer networks. Some single-flagellated bacteria swim with two distinct modes, one
with its flagellum extended away from its body and another with its flagellum wrapped
around it. The wrapped mode has been observed when the bacteria swim under tight
confinements or in highly viscous polymeric melts. In this study we investigate the
hydrodynamics of these two modes inside a circular pipe. We find that the wrap mode is
slower than the extended mode in bulk but more efficient under strong confinement due
to a hydrodynamic increased of its flagellum translation-rotation coupling
Swimming Efficiently by Wrapping
Single flagellated bacteria are ubiquitous in nature. They exhibit various
swimming modes using their flagella to explore complex surroundings such as
soil and porous polymer networks. Some single-flagellated bacteria swim with
two distinct modes, one with its flagellum extended away from its body and
another with its flagellum wrapped around it. The wrapped mode has been
observed when the bacteria swim under tight confinements or in highly viscous
polymeric melts. In this study we investigate the hydrodynamics of these two
modes inside a circular pipe. We find that the wrap mode is slower than the
extended mode in bulk but more efficient under strong confinement due to a
hydrodynamic increased of its flagellum translation-rotation coupling.Comment: 10 pages, 5 figure
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
Mapping flagellated swimmers to surface-slip driven swimmers.
Flagellated microswimmers are ubiquitous in natural habitats. Understanding the hydrodynamic behavior of these cells is of paramount interest, owing to their applications in
bio-medical engineering and disease spreading. Since the last two decades, computational
efforts have been continuously improved to accurately capture the complex hydrodynamic
behavior of these model systems. However, modeling the dynamics of such swimmers with
fine details is computationally expensive due to the large number of unknowns and the small
time-steps required to solve the equations. In this work we propose a method to map fully
resolved flagellated microswimmers to coarse, active slip driven swimmers which can be simulated at a reduced computational cost. Using the new method, the slip driven swimmers
move with the same velocity, to machine precision, as the flagellated swimmers and generate a similar flow field with a controlled accuracy. The method is validated for swimming
patterns near a no-slip boundary, interactions between swimmers and scattering with large
obstacles
Computational modeling of passive transport of functionalized nanoparticles
Functionalized nanoparticles (NPs) are complex objects present in a variety of systems ranging from synthetic grafted nanoparticles to viruses. The morphology and number of the decorating groups can vary widely between systems. Thus, the modeling of functionalized NPs typically considers simplified spherical objects as a first-order approximation. At the nanoscale label, complex hydrodynamic interactions are expected to emerge as the morphological features of the particles change, and they can be further amplified when the NPs are confined or near walls. Direct estimation of these variations can be inferred via diffusion coefficients of the NPs. However, the evaluation of the coefficients requires an improved representation of the NPs morphology to reproduce important features hidden by simplified spherical models. Here, we characterize the passive transport of free and confined functionalized nanoparticles using the Rigid Multi-Blob (RMB) method. The main advantage of RMB is its versatility to approximate the mobility of complex structures at the nanoscale with significant accuracy and reduced computational cost. In particular, we investigate the effect of functional groups' distribution, size, and morphology over nanoparticle translational and rotational diffusion. We identify that the presence of functional groups significantly affects the rotational diffusion of the nanoparticles; moreover, the morphology of the groups and number induce characteristic mobility reduction compared to non-functionalized nanoparticles. Confined NPs also evidenced important alterations in their diffusivity, with distinctive signatures in the off-diagonal contributions of the rotational diffusion. These results can be exploited in various applications, including biomedical, polymer nanocomposite fabrication, drug delivery, and imaging
- …