A numerical method for suspensions of articulated bodies in viscous flows
arXiv:2107.10800 · doi:10.1016/j.jcp.2022.111365
Abstract
An articulated body is defined as a finite number of rigid bodies connected by a set of arbitrary constraints that limit the relative motion between pairs of bodies. Such a general definition encompasses a wide variety of situations in the microscopic world, from bacteria to synthetic micro-swimmers, but it is also encountered when discretizing inextensible bodies, such as filaments or membranes. Simulating suspensions of such articulated bodies requires to solve the hydrodynamic interactions between large collections of objects of arbitrary shape while satisfying the multiple constraints that connect them. Two main challenges arise in this task: limiting the cost of the hydrodynamic solves, and enforcing the constraints within machine precision at each time-step. To address these challenges we propose a formalism that combines the body mobility problem in Stokes flow with a velocity formulation of the constraints, resulting in a mixed mobility-resistance problem. While resistance problems are known to scale poorly with the particle number, our preconditioned iterative solver is not sensitive to the system size. Constraint violations, e.g. due to discrete time-integration errors, are prevented by correcting the particles' positions and orientations at the end of each time-step. Our correction procedure, based on a nonlinear minimisation algorithm, is negligible in terms of computational cost and preserves the accuracy of the time-integration scheme. We showcase the robustness and scalability of our method by exploring the locomotion modes of a model microswimmer inspired by the diatom colony Bacillaria Paxillifer, and by simulating large suspensions of bacteria interacting near a no-slip boundary. Finally, we provide a Python implementation of our framework in a collaborative publicly available code.
arXiv comment35 pages, 11 figures; discussion expanded, typos fixed and minor updates on figures
References in corpus (27)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- The hydrodynamics of swimming microorganisms
- Swimming in circles: Motion of bacteria near solid boundaries
- Hydrodynamic attraction of swimming microorganisms by surfaces
- Hydrodynamics of self-propulsion near a boundary: predictions and accuracy of far-field approximations
- Fluid Flows Created by Swimming Bacteria Drive Self-Organization in Confined Suspensions
- Numerical Study of a Microscopic Artificial Swimmer
- Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors
- Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations
- Hydrodynamics of Suspensions of Passive and Active Rigid Particles: A Rigid Multiblob Approach
- Colloidomers: freely-jointed polymers made of droplets
- Brownian Dynamics of Confined Rigid Bodies
- A fast platform for simulating flexible fiber suspensions applied to cell mechanics
- Active turbulence in microswimmer suspensions -- the role of active hydrodynamic stress and volume exclusion
- A general formulation of Bead Models applied to flexible fibers and active filaments at low Reynolds number
- Relating microswimmer synthesis to rheotactic tunability
- Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions
- From flagellar undulations to collective motion: predicting the dynamics of sperm suspensions
- Methods for suspensions of passive and active filaments
- A regularised singularity approach to phoretic problems
- Wall entrapment of peritrichous bacteria: A mesoscale hydrodynamics simulation study
- Boundary behaviours of Leishmania mexicana: a hydrodynamic simulation study
- Buckling, Crumpling, and Tumbling of Semiflexible Sheets in Simple Shear Flow
- A fluctuating boundary integral method for Brownian suspensions
- Universal image systems for non-periodic and periodic Stokes flows above a no-slip wall
- A generalised drift-correcting time integration scheme for Brownian suspensions of rigid particles with arbitrary shape
- Stability and dynamics of anisotropically-tumbling chemotactic swimmers