Non-singular orbital elements for special perturbations in the two-body problem

Seven spatial elements and a time element are proposed as the state variables of a new special perturbation method for the two-body problem. The new elements hold for zero eccentricity and inclination and for negative values of the total energy. They are developed by combining a spatial transformation into projective coordinates (as in the Burdet–Ferrándiz regularization) with a time transformation in which the exponent of the orbital radius is equal to one instead of two (as commonly done in the literature). By following this approach, we discover a new linearization of the two-body problem, from which the orbital elements can be generated by the variation of parameters method. The geometrical significance of the spatial quantities is revealed by a new intermediate frame which differs from a local vertical local horizontal frame by one rotation in the instantaneous orbital plane. Four elements parametrize the attitude in space of this frame, which in turn defines the orientation of the orbital plane and fixes the departure direction for the longitude of the propagated body. The remaining three elements determine the motion along the radial unit vector and the orbital longitude. The performance of the method, tested using a series of benchmark orbit propagation scenarios, is extremely good when compared to several regularized formulations, some of which have been modified and improved here for the first time.