GRSS Propagator Module (grss.prop)

The propagator within GRSS is based on the IAS15 algorithm [1], which is a 15th-order integrator based on Gauss-Radau quadrature. The algorithm is adaptive, meaning that the step size is adjusted to ensure that the error is below a certain threshold. However, the integrator can also be run in fixed-step mode, where a user-defined fixed time step is used until the end of the integration.

The force model includes the following effects:

1. Newtonian point-mass gravity

2. Einstein-Infeld-Hoffmann (EIH) formulation of point-mass relativistic gravity [2] [3]

3. J₂ zonal harmonic effects from the Sun and the Earth

4. A₁, A₂, A₃ nongravitational acceleration model [4]

In addition to any integrated bodies defined by the user, the default force model includes effects from the Sun, the planets, the Moon, Pluto, and the Big 16 main-belt asteroids. The states of these bodies are read in from either the JPL DE431 or DE441 ephemerides, based on the user’s choice. The propagator can also interpolate any of the integrated bodies’ states to any time within the integration (or within a margin around the integration time bounds). This interpolation is done via interpolation polynomials based on the Newton divided difference method.

The propagator module also has the ability to calculate apparent states of any integrated body. This is done through two methods. The first is the simpler of the two, and uses a simple approximation of the light-time between an observer and the integrated body. The second method is more accurate, and uses an iterative method to calculate the light-time between an observer and the integrated body. The second method is slightly more computationally expensive since it requires a few extra interpolations.

This ability to calculate apparent states is used to calculate optical observables (Right Ascension and Declination). Furthermore, the propagator suite can also calculate range observables (delay and doppler). The delay measurement is trivially returned from the light-time calculation. The doppler measurable, however, is calculated using a dedicated doppler measurement model. This entire codebase is written in C++ and a Python binding is using pybind11.

References

[1]

Rein, H. and Spiegel D.S., “IAS15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits”, Monthly Notices of the Royal Astronomical Society, Volume 446, Issue 2, Pages 1424-1437, 2015. https://doi.org/10.1093/mnras/stu2164.

[2]

Einstein, A., Infeld, L., and Hoffmann, B., “The Gravitational Equations and the Problem of Motion”, Annals of Mathematics, Volume 39, Issue 1, Pages 65-100, 1938. https://doi.org/10.2307/1968714.

[3]

Moyer, T.D., “Mathematical formulation of the Double Precision Orbit Determination Program (DPODP)”, Jet Propulsion Laboratory Technical Report 32-1527, 1971. https://ntrs.nasa.gov/citations/19710017134.

[4]

Marsden, B.G., Sekanina, Z., and Yeomans, D.K., “Comets and nongravitational forces. V”, The Astronomical Journal, Volume 78, Issue 2, Pages 211-225, 1973. https://doi.org/10.1086/111402.