Regular Article
A Spectral Method for Polar Coordinates

https://doi.org/10.1006/jcph.1995.1171Get rights and content

Abstract

A new set of polynomial functions that can be used in spectral expansions of C functions in polar coordinates (r, φ) is defined by a singular Sturm-Liouville equation. With the use of the basis functions, the spectral representations remain analytic at the pole despite the coordinate singularity because the pole condition is exactly satisfied at the origin for ag azimuthal modes, not just a few of the gravest modes (which is the usual case). Based on recurrence relations, fast and stable numerical operators for 1/r, r(d/dr), the Laplacian and Helmholtz operators and their inverses are developed. Although the spacings in the azimuthal direction of the collocation points near the origin are small (i.e., α 1/M2, where M is the number of radial modes), the explicit numerical method for Euler's equation is not stiff at the origin. Namely, the CFL number σ is O(1) where the grid size in is defined as π/M (i.e., the maximum allowable timestep is proportional to 1/M, not 1/M2).

References (0)

Cited by (93)

  • Gyroscopic polynomials

    2023, Journal of Computational Physics
  • A gyroscopic polynomial basis in the sphere

    2022, Journal of Computational Physics
  • Tensor calculus in spherical coordinates using Jacobi polynomials. Part-I: Mathematical analysis and derivations

    2019, Journal of Computational Physics: X
    Citation Excerpt :

    We also show how our basis sets map to other past works. Some authors use rescaled Jacobi polynomials in one form or another, but do not allow the parameters to change as derivatives act [3,29,32,51,64]. In §6 we provide conclusions for this Part-I of our two-part series.

View all citing articles on Scopus
View full text