Fast and High Accuracy Multigrid Solution of the
Three Dimensional Poisson Equations

Jun Zhang
Department of Computer Science and Engineering
University of Minnesota
Minneapolis, MN 55455, USA


We employ a fourth-order compact finite difference scheme (FOS) with the multigrid algorithm to solve the three dimensional Poisson equation. We test the influence of different orderings of the grid space and different grid-transfer operators on the convergence and efficiency of our high accuracy algorithm. Fourier smoothing analysis is conducted to show that FOS has a smaller smoothing factor than the traditional second-order central difference scheme (CDS). A new method of Fourier smoothing analysis is proposed for the partially decoupled red-black Gauss-Seidel relaxation with FOS. Numerical results are given to compare the computed accuracy and the computational efficiency of FOS with multigrid against CDS with multigrid.

Key words: Poisson equation, multigrid method, fourth-order compact scheme, Fourier smoothing analysis.

Download the compressed postscript file, or the PDF file 3dpoisson.pdf.gz.
This paper has been published in Journal of Computational Physics, Vol. 143, 449--461 (1998).