Journal article
Symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems
SIAM Journal on Numerical Analysis, Vol.33(1), pp.368-387
1996
DOI: 10.1137/0733019
Abstract
This article deals with the numerical treatment of Hamiltonian systems with holonomic constraints. A class of partitioned Runge-Kutta methods, consisting of the couples of s-stage Lobatto IIIA and Lobatto IIIB methods, has been discovered to solve these problems efficiently. These methods are symplectic, preserve all underlying constraints, and are superconvergent with order 2s - 2. For separable Hamiltonians of the form H(q, p) = 1/2 pTM-1p + U(q) the Rattle algorithm based on the Verlet method was up to now the only known symplectic method preserving the constraints. In fact this method turns out to be equivalent to the 2-stage Lobatto IIIA-IIIB method of order 2. Numerical examples have been performed which illustrate the theoretical results.
Details
- Title: Subtitle
- Symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems
- Creators
- Laurent Jay
- Resource Type
- Journal article
- Publication Details
- SIAM Journal on Numerical Analysis, Vol.33(1), pp.368-387
- DOI
- 10.1137/0733019
- ISSN
- 0036-1429
- eISSN
- 1095-7170
- Language
- English
- Date published
- 1996
- Academic Unit
- Mathematics
- Record Identifier
- 9983985932302771
Metrics
20 Record Views