-
Notifications
You must be signed in to change notification settings - Fork 69
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Bizarre trajectories in frames tied to Pol #1741
Comments
Update: It appears that the state of the symmetric linear multistep integrator gets somehow damaged by close encounters between Bop and Pol. @eggrobin wrote a test (#1744) that demonstrates that the local error grows by steps each time a close encounter occurs. Around day 409 of year 4, for instance, the two bodies get within a few thousand kilometres of each other and the local error grows from 1 to 500 metres per day. In Cohen, where the velocities were computed using a backward difference method, the position errors resulted in amplified velocity errors, thus the picture above. After #1734 the velocities only use one position difference together with the accelerations, and the errors are much less visible (but still present). We are still trying to understand better how the integrator gets damaged, but this is clearly tied to the system being on the edge of instability. We are also trying to understand its Ляпунов exponent and the backward error. A possible fix will be to stop using a symmetric linear multistep integrator for the KSP system, and go for a symplectic Runge-Kutta-Nyström integrator. The errors will be larger but not subject to catastrophic growth. For the real solar system, the symmetric linear multistep integrator works just fine because that system doesn't have close encounters (good news everyone!). |
Reported by Eriksonn on the forum:
There is a rather embarrassing bug in Cohen (and, we suspect, Clifford and Christoffel) where trajectories drawn in a reference frame involving both Pol and Jool (either Pol-Centred-Jool-Aligned or Jool-Pol-Barycentric) appears horribly wiggly:
Interestingly enough, this is already fixed by #1734. More interestingly, we don't have a real understanding of what's happening here other than "something is ill-conditioned". Clearly the backward difference computation is part of the problem, but there may be deeper issues.
The text was updated successfully, but these errors were encountered: