论文信息 - Comparison of splitting algorithms for the rigid body

Comparison of splitting algorithms for the rigid body

We compare several different second-order splitting algorithms for the asymmetric rigid body, with the aim of determining which one produces the smallest energy error for a given rigid body, namely, for given moments of inertia. The investigation is based on the analysis of the dominant term of the modified Hamiltonian and indicates that different algorithms can produce energy errors which differ by several orders of magnitude. As a byproduct of this analysis we remark that, for the special case of a flat rigid body with moments of inertia proportional to (1, 0.75, 0.25), one of the considered algorithms is in fact of order four.

Francesco Fassò

[1] R. McLachlan,et al. The accuracy of symplectic integrators , 1992 .

[2] S. Reich. Momentum conserving symplectic integrators , 1994 .

[3] Francesco Fassò,et al. A Changing-Chart Symplectic Algorithm for Rigid Bodies and Other Hamiltonian Systems on Manifolds , 2001, SIAM J. Sci. Comput..

[4] B. Leimkuhler,et al. Symplectic splitting methods for rigid body molecular dynamics , 1997 .

[5] R. McLachlan,et al. Explicit Lie-Poisson integration and the Euler equations. , 1993, Physical review letters.

[6] E. Leimanis. The General Problem of the Motion of Coupled Rigid Bodies about a Fixed Point , 1965 .

[7] K. Brown,et al. Graduate Texts in Mathematics , 1982 .

[8] J. Marsden,et al. Introduction to mechanics and symmetry , 1994 .

[9] Jack Wisdom,et al. Lie-Poisson integrators for rigid body dynamics in the solar system , 1994 .

[10] Ralph Abraham,et al. Foundations Of Mechanics , 2019 .