Validated forward integration scheme for parabolic PDEs via Chebyshev series

Abstract

In this paper we introduce a new approach to compute rigorously solutions of Cauchy problems for a class of semi-linear parabolic partial differential equations. Expanding solutions with Chebyshev series in time and Fourier series in space, we introduce a zero finding problem $F(a)=0$ on a Banach algebra $X$ of Fourier-Chebyshev sequences, whose solution solves the Cauchy problem. The challenge lies in the fact that the linear part $\mathcal{L}=DF(0)$ has an infinite block diagonal structure with blocks becoming less and less diagonal dominant at infinity. We introduce analytic estimates to show that $\mathcal{L}$ is an invertible linear operator on $X$, and we obtain explicit, rigorous and computable bounds for the operator norm $\Vert {\mathcal{L}}^{-1}{\Vert }_{B(X)}$. These bounds are then used to verify the hypotheses of a Newton-Kantorovich type argument which shows that the (Newton-like) operator $\mathcal{T}(a)=a-{\mathcal{L}}^{-1}F(a)$ is a contraction on a small ball centered at a numerical approximation of the Cauchy problem. The contraction mapping theorem yields a fixed point which corresponds to a classical (strong) solution of the Cauchy problem. The approach is simple to implement, numerically stable and is applicable to a class of PDE models, which include for instance Fisher's equation and the Swift-Hohenberg equation. We apply our approach to each of these models.

References in corpus (13)

• Modeling Elasticity in Crystal Growth
• Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: an overview
• Hopf's last hope: spatiotemporal chaos in terms of unstable recurrent patterns
• Unstable recurrent patterns in Kuramoto-Sivashinsky dynamics
• Symbolic dynamics for Kuramoto-Sivashinsky PDE on the line --- a computer-assisted proof
• Computer-assisted proof of heteroclinic connections in the one-dimensional Ohta-Kawasaki model
• Existence of globally attracting solutions for one-dimensional viscous Burgers equation with nonautonomous forcing - a computer assisted proof
• A new formulation for the numerical proof of the existence of solutions to elliptic problems
• Singularities and heteroclinic connections in complex-valued evolutionary equations with a quadratic nonlinearity
• Constructive a priori error estimates for a full discrete approximation of periodic solutions for the heat equation
• Rigorous numerics for nonlinear heat equations in the complex plane of time
• Accurate method of verified computing for solutions of semilinear heat equations
• Spontaneous periodic orbits in the Navier-Stokes flow

Cited by in corpus (2)

• Stationary non-radial localized patterns in the planar Swift-Hohenberg PDE: constructive proofs of existence
• A rigorous integrator and global existence for higher-dimensional semilinear parabolic PDEs via semigroup theory