Discrepancy in osculating elements between DSST and numerical propagators with 2x0 gravity field
I've noticed a discrepancy in osculating elements when comparing orbits between DSST and numerical propagators.
Configuration:
- LEO satellite
- Gravity field 2x0 (but an higher degree x order gives the same results)
- Zonal and Tesseral DSST forces (with or without J2² mean effect)
- Propagation on 50 orbits
Result:
- Discrepancy between osculating elements on all 6 orbital parameters
- Ex. bias of about 20m on semi-major axis and parabolic evolution of inclination leading to a mean error of 0.4 mdeg after 50 orbits → See figure below