×

A simple compact fourth-order Poisson solver on polar geometry. (English) Zbl 1016.65093

Summary: We present a simple and efficient compact fourth-order Poisson solver in polar coordinates. This solver relies on the truncated Fourier series expansion, where the differential equations of the Fourier coefficients are solved by the compact fourth-order finite difference scheme. By shifting a grid a half mesh away from the origin and incorporating the symmetry constraint of Fourier coefficients, we can easily handle coordinate singularities without pole conditions. The numerical evidence confirms fourth-order accuracy for the problem on an annulus and third-order accuracy for the problem on a disk. In addition, a simple and comparably accurate approximation for the derivatives of the solution is also presented.

MSC:

65N35 Spectral, collocation and related methods for boundary value problems involving PDEs
65N06 Finite difference methods for boundary value problems involving PDEs
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
65T50 Numerical methods for discrete and fast Fourier transforms

Software:

FISHPAK

References:

[1] J. Adams, P. Swarztrauber, and, R. Sweet, Fishpack—A Package of Fortran Subprograms for the Solution of Separable Elliptic Partial Differential Equationshttp://www.netlib.org/fishpack; J. Adams, P. Swarztrauber, and, R. Sweet, Fishpack—A Package of Fortran Subprograms for the Solution of Separable Elliptic Partial Differential Equationshttp://www.netlib.org/fishpack
[2] Borges, L.; Daripa, P., A fast parallel algorithm for the Poisson equation on a disk, J. Comput. Phys., 169, 151 (2001) · Zbl 0980.65133
[3] Collatz, L., The Numerical Treatment of Differential Equation (1960), Springer-Verlag: Springer-Verlag New York · Zbl 0086.32601
[4] Gupta, M. M., A fourth-order Poisson solver, J. Comput. Phys., 55, 166 (1984) · Zbl 0548.65075
[5] Iyengar, S. R.K.; Manohar, R., High order difference methods for heat equation in polar cylindrical coordinates, J. Comput. Phys., 77, 425 (1988) · Zbl 0654.65069
[6] Jain, M. K.; Jain, R. K.; Krishna, M., A fourth-order difference scheme for quasilinear Poisson equation in polar co-ordinates, Commun. Numer. Methods Eng., 10, 791 (1994) · Zbl 0820.65064
[7] Keller, H. B., Numerical Solution of Two Point Boundary Value Problems (1976), Soc. for Indr. & Appl. Math: Soc. for Indr. & Appl. Math Philadelphia
[8] Lai, M.-C.; Wang, W.-C., Fast direct solvers for Poisson equation on 2D polar and spherical geometries, Numer. Methods Partial Differential Eq., 18, 56 (2002) · Zbl 0993.65114
[9] Mohseni, K.; Colonius, T., Numerical treatment of polar coordinate singularities, J. Comput. Phys., 157, 787 (2000) · Zbl 0981.76075
[10] Mittal, R. C.; Gahlaut, S., High-order finite-difference schemes to solve Poisson’s equation in polar coordinates, IMA J. Numer. Anal., 11, 261 (1991) · Zbl 0732.65097
[11] Strikwerda, J. C., Finite Difference Schemes and Partial Differential Equations (1989), Wadsworth & Brooks/Cole: Wadsworth & Brooks/Cole Belmont · Zbl 0681.65064
[12] Zhuang, Y.; Sun, X.-H., A high-order fast direct solver for singular Poisson equations, J. Comput. Phys., 171, 79 (2001) · Zbl 0990.65120
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.