On higher-order semi-explicit symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems
algorithmsmotion of the double pendulumnumerical resultssemi-explicit symplectic partitioned Runge-Kutta methodsseparable Hamiltonian systems
Dynamical aspects of finite-dimensional Hamiltonian and Lagrangian systems (37J99) Numerical methods for initial value problems involving ordinary differential equations (65L05) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Stability and convergence of numerical methods for ordinary differential equations (65L20) Error bounds for numerical methods for ordinary differential equations (65L70)
The author introduces a new family of partitioned Runge-Kutta algorithms (PRK) for approximation to the solution of separable Hamiltonian systems with holonomic constraints. The new algorithms are based on modification of explicit PRK methods for unconstrained systems in the form given by \textit{D. Okunbor} and \textit{R. D. Skeel} [Math. Comput. 59, No. 200, 439-455 (1992; Zbl 0765.65079)]. They are semi-explicit and in applications require less computation than PRK programs in current use. An analysis is made of the order of error of the algorithms based on a more general study of the backward error for numerical integrators defined on a \(C^\infty\) manifold \(\mathcal M\) contained in \(\mathbb{R}^m\). It is proved that, except for orders one and two, the order of error in the constrained system may be less than the order of error of the corresponding unconstrained system. The author describes two extensions of the program which are designed so that the order of error of the constrained program is the same as the order of error of the unconstrained program. The first of these requires a special choice of coefficients in the approximating routine, and the second is based on composition of PRK algorithms which have error of first or second order. Numerical results are presented which result from applying the three programs to the equations for the motion of the double pendulum in a gravitational field. The effective order of error is found to be similar to error constants for unconstrained explicit PRK methods discussed by \textit{R. I. McLachlan} [SIAM J. Sci. Comput. 16, No. 1, 151-168 (1995; Zbl 0821.65048)].
- Symplectic Partitioned Runge–Kutta Methods for Constrained Hamiltonian Systems
- Global modified Hamiltonian for constrained symplectic integrators
- Specialized Partitioned Additive Runge–Kutta Methods for Systems of Overdetermined DAEs with Holonomic Constraints
- Symmetrized half-explicit methods for constrained mechanical systems
- Energy and Quadratic Invariants Preserving Methods for Hamiltonian Systems With Holonomic Constraints
- Symplecticness conditions of some low order partitioned methods for non-autonomous Hamiltonian systems
- Partitioned Runge-Kutta Methods for Separable Hamiltonian Problems
- Symplectic trigonometrically fitted partitioned Runge-Kutta methods
- Structure Preservation for Constrained Dynamics with Super Partitioned Additive Runge--Kutta Methods
- scientific article; zbMATH DE number 1437533
- Global modified Hamiltonian for constrained symplectic integrators
- Line integral solution of Hamiltonian systems with holonomic constraints
- Projected Runge-Kutta methods for constrained Hamiltonian systems
- Specialized Partitioned Additive Runge–Kutta Methods for Systems of Overdetermined DAEs with Holonomic Constraints
- Structure Preservation for Constrained Dynamics with Super Partitioned Additive Runge--Kutta Methods
- 5. Efficient Higher Order Time Discretization Schemes For Hamilton–Jacobi–Bellman Equations Based On Diagonally Implicit Symplectic Runge–Kutta Methods
- Numerical integration of Hamiltonian problems by G-symplectic methods
- Parametric symplectic partitioned Runge-Kutta methods with energy-preserving properties for Hamiltonian systems
- Symplectic Partitioned Runge–Kutta Methods for Constrained Hamiltonian Systems
- Partitioned Runge-Kutta Methods for Separable Hamiltonian Problems
- Energy and Quadratic Invariants Preserving Methods for Hamiltonian Systems With Holonomic Constraints
- Arbitrary-order trigonometric Fourier collocation methods for multi-frequency oscillatory systems
- Arbitrarily high-order energy-conserving methods for Hamiltonian problems with quadratic holonomic constraints
This page was built for publication: On higher-order semi-explicit symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1358154)