Recent advances about the rigorous integration of parabolic PDEs via fully spectral Fourier-Chebyshev expansions
Matthieu Cadiot, Jean-Philippe Lessard
TL;DR
This work develops a rigorous, computer-assisted framework for solving semilinear parabolic PDEs using fully spectral Fourier-Chebyshev expansions. By recasting the Cauchy problem as a zero of a Fourier-Chebyshev coefficient map $F$, it combines an explicit decay estimate for the linear inverse with a Newton-Kantorovich scheme and a carefully constructed approximate inverse to obtain existence with explicit error bounds and the ability to use larger time steps. The approach is then demonstrated on the 2D Navier–Stokes equations (including a global existence result for a nontrivial initial condition) and the Swift–Hohenberg equation, with explicit computable bounds and numerical verifications provided; the techniques are supported by detailed tail estimates and a rigorous analysis of the infinite-dimensional operators. Overall, the paper provides a robust, scalable method for guaranteed existence and error control in high-dimensional dissipative PDEs via fully spectral methods, accompanied by publicly available code.
Abstract
This paper presents a novel approach to rigorously solving initial value problems for semilinear parabolic partial differential equations (PDEs) using fully spectral Fourier-Chebyshev expansions. By reformulating the PDE as a system of nonlinear ordinary differential equations and leveraging Chebyshev series in time, we reduce the problem to a zero-finding task for Fourier-Chebyshev coefficients. A key theoretical contribution is the derivation of an explicit decay estimate for the inverse of the linear part of the PDE, enabling larger time steps. This allows the construction of an approximate inverse for the Fréchet derivative and the application of a Newton-Kantorovich theorem to establish solution existence within explicit error bounds. Building on prior work, our method is extended to more complex partial differential equations, including the 2D Navier-Stokes equations, for which we establish global existence of the solution of the IVP for a given nontrivial initial condition.
