Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)



Specialized Runge-Kutta methods for index $ 2$ differential-algebraic equations

Author: Laurent O. Jay
Journal: Math. Comp. 75 (2006), 641-654
MSC (2000): Primary 65L05, 65L06, 65L80
Published electronically: December 19, 2005
MathSciNet review: 2196984
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: We consider the numerical solution of systems of semi-explicit index $ 2$ differential-algebraic equations (DAEs) by methods based on Runge-Kutta (RK) coefficients. For nonstiffly accurate RK coefficients, such as Gauss and Radau IA coefficients, the standard application of implicit RK methods is generally not superconvergent. To reestablish superconvergence projected RK methods and partitioned RK methods have been proposed. In this paper we propose a simple alternative which does not require any extra projection step and does not use any additional internal stage. Moreover, symmetry of Gauss methods is preserved. The main idea is to replace the satisfaction of the constraints at the internal stages in the standard definition by enforcing specific linear combinations of the constraints at the numerical solution and at the internal stages to vanish. We call these methods specialized Runge-Kutta methods for index $ 2$ DAEs (SRK-DAE$ 2$).

References [Enhancements On Off] (What's this?)

  • 1. U. Ascher and L. R. Petzold, Projected implicit Runge-Kutta methods for differential-algebraic equations, SIAM J. Numer. Anal. 28 (1991), 1097-1120. MR 1111456 (92f:65082)
  • 2. K. E. Brenan, S. L. Campbell, and L. R. Petzold, Numerical solution of initial-value problems in differential-algebraic equations, SIAM Classics in Appl. Math., SIAM, Philadelphia, Second Edition, 1996.MR 1363258 (96h:65083)
  • 3. R. P. K. Chan, P. Chartier, and A. Murua, Post-projected Runge-Kutta methods for index-$ 2$ differential-algebraic equations, Appl. Numer. Math. 42 (2002), 77-94. MR 1921330 (2003f:65124)
  • 4. -, Reversible methods of Runge-Kutta type for index-$ 2$ differential-algebraic equations, Numer. Math. 97 (2004), 427-440.MR 2059464 (2005b:65088)
  • 5. E. Hairer and L. O. Jay, Implicit Runge-Kutta methods for higher index differential-algebraic systems, WSSIAA Contributions in numerical mathematics, 2 (1993), 213-224. MR 1299761 (95g:65099)
  • 6. E. Hairer, Ch. Lubich, and M. Roche, The numerical solution of differential-algebraic systems by Runge-Kutta methods, Lect. Notes in Math., vol. 1409, Springer, Berlin, 1989. MR 1027594 (91a:65178)
  • 7. E. Hairer and G. Wanner, Solving ordinary differential equations ii. stiff and differential-algebraic problems, Comput. Math., vol. 14, Springer, Berlin, Second Revised Edition, 1996. MR 1439506 (97m:65007)
  • 8. L. O. Jay, Convergence of a class of Runge-Kutta methods for differential-algebraic systems of index $ 2$, BIT 33 (1993), 137-150.MR 1326008 (96a:65106)
  • 9. -, Solution of index $ 2$ implicit differential-algebraic equations by Lobatto Runge-Kutta methods, BIT 43 (2003), 91-104.MR 1981642 (2004g:65100)
  • 10. Ch. Lubich, On projected Runge-Kutta methods for differential-algebraic equations, BIT 31 (1991), 545-550. MR 1127491 (92h:65109)
  • 11. A. Murua, Partitioned Runge-Kutta methods for semi-explicit differential-algebraic systems of index $ 2$, Tech. Report EHU-KZAA-IKT-196, Univ. of the Basque country, 1996.

Similar Articles

Retrieve articles in Mathematics of Computation with MSC (2000): 65L05, 65L06, 65L80

Retrieve articles in all journals with MSC (2000): 65L05, 65L06, 65L80

Additional Information

Laurent O. Jay
Affiliation: Department of Mathematics, 14 MacLean Hall, The University of Iowa, Iowa City, Iowa 52242-1419
Email: E-mail address:

Keywords: Differential-algebraic equations, index $2$, Runge-Kutta methods
Received by editor(s): January 15, 2004
Received by editor(s) in revised form: January 26, 2005
Published electronically: December 19, 2005
Additional Notes: This material is based upon work supported by the National Science Foundation under Grant No. 9983708.
Article copyright: © Copyright 2005 American Mathematical Society
The copyright for this article reverts to public domain 28 years after publication.

American Mathematical Society