We introduce a new class of parametrized structure-preserving partitioned Runge-Kutta (α-PRK) methods for Hamiltonian systems with holonomic constraints. The methods are symplectic for any fixed scalar parameter α, and are reduced to the usual symplectic PRK methods like Shake-Rattle method or PRK schemes based on Lobatto IIIA-IIIB pairs when α = 0. We provide a new variational formulation for symplectic PRK schemes and use it to prove that the α-PRK methods can preserve the quadratic invariants for Hamiltonian systems subject to holonomic constraints. Meanwhile, for any given consistent initial values (p0, q0) and small step size h > 0, it is proved that there exists α * = α(h, p0, q0) such that the Hamiltonian energy can also be exactly preserved at each step. Based on this, we propose some energy and quadratic invariants preserving α-PRK methods. These α-PRK methods are shown to have the same convergence rate as the usual PRK methods and perform very well in various numerical experiments.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.