On the Extension of Adams--Bashforth--Moulton Methods for Numerical Integration of Delay Differential Equations and Application to the Moon's Orbit
This work addresses the need for precise two-directional numerical integration of delay differential equations in celestial mechanics, but the contribution is incremental.
The authors extend Adams-Bashforth-Moulton methods to integrate delay differential equations backwards in time, enabling forward and backward numerical integration of the Moon's orbit with tidal time delays. They demonstrate the method's applicability but provide no concrete numerical results.
One of the problems arising in modern celestial mechanics is the need of precise numerical integration of dynamical equations of motion of the Moon. The action of tidal forces is modeled with a time delay and the motion of the Moon is therefore described by a functional differential equation (FDE) called delay differential equation (DDE). Numerical integration of the orbit is normally being performed in both directions (forwards and backwards in time) starting from some epoch (moment in time). While the theory of normal forwards-in-time numerical integration of DDEs is developed and well-known, integrating a DDE backwards in time is equivalent to solving a different kind of FDE called advanced differential equation, where the derivative of the function depends on not yet known future states of the function. We examine a modification of Adams--Bashforth--Moulton method allowing to perform integration of the Moon's DDE forwards and backwards in time and the results of such integration.