Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4967236 | Journal of Computational Physics | 2017 | 27 Pages |
Abstract
We present a numerical spectral method to solve systems of differential equations on an infinite interval yâ(ââ,â) in presence of linear differential operators of the form Q(y)(â/ây)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 TBn(y) and SBn(y) preserves the sparsity of their discretization. This paves the way for fast O(Nlnâ¡N) and spectrally accurate mixed implicit-explicit time-marching of sets of linear and nonlinear equations in unbounded geometries.
Keywords
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Benjamin Miquel, Keith Julien,