The Accurate Solution of Poissons Equation by Expansion in Chebyshev Polynomials-学术期刊.pdf

JOURNAL OF COMPUTATIONAL PHYSICS 30, 167-180 (1979) The Accurate Solution of Poisson’s Equation by Expansion in Chebyshev Polynomials DALE B. HAIDVOGEL Center for Earth and Planetary Physics, Pierce Hall, Harvard University, Cambridge, Massachusetts 02138 AND THOMAS ZANG ICASE, NASA Langley Research Center, Hampton, Virginia 23665 Received September 14, 1977; revised January 30, 1978 A Chebyshev expansion technique is applied to Poisson’s equation on a square with homogeneous Dirichlet boundary conditions. The spectral equations are solved in two ways-by alternating direction and by matrix diagonalization methods. Solutions are sought to both oscillatory and mildly singular problems. The accuracy and efficiency of the Chebyshev approach compare favorably with those of standard second- and fourth-order finite-difference methods. 1. INTR~DUCTJ~N Many physical problems require a numerical solution of the two-dimensional Poisson’s equation. On certain domains, such as rectangles, the so-called fast Poisson solvers that are now widely available (cf. the reviews by Dorr [l] and Swarztrauber [2]) provide a rapid solution of the standard five-point difference approximation to the partial differential equation. However, the resolution of these methods