An integral equation method for the simulation of doubly-periodic suspensions of rigid bodies in a shearing viscous flow
From MaRDI portal
Publication:6330809
DOI10.1016/J.JCP.2020.109809arXiv1912.04501MaRDI QIDQ6330809FDOQ6330809
Authors: Jun Wang, Ehssan Nazockdast, Alex H. Barnett 
Publication date: 10 December 2019
Abstract: With rheology applications in mind, we present a fast solver for the time-dependent effective viscosity of an infinite lattice containing one or more neutrally buoyant smooth rigid particles per unit cell, in a two-dimensional Stokes fluid with given shear rate. At each time, the mobility problem is reformulated as a 2nd-kind boundary integral equation, then discretized to spectral accuracy by the Nystrom method and solved iteratively, giving typically 10 digits of accuracy. Its solution controls the evolution of particle locations and angles in a first-order system of ordinary differential equations. The formulation is placed on a rigorous footing by defining a generalized periodic Green's function for the skew lattice. Numerically, the periodized integral operator is split into a near image sum|applied in linear time via the fast multipole method|plus a correction field solved cheaply via proxy Stokeslets. We use barycentric quadratures to evaluate particle interactions and velocity fields accurately, even at distances much closer than the node spacing. Using first-order time-stepping we simulate, eg, 25 ellipses per unit cell to 3-digit accuracy on a desktop in 1 hour per shear time. Our examples show equilibration at long times, force chains, and two types of blow-ups (jamming) whose power laws match lubrication theory asymptotics.
This page was built for publication: An integral equation method for the simulation of doubly-periodic suspensions of rigid bodies in a shearing viscous flow
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6330809)
