ID: 570
MS01 Recent Advances in Fast Algorithms and Integral Equation Methods

Fast and accurate simulation of close-to-touching discs in 2D Stokes flow

Dhairya Malhotra1, Mariana Martinez Aguilar2, Dan Fortunato1

1 Flatiron Institute, United States of America; 2 EPFL, Switzerland; dmalhotra@flatironinstitute.org

We introduce a high-order boundary integral method for simulating dense
suspensions of rigid discs in a Stokesian fluid in 2D. Our method efficiently
handles close-to-touching interactions down to distances of 1e-10 with only a
coarse discretization of the boundary. Additionally, we present a
preconditioner that significantly reduces the number of GMRES iterations
required for solving the Stokes mobility problem at each time step. Coupled
with high-order, adaptive time-stepping using spectral deferred correction, we
are able to take larger time steps, mitigating the temporal stiffness resulting
from close-to-touching interactions.