Spectrally accurate fast summation for periodic Stokes potentials

A spectrally accurate method for the fast evaluation of N-particle sums of the periodic Stokeslet is presented. Two different decomposition methods, leading to one sum in real space and one in reciprocal space, are considered. An FFT based method is applied to the reciprocal part of the sum, invoking the equivalence of multiplications in reciprocal space to convolutions in real space, thus using convolutions with a Gaussian function to place the point sources on a grid. Due to the spectral accuracy of the method, the grid size needed is low and also in practice, for a fixed domain size, independent of N. The leading cost, which is linear in N, arises from the to-grid and from-grid operations. Combining this FFT based method for the reciprocal sum with the direct evaluation of the real space sum, a spectrally accurate algorithm with a total complexity of O(NlogN) is obtained. This has been shown numerically as the system is scaled up at constant density.

[1]  John Strain,et al.  Locally-corrected spectral methods and overdetermined elliptic systems , 2007, J. Comput. Phys..

[2]  Lexing Ying,et al.  A high-order 3D boundary integral equation solver for elliptic PDEs in smooth domains , 2006, J. Comput. Phys..

[3]  Leslie Greengard,et al.  Accelerating the Nonuniform Fast Fourier Transform , 2004, SIAM Rev..

[4]  P. P. Ewald Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[5]  Jingfang Huang,et al.  A parallel fast multipole accelerated integral equation scheme for 3D Stokes equations , 2007 .

[6]  Leslie Greengard,et al.  Integral equation methods for Stokes flow in doubly-periodic domains , 2004 .

[7]  C. Pozrikidis Boundary Integral and Singularity Methods for Linearized Viscous Flow: Index , 1992 .

[8]  Eric F Darve,et al.  A smooth particle-mesh Ewald algorithm for Stokes suspension simulations: The sedimentation of fibers , 2005 .

[9]  R. Dror,et al.  Gaussian split Ewald: A fast Ewald mesh method for molecular simulation. , 2005, The Journal of chemical physics.

[10]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[11]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[12]  C. Beenakker Ewald sum of the Rotne-Prager tensor , 1986 .

[13]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[14]  C. Pozrikidis,et al.  Computation of periodic Green's functions of Stokes flow , 1996 .

[15]  J. Perram,et al.  Cutoff Errors in the Ewald Summation Formulae for Point Charge Systems , 1992 .

[16]  John Strain,et al.  A geometric nonuniform fast Fourier transform , 2009, J. Comput. Phys..

[17]  Christian Holm,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998 .

[18]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[19]  Vladimir Rokhlin,et al.  Fast Fourier Transforms for Nonequispaced Data , 1993, SIAM J. Sci. Comput..

[20]  J. Thomas Beale,et al.  Locally corrected semi-Lagrangian methods for Stokes flow with moving elastic interfaces , 2008, J. Comput. Phys..

[21]  John F. Brady,et al.  Accelerated Stokesian Dynamics simulations , 2001, Journal of Fluid Mechanics.

[22]  Leslie Greengard,et al.  A fast multipole method for the three-dimensional Stokes equations , 2008, J. Comput. Phys..

[23]  H. Hasimoto On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres , 1959, Journal of Fluid Mechanics.

[24]  J. Strain Fast potential theory. II: Layer potentials and discrete sums , 1992 .