Kernel Aggregated Fast Multipole Method: Efficient summation of Laplace and Stokes kernel functions

Abstract

Many different simulation methods for Stokes flow problems involve a common computationally intense task -- the summation of a kernel function over $O(N^2)$ pairs of points. One popular technique is the Kernel Independent Fast Multipole Method (KIFMM), which constructs a spatial adaptive octree for all points and places a small number of equivalent multipole and local equivalent points around each octree box, and completes the kernel sum with $O(N)$ cost, using these equivalent points. Simpler kernels can be used between these equivalent points to improve the efficiency of KIFMM. Here we present further extensions and applications to this idea, to enable efficient summations and flexible boundary conditions for various kernels. We call our method the Kernel Aggregated Fast Multipole Method (KAFMM), because it uses different kernel functions at different stages of octree traversal. We have implemented our method as an open-source software library STKFMM based on the high performance library PVFMM, with support for Laplace kernels, the Stokeslet, regularized Stokeslet, Rotne-Prager-Yamakawa (RPY) tensor, and the Stokes double-layer and traction operators. Open and periodic boundary conditions are supported for all kernels, and the no-slip wall boundary condition is supported for the Stokeslet and RPY tensor. The package is designed to be ready-to-use as well as being readily extensible to additional kernels.

References in corpus (12)

• Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors
• Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions
• Fast and spectrally accurate Ewald summation for 2-periodic electrostatic systems
• Regularized Stokeslet segments
• A scalable computational platform for particulate Stokes suspensions
• Fast Ewald summation for free-space Stokes potentials
• Flexibly imposing periodicity in kernel independent FMM: A Multipole-To-Local operator approach
• Complex dynamics of long, flexible fibers in shear
• Universal image systems for non-periodic and periodic Stokes flows above a no-slip wall
• The Spectral Ewald method for singly periodic domains
• RPYFMM: Parallel Adaptive Fast Multipole Method for Rotne-Prager-Yamakawa Tensor in Biomolecular Hydrodynamics Simulations
• Fast Ewald summation for Green's functions of Stokes flow in a half-space

Cited by in corpus (6)

• A multiscale biophysical model gives quantized metachronal waves in a lattice of cilia
• Hydrochemical interactions of phoretic particles: a regularized multipole framework
• Computing hydrodynamic interactions in confined doubly-periodic geometries in linear time
• Fast Ewald summation for Stokes flow with arbitrary periodicity
• Efficient Convergent Boundary Integral Methods for Slender Bodies
• Rheology of moderated dilute suspensions of star colloids: the shape factor