Accurate computation of highly eccentric satellite orbits

For computing highly eccentric (e≳0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen.