A local corrections algorithm for solving Poisson’s equation in three dimensions

We present a second-order accurate algorithm for solving the free-space Poisson's equation on a locally-refined nested grid hierarchy in three dimensions. Our approach is based on linear superposition of local convolutions of localized charge distributions, with the nonlocal coupling represented on coarser grids. There presentation of the nonlocal coupling on the local solutions is based on Anderson's Method of Local Corrections and does not require iteration between different resolutions. A distributed-memory parallel implementation of this method is observed to have a computational cost per grid point less than three times that of a standard FFT-based method on a uniform grid of the same resolution, and scales well up to 1024 processors.

[1]  O. D. Kellogg.  Foundations of potential theory , 1934 .

[2]  P. P. Starling.  The numerical solution of Laplace's equation , 1963 .

[3]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[4]  K. Lackner,et al.  Computation of ideal MHD equilibria , 1976 .

[5]  R. A. James,et al.  The solution of Poisson''s equation for isolated source distributions , 1977 .

[6]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[7]  C. R. Anderson.  A method of local corrections for computing the velocity field due to a distribution of vortex blobs , 1986 .

[8]  Achi Brandt,et al.  Local mesh refinement multilevel techniques , 1987 .

[9]  Scott B. Baden,et al.  Programming Abstractions for Run-Time Partitioning of Scientific Continuum Calculations Running on Multiprocessors , 1987, PPSC.

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

[11]  P. Colella,et al.  A Fast Adaptive Vortex Method in Three Dimensions , 1994 .

[12]  L. Greengard,et al.  A Direct Adaptive Poisson Solver of Arbitrary Order Accuracy , 1996 .

[13]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[14]  Steven G. Johnson,et al.  FFTW: an adaptive software architecture for the FFT , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[15]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[16]  P. Colella,et al.  A Finite Difference Domain Decomposition Method Using Local Corrections for the Solution of Poisson's Equation , 1999 .

[17]  Leslie Greengard,et al.  A New Fast-Multipole Accelerated Poisson Solver in Two Dimensions , 2001, SIAM J. Sci. Comput..

[18]  P. Colella,et al.  A node-centered local refinement algorithm for Poisson's equation in complex geometries , 2004 .

[19]  Jingfang Huang,et al.  An adaptive fast solver for the modified Helmholtz equation in two dimensions , 2006 .