Fast Finite Difference Solutions of the Three Dimensional Poisson’s Equation in Cylindrical Coordinates

Alemayehu Shiferaw, R. C. Mittal

Department of Mathematics, Indian Institute of Technology, Roorkee, India.

DOI: 10.4236/ajcm.2013.34045                

Abstract

In this work, the three-dimensional Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly, by extending the method of Hockney. The Poisson equation is approximated by second-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results.

Keywords

Poisson’s Equation; Hockney’s Method; Thomas Algorithm

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	E. H Chao, S. F. Paul, R. C. Davidson and K. S. Fine, “Direct Numerical Solution of Poisson’s Equation in Cylindrical (r, z) Coordinates,” Princeton University, Princeton, 1997.
[2]	H. Chen, Y. Su and B. D. Shizgal, “A Direct Spectral Collocation Poisson Solver in Polar and Cylindrical Coordinates,” Journal of Computational Physics, Vol. 160, No. 2, 2000, pp. 453-469. http://dx.doi.org/10.1006/jcph.2000.6461
[3]	I. Christopher, G. Knorr, M. Shoucri and P. Bertrand, “Solution of the Poisson Equation in an Annulus,” Journal of Computational Physics, Vol. 131, No. 2, 1997, pp. 323-326. http://dx.doi.org/10.1006/jcph.1996.5598
[4]	J. C. Kalita and K. K. Ray, “A Transformation-Free HOC Scheme for Incompressible Viscous Flows Past an Impulsively Started Circular Cylinder,” Journal of Computational Physics, Vol. 228, No. 14, 2009, pp. 5207-5236. http://dx.doi.org/10.1016/j.jcp.2009.04.016
[5]	M. C. Lai and W. C. Wang, “Fast Direct Solvers for Poisson Equation on 2D Polar and Spherical Geometries,” Numerical Methods for Partial Differential Equations, Vol. 18, No. 1, 2002, pp. 56-68.
[6]	P. N. Swarztrauber and R. A. Sweet, “The Direct Solution of the Discrete Poisson Equation on a Disk,” SIAM Journal on Numerical Analysis, Vol. 10, No. 5, 1973, pp. 900-907.
[7]	R. C. Mittal and S. Gahlaut, “High Order Finite Difference Schemes to Solve Poisson’s Equation in Cylindrical Symmetry,” Communications in Applied Numerical Methods, Vol. 3, 1987, pp. 457-461.
[8]	R. C. Mittal and S. Gahlaut, “High-Order Finite Differences Schemes to Solve Poisson’s Equation in Polar Coordinates,” IMA Journal of Numerical Analysis, Vol. 11, No. 2, 1991, pp. 261-270. http://dx.doi.org/10.1093/imanum/11.2.261
[9]	C. S. Tan, “Accurate Solution of three Dimensional Poisson’s Equation in Cylindrical Coordinate by Expansion in Chebyshev Polynomials,” Journal of Computational Physics, Vol. 59, No. 1, 1985, pp. 81-95. http://dx.doi.org/10.1016/0021-9991(85)90108-1
[10]	S. R. K. Iyengar and R. Manohar, “High Order Difference Methods for Heat Equation in Polar Cylindrical Polar Cylindrical Coordinates,” Journal of Computational Physics, Vol. 77, No. 2, 1988, pp. 425-438. http://dx.doi.org/10.1016/0021-9991(88)90176-3
[11]	S. R. K. Iyengar and A. Goyal, “A Note on Multigrid for the Three-Dimensional Poisson Equation in Cylindrical Coordinates,” Journal of Computational and Applied Mathematics, Vol. 33, No. 2, 1990, pp. 163-169. http://dx.doi.org/10.1016/0377-0427(90)90366-8
[12]	M. C. Lai and J. M. Tseng, “A Formally Fourth-Order Accurate Compact Scheme for 3D Poisson Equation in Cylindrical and Spherical Coordinates,” Journal of Computational and Applied Mathematics, Vol. 201, No. 1, 2007, pp. 175-181. http://dx.doi.org/10.1016/j.cam.2006.02.011
[13]	J. Xu, P. N. Ostroumov and J. Nolen, “A Parallel 3D Poisson Solver for Space Charge Simulation in Cylindrical Coordinates,” Computer Physics Communications, Vol. 178, No. 4, 2008, pp. 290-300. http://dx.doi.org/10.1016/j.cpc.2007.09.008
[14]	M. A. Malcolm and J. Palmer, “A Fast Method for Solving a Class of Tri-Diagonal Linear Systems,” Communications of the ACM, Vol. 17, No. 1, 1974, pp. 14-17. http://dx.doi.org/10.1145/360767.360777
[15]	R. W. Hockney, “A Fast Direct Solution of Poisson Equation Using Fourier Analysis,” JACM, Vol. 12, No. 1, 1965, pp. 95-113. http://dx.doi.org/10.1145/321250.321259