Strip Decomposition Parallelization of Fast Direct Poisson Solver on a 3D Cartesian Staggered Grid

A strip domain decomposition parallel algorithm for fast direct Poisson solver is presented on a 3D Cartesian staggered grid. The parallel algorithm follows the principles of sequential algorithm for fast direct Poisson solver. Both Dirichlet and Neumann boundary conditions are addressed. Several test cases are likewise addressed in order to shed light on accuracy and efficiency in the strip domain parallelization algorithm. Actually the current implementation shows a very high efficiency when dealing with a large grid mesh up to 9 3.6 10 × under massive parallel approach, which explicitly demonstrates that the proposed algorithm is ready for massive parallel computing. Keywords—Strip-decomposition, parallelization, fast direct poisson solver.

[1]  A. E. Fuhs,et al.  The Fourier and Cyclic Reduction Methods for Solving Poisson ’ s Equation , 1996 .

[2]  Gunilla Sköllermo A Fourier method for the numerical solution of Poisson's equation , 1975 .

[3]  Tsutomu Hoshino,et al.  Parallel Poisson Solver FAGECR-Implementation and Performance Evaluation on PAX Computer , 1989 .

[4]  A Multidomain Spectral Method for Scalar and Vectorial Poisson Equations with Noncompact Sources , 2000, gr-qc/0003072.

[5]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[6]  Elias Balaras,et al.  An embedded-boundary formulation for large-eddy simulation of turbulent flows interacting with moving boundaries , 2006, J. Comput. Phys..

[7]  Roger W. Hockney,et al.  A Fast Direct Solution of Poisson's Equation Using Fourier Analysis , 1965, JACM.

[8]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[9]  Prabir Daripa,et al.  A fast parallel algorithm for the Poisson equation on a disk , 2001 .

[10]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[11]  June-Yub Lee,et al.  A parallel Poisson solver using the fast multipole method on networks of workstations , 1998 .

[12]  Paul N. Swarztrauber,et al.  The Vector Multiprocessor , 2000, Int. J. High Speed Comput..

[13]  Rainald Löhner,et al.  Improving the speed and accuracy of projection-type incompressible flow solvers , 2006 .

[14]  S. Ghanemi,et al.  A Domain Decomposition Method for Helmholtz Scattering Problems , 1997 .

[15]  Paul N. Swarztrauber,et al.  Multiprocessor FFTs , 1987, Parallel Comput..

[16]  William L. Briggs,et al.  Multiprocessor FFT Methods , 1985, PPSC.

[17]  Tony F. Chan,et al.  A domain-decomposed fast Poisson solver on a rectangle , 1985, PPSC.

[18]  Amir Averbuch,et al.  A Fast Spectral Solver for a 3D Helmholtz Equation , 1999, SIAM J. Sci. Comput..

[19]  Gene H. Golub,et al.  A Fast Poisson Solver for the Finite Difference Solution of the Incompressible Navier-Stokes Equations , 1998, SIAM J. Sci. Comput..

[20]  P. Swarztrauber,et al.  SPHEREPACK 3.0: A Model Development Facility , 1999 .

[21]  Clive Temperton On the FACR( l) algorithm for the discrete Poisson equation , 1980 .

[22]  Amir Averbuch,et al.  Parallelizing implicit algorithms for time-dependent problems by parabolic domain decomposition , 1993 .

[23]  Ulrich Schumann,et al.  Parallel solution of tridiagonal systems for the Poisson equation , 1995 .

[24]  L. Greengard,et al.  A Fast Poisson Solver for Complex Geometries , 1995 .

[25]  Paul N. Swarztrauber,et al.  Vector and parallel methods for the direct solution of Poisson's equation , 1989 .

[26]  E. Balaras Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations , 2004 .

[27]  R. Sweet,et al.  Fast Fourier transforms for direct solution of Poisson's equation with staggered boundary conditions , 1988 .

[28]  Peter S. Pacheco Parallel programming with MPI , 1996 .