Jump to content

Deriving the Equation for Burnout Height in KSP


bowtiesRcool

Recommended Posts

Hi all,

I have been trying to work on an equation I can use to calculate the burnout height of a single stage rocket in vertical flight from the KSP Launchpad where the forces of Drag and Gravity are both variable.

I start with:

SUM(F) = Fthrust + Fdrag +Fgravity

where,

Fdrag = (0.004892*Cd * Po)mrv2 e-h/H

Fgravity = G*Mmrocket/(r+h)2

Fthrust = some constant

I attempted to modify the derivation of Tsiolkovski's Ideal Rocket Equation to account for the forces of Drag and Gravity where both Fg and Fd are variable.

I have followed the method for deriving the Rocket Equation from MIT's website (here) and ended up with this:

dv = [Fd/mr - Fg/mr]dt - (vedmr)/mr

However, the math required to make this equation easily integratabtle is utterly baffling me because of the number of variables present in the function. Does anyone know how to simplify this?

Edited by bowtiesRcool
corrected equation (vedmr) should be divided by mr
Link to comment
Share on other sites

Tsiolkovsky's famous equation is only really meant for a rocket already in deep space. Where forces like drag and gravity are variable with velocity and altitude, there isn't really any one-step equation that can be used to solve for the altitude. You would have to do a numerical solution using time-steps. I guess you could try to do this with sufficiently steps and then try to curve fit the entire thing, but I don't think this is what you were going for, was it?

Link to comment
Share on other sites

Hi all,

I have been trying to work on an equation I can use to calculate the burnout height of a single stage rocket in vertical flight from the KSP Launchpad where the forces of Drag and Gravity are both variable.

I start with:

SUM(F) = Fthrust + Fdrag +Fgravity

where,

Fdrag = (0.004892*Cd * Po)mrv2 e-h/H

Fgravity = G*Mmrocket/(r+h)2

Fthrust = some constant

I attempted to modify the derivation of Tsiolkovski's Ideal Rocket Equation to account for the forces of Drag and Gravity where both Fg and Fd are variable.

I have followed the method for deriving the Rocket Equation from MIT's website (here) and ended up with this:

dv = [Fd/mr - Fg/mr]dt - (vedmr)/mr

However, the math required to make this equation easily integratabtle is utterly baffling me because of the number of variables present in the function. Does anyone know how to simplify this?

Okay, all the force equations look mostly accurate, but thrust is actually a function of height, and the mass of the vehicle is Minitial - t * mass flow rate.

I think you'll get something absurd like d2h/dt2 = (thrust as a function of altitude)/(Minitial-propellant mass flow rate * t) - (some constant) * (dh/dt)2 * e-h/H / (Minitial-propellant mass flow rate * t) - G*Mplanet/(r+h)2.

This is a nightmarish differential equation, including one term with an awkward function of h(t) and t, one term with a nonlinear function of h(t), and one term with an awkward function of h(t) and t multiplied by a nonlinear function of dh/dt. I kind of doubt this is analytically solvable.

Link to comment
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now
×
×
  • Create New...