On the construction of geometric integrators in the RKMK class
نویسنده
چکیده
We consider the construction of geometric integrators in the class of RKMK methods. Any di erential equation in the form of an in nitesimal generator on a homogeneous space is shown to be locally equivalent to a di erential equation on the Lie algebra corresponding to the Lie group acting on the homogenous space. This way we obtain a distinction between the coordinate-free phrasing of the di erential equation and the local coordinates used. In this paper we study methods based on arbitrary local coordinates on the Lie group manifold. By choosing the coordinates to be canonical coordinates of the rst kind we obtain the original method of Munthe-Kaas [14]. Methods similar to the RKMK method are developed based on the di erent coordinatizations of the Lie group manifold, given by the Cayley transform, diagonal Pad e approximants of the exponential map, canonical coordinates of the second kind, etc. Some numerical experiments are also given.
منابع مشابه
Norges Teknisk-naturvitenskapelige Universitet Cost Efficient Lie Group Integrators in the Rkmk Class Cost Efficient Lie Group Integrators in the Rkmk Class
In this work a systematic procedure is implemented in order to minimise the computational cost of the Runge–Kutta–Munthe-Kaas (RKMK) class of Lie-group solvers. The process consists of the application of a linear transformation to the stages of the method and the analysis of a graded free Lie algebra to reduce the number of commutators involved. We consider here RKMK integration methods up to o...
متن کاملNumerical Integration of Lie-Poisson Systems While Preserving Coadjoint Orbits and Energy
In this paper we apply geometric integrators of the RKMK type to the problem of integrating Lie– Poisson systems numerically. By using the coadjoint action of the Lie group G on the dual Lie algebra g∗ to advance the numerical flow, we devise methods of arbitrary order that automatically stay on the coadjoint orbits. First integrals known as Casimirs are retained to machine accuracy by the nume...
متن کاملVariational Integrators for Underactuated Mechanical Control Systems with Symmetries
Optimal control problems for underactuated mechanical systems can be seen as a higher-order variational problem subject to higher-order constraints (that is, when the Lagrangian function and the constraints depend on higher-order derivatives such as the acceleration, jerk or jounces). In this paper we discuss the variational formalism for the class of underactuated mechanical control systems wh...
متن کاملConstructing Equivalence-preserving Dirac Variational Integrators with Forces
The dynamical motion of mechanical systems possesses underlying geometric structures, and preserving these structures in numerical integration improves the qualitative accuracy and reduces the long-time error of the simulation. For a single mechanical system, structure preservation can be achieved by adopting the variational integrator construction. This construction has been generalized to mor...
متن کاملHamilton-Pontryagin Integrators on Lie Groups
In this thesis structure-preserving time integrators for mechanical systems whose configuration space is a Lie Group are derived from a Hamilton-Pontryagin (HP) variational principle. In addition to its attractive properties for degenerate mechanical systems, the HP viewpoint also affords a practical way to design discrete Lagrangians, which are the cornerstone of variational integration theory...
متن کامل