- #1
MattRob
- 211
- 29
Sooo, this is something I'm really happy I figured out, buuuut, I want to make sure it's correct.
I guess the context doesn't really matter, so you can skip the slashed part if you like. Honestly I'm not sure what forum this belongs in (Mathematics? Topology? There's no "geometry" forum, heh), but the application of it is astronomical and has to do with orbits and launch windows, so I figured there'd be some expertise, here, at least.
Context/Application:
///
I got the full realism overhaul on Kerbal Space Program, and I'm planning a mission that involves a docked two-craft assembly inserting into a low orbit over Io, undocking, the lander landing then launching again to conduct a rendezvous with the orbiter, a la' Apollo. The only problem is, I can guarantee I won't be in an equatorial orbit, and I want to be a stickler for realism and get back up before radiation doses build up too much. Issue is, Io spins, so I'm no longer under the orbiter's orbit if I stay too long. But I've started to realize, that the launch windows to get into the correct orbital plane are a bit complicated to figure out, and that I could use my latitude to pick how long until that launch window opens again.
So I'm trying to figure out how to determine that launch window. If I can find this angle, then it's really simple math based on how quickly Io spins, to pick my latitude for a certain surface stay time.
///
The problem itself:
I've got an orbit inclined by a certain amount [itex]I[/itex]. I need to find the longitude/angle in-between the two points that this orbit intersects at any particular latitude. I'll call this angle [itex]\theta[/itex].
First step, to simplify the problem, I assumed that I could stretch out the sphere to a cylinder while preserving the relationship of theta to other relevant features (ie, inclination and latitude), so long as the orbit is stretched as well.
Next, I realized that viewed from above, the orbit and the planet are a circle (this is why the cylindrical projecting happened - if I'd viewed this from above without it, then the orbit would be an ellipse and this would be much more difficult).
And the angle suddenly became, instead of a complex problem in 3d space with a globe and an inclined circle, a simple problem of a 2d circle. "I know the math for this!"
If I use the center of the circle as the origin, then [itex]\phi = acos(\frac{x}{r})[/itex], and [itex]\theta = 2\phi[/itex] . Given an x-coordinate, I could find [itex]\theta[/itex]! So how to find that?
Well, there's an illustration here I haven't covered yet. Viewed on the X-Z plane, the relationship in-between [itex]z[/itex] and [itex]x[/itex] is obvious.
[itex]tan(I) = \frac{z}{x}[/itex], [itex]x = \frac{z}{tan(I)}[/itex]
Subbing back in,
[itex]\frac{1}{2}\theta = \phi = acos(\frac{z}{r tan(I)})[/itex]
And for z, I simply had to remember that I'm working with a cylinder instead of a sphere. My vertical distance from the equatorial plane ([itex]z[/itex]) is simply given by my radius and my latitude ([itex]\varphi[/itex]) :
*illustrating how despite the fact I've moved onto a cylinder, rather than a sphere, "latitude" still retains meaning, and how it carries over and its trigonometric relations.
[itex]tan(\varphi) = \frac{z}{r}[/itex]
Which leads to:
[itex]\theta = 2acos(\frac{tan(\varphi)}{tan(I)})[/itex]
Where [itex]I[/itex] is the inclination of the orbit to the equator, and [itex]\varphi[/itex] is the latitude.
So, did I get it right? I've shown that this is the case for a cylinder, but I'm wondering if my initial assumption about warping the sphere into a cylinder holds.
I'm just starting my Sophomore year as an undergraduate for "Physics-Astronomy," and I can't help but shake the feeling that this is a little bit what it must feel like to do original research, except it would be far more complex and something nobody's done before :p
A ton of fun, nonetheless. Note, I titled the thread more generally (and only after I wrote the main body of text) because I'm going to be doing a lot of these little things, all generally about orbital mechanics. I'm encountering a huge number of these problems as I plan this mission in Kerbal Space Program (gravity assists, finding relative inclination in-between two orbits given their LAN and [itex]I[/itex] to a parent body, how to find out when to launch to enter the ecliptic plane from a launch site at [itex]\varphi[/itex] latitude, etc.), and I figured I'd might as well keep it in one thread.
I guess the context doesn't really matter, so you can skip the slashed part if you like. Honestly I'm not sure what forum this belongs in (Mathematics? Topology? There's no "geometry" forum, heh), but the application of it is astronomical and has to do with orbits and launch windows, so I figured there'd be some expertise, here, at least.
Context/Application:
///
I got the full realism overhaul on Kerbal Space Program, and I'm planning a mission that involves a docked two-craft assembly inserting into a low orbit over Io, undocking, the lander landing then launching again to conduct a rendezvous with the orbiter, a la' Apollo. The only problem is, I can guarantee I won't be in an equatorial orbit, and I want to be a stickler for realism and get back up before radiation doses build up too much. Issue is, Io spins, so I'm no longer under the orbiter's orbit if I stay too long. But I've started to realize, that the launch windows to get into the correct orbital plane are a bit complicated to figure out, and that I could use my latitude to pick how long until that launch window opens again.
So I'm trying to figure out how to determine that launch window. If I can find this angle, then it's really simple math based on how quickly Io spins, to pick my latitude for a certain surface stay time.
///
The problem itself:
I've got an orbit inclined by a certain amount [itex]I[/itex]. I need to find the longitude/angle in-between the two points that this orbit intersects at any particular latitude. I'll call this angle [itex]\theta[/itex].
First step, to simplify the problem, I assumed that I could stretch out the sphere to a cylinder while preserving the relationship of theta to other relevant features (ie, inclination and latitude), so long as the orbit is stretched as well.
Next, I realized that viewed from above, the orbit and the planet are a circle (this is why the cylindrical projecting happened - if I'd viewed this from above without it, then the orbit would be an ellipse and this would be much more difficult).
And the angle suddenly became, instead of a complex problem in 3d space with a globe and an inclined circle, a simple problem of a 2d circle. "I know the math for this!"
If I use the center of the circle as the origin, then [itex]\phi = acos(\frac{x}{r})[/itex], and [itex]\theta = 2\phi[/itex] . Given an x-coordinate, I could find [itex]\theta[/itex]! So how to find that?
Well, there's an illustration here I haven't covered yet. Viewed on the X-Z plane, the relationship in-between [itex]z[/itex] and [itex]x[/itex] is obvious.
[itex]tan(I) = \frac{z}{x}[/itex], [itex]x = \frac{z}{tan(I)}[/itex]
Subbing back in,
[itex]\frac{1}{2}\theta = \phi = acos(\frac{z}{r tan(I)})[/itex]
And for z, I simply had to remember that I'm working with a cylinder instead of a sphere. My vertical distance from the equatorial plane ([itex]z[/itex]) is simply given by my radius and my latitude ([itex]\varphi[/itex]) :
*illustrating how despite the fact I've moved onto a cylinder, rather than a sphere, "latitude" still retains meaning, and how it carries over and its trigonometric relations.
[itex]tan(\varphi) = \frac{z}{r}[/itex]
Which leads to:
[itex]\theta = 2acos(\frac{tan(\varphi)}{tan(I)})[/itex]
Where [itex]I[/itex] is the inclination of the orbit to the equator, and [itex]\varphi[/itex] is the latitude.
So, did I get it right? I've shown that this is the case for a cylinder, but I'm wondering if my initial assumption about warping the sphere into a cylinder holds.
I'm just starting my Sophomore year as an undergraduate for "Physics-Astronomy," and I can't help but shake the feeling that this is a little bit what it must feel like to do original research, except it would be far more complex and something nobody's done before :p
A ton of fun, nonetheless. Note, I titled the thread more generally (and only after I wrote the main body of text) because I'm going to be doing a lot of these little things, all generally about orbital mechanics. I'm encountering a huge number of these problems as I plan this mission in Kerbal Space Program (gravity assists, finding relative inclination in-between two orbits given their LAN and [itex]I[/itex] to a parent body, how to find out when to launch to enter the ecliptic plane from a launch site at [itex]\varphi[/itex] latitude, etc.), and I figured I'd might as well keep it in one thread.
Last edited: