An improved approximate Newton method for implicit Runge-Kutta formulas

Abstract Implicit Runge–Kutta (IRK) methods (such as the s -stage Radau IIA method with s = 3 , 5 , or 7) for solving stiff ordinary differential equation systems have excellent stability properties and high solution accuracy orders, but their high computing costs in solving their nonlinear stage equations have seriously limited their applications to large scale problems. To reduce such a cost, several approximate Newton algorithms were developed, including a commonly used one called the simplified Newton method. In this paper, a new approximate Jacobian matrix and two new test rules for controlling the updating of approximate Jacobian matrices are proposed, yielding an improved approximate Newton method. Theoretical and numerical analysis show that the improved approximate Newton method can significantly improve the convergence and performance of the simplified Newton method.

[1]  Laurent O. Jay,et al.  Inexact Simplified Newton Iterations for Implicit Runge-Kutta Methods , 2000, SIAM J. Numer. Anal..

[2]  L. M. Skvortsov.  An efficient scheme for the implementation of implicit Runge-Kutta methods , 2008 .

[3]  Kjell Gustafsson,et al.  Control Strategies for the Iterative Solution of Nonlinear Equations in ODE Solvers , 1997, SIAM J. Sci. Comput..

[4]  E. Hairer,et al.  Stiff differential equations solved by Radau methods , 1999 .

[5]  John C. Butcher,et al.  On the implementation of implicit Runge-Kutta methods , 1976 .

[6]  R. Alexander.  Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[7]  Juan I. Montijano,et al.  Implementation of high-order implicit runge-kutta methods , 2001 .

[8]  Juan I. Montijano,et al.  Variable-order starting algorithms for implicit Runge-Kutta methods on stiff problems , 2002 .

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

[10]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[11]  M. N. Spijker,et al.  The error committed by stopping the Newton iteration in the numerical solution of stiff initial value problems , 1994 .

[12]  Hans Olsson,et al.  Stage Value Predictors and Efficient Newton Iterations in Implicit Runge-Kutta Methods , 1998, SIAM J. Sci. Comput..

[13]  S. González-Pinto,et al.  Speeding up Netwton-type iterations for stiff problems , 2005 .

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

[15]  Theodore A. Bickart,et al.  An Efficient Solution Process for Implicit Runge–Kutta Methods , 1977 .

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

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

[18]  L. Shampine.  Implementation of Implicit Formulas for the Solution of ODEs , 1980 .