Abstract
In this article a broad class of systems of implicit differential–algebraic equations (DAEs) is considered, including the equations of mechanical systems with holonomic and nonholonomic constraints. Solutions to these DAEs can be approximated numerically by applying a class of super partitioned additive Runge–Kutta (SPARK) methods. Several properties of the SPARK coefficients, satisfied by the family of Lobatto IIIA-B-C-C*-D coefficients, are crucial to deal properly with the presence of constraints and algebraic variables. A main difficulty for an efficient implementation of these methods lies in the numerical solution of the resulting systems of nonlinear equations. Inexact modified Newton iterations can be used to solve these systems. Linear systems of the modified Newton method can be solved approximately with a preconditioned linear iterative method. Preconditioners can be obtained after certain transformations to the systems of nonlinear and linear equations. These transformations rely heavily on specific properties of the SPARK coefficients. A new truly parallelizable preconditioner is presented.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
K.E. Brenan, S.L. Campbell and L.R. Petzold, Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, SIAM Classics in Applied Mathematics, 2nd ed. (SIAM, Philadelphia, PA, 1996).
P.N. Brown, A.C. Hindmarsh and L.R. Petzold, Using Krylov methods in the solution of large-scale differential-algebraic systems, SIAM J. Sci. Comput. 15 (1994) 1467–1488.
G. Golub and C.F. van Loan, Matrix Computations, Johns Hopkins Studies in the Mathematical Science, 3rd ed. (Johns Hopkins Univ. Press, Baltimore/London, 1996).
A. Greenbaum, Iterative Methods for Solving Linear Systems, Frontiers in Applied Mathematics, Vol. 17 (SIAM, Philadelphia, PA, 1997).
E. Hairer, Ch. Lubich and M. Roche, The Numerical Solution of Differential-Algebraic Systems by Runge-Kutta Methods, Lecture Notes in Mathematics, Vol. 1409 (Springer, Berlin, 1989).
E. Hairer and G. Wanner, Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems, Comput. Math., Vol. 14, 2nd revised ed. (Springer, Berlin, 1996).
E.J. Haug, Computer Aided Kinematics and Dynamics of Mechanical Systems, Vol. I: Basic Methods (Allyn and Bacon, Boston, USA, 1989).
L.O. Jay, Runge-Kutta type methods for index three differential-algebraic equations with applications to Hamiltonian systems, Ph.D. thesis, Department ofMathematics, University of Geneva, Switzerland (1994).
L.O. Jay, Structure preservation for constrained dynamics with super partitioned additive Runge-Kutta methods, SIAM J. Sci. Comput. 20 (1998) 416–446.
L.O. Jay, A parallelizable preconditioner for the iterative solution of implicit Runge-Kutta type methods, J. Comput. Appl. Math. 111 (1999) 63–76.
L.O. Jay, Inexact simplified Newton iterations for implicit Runge-Kutta methods, SIAM J. Numer. Anal. 38 (2000) 1369–1388.
L.O. Jay, Preconditioning and parallel implementation of implicit Runge-Kutta methods, Technical Report, Dept. of Math., University of Iowa, USA (2001) in preparation.
P.J. Rabier and W.C. Rheinboldt, Nonholonomic Motion of Rigid Mechanical Systems from a DAE Viewpoint (SIAM, Philadelphia, PA, 2000).
Y. Saad, Iterative Methods for Sparse Linear Systems (PWS, Boston, MA, 1996).
Y. Saad and M.H. Schultz, GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 7 (1986) 856–869.
W.O. Schiehlen, ed., Multibody Systems Handbook (Springer, Berlin, 1990).
W.O. Schiehlen, ed., Advanced Multibody System Dynamics, Simulation and Software Tools (Kluwer Academic, London, 1993).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Jay, L.O. Iterative Solution of SPARK Methods Applied to DAEs. Numerical Algorithms 31, 171–191 (2002). https://doi.org/10.1023/A:1021116800010
Issue Date:
DOI: https://doi.org/10.1023/A:1021116800010