High order operator splitting methods based on an integral deferred correction framework

Integral deferred correction (IDC) methods have been shown to be an efficient way to achieve arbitrary high order accuracy and possess good stability properties. In this paper, we construct high order operator splitting schemes using the IDC procedure to solve initial value problems (IVPs). We present analysis to show that the IDC methods can correct for both the splitting and numerical errors, lifting the order of accuracy by r with each correction, where r is the order of accuracy of the method used to solve the correction equation. We further apply this framework to solve partial differential equations (PDEs). Numerical examples in two dimensions of linear and nonlinear initial–boundary value problems are presented to demonstrate the performance of the proposed IDC approach.