Overcoming order reduction in diffusion-reaction splitting. Part 1: Dirichlet boundary conditions

Abstract

For diffusion-reaction equations employing a splitting procedure is attractive as it reduces the computational demand and facilitates a parallel implementation. Moreover, it opens up the possibility to construct second-order integrators that preserve positivity. However, for boundary conditions that are neither periodic nor of homogeneous Dirichlet type order reduction limits its usefulness. In the situation described the Strang splitting procedure is not more accurate than Lie splitting. In this paper, we propose a splitting procedure that, while retaining all the favorable properties of the original method, does not suffer from order reduction. We demonstrate our results by conducting numerical simulations in one and two space dimensions with inhomogeneous and time dependent Dirichlet boundary conditions. In addition, a mathematical rigorous convergence analysis is conducted that confirms the results observed in the numerical simulations.

Cited by in corpus (9)

• Stochastic Simulation of Reaction-Diffusion Systems: A Fluctuating-Hydrodynamics Approach
• Overcoming order reduction in diffusion-reaction splitting. Part 2: oblique boundary conditions
• A $\mu$-mode integrator for solving evolution equations in Kronecker form
• Efficient boundary corrected Strang splitting
• A comparison of boundary correction methods for Strang splitting
• Splitting physics-informed neural networks for inferring the dynamics of integer- and fractional-order neuron models
• An adaptive low-rank splitting approach for the extended Fisher--Kolmogorov equation
• Avoiding order reduction when integrating diffusion-reaction boundary value problems with exponential splitting methods
• Nonsmooth data error estimates for exponential Runge-Kutta methods and applications to split exponential integrators