Abstract
The numerical integration of stiff mechanical systems is studied in which a strong potential forces the motion to remain close to a manifold. The equations of motion are written as a singular singular perturbation problem with a small stiffness parameter ɛ. Smooth solutions of such systems are characterized, in distinction to highly oscillatory general solutions. Implicit Runge-Kutta methods using step sizes larger than ɛ are shown to approximate smooth solutions, and precise error estimates are derived. As ɛ → 0, Runge-Kutta solutions of the stiff system converge to Runge-Kutta solutions of the associated constrained system formulated as a differential-algebraic equation of index 3. Standard software for stiff initial-value problems does not work satisfactorily on the stiff systems considered here. The reasons for this failure are explained, and remedies are proposed.
Similar content being viewed by others
References
R. Abraham and J. E. Marsden,Foundations of Mechanics. Benjamin/Cumming Publ. Co., Reading, Mass., 1978.
K. E. Brenan, S. L. Campbell and L. R. Petzold,Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations. North-Holland, New York 1989.
J. C. Butcher,On the implementation of implicit Runge-Kutta methods. BIT16, 237–240 (1976).
M. Crouzeix and J. Rappaz,On Numerical Approximation in Bifurcation Theory. Masson, Paris, and Springer-Verlag, Berlin 1990.
E. Eich and M. Hanke,Regularization methods for constrained mechanical multibody systems. Preprint, Humboldt-Universität zu Berlin 1991.
E. Hairer, Ch. Lubich and M. Roche,Error of Runge-Kutta methods for stiff problems studied via differential algebraic equations. BIT28, 678–700 (1988).
E. Hairer, Ch. Lubich and M. Roche,The Numerical Solution of Differential-Algrbraic Systems by Runge-Kutta Methods. Springer LNM 1409, 1989.
E. Hairer and G. Wanner,Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems. Springer-Verlag, Berlin 1991.
L. Jay,Collocation methods for differential-algebraic equations of index 3. Preprint, Section de Math., Univ. de Genève 1992.
L. Kalachev and R. E. O'Malley, Jr.,Regularization of nonlinear differential-algebraic equations. Report 92-4, Dept. Appl. Math., Univ. of Washington 1992.
T. Kato,Perturbation Theory for Linear Operators. 2nd Ed., Grundlehren math. Wiss. 132, Springer-Verlag, Berlin 1980.
U. Kirchgraber, F. Lasagni, K. Nipp and D. Stoffer,On the application of invariant manifold theory, in particular to numerical analysis. Int. Ser. Num. Math.97, 189–197 (1991).
M. Knorrenschild,Regularisierung von differentiell-algebraischen Systemen-Theoretische und numerische Aspekte. Doctoral dissertation, RWTH Aachen 1988.
P. Lötstedt,On a penalty function method for the simulation of mechanical systems subject to constraints. Report TRITA-NA-7919, Royal Inst. of Technology, Stockholm 1979.
D. Morgenstern and I. Szabó,Vorlesungen über theoretische Mechanik. Grundlehren math. Wiss. 112, Springer-Verlag, Berlin 1961.
K. Nipp and D. Stoffer,Attractive invariant manifolds for maps: Existence, smoothness, and continuous dependence on the map. Report 92-11, SAM, ETH Zürich 1992.
R. E. O'Malley, Jr.,Singular Perturbation Methods for Ordinary Differential Equations. Springer-Verlag, New York 1991.
Author information
Authors and Affiliations
Additional information
This work was supported in part by the Austrian Science Foundation, grant P8443-PHY.
Rights and permissions
About this article
Cite this article
Lubich, C. Integration of stiff mechanical systems by Runge-Kutta methods. Z. angew. Math. Phys. 44, 1022–1053 (1993). https://doi.org/10.1007/BF00942763
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF00942763