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

DOI10.1007/S10444-021-09896-1zbMATH Open1492.65348arXiv2010.15155OpenAlexW3197771639MaRDI QIDQ2230587FDOQ2230587

Authors: Wen Yan, Robert Blackwell 

Publication date: 24 September 2021

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

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.

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

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

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)

Cites Work

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

Cited In (5)

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

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)