- #1
mistergrinch
- 44
- 0
Hello, I'm trying to solve the rocket equation in a nonuniform gravitational field. The standard rocket equation at a distance x from the center of a planet of mass M gives:
m*d^2x/dt^2 + v_e*dm/dt = -GmM/x^2
where v_e is the relative exhaust velocity.
Multiplying by dx/dt and assuming dm/dt=a (const.), this can be written as:
d/dt( 1/2m(dx/dt)^2 - GMm/x + v_e*a*x ) = 0 =>
1/2m(dx/dt)^2 - GMm/x + v_e*a*x = K (const.)
I can then separate variables and integrate, which gives me a very ugly expression for t in terms of x.
Another approach is to write the original equation as:
dv/dt + 1/m*dm/t*v_e = -GM/x^2 which can be integrated to give
v = v_e * ln(m0/m) - GM * integral_0_t{1/x^2 * dt}
I.e. the rocket eqn. modified by an integral of the gravitational force. But I have no idea what to do with this integral so I don't see how this helps me!
I was wondering if anyone has a better idea how to solve this problem. How do they solve problems like this at NASA, or are the rocket burns short enough that they just assume F_gravity is constant? This seems like a useful problem to be able to solve -- I saw a reference at google saying Tsiolkovsky was the first to solve it, but I can't find any references to any solutions. Does anyone have any ideas?
m*d^2x/dt^2 + v_e*dm/dt = -GmM/x^2
where v_e is the relative exhaust velocity.
Multiplying by dx/dt and assuming dm/dt=a (const.), this can be written as:
d/dt( 1/2m(dx/dt)^2 - GMm/x + v_e*a*x ) = 0 =>
1/2m(dx/dt)^2 - GMm/x + v_e*a*x = K (const.)
I can then separate variables and integrate, which gives me a very ugly expression for t in terms of x.
Another approach is to write the original equation as:
dv/dt + 1/m*dm/t*v_e = -GM/x^2 which can be integrated to give
v = v_e * ln(m0/m) - GM * integral_0_t{1/x^2 * dt}
I.e. the rocket eqn. modified by an integral of the gravitational force. But I have no idea what to do with this integral so I don't see how this helps me!
I was wondering if anyone has a better idea how to solve this problem. How do they solve problems like this at NASA, or are the rocket burns short enough that they just assume F_gravity is constant? This seems like a useful problem to be able to solve -- I saw a reference at google saying Tsiolkovsky was the first to solve it, but I can't find any references to any solutions. Does anyone have any ideas?