Fast Ewald summation for Stokesian particle suspensions
2014 (English)In: International Journal for Numerical Methods in Fluids, ISSN 0271-2091, E-ISSN 1097-0363, Vol. 76, no 10, 669-698 p.Article in journal (Refereed) Published
We present a numerical method for suspensions of spheroids of arbitrary aspect ratio, which sediment under gravity. The method is based on a periodized boundary integral formulation using the Stokes double layer potential. The resulting discrete system is solved iteratively using generalized minimal residual accelerated by the spectral Ewald method, which reduces the computational complexity to O(N log N), where N is the number of points used to discretize the particle surfaces. We develop predictive error estimates, which can be used to optimize the choice of parameters in the Ewald summation. Numerical tests show that the method is well conditioned and provides good accuracy when validated against reference solutions.
Place, publisher, year, edition, pages
John Wiley & Sons, 2014. Vol. 76, no 10, 669-698 p.
viscous flows, integral equations, error estimation, microfluidics, multibody dynamics, spectral, double layer, boundary integral, ewald summation
Fluid Mechanics and Acoustics Computational Mathematics
Research subject Applied and Computational Mathematics
IdentifiersURN: urn:nbn:se:kth:diva-116383DOI: 10.1002/fld.3953ISI: 000344349000004OAI: oai:DiVA.org:kth-116383DiVA: diva2:589279
FunderSwedish Research Council, 2011-3178
QC 201411192013-01-172013-01-172016-04-27Bibliographically approved