DeepAI AI Chat
Log In Sign Up

A scalable computational platform for particulate Stokes suspensions

by   Wen Yan, et al.
Simons Foundation

We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle collisions. This algorithm extends the well-known complementarity method for non-smooth multi-body dynamics to resolve collisions in dense rigid body suspensions. This approach formulates the collision resolution problem as a linear complementarity problem with geometric `non-overlapping' constraints imposed at each timestep. It is then reformulated as a constrained quadratic programming problem and the Barzilai-Borwein projected gradient descent method is applied for its solution. This framework is designed to be applicable for any convex particle shape, e.g., spheres and spherocylinders, and applicable to any Stokes mobility solver, including the Rotne-Prager-Yamakawa approximation, Stokesian Dynamics, and PDE solvers (e.g., boundary integral and immersed boundary methods). In particular, this method imposes Newton's Third Law and records the entire contact network. Further, we describe a fast, parallel, and spectrally-accurate boundary integral method tailored for spherical particles, capable of resolving lubrication effects. We show weak and strong parallel scalings up to 8× 10^4 particles with approximately 4× 10^7 degrees of freedom on 1792 cores. We demonstrate the versatility of this framework with several examples, including sedimentation of particle clusters, and active matter systems composed of ensembles of particles driven to rotate.


page 4

page 18

page 20

page 21

page 24

page 26

page 29

page 30


Fast and accurate solvers for simulating Janus particle suspensions in Stokes flow

We present a novel computational framework for simulating suspensions of...

A Barrier Method for Contact Avoiding Particles in Stokes Flow

Rigid particles in a Stokesian fluid can physically not overlap, as a th...

Boundary layer in linear viscoelasticity

It is well known that a boundary layer develops along an infinite plate ...

Scalable Simulation of Realistic Volume Fraction Red Blood Cell Flows through Vascular Networks

High-resolution blood flow simulations have potential for developing bet...

An integral model based on slender body theory, with applications to curved rigid fibers

We propose a novel integral model describing the motion of curved slende...

A parallel-GPU code for asteroid aggregation problems with angular particles

The paper presents a numerical implementation of the gravitational N-bod...