Fast Ewald summation for Stokesian particle suspensions

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

Ämnesord

TEKNIK OCH TEKNOLOGIER  -- Maskinteknik -- Strömningsmekanik och akustik (hsv//swe)

ENGINEERING AND TECHNOLOGY  -- Mechanical Engineering -- Fluid Mechanics and Acoustics (hsv//eng)

NATURVETENSKAP  -- Matematik -- Beräkningsmatematik (hsv//swe)

NATURAL SCIENCES  -- Mathematics -- Computational Mathematics (hsv//eng)

Nyckelord

viscous flows

integral equations

error estimation

microfluidics

multibody dynamics

spectral

double layer

boundary integral

ewald summation

Tillämpad matematik och beräkningsmatematik

Applied and Computational Mathematics

Publikations- och innehållstyp

ref (ämneskategori)

art (ämneskategori)