Almost symplectic Runge-Kutta schemes for Hamiltonian systems

JOURNAL OF COMPUTATIONAL PHYSICS

卷:203

Article

Tan, XB

关键词: geometric integrators;    Hamiltonian structure;    symplectic Runge-Kutta methods;    pseudo-symplecticity;    fixed-point iteration;    Newton's method;    convergence;   

DOI  :  10.1016/j.jcp.2004.08.012

来源: Elsevier

【 摘 要 】

Symplectic Runge-Kutta schemes for the integration of general Hamiltonian systems are implicit. In practice, one has to solve the implicit algebraic equations using some iterative approximation method, in which case the resulting integration scheme is no longer symplectic. In this paper, the preservation of the symplectic structure is analyzed under two popular approximation schemes, fixed-point iteration and Newton's method, respectively. Error bounds for the symplectic structure are established when N fixed-point iterations or N iterations of Newton's method are used. The implications of these results for the implementation of symplectic methods are discussed and then explored through numerical examples. Numerical comparisons with non-symplectic Runge-Kutta methods and pseudo-symplectic methods are also presented. (C) 2004 Elsevier Inc. All rights reserved.

【 授权许可】

Free