Advanced Search

Query:

Fill in the form and click »Search«...

Format:

Display: entries per page entries

Zbl 1232.65177
Leok, Melvin; Zhang, Jingjing
Discrete Hamiltonian variational integrators. (English)
[J] IMA J. Numer. Anal. 31, No. 4, 1497-1532 (2011). ISSN 0272-4979; ISSN 1464-3642/e

This paper is concerned with the study of some properties of a discrete Hamiltonian dynamics corresponding to a continuous Hamiltonian system with a Type II generating flow. Recall that an $n$-dimensional Hamiltonian system defined by $H = H(p, q)$, $q=(q^1, q^2, \dots , q^n)^T$, $p=(p_1,p_2, \dots ,p_n)^T$, is said of Type II generating flow if the solutions $q=q(t)$, $p=p(t)$, $t \in [0,T]$ make stationary the functional $F[q(\cdot), p(\cdot)]\equiv p(T) q(T) - \int_0^T [ p(t) \dot{q}(t) - H( q(t), p(t) )] dt$ in the set of smooth functions with the boundary conditions $\delta q(0)=0$, $\delta p(T)=0$. Here due to the relation between the Type II continuous and discrete Hamiltonian dynamics through a variational formulation the authors are able to derive discretizations of Hamiltonian systems that retain the essential properties of the flow. Furthermore, when the Hamiltonian is hyperregular it is possible through the Legendre transformation to translate this correspondence to the Lagrangian formulation. \par The paper is organized as follows: In Section 2 the variational formulation of Type II continuous systems is revised together the generating functions. Next, in Section 3, for a given partition $ 0 =t_0< t_1< \cdots <t_N=T$ of $[0,T]$ and a discrete curve $ \{ (q_k, p_k)\}_{k=0}^N$ associated to this partition, a discrete Type II Hamilton's variational principle associated to a discrete functional $F_d$ of $F$ is presented. The discrete right Hamilton's equations and the discrete generating function are given. The theory is clearly summarized at the end of Section 3. In Section 4, Galerkin type discretizations are considered for the practical solution of these discrete Hamiltonians. For a given set of basis functions $\{\psi_i(\tau)\}_{i=1}^s$, $ \tau \in [0,1]$ and quadrature nodes $\{c_i\}_{i=1}^s$ the authors give a systematic procedure to construct a Galerkin variational integrator, presenting the Lagrangian point of view when the Hamiltonian is hyperregular. It is shown that a special class of Galerkin variational integrators can be implemented as symplectic partitioned Runge-Kutta methods (SPRK). Some simple examples are presented to illustrate the theory. In Section 5 the preservation properties of the discrete Hamiltonian are considered. Finally, in Section 6, the results of some numerical experiments are given to illustrate the preservation properties of some simple Type II Hamiltonians.
[Manuel Calvo (Zaragoza)]

MSC 2000:: *65P10 Hamiltonian systems including symplectic integrators
37M15 Symplectic integrators
65L06 Multistep, Runge-Kutta, and extrapolation methods
65L60 Finite numerical methods for ODE
37J05 Relations with symplectic geometry and topology
37J45 Periodic, homoclinic and heteroclinic orbits, etc.

Keywords: Hamiltonian mechanics; variational integrators; symplectic integrators; preserving properties in Hamiltonian systems; discrete Hamiltonian dynamics; type II generating flow; Legendre transformation; Lagrangian formulation; generating functions; Galerkin type discretizations; Galerkin variational integrators; Runge-Kutta methods; numerical experiments

Comment on this Item PDF MathML XML ASCII DVI PS BibTeX Online Ordering Article Journal

	Scientific prize winners of the ICM 2010
	Overhang
	Lie groups, physics and geometry. An introduction for physicists, engineers and chemists.

Advanced Search

Zentralblatt MATH Berlin [Germany]