Tensor calculus in polar coordinates using Jacobi polynomials

Abstract

Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk. By construction, the bases satisfy regularity conditions at r=0 for any tensorial field. The coordinate singularity in a disk is a prototypical case for many coordinate singularities. The work presented here extends to other geometries. The operators represent covariant derivatives, multiplication by azimuthally symmetric functions, and the tensorial relationship between fields. These arise naturally from relations between classical orthogonal polynomials, and form a Heisenberg algebra. Other past work uses more specific polynomial bases for solving equations in polar coordinates. The main innovation in this paper is to use a larger set of possible bases to achieve maximum bandedness of linear operations. We provide a series of applications of the methods, illustrating their ease-of-use and accuracy.

48 pages, 8 figures. Accepted for publication in the Journal of Computational Physics

References in corpus (1)

• Orthogonal polynomials of several variables

Cited by in corpus (22)

• Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods
• Computation of Power Law Equilibrium Measures on Balls of Arbitrary Dimension
• On symmetrizing the ultraspherical spectral method for self-adjoint problems
• Tensor calculus in spherical coordinates using Jacobi polynomials. Part-I: Mathematical analysis and derivations
• A gyroscopic polynomial basis in the sphere
• Recurrence relations for orthogonal polynomials on a triangle
• A static memory sparse spectral method for time-fractional PDEs
• Gyroscopic polynomials
• Hemispheric handedness in the Galactic synchrotron polarization foreground
• Orthogonal signed-distance coordinates and vector calculus near evolving curves and surfaces
• Linear stability analysis of wake vortices by a spectral method using mapped Legendre functions
• A sparse spectral method for Volterra integral equations using orthogonal polynomials on the triangle
• Rapidly Rotating Wall-Mode Convection
• Orthogonal polynomial projection error measured in Sobolev norms in the unit disk
• A self-consistent numerical model of internal wave-induced mean flow oscillations in polar geometry
• Non-homogeneous wave equation on a cone
• Sparse spectral and p-finite element methods for partial differential equations on disk slices and trapeziums
• Tensor calculus in spherical coordinates using Jacobi polynomials, Part-II: Implementation and Examples
• A fully diagonalized spectral method on the unit ball
• Rational function approximation with normalized positive denominators
• A fast and spectrally convergent algorithm for rational-order fractional integral and differential equations
• Computing the Dirichlet-Neumann Operator on a Cylinder