A Barrier Method for Contact Avoiding Particles in Stokes Flow

01/04/2023
by   Anna Broms, et al.
0

Rigid particles in a Stokesian fluid can physically not overlap, as a thin layer of fluid always separates a particle pair, exerting increasingly strong repulsive forces on the bodies for decreasing separations. Numerically, resolving these lubrication forces comes at an intractably large cost even for moderate system sizes. Hence, it can typically not be guaranteed that particle collisions and overlaps do not occur in a dynamic simulation, independently of the choice of method to solve the Stokes equations. In this work, non-overlap constraints, in terms of the Euclidean distance between boundary points on the particles, are represented via a barrier energy. We solve for the minimum magnitudes of repelling contact forces between any particle pair in contact to correct for overlaps by enforcing a zero barrier energy at the next time level, given a contact-free configuration at a previous instance in time. The method is tested using a multiblob method to solve the mobility problem in Stokes flow applied to suspensions of spheres, rods and boomerang shaped particles. Collision free configurations are obtained at all instances in time. The effect of the contact forces on the collective order of a set of rods in a background flow that naturally promote particle interactions is also illustrated.

READ FULL TEXT
research
02/09/2021

Dynamic Mode Decomposition of inertial particle caustics in Taylor-Green flow

Inertial particles advected by a background flow can show complex struct...
research
09/14/2019

A scalable computational platform for particulate Stokes suspensions

We describe a computational framework for simulating suspensions of rigi...
research
07/22/2022

A Locally Corrected Multiblob Method with Hydrodynamically Matched Grids for the Stokes Mobility Problem

Inexpensive numerical methods are key to enable simulations of systems o...
research
04/19/2022

ParticLS: Object-oriented software for discrete element methods and peridynamics

ParticLS (Particle Level Sets) is a software library that implements the...
research
02/23/2021

Using a deep neural network to predict the motion of under-resolved triangular rigid bodies in an incompressible flow

We consider non-spherical rigid body particles in an incompressible flui...
research
12/10/2019

An integral equation method for the simulation of doubly-periodic suspensions of rigid bodies in a shearing viscous flow

With rheology applications in mind, we present a fast solver for the tim...
research
09/23/2019

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

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

Please sign up or login with your details

Forgot password? Click here to reset