A novel analytic continuation power series solution for the perturbed two-body problem

Inspired by the original developments of recursive power series by means of Lagrange invariants for the classical two-body problem, a new analytic continuation algorithm is presented and studied. The method utilizes kinematic transformation scalar variables differentiated to arbitrary order to generate the required power series coefficients. The present formulation is extended to accommodate the spherical harmonics gravity potential model. The scalar variable transformation essentially eliminates any divisions in the analytic continuation and introduces a set of variables that are closed with respect to differentiation, allowing for arbitrary-order time derivatives to be computed recursively. Leibniz product rule is used to produce the needed arbitrary-order expansion variables. With arbitrary-order time derivatives available, Taylor series-based analytic continuation is applied to propagate the position and velocity vectors for the nonlinear two-body problem. This foundational method has been extended to also demonstrate an effective variable step size control for the Taylor series expansion. The analytic power series approach is demonstrated using trajectory calculations for the main problem in satellite orbit mechanics including high-order spherical harmonics gravity perturbation terms. Numerical results are presented to demonstrate the high accuracy and computational efficiency of the produced solutions. It is shown that the present method is highly accurate for all types of studied orbits achieving 12-16 digits of accuracy (the extent of double precision). While this double-precision accuracy exceeds typical engineering accuracy, the results address the precision versus computational cost issue and also implicitly demonstrate the process to optimize efficiency for any desired accuracy. We comment on the shortcomings of existing power series-based general numerical solver to highlight the benefits of the present algorithm, directly tailored for solving astrodynamics problems. Such efficient low-cost algorithms are highly needed in long-term propagation of cataloged RSOs for space situational awareness applications. The present analytic continuation algorithm is very simple to implement and efficiently provides highly accurate results for orbit propagation problems. The methodology is also extendable to consider a wide variety of perturbations, such as third body, atmospheric drag and solar radiation pressure.