Numerical propagation of trajectories of Earth-orbiting spacecraft
Introduction
In this project we compute the compare the errors obtained from propagating the trajectories of different satellites with to models: our model and the SGP4 model. Our model consist in solving the following differential equation:
Here, $\mathbf{a}_{\mathrm{GP}}=(\ddot{x},\ddot{y},\ddot{z})$ is the acceleration caused by the geopotential (expanded in spherical harmonics) and $\ddot{x}$, $\ddot{y}$ and $\ddot{z}$ are given in Cunningham recursions; $\mathbf{a}_{\mathrm{drag}}$ is the acceleration caused by the atmospheric drag; $\mathbf{a}_{\mathrm{rad}}$ is the acceleration caused by the solar radiation pressure; $\mathbf{a}_{\mathrm{sun}}$ is the acceleration caused by the Sun; and $\mathbf{a}_{\mathrm{moon}}$ is the acceleration caused by the Moon. The coefficients $\delta_{i}\in\{0,1\}$ are used to enable and disable the different perturbations.
The initial conditions of the initial value problem will be the position and velocity obtained from the TLE. In order to solve this system of 6 differential equations, we have opted to use the Runge-Kutta-Fehlberg method of order 7(8).
Requirements
gcc and g++ compilers
make
python3
gnuplot
Compilation and execution of the code
Download the code:
[email protected]:victorballester7/final-bachelor-thesis.git
cd final-bachelor-thesis
chmod +x execute.sh
Here <t/f> are flags to enable or disable the different pertubations and the argument <tle/sgp4> is used to compare the errors with the TLE data set or with the SGP4 propagator. The argument `
The available satellites names are:
Low Earth Orbit:
HUBBLE
IRIDIUM_71
ISS
NUTSAT
STARLINK_1007
Medium Earth Orbit:
GALILEO_20
GLONASS
NAVSTAR_61
SIRIUS_3
Geostationary Orbit:
GALAXY_16
GALAXY_18
GALAXY_19
TDRS-3
TDRS-5
High Earth Orbit:
TESS
Example:
./execute.sh ISS f t t f f t sgp4 10
This command computes the errors of the ISS satellite with the SGP4 propagator and with the following perturbations enabled: Sun, Moon and atmospheric drag.
Results
We now expose the main results of this work:
Comparison of the errors for LEO satellites
ISS position error when considering the Earth as a point mass or as a non-homogeneous spherical distribution of mass (with the geopotential model).
Propagation of the ISS satellite when considering only the geopotential model for the Earth and the atmospheric drag.
Comparison of the errors for MEO satellites
Propagation of the Sirius-3 satellite considering the perturbations from the Moon, the Sun and the solar radiation pressure.
Propagation of the Galileo-20 satellite considering the perturbations from the Moon, the Sun and the solar radiation pressure.
Comparison of the errors for GEO satellites
Propagation of the TDRS-3 satellite considering the perturbations from the Moon, the Sun and the solar radiation pressure.