Integral deferred correction methods constructed with high order Runge-Kutta integrators

Spectral deferred correction (SDC) methods for solving ordinary differential equations (ODEs) were introduced by Dutt, Greengard and Rokhlin (2000). It was shown in that paper that SDC methods can achieve arbitrary high order accuracy and possess nice stability properties. Their SDC methods are constructed with low order integrators, such as forward Euler or backward Euler, and are able to handle stiff and non-stiff terms in the ODEs. In this paper, we use high order Runge-Kutta (RK) integrators to construct a family of related methods, which we refer to as integral deferred correction (IDC) methods. The distribution of quadrature nodes is assumed to be uniform, and the corresponding local error analysis is given. The smoothness of the error vector associated with an IDC method, measured by the discrete Sobolev norm, is a crucial tool in our analysis. The expected order of accuracy is demonstrated through several numerical examples. Superior numerical stability and accuracy regions are observed when high order RK integrators are used to construct IDC methods.

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

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

[3]  Yinhua Xia,et al.  Efficient time discretization for local discontinuous Galerkin methods , 2007 .

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

[5]  Anders C. Hansen,et al.  On the order of deferred correction , 2011 .

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

[7]  Winfried Auzinger,et al.  Modified Defect Correction Algorithms for ODEs. Part I: General Theory , 2004, Numerical Algorithms.

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

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

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

[11]  Michael L. Minion,et al.  Semi-implicit projection methods for incompressible flow based on spectral deferred corrections , 2004 .

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

[13]  Robert D. Skeel,et al.  A Theoretical Framework for Proving Accuracy Results for Deferred Corrections , 1982 .

[14]  Michael L. Minion,et al.  Conservative multi-implicit spectral deferred correction methods for reacting gas dynamics , 2004 .

[15]  Anders C. Hansen,et al.  Convergence , 1992, Journal of Political Economy.

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

[17]  Thomas Hagstrom,et al.  On the spectral deferred correction of splitting methods for initial value problems , 2006 .

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