Contact-aware simulations of particulate Stokesian suspensions

Libin Lu, Abtin Rahimian, Denis Zorin

Research output: Contribution to journalArticle

Abstract

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.

Original languageEnglish (US)
Pages (from-to)160-182
Number of pages23
JournalJournal of Computational Physics
Volume347
DOIs
StatePublished - Oct 15 2017

Fingerprint

particulates
Particle interactions
formulations
simulation
Fluids
intersections
Costs
interference
fluids
configurations
interactions

Keywords

  • Boundary integral
  • Complex fluids
  • Constraint-based collision handling
  • High volume fraction flow
  • Particulate Stokes flow
  • SDC time stepping

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

Contact-aware simulations of particulate Stokesian suspensions. / Lu, Libin; Rahimian, Abtin; Zorin, Denis.

In: Journal of Computational Physics, Vol. 347, 15.10.2017, p. 160-182.

Research output: Contribution to journalArticle

@article{f8802d6259744a8daece16a93151bb62,
title = "Contact-aware simulations of particulate Stokesian suspensions",
abstract = "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.",
keywords = "Boundary integral, Complex fluids, Constraint-based collision handling, High volume fraction flow, Particulate Stokes flow, SDC time stepping",
author = "Libin Lu and Abtin Rahimian and Denis Zorin",
year = "2017",
month = "10",
day = "15",
doi = "10.1016/j.jcp.2017.06.039",
language = "English (US)",
volume = "347",
pages = "160--182",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Contact-aware simulations of particulate Stokesian suspensions

AU - Lu, Libin

AU - Rahimian, Abtin

AU - Zorin, Denis

PY - 2017/10/15

Y1 - 2017/10/15

N2 - 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.

AB - 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.

KW - Boundary integral

KW - Complex fluids

KW - Constraint-based collision handling

KW - High volume fraction flow

KW - Particulate Stokes flow

KW - SDC time stepping

UR - http://www.scopus.com/inward/record.url?scp=85022227390&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85022227390&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2017.06.039

DO - 10.1016/j.jcp.2017.06.039

M3 - Article

AN - SCOPUS:85022227390

VL - 347

SP - 160

EP - 182

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -