Fast Ewald summation for Stokes flow with arbitrary periodicity
arXiv:2210.01255 · doi:10.1016/j.jcp.2023.112473
Abstract
A fast and spectrally accurate Ewald summation method for the evaluation of stokeslet, stresslet and rotlet potentials of three-dimensional Stokes flow is presented. This work extends the previously developed Spectral Ewald method for Stokes flow to periodic boundary conditions in any number (three, two, one, or none) of the spatial directions, in a unified framework. The periodic potential is split into a short-range and a long-range part, where the latter is treated in Fourier space using the fast Fourier transform. A crucial component of the method is the modified kernels used to treat singular integration. We derive new modified kernels, and new improved truncation error estimates for the stokeslet and stresslet. An automated procedure for selecting parameters based on a given error tolerance is designed and tested. Analytical formulas for validation in the doubly and singly periodic cases are presented. We show that the computational time of the method scales like O(N log N) for N sources and targets, and investigate how the time depends on the error tolerance and window function, i.e. the function used to smoothly spread irregular point data to a uniform grid. The method is fastest in the fully periodic case, while the run time in the free-space case is around three times as large. Furthermore, the highest efficiency is reached when applying the method to a uniform source distribution in a primary cell with low aspect ratio. The work presented in this paper enables efficient and accurate simulations of three-dimensional Stokes flow with arbitrary periodicity using e.g. boundary integral and potential methods.
arXiv comment54 pages, 15 figures
References in corpus (21)
- Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations
- Fast convolution with free-space Green's functions
- Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions
- Fast and spectrally accurate Ewald summation for 2-periodic electrostatic systems
- Large Scale Brownian Dynamics of Confined Suspensions of Rigid Particles
- Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations with Stresslet Constraints
- Fast Ewald summation for free-space Stokes potentials
- An integral-based spectral method for inextensible slender fibers in Stokes flow
- Flexibly imposing periodicity in kernel independent FMM: A Multipole-To-Local operator approach
- A fast spectral method for electrostatics in doubly-periodic slit channels
- A fluctuating boundary integral method for Brownian suspensions
- Universal image systems for non-periodic and periodic Stokes flows above a no-slip wall
- Computing hydrodynamic interactions in confined doubly-periodic geometries in linear time
- A Discrete Ion Stochastic Continuum Overdamped Solvent Algorithm for Modeling Electrolytes
- The Spectral Ewald method for singly periodic domains
- Kernel Aggregated Fast Multipole Method: Efficient summation of Laplace and Stokes kernel functions
- An integral equation method for closely interacting surfactant-covered droplets in wall-confined Stokes flow
- Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries
- Fast Ewald summation for Green's functions of Stokes flow in a half-space
- Accelerating the calculation of dipolar interactions in particle based simulations with open boundary conditions by means of the P2NFFT method
- Fast Ewald summation for electrostatic potentials with arbitrary periodicity
Cited by in corpus (4)
- Computing hydrodynamic interactions in confined doubly-periodic geometries in linear time
- Fast Fourier Transform periodic interpolation method for superposition sums in a periodic unit cell
- A Simplified Fast Multipole Method Based on Strong Recursive Skeletonization
- Ewald summing irreducible components of flow around active particles