Efficient integration of large stiff systems of ODEs with exponential propagation iterative (EPI) methods

A new class of exponential propagation techniques which we call exponential propagation iterative (EPI) methods is introduced in this paper. It is demonstrated how for large stiff systems these schemes provide an efficient alternative to standard integrators for computing solutions over long time intervals. The EPI methods are constructed by reformulating the integral form of a solution to a nonlinear autonomous system of ODEs as an expansion in terms of products between special functions of matrices and vectors that can be efficiently approximated using Krylov subspace projections. The methodology for constructing EPI schemes is presented and their performance is illustrated using numerical examples and comparisons with standard explicit and implicit integrators. The history of the exponential propagation type integrators and their connection with EPI schemes are also discussed.

[1]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[2]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[3]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[4]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[5]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[6]  G. Nicolis,et al.  Chemical instabilities and sustained oscillations. , 1971, Journal of theoretical biology.

[7]  J. NAGUMOt,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 2006 .

[8]  P. Bellan,et al.  Three-dimensional Model of the Structure and Evolution of Coronal Mass Ejections , 2002 .

[9]  W. Arnoldi.  The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[10]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[11]  H. V. D. Vorst,et al.  An iterative solution method for solving f ( A ) x = b , using Krylov subspace information obtained for the symmetric positive definite matrix A , 1987 .

[12]  Andr'e Nauts,et al.  New Approach to Many-State Quantum Dynamics: The Recursive-Residue-Generation Method , 1983 .

[13]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[14]  Mayya Tokman,et al.  Magnetohydrodynamic modeling of solar magnetic arcades using exponential propagation methods , 2001 .

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

[16]  S. Thomas,et al.  WEIGHTED QUADRATURE IN KRYLOV , 2022 .

[17]  S. Krogstad.  Generalized integrating factor methods for stiff PDEs , 2005 .

[18]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[19]  J. M. Keiser,et al.  A New Class of Time Discretization Schemes for the Solution of Nonlinear PDEs , 1998 .

[20]  C. H. WADDINGTON,et al.  Towards a Theoretical Biology , 1968, Nature.

[21]  Homer F. Walker,et al.  Choosing the Forcing Terms in an Inexact Newton Method , 1996, SIAM J. Sci. Comput..

[22]  David A. Pope.  An exponential method of numerical integration of ordinary differential equations , 1963, CACM.

[23]  R. Dembo,et al.  INEXACT NEWTON METHODS , 1982 .

[24]  C. W. Gear,et al.  Iterative Solution of Linear Equations in ODE Codes , 1983 .

[25]  J. D. Lawson.  Generalized Runge-Kutta Processes for Stable Systems with Large Lipschitz Constants , 1967 .

[26]  M. Hochbruck,et al.  Exponential Runge--Kutta methods for parabolic problems , 2005 .

[27]  I. N. Sneddon,et al.  The Solution of Ordinary Differential Equations , 1987 .

[28]  T. Park,et al.  Unitary quantum time evolution by iterative Lanczos reduction , 1986 .

[29]  S. P. Nørsett.  An A-stable modification of the Adams-Bashforth methods , 1969 .

[30]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[31]  Yousef Saad,et al.  Efficient Solution of Parabolic Equations by Krylov Approximation Methods , 1992, SIAM J. Sci. Comput..

[32]  L. Tuckerman,et al.  A method for exponential propagation of large systems of stiff nonlinear differential equations , 1989 .

[33]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[34]  J. Butcher.  Coefficients for the study of Runge-Kutta integration processes , 1963, Journal of the Australian Mathematical Society.

[35]  L. Knizhnerman.  Calculation of functions of unsymmetric matrices using Arnoldi's method , 1991 .

[36]  H. Schwan,et al.  Biological Engineering , 1970 .

[37]  Lloyd N. Trefethen,et al.  Fourth-Order Time-Stepping for Stiff PDEs , 2005, SIAM J. Sci. Comput..

[38]  E. C. Zeeman,et al.  Differential equations for the heartbeat and nerve impulse , 1973 .