I am using GMAT to propagate satellite motion. At first, I excluded all the propagators besides harmonics(of second degree and 0th order, as shown in the picture):
Here are the results(x,y,z
coordinates) of propagating with the time interval of 10 days from the 668 km height:
1207.4141854713, 193.09221015360, 6947.8246951043
Then I added the solar radiation pressure(Spherical model):
Here are satellite coordinates after that:
1207.4150240009, 193.09806743665, 6947.8324711391
After that I excluded the radiation and added relativistic effects, as demonstrated here:
New coordinates are provided below:
1207.4097579535, 193.11627451516, 6947.8248222515
The difference between method was estimated by calculating the distance:
$$r_{crad}=\sqrt{(x_{c}-x_{rad})^2+(y_{c}-y_{rad})^2+(z_{c}-z_{rad})^2},$$ $$r_{crel}=\sqrt{(x_{c}-x_{rel})^2+(y_{c}-y_{rel})^2+(z_{c}-z_{rel})^2},$$
where $x_{c},y_{c},z_{c}$ are "clear" coordinates of satellite(i.e. without relativity or radiation), $x_{rad},y_{rad},z_{rad}$ and $x_{rel},y_{rel},z_{rel}$ are "radiation" and "relativistic" ones, respectively. The results of difference estimation: $r_{crad}=$9.8 m,$r_{rel}=$24 m. That means that relativistic factors affect more than solar radiation. I would like to know, does it make sense? And what can be a reason of it, if yes? Can this reason be the height of the satellite(the value of 668 km means that the satellite is in termosphere)?
-976.3107644649057, -4835.627052558522,-5031.728586125443,-0.7944031487816871,5.474532271429767,-5.094496750907486
. Dry mass is equal to 850 kg, the coefficient of relativity - 1.8, SRP area - 1 $m^2$ while drag area is 15 $m^2$(I understand that this value makes little sense but it was given by GMAT, and I will experiment on its changing). $\endgroup$ – Alex Johnson Aug 16 '18 at 2:34