Hybrid Chebyshev function bases for sparse spectral methods in parity-mixed PDEs on an infinite domain

Abstract

We present a numerical spectral method to solve systems of differential equations on an infinite interval $y\in (-\infty ,\infty )$ in presence of linear differential operators of the form $Q(y){\left(\partial /{\partial }_y\right)}^b$ (where $Q(y)$ is a rational fraction and $b$ a positive integer). Even when these operators are not parity-preserving, we demonstrate how a mixed expansion in interleaved Chebyshev rational functions $T{B}_n(y)$ and $S{B}_n(y)$ preserves the sparsity of their discretization. This paves the way for fast $O(N\ln N)$ and spectrally accurate mixed implicit-explicit time-marching of sets of linear and nonlinear equations in unbounded geometries.

26 pages, 11 figures, submitted for publication

Cited by in corpus (2)

• Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods
• Equatorially Trapped Convection in a Rapidly Rotating Shallow Shell