- #1
TheHarvesteR
- 14
- 0
Hi, first post here at PF :)
I have a problem here regarding orbit propagation. Basically my situation is as follows:
I have coded a system that can track the orbital parameters from an object in a simulated orbit (Basic rigidbody physics). The system takes the state vectors of the orbiting rigidbody, and from those it computes the keplerian parameters. So far, so good.
Then, the system allows removing this tracked object from the physics simulation, to have it's orbit propagated using the parameters computed earlier (or any arbitrary parameters at this point). This also works well, for closed orbits.
My problem now is figuring out how to numerically propagate parabolic and hyperbolic trajectories. The current system fails with those, since several of the calculated parameters result in NaN when eccentricity > 1.
What happens is since I'm getting all this from state vectors, each parameter has a few dependencies on other parameters calculated earlier, and if a single one fails to produce a valid result, all others that depend on it will also fail.
More specifically, I'm getting a NaN for eccentric and mean anomalies on hyperbolic/parabolic trajectories... I don't know if this is correct or not. But these values are very much necessary in my current system to get a position as a function of time.
EDIT: Well, I've continued trying different things, and I now have what seems to be valid mean and eccentric anomaly values.
The problem is that for the resulting eccentric anomaly, my true anomaly is always NaN for hyperbolic orbits, and PI for parabolic orbits.
I've discovered that there is a different form of Kepler's equation for hyperbolic orbits, so I added a new solver for that (which hopefully is doing things right). But still no luck with that true anomaly.
What I ultimately need is to find out how to determine the position of an object in a hyperbolic trajectory as a function of time.
Oh and one more question: I'm somewhat confused as to whether my semimajor axis should be positive or negative on hyp orbits. Right now I'm doing tests with positive SMA values, but the orbit tracking code (which calculates parameters from state vectors) gives out a negative SMA. Which is more correct here, positive or negative?
Any help at this point would be immensely appreciated. :)
Thanks in advance,
Cheers
I have a problem here regarding orbit propagation. Basically my situation is as follows:
I have coded a system that can track the orbital parameters from an object in a simulated orbit (Basic rigidbody physics). The system takes the state vectors of the orbiting rigidbody, and from those it computes the keplerian parameters. So far, so good.
Then, the system allows removing this tracked object from the physics simulation, to have it's orbit propagated using the parameters computed earlier (or any arbitrary parameters at this point). This also works well, for closed orbits.
My problem now is figuring out how to numerically propagate parabolic and hyperbolic trajectories. The current system fails with those, since several of the calculated parameters result in NaN when eccentricity > 1.
What happens is since I'm getting all this from state vectors, each parameter has a few dependencies on other parameters calculated earlier, and if a single one fails to produce a valid result, all others that depend on it will also fail.
More specifically, I'm getting a NaN for eccentric and mean anomalies on hyperbolic/parabolic trajectories... I don't know if this is correct or not. But these values are very much necessary in my current system to get a position as a function of time.
EDIT: Well, I've continued trying different things, and I now have what seems to be valid mean and eccentric anomaly values.
The problem is that for the resulting eccentric anomaly, my true anomaly is always NaN for hyperbolic orbits, and PI for parabolic orbits.
I've discovered that there is a different form of Kepler's equation for hyperbolic orbits, so I added a new solver for that (which hopefully is doing things right). But still no luck with that true anomaly.
What I ultimately need is to find out how to determine the position of an object in a hyperbolic trajectory as a function of time.
Oh and one more question: I'm somewhat confused as to whether my semimajor axis should be positive or negative on hyp orbits. Right now I'm doing tests with positive SMA values, but the orbit tracking code (which calculates parameters from state vectors) gives out a negative SMA. Which is more correct here, positive or negative?
Any help at this point would be immensely appreciated. :)
Thanks in advance,
Cheers
Last edited: