Python-JAX-based Fast Stokesian Dynamics

Abstract

Stokesian Dynamics (SD) is a powerful computational framework for simulating the motion of particles in a viscous Newtonian fluid under Stokes-flow conditions. Traditional SD implementations can be computationally expensive as they rely on the inversion of large mobility matrices to determine hydrodynamic interactions. Recently, however, the simulation of thermalized systems with large numbers of particles has become feasible [Fiore and Swan, J. Fluid. Mech. $\text{878}$, 544 (2019)]. Their ``fast Stokesian dynamics'' (FSD) method leverages a saddle-point formulation to ensure overall scaling of the algorithm that is linear in the number of particles $\mathcal{O}(N)$; performance relies on dedicated graphics-processing-unit computing. Here, we present a different route toward implementing FSD, which instead leverages the Just-in-Time (JIT) compilation capabilities of Google JAX. We refer to this implementation as JFSD and perform benchmarks on it to verify that it has the right scaling and is sufficiently fast by the standards of modern computational physics. In addition, we provide a series of physical test cases that help ensure accuracy and robustness, as the code undergoes further development. Thus, JFSD is ready to facilitate the study of hydrodynamic effects in particle suspensions across the domains of soft, active, and granular matter.

References in corpus (11)

• Lattice Boltzmann simulations of soft matter systems
• Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids
• Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations
• Efficient mesoscale hydrodynamics: multiparticle collision dynamics with massively parallel GPU acceleration
• Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions
• Probing colloidal gels at multiple lengthscales: the role of hydrodynamics
• Hydrodynamics strongly affect the dynamics of colloidal gelation but not gel structure
• Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations with Stresslet Constraints
• Active Stokesian Dynamics
• Modelling a Hydrodynamic Instability in Freely Settling Colloidal Gels
• Hydrodynamic Stability Criterion for Colloidal Gelation under Gravity