Kernel aggregated fast multipole method. Efficient summation of Laplace and Stokes kernel functions

DOI10.1007/S10444-021-09896-1MaRDI QIDQ2230587zbMATH OpenOpenAlexFDO

Authors Wen Yan, Robert Blackwell

Publication date 24 September 2021

Published in Advances in Computational Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2010.15155

zbMATH Keywords

periodic boundary conditionsfast multipole methodstokeslet

Mathematics Subject Classification ID

Numerical summation of series (65B10) Complexity and performance of numerical algorithms (65Y20) Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50) Stokes and related (Oseen, etc.) flows (76D07) Software, source code, etc. for problems pertaining to numerical analysis (65-04) Numerical methods for partial differential equations, boundary value problems (65N99)

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.

Recommendations

• Flexibly imposing periodicity in kernel independent FMM: a multipole-to-local operator approach
• A kernel-independent adaptive fast multipole algorithm in two and three dimensions
• A fast multipole method for the three-dimensional Stokes equations
• Fast Ewald summation for free-space Stokes potentials
• A Fourier-series-based kernel-independent fast multipole method

Cites work

• A fast algorithm for particle simulations
• A fast multipole method for the Rotne-Prager-Yamakawa tensor and its applications
• A kernel-independent adaptive fast multipole algorithm in two and three dimensions
• A kernel-independent treecode for general Rotne-Prager-Yamakawa tensor
• A scalable computational platform for particulate Stokes suspensions
• Fast Ewald summation for Green's functions of Stokes flow in a half-space
• Fast Ewald summation for free-space Stokes potentials
• Flexibly imposing periodicity in kernel independent FMM: a multipole-to-local operator approach
• Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors
• Modeling the dynamics of an elastic rod with intrinsic curvature and twist using a regularized Stokes formulation
• PVFMM: A parallel kernel independent FMM for particle and volume potentials
• RPYFMM: parallel adaptive fast multipole method for Rotne-Prager-Yamakawa tensor in biomolecular hydrodynamics simulations
• Regularized Stokeslet segments
• Rotne-Prager-Yamakawa approximation for different-sized particles in application to macromolecular bead models
• Spectral Ewald acceleration of Stokesian dynamics for polydisperse suspensions
• The Ewald sums for singly, doubly and triply periodic electrostatic systems
• The method of regularized Stokeslets in three dimensions:  Analysis, validation, and application to helical swimming
• The spectral Ewald method for singly periodic domains
• Universal image systems for non-periodic and periodic Stokes flows above a no-slip wall

Cited in

(5)

• Fast orthogonalization to the kernel of the discrete gradient operator with application to Stokes problem
• Fast Ewald summation for Stokes flow with arbitrary periodicity
• STKFMM
• Flexibly imposing periodicity in kernel independent FMM: a multipole-to-local operator approach
• A fast multipole method for the Rotne-Prager-Yamakawa tensor and its applications

Describes a project that uses

Uses Software

• pvfmm
• RPYFMM

This page was built for publication: Kernel aggregated fast multipole method. Efficient summation of Laplace and Stokes kernel functions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2230587)