Polyhedral Potential and Variational Integrator Computation of the Full Two Body Problem

Optimal control problems are formulated and efficient computational procedures are proposed for combined orbital and rotational maneuvers of a rigid body in three dimensions. The rigid body is assumed to act under the influence of forces and moments that arise from a potential and from control forces and moments. The key features of this paper are its use of computational procedures that are guaranteed to preserve the geometry of the optimal solutions. The theoretical basis f...

This paper presents a counterexample to the conjecture that the semi-explicit Lie-Newmark algorithm is variational. As a consequence the Lie-Newmark method is not well-suited for long-time simulation of rigid body-type mechanical systems. The counterexample consists of a single rigid body in a static potential field, and can serve as a test of the variational nature of other rigid-body integrators.

In this paper the exact analytical solution of the motion of a rigid body with arbitrary mass distribution is derived in the absence of forces or torques. The resulting expressions are cast into a form where the dependence of the motion on initial conditions is explicit and the equations governing the orientation of the body involve only real numbers. Based on these results, an efficient method to calculate the location and orientation of the rigid body at arbitrary times is ...

We construct a highly-symmetric periodic orbit of six bodies in three dimensions. In this orbit, binary collisions occur at the origin in a regular periodic fashion, rotating between pairs of bodies located on the coordinate axes. Regularization of the collisions in the orbit is achieved by an extension of the Levi-Civita method. Initial conditions for the orbit are found numerically. In contrast to an earlier periodic collision-based orbit in three dimensions, this orbit is ...

In this paper we present \texttt{SymOrb.jl}, a software which combines group representation theory and variational methods to provide numerical solutions of singular dynamical systems of paramount relevance in Celestial Mechanics and other interacting particles models. Among all, it prepares for large-scale search of symmetric periodic orbits for the classical $n$-body problem and their classification, paving the way towards a computational validation of Poincar\'e conjecture...

The three-body problem is a fundamental long-standing open problem, with applications in all branches of physics, including astrophysics, nuclear physics and particle physics. In general, conserved quantities allow to reduce the formulation of a mechanical problem to fewer degrees of freedom, a process known as dynamical reduction. However, extant reductions are either non-general, or hide the problem's symmetry or include unexplained definitions. This paper presents a dynami...

We develop an explicit, second-order, variational time integrator for full body dynamics that preserves the momenta of the continuous dynamics, such as linear and angular momenta, and exhibits near-conservation of total energy over exponentially long times. In order to achieve these properties, we parametrize the space of rotations using exponential local coordinates represented by a rescaled form of the Rodrigues rotation vector and we systematically derive the time integrat...

This letter introduces an advanced novel theory for calculating non-linear Newtonian hydrostatic perturbations in the density, shape, and gravitational field of fluid stars and planets subjected to external tidal and rotational forces. The theory employs a Lie group approach using exponential mappings to derive exact differential equations for large gravitational field perturbations and the shape function, which describes the finite deformation of the body's figure. This appr...

This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the algorithm to handle individual time steps; this introduces fifth-order errors in angular momentum conservation and symplecticity. We show that using adaptive block power of two timesteps does not increase the error in symplecticity. In contr...

Elegant integration schemes of second and fourth order for simulations of rigid body systems are presented which treat translational and rotational motion on the same footing. This is made possible by a recent implementation of the exact solution of free rigid body motion. The two schemes are time-reversible, symplectic, and exactly respect conservation principles for both the total linear and angular momentum vectors. Simulations of simple test systems show that the second o...