On the convergence of spectral deferred correction methods

In this work we analyze the convergence properties of the Spectral Deferred Correction (SDC) method originally proposed by Dutt et al. [BIT, 40 (2000), pp. 241--266]. The framework for this high-order ordinary differential equation (ODE) solver is typically described wherein a low-order approximation (such as forward or backward Euler) is lifted to higher order accuracy by applying the same low-order method to an error equation and then adding in the resulting defect to correct the solution. Our focus is not on solving the error equation to increase the order of accuracy, but on rewriting the solver as an iterative Picard integral equation solver. In doing so, our chief finding is that it is not the low-order solver that picks up the order of accuracy with each correction, but it is the underlying quadrature rule of the right hand side function that is solely responsible for picking up additional orders of accuracy. Our proofs point to a total of three sources of errors that SDC methods carry: the error at the current time point, the error from the previous iterate, and the numerical integration error that comes from the total number of quadrature nodes used for integration. The second of these two sources of errors is what separates SDC methods from Picard integral equation methods; our findings indicate that as long as difference between the current and previous iterate always gets multiplied by at least a constant multiple of the time step size, then high-order accuracy can be found even if the underlying "solver" is inconsistent the underlying ODE. From this vantage, we solidify the prospects of extending spectral deferred correction methods to a larger class of solvers to which we present some examples.

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

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

[3]  Robert Speck,et al.  Spectral deferred corrections with fast-wave slow-wave splitting , 2016, SIAM J. Sci. Comput..

[4]  Martin Weiser,et al.  Faster SDC convergence on non-equidistant grids by DIRK sweeps , 2015 .

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

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

[7]  T. Jahnke.  Geometric Numerical Integration , 2013 .

[8]  J. Kuntzmann,et al.  Neuere Entwicklungen der Methode von Runge und Kutta , 1961 .

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

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

[11]  Michael L. Minion,et al.  TOWARD AN EFFICIENT PARALLEL IN TIME METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS , 2012 .

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

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

[14]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[15]  Max Duarte,et al.  High order schemes based on operator splitting and deferred corrections for stiff time dependent PDEs , 2014, 1407.0195.

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

[17]  Yuan Liu,et al.  High order operator splitting methods based on an integral deferred correction framework , 2014, J. Comput. Phys..

[18]  P. Alam.  ‘S’ , 2021, Composites Engineering: An A–Z Guide.

[19]  Jingfang Huang,et al.  A Numerical Framework for Integrating Deferred Correction Methods to Solve High Order Collocation Formulations of ODEs , 2015, Journal of Scientific Computing.

[20]  Wei Guo,et al.  A high order time splitting method based on integral deferred correction for semi-Lagrangian Vlasov simulations , 2014, J. Comput. Phys..

[21]  István Faragó,et al.  Note on the Convergence of the Implicit Euler Method , 2012, NAA.

[22]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

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

[24]  Colin B. Macdonald,et al.  Parallel High-Order Integrators , 2010, SIAM J. Sci. Comput..

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

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

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

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

[29]  L. Einkemmer.  Structure preserving numerical methods for the Vlasov equation , 2016, 1604.02616.

[30]  Per-Olof Persson,et al.  Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations , 2017, J. Comput. Phys..

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

[32]  Michael L. Minion,et al.  A HYBRID PARAREAL SPECTRAL DEFERRED CORRECTIONS METHOD , 2010 .

[33]  Anita T. Layton,et al.  On the efficiency of spectral deferred correction methods for time-dependent partial differential equations , 2009 .

[34]  Andrew J. Christlieb,et al.  Implicit Parallel Time Integrators , 2011, J. Sci. Comput..

[35]  Colin B. Macdonald,et al.  Revisionist integral deferred correction with adaptive step-size control , 2013, 1310.6331.

[36]  Tommaso Buvoli.  A Class of Exponential Integrators Based on Spectral Deferred Correction , 2014 .

[37]  P. Zadunaisky.  On the estimation of errors propagated in the numerical integration of ordinary differential equations , 1976 .

[38]  S. Kadioglu,et al.  An Essentially Non-Oscillatory Spectral Deferred Correction Method for Conservation Laws , 2016 .

[39]  Tao Tang,et al.  High-Order Convergence of Spectral Deferred Correction Methods on General Quadrature Nodes , 2013, J. Sci. Comput..

[40]  Rolf Krause,et al.  A multi-level spectral deferred correction method , 2013, BIT Numerical Mathematics.

[41]  Katherine J. Evans,et al.  A Spectral Deferred Correction Method Applied to the Shallow Water Equations on a Sphere , 2013 .

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