Publication Date:
2021-01-21
Description:
A Hamiltonian system subject to smooth constraints can typically be viewed as a Hamiltonian system on a manifold. Numerical computations, however, must be performed in $ R^n$. In this paper, canonical transformations from ``Hamiltonian differential--algebraic equations'' to ODEs in Euclidean space are considered. In \S2, canonical parameterizations or local charts are developed and it is shown how these can be computed in a practical framework. In \S3 we consider the construction of unconstrained Hamiltonian ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an integral invariant and whose flow reduces to the flow of the constrained system along the manifold. It is shown that certain of these unconstrained Hamiltonian systems force Lyapunov stability of the constraint--invariants, while others lead to an unstable invariant. In \S4, we compare various projection techniques which might be incorporated to better insure preservation of the constraint--invariants in the context of numerical discretization. Numerical experiments illustrate the degree to which the constraint and symplectic invariants are maintained under discretization of various formulations. {\bf Keywords:} differential--algebraic equations, Hamiltonian systems, canonical discretization schemes. {\bf AMS(MOS):} subject classification 65L05.
Keywords:
ddc:000
Language:
English
Type:
reportzib
,
doc-type:preprint
Format:
application/pdf
