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.

[1]  K. Bohmer.  Defect Correction Methods: Theory and Applications , 1984 .

[2]  J. Brandts.  [Review of: W. Hundsdorfer, J.G. Verwer (2003) Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations] , 2006 .

[3]  G. Quispel,et al.  Acta Numerica 2002: Splitting methods , 2002 .

[4]  Jim Douglas,et al.  Numerical solution of two‐dimensional heat‐flow problems , 1955 .

[5]  Jingfang Huang,et al.  Accelerating the convergence of spectral deferred correction methods , 2006, J. Comput. Phys..

[6]  Michael L. Minion,et al.  Implications of the choice of predictors for semi-implicit Picard integral deferred correction methods , 2007 .

[7]  Benjamin W. Ong,et al.  COMMENTS ON HIGH-ORDER INTEGRATORS EMBEDDED WITHIN INTEGRAL DEFERRED CORRECTION METHODS , 2009 .

[8]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[9]  Kaitai Li,et al.  A third accurate operator splitting method , 2011, Math. Comput. Model..

[10]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[11]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[12]  G. Quispel,et al.  Geometric integrators for ODEs , 2006 .

[13]  Gilbert Strang.  Approximating semigroups and the consistency of difference schemes , 1969 .

[14]  M. Minion.  Semi-implicit spectral deferred correction methods for ordinary differential equations , 2003 .

[15]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[16]  Klaus Böhmer,et al.  Defect Correction Methods , 1984, Computing Supplementum.

[17]  Anita T. Layton,et al.  On the choice of correctors for semi-implicit Picard deferred correction methods , 2008 .

[18]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[19]  Michael L. Minion,et al.  Implications of the Choice of Quadrature Nodes for Picard Integral Deferred Corrections Methods for Ordinary Differential Equations , 2005 .

[20]  Jim Douglas.  On Accuracy of Alternating Direction Implicit Methods for Parabolic Equations , 1999 .

[21]  A. Ostermann,et al.  High order splitting methods for analytic semigroups exist , 2009 .

[22]  Andrew J. Christlieb,et al.  Integral deferred correction methods constructed with high order Runge-Kutta integrators , 2009, Math. Comput..

[23]  Gurij Ivanovich Marchuk,et al.  Some application of splitting-up methods to the solution of mathematical physics problems , 1968 .

[24]  Mengping Zhang,et al.  STRONG STABILITY PRESERVING PROPERTY OF THE DEFERRED CORRECTION TIME DISCRETIZATION , 2008 .

[25]  Jingfang Huang,et al.  Arbitrary order Krylov deferred correction methods for differential algebraic equations , 2007, J. Comput. Phys..

[26]  Vladimir Rokhlin,et al.  A Highly Accurate Solver for Stiff Ordinary Differential Equations , 2012, SIAM J. Sci. Comput..

[27]  J. Schnakenberg,et al.  Simple chemical reaction systems with limit cycle behaviour. , 1979, Journal of theoretical biology.

[28]  J. J. Douglas.  On the Numerical Integration of $\frac{\partial ^2 u}{\partial x^2 } + \frac{\partial ^2 u}{\partial y^2 } = \frac{\partial u}{\partial t}$ by Implicit Methods , 1955 .

[29]  Paul C. Fife,et al.  Mathematical Aspects of Reacting and Diffusing Systems , 1979 .

[30]  Mechthild Thalhammer,et al.  High-Order Exponential Operator Splitting Methods for Time-Dependent Schrödinger Equations , 2008, SIAM J. Numer. Anal..

[31]  J. Geiser,et al.  Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives , 2007 .

[32]  Benjamin W. Ong,et al.  Semi-implicit integral deferred correction constructed with additive Runge-Kutta methods , 2011 .

[33]  A. Bourlioux,et al.  High-order multi-implicit spectral deferred correction methods for problems of reactive flow , 2003 .

[34]  Bernie D. Shizgal,et al.  Pseudospectral method of solution of the Fitzhugh-Nagumo equation , 2009, Math. Comput. Simul..

[35]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .