Error Inhibiting Schemes for Initial Boundary Value Heat Equation

by   Adi Ditkowski, et al.

In this paper, we elaborate the analysis of some of the schemes which were presented in arXiv:1711.07926 for the heat equation with periodic boundary conditions. We adopt this methodology to derive finite-difference schemes for heat equation with Dirichlet and Neumann boundary conditions, whose convergence rates are higher than their truncation errors. We call these schemes error inhibiting schemes. When constructing a semi-discrete approximation to a partial differential equation (PDE), a discretization of the spatial operator has to be derived. For stability, this discrete operator must be semi-bounded. Under this semi-boundness, the Lax-Richtmyer equivalence theorem assures that the scheme converges at most, as the order of the error will be at most of the order of the truncation error. Usually, the error is in indeed of the order of the truncation error. In this paper, we demonstrate that schemes can be constructed such that their errors are smaller than their truncation errors. This property can enable us design schemes which are more efficient than standard schemes.


Dissipation-preserving discretization of the Cahn–Hilliard equation with dynamic boundary conditions

This paper deals with time stepping schemes for the Cahn–Hilliard equati...

Stability and error analysis of a class of high-order IMEX schemes for Navier-stokes equations with periodic boundary conditions

We construct high-order semi-discrete-in-time and fully discrete (with F...

Weak intermittency and second moment bound of a fully discrete scheme for stochastic heat equation

In this paper, we first prove the weak intermittency, and in particular ...

On the Formalization of the Heat Conduction Problem in HOL

Partial Differential Equations (PDEs) are widely used for modeling the p...

Effects of round-to-nearest and stochastic rounding in the numerical solution of the heat equation in low precision

Motivated by the advent of machine learning, the last few years saw the ...

Algorithm for numerical solutions to the kinetic equation of a spatial population dynamics model with coalescence and repulsive jumps

An algorithm is proposed for finding numerical solutions of a kinetic eq...