A fast integral equation method for solid particles in viscous flow using quadrature by expansion
2016 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 326, 420-445 p.Article in journal (Refereed) Published
Boundary integral methods are advantageous when simulating viscous flow around rigid particles, due to the reduction in number of unknowns and straightforward handling of the geometry. In this work we present a fast and accurate framework for simulating spheroids in periodic Stokes flow, which is based on the completed double layer boundary integral formulation. The framework implements a new method known as quadrature by expansion (QBX), which uses surrogate local expansions of the layer potential to evaluate it to very high accuracy both on and off the particle surfaces. This quadrature method is accelerated through a newly developed precomputation scheme. The long range interactions are computed using the spectral Ewald (SE) fast summation method, which after integration with QBX allows the resulting system to be solved in M log M time, where M is the number of particles. This framework is suitable for simulations of large particle systems, and can be used for studying e.g. porous media models.
Place, publisher, year, edition, pages
Academic Press, 2016. Vol. 326, 420-445 p.
Viscous flow, Stokes equations, Boundary integral methods, Quadrature by expansion, Fast Ewald summation
IdentifiersURN: urn:nbn:se:kth:diva-196589DOI: 10.1016/j.jcp.2016.09.006ISI: 000386067400023ScopusID: 2-s2.0-84988423251OAI: oai:DiVA.org:kth-196589DiVA: diva2:1047630
FunderSwedish Research Council, 2011-3178Göran Gustafsson Foundation for Research in Natural Sciences and MedicineSwedish e‐Science Research Center
QC 201611182016-11-182016-11-172016-11-18Bibliographically approved