Arbitrary order Krylov deferred correction methods for differential algebraic equations

In this paper, a new framework for the construction of accurate and efficient numerical methods for differential algebraic equation (DAE) initial value problems is presented. The methods are based on applying spectral deferred correction techniques as preconditioners to a Picard integral collocation formulation for the solution. The resulting preconditioned nonlinear system is solved using Newton-Krylov schemes such as the Newton-GMRES method. Least squares based orthogonal polynomial approximations are computed using Gaussian type quadratures, and spectral integration is used to avoid the numerically unstable differentiation operator. The resulting Krylov deferred correction (KDC) methods are of arbitrary order of accuracy and very stable. Preliminary results show that these new methods are very competitive with existing DAE solvers, particularly when high precision is desired.

[1]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[2]  D. Gottlieb,et al.  Numerical analysis of spectral methods , 1977 .

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

[4]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[5]  P. E. Zadunaisky.  A method for the estimation of errors propagated in the numerical solution of a system of ordinary differential equations , 1966 .

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

[7]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[8]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[9]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[10]  K. Chen,et al.  Matrix preconditioning techniques and applications , 2005 .

[11]  Prashanth K. Vijalapura,et al.  Fractional step methods for index-1 differential-algebraic equations , 2005 .

[12]  Ernst Hairer,et al.  The numerical solution of differential-algebraic systems by Runge-Kutta methods , 1989 .

[13]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[14]  L. Trefethen,et al.  An Instability Phenomenon in Spectral Methods , 1987 .

[15]  V. Rokhlin,et al.  Fast algorithms for polynomial interpolation, integration, and differentiation , 1996 .

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

[17]  Aaditya V. Rangan.  Ph.D. thesis: Adaptive solvers for partial differential and differential-algebraic equations , 2003 .

[18]  Leslie Greengard,et al.  Spectral integration and two-point boundary value problems , 1991 .

[19]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[20]  J. Damasceno,et al.  Numerical simulation of two-dimensional flows over a circular cylinder using the immersed boundary method , 2003 .

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

[22]  Petter Kolm,et al.  Quadruple and octuple layer potentials in two dimensions , 2000 .

[23]  Winfried Auzinger,et al.  Modified defect correction algorithms for ODEs. Part II: Stiff initial value problems , 2005, Numerical Algorithms.

[24]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[25]  J. M. Sanz-Serna,et al.  Convergence and order reduction of Runge-Kutta schemes applied to evolutionary problems in partial differential equations , 1987 .

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

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

[28]  Bradley K. Alpert,et al.  A Fast Algorithm for the Evaluation of Legendre Expansions , 1991, SIAM J. Sci. Comput..

[29]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[30]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

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

[32]  Victor Pereyna.  Iterated deferred corrections for nonlinear boundary value problems , 1968 .

[33]  Mari Paz Calvo,et al.  Avoiding the order reduction of Runge-Kutta methods for linear initial boundary value problems , 2002, Math. Comput..

[34]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

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

[36]  C. Kelley.  Solving Nonlinear Equations with Newton's Method , 1987 .

[37]  L. Petzold.  A description of dassl: a differential/algebraic system solver , 1982 .