Energy and quadratic invariants preserving methods for Hamiltonian systems with holonomic constraints

DOI10.4208/JCM.2106-M2020-0205zbMATH Open1515.65312arXiv2007.06338MaRDI QIDQ6344965FDOQ6344965

Authors: Lei Li, Dongling Wang 

Publication date: 13 July 2020

Abstract: We introduce a new class of parametricization structure-preserving partitioned Runge-Kutta ($alpha$-PRK) methods for Hamiltonian systems with holonomic constraints. When the scalar parameter $alpha=0$, the methods are reduced to the usual symplectic PRK methods like Shake-Rattle method or PRK schemes based on Lobatto IIIA-IIIB pairs, which can preserve all the quadratic invariants and the constraints. When $alphaeq0$, the methods are also shown to preserve all the quadratic invariants and the constraints manifold exactly. At the same time, for any given consistent initial values $(p_0,q_0)$ and small step size $h>0$, it is proved that there exists $alph{a}^{*}=alpha(h,p_0,q_0)$ such that the Hamiltonian energy can also be exactly preserved at each step. We provide a new variational formulation for symplectic PRK schemes and use it to prove that the parametrized PRK methods can preserve the quadratic invariants for Hamiltonian systems subject to holonomic constraints. The parametric $alpha$-PRK methods are shown to have the same convergence rate as the usual PRK methods and perform very well in various numerical experiments.

Mathematics Subject Classification ID

Numerical methods for Hamiltonian systems including symplectic integrators (65P10) Numerical methods for initial value problems involving ordinary differential equations (65L05)

This page was built for publication: Energy and quadratic invariants preserving methods for Hamiltonian systems with holonomic constraints

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6344965)