Fast Ewald summation for free-space Stokes potentials

DOI10.1186/S40687-016-0092-7zbMATH Open1360.65295arXiv1607.04808OpenAlexW3103605810WikidataQ59611925 ScholiaQ59611925MaRDI QIDQ516003FDOQ516003

Authors: Ludvig af Klinteberg, Davoud Saffar Shamshirgar, Anna-Karin Tornberg 

Publication date: 20 March 2017

Published in: Research in the Mathematical Sciences (Search for Journal in Brave)

Abstract: We present a spectrally accurate method for the rapid evaluation of free-space Stokes potentials, i.e. sums involving a large number of free space Green's functions. We consider sums involving stokeslets, stresslets and rotlets that appear in boundary integral methods and potential methods for solving Stokes equations. The method combines the framework of the Spectral Ewald method for periodic problems, with a very recent approach to solving the free-space harmonic and biharmonic equations using fast Fourier transforms (FFTs) on a uniform grid. Convolution with a truncated Gaussian function is used to place point sources on a grid. With precomputation of a scalar grid quantity that does not depend on these sources, the amount of oversampling of the grids with Gaussians can be kept at a factor of two, the minimum for aperiodic convolutions by FFTs. The resulting algorithm has a computational complexity of O(N log N) for problems with N sources and targets. Comparison is made with a fast multipole method (FMM) to show that the performance of the new method is competitive.

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

Recommendations

• Spectrally accurate fast summation for periodic Stokes potentials
• FFT based spectral Ewald methods as an alternative to fast multipole methods
• Fast Ewald summation for Green's functions of Stokes flow in a half-space
• A fast multipole method for the three-dimensional Stokes equations
• Fast potential theory. II: Layer potentials and discrete sums

zbMATH Keywords

computational complexityconvolutionalgorithmStokes equationsfast Fourier transformsfast multipole methodboundary integral methodsfree space Green's functionsfree-space Stokes potentialsharmonic and biharmonic equationspotential methodsspectral Ewald method

Mathematics Subject Classification ID

Complexity and performance of numerical algorithms (65Y20) Navier-Stokes equations (35Q30) Boundary element methods for boundary value problems involving PDEs (65N38) Numerical methods for discrete and fast Fourier transforms (65T50)

Cites Work

• A fast algorithm for particle simulations
• A kernel-independent adaptive fast multipole algorithm in two and three dimensions
• Fast multipole method for the biharmonic equation in three dimensions
• Spectral accuracy in fast Ewald-based methods for particle simulations
• Title not available (Why is that?)
• Accelerating the Nonuniform Fast Fourier Transform
• A fast adaptive multipole algorithm in three dimensions
• Computation of periodic Green's functions of Stokes flow
• The type 3 nonuniform FFT and its applications
• On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres
• A smooth particle-mesh Ewald algorithm for Stokes suspension simulations: The sedimentation of fibers
• Spectrally accurate fast summation for periodic Stokes potentials
• A fast multipole method for the three-dimensional Stokes equations
• Completed double layer boundary element method for periodic suspensions
• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• Fast convolution with free-space Green's functions
• A parallel fast multipole accelerated integral equation scheme for 3D Stokes equations
• Title not available (Why is that?)
• Spectral Ewald acceleration of Stokesian dynamics for polydisperse suspensions

Cited In (25)

• The spectral Ewald method for singly periodic domains
• A 3D boundary integral method for the electrohydrodynamics of surfactant-covered drops
• A kernel-independent treecode for general Rotne-Prager-Yamakawa tensor
• Kernel aggregated fast multipole method. Efficient summation of Laplace and Stokes kernel functions
• An integral equation method for the advection-diffusion equation on time-dependent domains in the plane
• A kernel-independent treecode based on barycentric Lagrange interpolation
• Fast Ewald summation for Green's functions of Stokes flow in a half-space
• Fast Ewald summation for Stokes flow with arbitrary periodicity
• A fast multipole method for the three-dimensional Stokes equations
• Fast Fourier transform periodic interpolation method for superposition sums in a periodic unit cell
• Fast and scalable evaluation of pairwise potentials
• FFT based spectral Ewald methods as an alternative to fast multipole methods
• The Effect of Global Smoothness on the Accuracy of Treecodes
• Fast convolution with free-space Green's functions
• An accurate integral equation method for Stokes flow with piecewise smooth boundaries
• A Treecode Algorithm for 3D Stokeslets and Stresslets
• NFFT based Ewald summation for electrostatic systems with charges and dipoles
• A fluctuating boundary integral method for Brownian suspensions
• Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme
• Spectrally accurate fast summation for periodic Stokes potentials
• A highly accurate boundary integral equation method for surfactant-laden drops in 3D
• Parallel accelerated Stokesian dynamics with Brownian motion
• Fast potential theory. II: Layer potentials and discrete sums
• A local target specific quadrature by expansion method for evaluation of layer potentials in 3D
• Partition of unity extension of functions on complex domains

Uses Software

• Spectral Ewald
• GitHub

This page was built for publication: Fast Ewald summation for free-space Stokes potentials

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