Contact-aware simulations of particulate Stokesian suspensions

• Published in: Journal of computational physics Vol. 347; pp. 160 - 182
• Main Authors: Lu, Libin, Rahimian, Abtin, Zorin, Denis
• Format: Journal Article
• Language: English
• Published: Cambridge Elsevier Inc 15.10.2017; Elsevier Science Ltd
• Subjects: Boundary conditions; Boundary integral; Complex fluids; Computational physics; Computer simulation; Constraint-based collision handling; Deformation; Discretization; Fluid mechanics; Formability; High volume fraction flow; Intersections; Navier-Stokes equations; Numerical analysis; Particulate Stokes flow; Robustness (mathematics); SDC time stepping; Stability; Studies; High volume fraction flow; Boundary integral; Constraint-based collision handling; Complex fluids; Particulate Stokes flow; SDC time stepping
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2017.06.039

We present an efficient, accurate, and robust method for simulation of dense suspensions of deformable and rigid particles immersed in Stokesian fluid in two dimensions. We use a well-established boundary integral formulation for the problem as the foundation of our approach. This type of formulation, with a high-order spatial discretization and an implicit and adaptive time discretization, have been shown to be able to handle complex interactions between particles with high accuracy. Yet, for dense suspensions, very small time-steps or expensive implicit solves as well as a large number of discretization points are required to avoid non-physical contact and intersections between particles, leading to infinite forces and numerical instability. Our method maintains the accuracy of previous methods at a significantly lower cost for dense suspensions. The key idea is to ensure interference-free configuration by introducing explicit contact constraints into the system. While such constraints are unnecessary in the formulation, in the discrete form of the problem, they make it possible to eliminate catastrophic loss of accuracy by preventing contact explicitly. Introducing contact constraints results in a significant increase in stable time-step size for explicit time-stepping, and a reduction in the number of points adequate for stability.