1,532 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 Lifting Relation from Macroscopic Variables to Mesoscopic Variables in Lattice Boltzmann Method: Derivation, Numerical Assessments and Coupling Computations Validation
In this paper, analytic relations between the macroscopic variables and the
mesoscopic variables are derived for lattice Boltzmann methods (LBM). The
analytic relations are achieved by two different methods for the exchange from
velocity fields of finite-type methods to the single particle distribution
functions of LBM. The numerical errors of reconstructing the single particle
distribution functions and the non-equilibrium distribution function by
macroscopic fields are investigated. Results show that their accuracy is better
than the existing ones. The proposed reconstruction operator has been used to
implement the coupling computations of LBM and macro-numerical methods of FVM.
The lid-driven cavity flow is chosen to carry out the coupling computations
based on the numerical strategies of domain decomposition methods (DDM). The
numerical results show that the proposed lifting relations are accurate and
robust
Finite Volume Streaming-based Lattice Boltzmann algorithm for fluid-dynamics simulations: a one-to-one accuracy and performance study
A new finite volume (FV) discretisation method for the Lattice Boltzmann (LB)
equation which combines high accuracy with limited computational cost is
presented. In order to assess the performance of the FV method we carry out a
systematic comparison, focused on accuracy and computational performances, with
the standard (ST) Lattice Boltzmann equation algorithm. To our
knowledge such a systematic comparison has never been previously reported. In
particular we aim at clarifying whether and in which conditions the proposed
algorithm, and more generally any FV algorithm, can be taken as the method of
choice in fluid-dynamics LB simulations. For this reason the comparative
analysis is further extended to the case of realistic flows, in particular
thermally driven flows in turbulent conditions. We report the first successful
simulation of high-Rayleigh number convective flow performed by a Lattice
Boltzmann FV based algorithm with wall grid refinement.Comment: 15 pages, 14 figures (discussion changes, improved figure
readability
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
Applications of Asymptotic Analysis
This workshop focused on asymptotic analysis and its fundamental role in the derivation and understanding of the nonlinear structure of mathematical models in various fields of applications, its impact on the development of new numerical methods and on other fields of applied mathematics such as shape optimization. This was complemented by a review as well as the presentation of some of the latest developments of singular perturbation methods
- …