A LOCAL-GLOBAL VERSION OF A STEPSIZE CONTROL FOR RUNGE-KUTTA METHODS

  • Kulikov, G.Yu (Department of Mathematical Cybernetics, Faculty of Mechanics and Mathematics, Ulyanovsk State University)
  • Published : 2000.05.01

Abstract

In this paper we develop a new procedure to control stepsize for Runge- Kutta methods applied to both ordinary differential equations and semi-explicit index 1 differential-algebraic equation In contrast to the standard approach, the error control mechanism presented here is based on monitoring and controlling both the local and global errors of Runge- Kutta formulas. As a result, Runge-Kutta methods with the local-global stepsize control solve differential of differential-algebraic equations with any prescribe accuracy (up to round-off errors)

Keywords

References

  1. BIT v.35 Convergence of multistep discretizations of DAEs C. Arevalo;G. Soderlind
  2. Mech. Struct. and Mach. v.23 no.2 On the error control for multistep methods applied to ODEs with invariants and DAEs in multibody dynamics E. Eich;C. Fuhrer;J. Yen
  3. SIAM J. Numer. Anal. v.21 ODE methods for the solution of differential/algebraic systems C.W. Gear;L.R. Petzold
  4. Solving ordinary differential equationⅠ:Nonstiff problems E. Hairer;S.P. Nørsett;G. Wanner
  5. Lecture Note in Math. 1409 The numerical solution of differential-algebraic systems by Runge-Kutta methods E. Hairer;Ch. Lubich;M. Roche
  6. Solving ordinary differential equations Ⅱ: Stiff and differential-algebraic problems E. Hairer;G. Wanner
  7. Applied Numerical Mathematics v.15 The use of Butcher series in the analysis of Newton-like iterations in Runge-Kutta formulas K.R. Jackson;A. Kværnø;S.P. Nørsett
  8. (in Russian) Vestnik Moskov. Univ. Ser. 1 Mat. Mekh. no.1 A method for the numerical solution of the autonomous Cauchy problem with an algebraic restriction of the phase variables G.Yu. Kulikov
  9. translation in Moscow Univ. Math. Bull. v.47 no.1 G.Yu. Kulikov
  10. (in Russian) Vestnik Moskov. Univ. Ser. 1 Mat. Mekh. no.3 The numerical solution of an autonomous Cauchy problem with an algebraic constraint of the phase variables (the nonsingular case) G.Yu. Kulikov
  11. translation in Moscow Univ. Math. Bull. v.48 no.3 G.Yu. Kulikov
  12. (in Russian) Zh. Vychisl. Mat. Mat. Fiz. v.33 no.4 The numerical solution of the autonomous Cauchy problem with an algebraic relation between the phase variables G.Yu. Kulikov
  13. translation in Comp. Maths Math. Phys. v.33 no.4 G.Yu. Kulikov
  14. Technical report 7/1996, IMM, Technical University of Denmark, Lyngby Convergence and implementation of implicit Runge-Kutta methods for DAEs G.Yu. Kulikov;P.G. Thomsen
  15. (in Russian) Zh. Vychisl. Mat. Mat. Fiz. v.36 no.8 Convergence theorems for iterative Runge-Kutta methods with a constant integration step G.Yu. Kulikov
  16. translation in Comp. Maths Math. Phys. v.36 no.8 G.Yu. Kulikov
  17. Korean J. Comput. & Appl. Math. v.4 no.2 Numerical methods solving the semi-explicit differential-algebraic equations by implicit multistep fixed stepsize methods G.Yu. Kulikov
  18. (in Russian) Vestnik Moskov. Univ. Ser. 1 Mat. Mekh. no.5 On the numerical solution of the Cauchy problem for a system of differential-algebraic equations by means of implicit Runge-Kutta methods with a variable integration step G.Yu. Kulikov
  19. translation in Moscow Univ. Math. Bull. v.52 no.5 G.Yu. Kulikov
  20. Russian J. Numer. Anal. Math. Modelling v.12 no.6 A theory of symmetric one-step methods for differential-algebraic equations G.Yu. Kulikov
  21. (in Russian) Zh. Vychisl. Mat. Fiz. v.38 no.1 Numerical solution of the Cauchy problem for a system of differential-algebraic equations with the use of implicit Runge-Kutta methods with a nontrivial predictor G.Yu. Kulikov
  22. translation in Comp. Maths Math. Phys. v.38 no.1 G.Yu. Kulikov
  23. Korean J. Comput. & Appl. v.6 no.3 A local-global stepsize control for multistep methods applied to semi-explicit index 1 differential-algebraic equations G.Yu. Kulikov
  24. Technical report 1/1992, Mathematical Sciences Div. The order of Runge-Kuta methods applied to semi-explicit DAEs of index 1, using Newton-type iterations to compute the internal stage values A. kværnø