A changing-chart symplectic algorithm for rigid bodies and other hamiltonian systems on manifolds

• Published in: SIAM journal on scientific computing 2002, Vol.23 (4), p.1189-1203
• Main Authors: BENETTIN, Giancarlo, CHERUBINI, Anna Maria, FASSO, Francesco
• Format: Article
• Language: English
• Subjects: Algorithms; Classical and quantum physics: mechanics and fields; Classical mechanics of discrete systems: general mathematical aspects; Energy conservation; Euclidean space; Exact sciences and technology; Methods; Physics
• ISSN: 1064-8275; 1095-7197
• DOI: 10.1137/S1064827500381720

We revive the elementary idea of constructing symplectic integrators for Hamiltonian flows on manifolds by covering the manifold with the charts of an atlas, implementing the algorithm in each chart (thus using coordinates) and switching among the charts whenever a coordinate singularity is approached. We show that this program can be implemented successfully by using a splitting algorithm if the Hamiltonian is the sum H1+H2 of two (or more) integrable Hamiltonians. Profiting from integrability, we compute exactly the flows of H1 and H2 in each chart and thus compute the splitting algorithm on the manifold by means of its representative in any chart. This produces a symplectic algorithm on the manifold which possesses an interpolating Hamiltonian, and hence it has excellent properties of conservation of energy. We exemplify the method for a point constrained to the sphere and for a symmetric rigid body under the influence of positional potential forces.