A higher order local linearization method for solving ordinary differential equations

The local linearization (LL) method for the integration of ordinary differential equations is an explicit one-step method that has a number of suitable dynamical properties. However, a major drawback of the LL integrator is that its order of convergence is only two. The present paper overcomes this limitation by introducing a new class of numerical integrators, called the LLT method, that is based on the addition of a correction term to the LL approximation. In this way an arbitrary order of convergence can be achieved while retaining the dynamic properties of the LL method. In particular, it is proved that the LLT method reproduces correctly the phase portrait of a dynamical system near hyperbolic stationary points to the order of convergence. The performance of the introduced method is further illustrated through computer simulations.

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

[2]  L. Dieci,et al.  Padé approximation for the exponential of a block triangular matrix , 2000 .

[3]  Antonella Zanna,et al.  Collocation and Relaxed Collocation for the Fer and the Magnus Expansions , 1999 .

[4]  Juan C. Jiménez,et al.  Local Linearization Method for Numerical Integration of Delay Differential Equations , 2006, SIAM J. Numer. Anal..

[5]  J. Butcher.  The Numerical Analysis of Ordinary Di erential Equa-tions , 1986 .

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

[7]  B. V. Pavlov,et al.  The method of local linearization in the numerical solution of stiff systems of ordinary differential equations , 1988 .

[8]  Fernando Casas,et al.  Improved High Order Integrators Based on the Magnus Expansion , 2000 .

[9]  Arieh Iserles.  Quadrature methods for stiff ordinary differential systems , 1981 .

[10]  Rolando J. Biscay,et al.  The Local Linearization Method for Numerical Integration of Random Differential Equations , 2005 .

[11]  Juan C. Jiménez,et al.  Rate of convergence of local linearization schemes for initial-value problems , 2005, Appl. Math. Comput..

[12]  Juan C. Jiménez,et al.  Dynamic properties of the local linearization method for initial-value problems , 2002, Appl. Math. Comput..

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

[14]  Juan C. Jiménez,et al.  Local Linearization-Runge Kutta (LLRK) Methods for Solving Ordinary Differential Equations , 2006, International Conference on Computational Science.

[15]  A. Iserles,et al.  Methods for the approximation of the matrix exponential in a Lie‐algebraic setting , 1999, math/9904122.

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

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

[18]  Marlis Hochbruck,et al.  Explicit Exponential Runge-Kutta Methods for Semilinear Parabolic Problems , 2005, SIAM J. Numer. Anal..

[19]  N. Higham.  The Scaling and Squaring Method for the Matrix Exponential Revisited , 2005, SIAM J. Matrix Anal. Appl..

[20]  Juan Carlos Jimenez,et al.  A simple algebraic expression to evaluate the local linearization schemes for stochastic differential equations , 2002, Appl. Math. Lett..

[21]  A. Iserles,et al.  On the Implementation of the Method of Magnus Series for Linear Differential Equations , 1999 .

[22]  J. Magnus,et al.  Matrix Differential Calculus with Applications in Statistics and Econometrics , 1991 .

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

[24]  John Carroll.  A Matricial Exponentially Fitted Scheme for the Numerical Solution of Stiff Initial-Value Problems , 1993 .

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

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

[27]  Arieh Iserles,et al.  Solving linear ordinary differential equations by exponentials of iterated commutators , 1984 .

[28]  James M. Bower,et al.  The Book of GENESIS , 1994, Springer New York.

[29]  W. Beyn.  On the Numerical Approximation of Phase Portraits Near Stationary Points , 1987 .

[30]  R. K. Jain.  Some A -Stable Methods for Stiff Ordinary Differential Equations , 1972 .

[31]  T. E. Hull,et al.  Comparing numerical methods for stiff systems of O.D.E:s , 1975 .

[32]  C. Loan.  Computing integrals involving the matrix exponential , 1978 .

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

[34]  G. Corliss,et al.  ATOMFT: solving ODEs and DAEs using Taylor series , 1994 .

[35]  L. M. Pedroso,et al.  Computing multiple integrals involving matrix exponentials , 2007 .

[36]  Ralph A. Willoughby,et al.  EFFICIENT INTEGRATION METHODS FOR STIFF SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS , 1970 .

[37]  L. M. Pedroso,et al.  Letter to the Editor: Computing multiple integrals involving matrix exponentials , 2008 .

[38]  T. Ozaki.  2 Non-linear time series models and dynamical systems , 1985 .

[39]  R. J. Biscay,et al.  Approximation of continuous time stochastic processes by the local linearization method revisited , 2002 .

[40]  J. C. Jimenez,et al.  Local Linearization method for the numerical solution of stochastic differential equations , 1996 .

[41]  J. Butcher.  The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

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

[43]  L. Perko.  Differential Equations and Dynamical Systems , 1991 .

[44]  V. Lakshmikantham,et al.  Stability Analysis of Nonlinear Systems , 1988 .

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

[46]  Juan I. Ramos,et al.  Piecewise-linearized methods for initial-value problems , 1997 .

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

[48]  Gene H. Golub,et al.  Matrix computations , 1983 .

[49]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .