A formally fourth-order accurate compact scheme for 3D Poisson equation in cylindrical coordinates

Ming-Chih Lai*, Jui Ming Tseng

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

11 Scopus citations

Abstract

In this paper, we extend our previous work (M.-C. Lai, A simple compact fourth-order Poisson solver on polar geometry, J. Comput. Phys. 182 (2002) 337-345) to 3D cases. More precisely, we present a spectral/finite difference scheme for Poisson equation in cylindrical coordinates. The scheme relies on the truncated Fourier series expansion, where the partial differential equations of Fourier coefficients are solved by a formally fourth-order accurate compact difference discretization. Here the formal fourth-order accuracy means that the scheme is exactly fourth-order accurate while the poles are excluded and is third-order accurate otherwise. Despite the degradation of one order of accuracy due to the presence of poles, the scheme handles the poles naturally; thus, no pole condition is needed. The resulting linear system is then solved by the Bi-CGSTAB method with the preconditioner arising from the second-order discretization which shows the scalability with the problem size.

Original languageEnglish
Pages (from-to)175-181
Number of pages7
JournalJournal of Computational and Applied Mathematics
Volume201
Issue number1
DOIs
StatePublished - 1 Apr 2007

Keywords

  • Bi-CGSTAB method
  • Cylindrical coordinates
  • Fast Fourier transform
  • Poisson equation
  • Symmetry constraint

Fingerprint Dive into the research topics of 'A formally fourth-order accurate compact scheme for 3D Poisson equation in cylindrical coordinates'. Together they form a unique fingerprint.

Cite this