Mathematical foundations of spectral methods for time-dependent PDEs

Abstract: The contention of this paper is that a spectral method for time-dependent PDEs is basically no more than a choice of an orthonormal basis of the underlying Hilbert space. This choice is governed by a long list of considerations: stability, speed of convergence, geometric numerical integration, fast approximation and efficient linear algebra. We subject different choices of orthonormal bases, focussing on the real line, to these considerations. While nothing is likely to improve upon a Fourier basis in the presence of periodic boundary conditions, the situation is considerably more interesting in other settings. We introduce two kinds of orthonormal bases, T-systems and W-systems, and investigate in detail their features. T-systems are designed to work with Cauchy boundary conditions, while W-systems are suited to zero Dirichlet boundary conditions.