Validated forward integration scheme for parabolic PDEs via Chebyshev series

01/03/2021
by   Jacek Cyranka, et al.
0

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 ℒ := 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 ℒ is a boundedly invertible linear operator on X, and we obtain explicit, rigorous and computable bounds for the operator norm ℒ^-1_B(X). These bounds are then used to verify the hypotheses of a Newton-Kantorovich type argument which shows that the (Newton-like) operator 𝒯(a) := a - ℒ^-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, the Kuramoto-Sivashinsky equation, the Swift-Hohenberg equation and the phase-field crystal (PFC) equation. We apply our approach to each of these models and report plausible experimental results, which motivate further research on the method.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/14/2023

Validated integration of semilinear parabolic PDEs

Integrating evolutionary partial differential equations (PDEs) is an ess...
research
10/02/2019

A new formulation for the numerical proof of the existence of solutions to elliptic problems

Infinite-dimensional Newton methods can be effectively used to derive nu...
research
01/24/2023

Koopman neural operator as a mesh-free solver of non-linear partial differential equations

The lacking of analytic solutions of diverse partial differential equati...
research
06/10/2021

Existence of Strong Solution for the Complexified Non-linear Poisson Boltzmann Equation

We prove the existence and uniqueness of the complexified Nonlinear Pois...
research
12/17/2020

Speed of convergence of Chernoff approximations for two model examples: heat equation and transport equation

Paul Chernoff in 1968 proposed his approach to approximations of one-par...
research
02/19/2020

On bounded mild solution for a class of semilinear stochastic partial differential equation driven by Levy and stable processes

We show the existence and uniqueness of bounded mild solutions for a cla...
research
06/16/2023

A Numerically Robust and Stable Time-Space Pseudospectral Approach for Generalized Burgers-Fisher Equation

In this article, we present the time-space Chebyshev pseudospectral meth...

Please sign up or login with your details

Forgot password? Click here to reset