A formal proof of the Lax equivalence theorem for finite difference schemes

03/25/2021
by   Mohit Tekriwal, et al.
0

The behavior of physical systems is typically modeled using differential equations which are too complex to solve analytically. In practical problems, these equations are discretized on a computational domain, and numerical solutions are computed. A numerical scheme is called convergent, if in the limit of infinitesimal discretization, the bounds on the discretization error is also infinitesimally small. The approximate solution converges to the "true solution" in this limit. The Lax equivalence theorem enables a proof of convergence given consistency and stability of the method. In this work, we formally prove the Lax equivalence theorem using the Coq Proof Assistant. We assume a continuous linear differential operator between complete normed spaces, and define an equivalent mapping in the discretized space. Given that the numerical method is consistent (i.e., the discretization error tends to zero as the discretization step tends to zero), and the method is stable (i.e., the error is uniformly bounded), we formally prove that the approximate solution converges to the true solution. We then demonstrate convergence of the difference scheme on an example problem by proving both its consistency and stability, and then applying the Lax equivalence theorem. In order to prove consistency, we use the Taylor-Lagrange theorem by formally showing that the discretization error is bounded above by the nth power of the discretization step, where n is the order of the truncated Taylor polynomial.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/11/2022

Formal verification of iterative convergence of numerical algorithms

Physical systems are usually modeled by differential equations, but solv...
research
07/28/2019

Equivalence of pth moment stability between stochastic differential delay equations and their numerical methods

In this paper, a general theorem on the equivalence of pth moment stabil...
research
03/24/2023

Higher order time discretization method for a class of semilinear stochastic partial differential equations with multiplicative noise

In this paper, we consider a new approach for semi-discretization in tim...
research
02/18/2022

Central limit theorem for full discretization of parabolic SPDE

In order to characterize the fluctuation between the ergodic limit and t...
research
11/29/2020

Error Inhibiting Methods for Finite Elements

Finite Difference methods (FD) are one of the oldest and simplest method...
research
06/21/2022

Efficient Numerical Schemes for Multidimensional Population Balance Models

Multidimensional population balance models (PBMs) describe chemical and ...

Please sign up or login with your details

Forgot password? Click here to reset