Perigee Reosculation fails for highly inclined circular orbits #2
Description
On specific orbits when reosculating, the method for finding the value of tau (secant method) will rarely be found to be outside of the possible range after a few iterations (in this case it was ~ -28). I’m still trying to find the exact reason that the reosculating is failing.
The ICs that this occurred on are:
r0 = [6678. 0. 0.]
v0 = [ 0. -0.89755918 7.67352051]
a=6678
e=0
i=1.6872355941793078
t0 =0, tf=54,000
It seems like after a few orbits the eccentricity grows large enough for the reosc code to trigger
Some additional info I’ve found. This only occurs with circular orbits, and occurs more quickly the more polar the orbit is. It also seems to always happen but might not trigger if the inclination is small. at 45 degrees it took more than 100 orbits to occur, and with almost any amount of eccentricity it does not occur.
I just ran a=6878 e=0.005 i=90 degrees for 1000 orbits without issue, when the same elements except 0 eccentricity would fail to osculate properly at ~5 orbits and 0.001 eccentricity would fail at ~10