Discrete Routh Reduction

This paper develops the theory of abelian Routh reduction for discrete mechanical systems and applies it to the variational integration of mechanical systems with abelian symmetry. The reduction of variational Runge-Kutta discretizations is considered, as well as the extent to which symmetry reduction and discretization commute. These reduced methods allow the direct simulation of dynamical features such as relative equilibria and relative periodic orbits that can be obscured or difficult to identify in the unreduced dynamics. The methods are demonstrated for the dynamics of an Earth orbiting satellite with a non-spherical $J_2$ correction, as well as the double spherical pendulum. The $J_2$ problem is interesting because in the unreduced picture, geometric phases inherent in the model and those due to numerical discretization can be hard to distinguish, but this issue does not appear in the reduced algorithm, where one can directly observe interesting dynamical structures in the reduced phase space (the cotangent bundle of shape space), in which the geometric phases have been removed. The main feature of the double spherical pendulum example is that it has a nontrivial magnetic term in its reduced symplectic form. Our method is still efficient as it can directly handle the essential non-canonical nature of the symplectic structure. In contrast, a traditional symplectic method for canonical systems could require repeated coordinate changes if one is evoking Darboux' theorem to transform the symplectic structure into canonical form, thereby incurring additional computational cost. Our method allows one to design reduced symplectic integrators in a natural way, despite the noncanonical nature of the symplectic structure.Comment: 24 pages, 7 figures, numerous minor improvements, references added, fixed typo

Variational and Geometric Structures of Discrete Dirac Mechanics

In this paper, we develop the theoretical foundations of discrete Dirac mechanics, that is, discrete mechanics of degenerate Lagrangian/Hamiltonian systems with constraints. We first construct discrete analogues of Tulczyjew's triple and induced Dirac structures by considering the geometry of symplectic maps and their associated generating functions. We demonstrate that this framework provides a means of deriving discrete Lagrange-Dirac and nonholonomic Hamiltonian systems. In particular, this yields nonholonomic Lagrangian and Hamiltonian integrators. We also introduce discrete Lagrange-d'Alembert-Pontryagin and Hamilton-d'Alembert variational principles, which provide an alternative derivation of the same set of integration algorithms. The paper provides a unified treatment of discrete Lagrangian and Hamiltonian mechanics in the more general setting of discrete Dirac mechanics, as well as a generalization of symplectic and Poisson integrators to the broader category of Dirac integrators.Comment: 26 pages; published online in Foundations of Computational Mathematics (2011