An initial–boundary corrected splitting method for diffusion-reaction problems

Strang splitting is a widely used second-order method for solving diffusion-reaction problems. However, its convergence order is often reduced to order $1$ for Dirichlet boundary conditions and to order $1.5$ for Neumann and Robin boundary conditions, leading to lower accuracy and reduced efficiency. In this paper, we consider a new splitting approach, called an initial-boundary corrected splitting, which avoids order reduction while improving computational efficiency for a wider range of applications. In contrast to the corrections proposed in the literature, it does not require the computation of correction terms that depend on the boundary conditions and boundary data. Through rigorous analytical convergence analysis and numerical experiments, we demonstrate the improved accuracy and performance of the proposed method.