Talk:PlanetPhysics/Single Stage Rocket Burnout Height

Original TeX Content from PlanetPhysics Archive
%%% This file is part of PlanetPhysics snapshot of 2011-09-01 %%% Primary Title: single stage rocket burnout height %%% Primary Category Code: 45.50.Dd %%% Filename: SingleStageRocketBurnoutHeight.tex %%% Version: 13 %%% Owner: bloftin %%% Author(s): bloftin %%% PlanetPhysics is released under the GNU Free Documentation License. %%% You should have received a file called fdl.txt along with this file. %%% If not, please write to gnu@gnu.org. \documentclass[12pt]{article} \pagestyle{empty} \setlength{\paperwidth}{8.5in} \setlength{\paperheight}{11in}

\setlength{\topmargin}{0.00in} \setlength{\headsep}{0.00in} \setlength{\headheight}{0.00in} \setlength{\evensidemargin}{0.00in} \setlength{\oddsidemargin}{0.00in} \setlength{\textwidth}{6.5in} \setlength{\textheight}{9.00in} \setlength{\voffset}{0.00in} \setlength{\hoffset}{0.00in} \setlength{\marginparwidth}{0.00in} \setlength{\marginparsep}{0.00in} \setlength{\parindent}{0.00in} \setlength{\parskip}{0.15in}

\usepackage{html}

% this is the default PlanetMath preamble. as your knowledge % of TeX increases, you will probably want to edit this, but % it should be fine as is for beginners.

% almost certainly you want these \usepackage{amssymb} \usepackage{amsmath} \usepackage{amsfonts}

% used for TeXing text within eps files %\usepackage{psfrag} % need this for including graphics (\includegraphics) %\usepackage{graphicx} % for neatly defining theorems and propositions %\usepackage{amsthm} % making logically defined graphics %\usepackage{xypic}

% there are many more packages, add them here as you need them

% define commands here

\begin{document}

Applying \htmladdnormallink{Newton's laws}{http://planetphysics.us/encyclopedia/Newtons3rdLaw.html} to rocket \htmladdnormallink{motion}{http://planetphysics.us/encyclopedia/CosmologicalConstant.html} is not only exciting, but also quite instructional. Problems involving rocket motion illustrate how to use Newton's 2nd law when \htmladdnormallink{mass}{http://planetphysics.us/encyclopedia/Mass.html} is not constant. Here we will look at how high a single stage rocket will go under the influence of gravity. To make the problem manageable, a few simplifying assumptions are made:

\begin{itemize} \item Motion in the y direction only \item Drag is neglected \item Constant burn rate for rocket \item \htmladdnormallink{force}{http://planetphysics.us/encyclopedia/Thrust.html} due to gravity is constant \item Rocket does not escape Earth's gravity \item Ideal rocket \item Lots of other miscellaneous terms \end{itemize}

To familiarize the reader with what is involved in this calculation, we will start with the answer and then derive the equation. So the max altitude the rocket will achieve is given by

\begin{equation} y_{max} = -\frac{g t^2_b}{2} + \frac{v_e m_f t_b}{m_i - m_f} ln \left (\frac{m_f}{m_i} \right ) + v_e t_b + \frac{v^2_b}{2 g} \end{equation}

Description of variables:

$t_b$: Time that the rocket is burning fuel (given) \\ $g$: \htmladdnormallink{acceleration}{http://planetphysics.us/encyclopedia/Acceleration.html} of Gravity (given) \\ $m_f$: Final mass of rocket after burn (given) \\ $m_i$: Initial mass of rocket (given) \\ $v_b$: \htmladdnormallink{velocity}{http://planetphysics.us/encyclopedia/Velocity.html} of rocket at burnout (calculated) \\ $v_e$: Exhaust velocity, velocity of the fuel as it is ejected out of the Rocket (given)\\

The problem is best approached by breaking it into two parts. First, we calculate the altitude that the rocket reaches when all its fuel is burned. After burnout, the rocket still climbs to a higher altitude until gravity finally brings its velocity to zero. Think of it like shooting a bullet into the sky, after the initial \htmladdnormallink{thrust}{http://planetphysics.us/encyclopedia/Thrust.html} of the gun, the bullet still goes higher (duh!). So the second calculation is to add the distance traveled after burnout.

The goal is to apply Newton's 2nd law, so let us start there. For our one dimensional case

\begin{equation} F = ma = -mg \end{equation}

Since we do not have constant mass throughout the rocket burn, we also have

\begin{equation} F = \frac{d}{dt}\left (m v \right) \end{equation}

Using the chain rule and setting (2) equal to (3)

$$-mg = \frac{dm}{dt} v_e + m \frac{dv}{dt} $$

multiply by dt \begin{equation} -mg dt = v_e dm + m dv \end{equation}

We still have three differentials, so we cannot directly integrate this equation. Since we are assuming a constant burn rate k and it is positive

$$ \frac{dm}{dt} = -k $$

so

$$ dt = -\frac{dm}{k} $$

plug this into (4) to get

$$ -\frac{mg}{k} dm = v_e dm + m dv $$

Separate variables to setup the integration

$$ \left ( \frac{mg}{k} - v_e \right ) dm = m dv $$

divide by m

$$ dv = \left ( \frac{g}{k} - \frac{v_e}{m} \right ) dm $$

integrate

$$ \int_{v_i}^{v} dv = \int_{m_i}^{m} \frac{g}{k} dm - \int_{m_i}^{m} \frac{v_e}{m} dm $$

The initial velocity of the rocket is zero, so carrying out the integration gives us the velocity at a given mass

$$ v(m) = \frac{g}{k} \left ( m - m_i \right) - v_e ln(m) + v_e ln(m_i) $$

simlify using properties of the log \htmladdnormallink{function}{http://planetphysics.us/encyclopedia/Bijective.html} \begin{equation} v(m) = \frac{g}{k} \left ( m - m_i \right) - v_e ln(\frac{m}{m_i}) \end{equation}

It is time to take care of the constant k.

$$ \frac{dm}{dt} = -k $$

rearrange

$$dm = -k dt $$

integrate

$$ \int_{m_i}^{m} dm = -k \int_{0}^{t} dt $$

$$ m - m_i = -kt $$

which gives

$$k = \frac{m_i - m}{t} $$

Note the sign, we have a positive k, since $m_i$ will always be bigger than m (the rocket is ejecting mass). Plug this into (5).

$$v(t) = \frac{g \left(m - m_i \right) t}{\left( m_i - m \right)} - v_e ln \left( \frac{m}{m_i} \right )$$

cancel terms to get

$$v(t) = -gt - v_e ln \left( \frac{m}{m_i} \right ) $$

While this equation gives us the velocity of the rocket at burnout, we also want altitude. As usual, integrate velocity to get \htmladdnormallink{position}{http://planetphysics.us/encyclopedia/Position.html} $$ \frac{dy}{dt} = v(t) $$

$$ \int_0^y dy = -g \int_0^{t_b} t dt - v_e ln \left( \frac{m}{m_i} \right) dt $$

To integrate the last term, we need to replace dt with dm, since m is a function of t

$$ dt = -\frac{dm}{k} $$

$$ \int_0^y dy = -g \int_0^{t_b} t dt + \frac{v_e}{k} \int_{m_i}^{m_f} ln \left( \frac{m}{m_i} \right) dm $$

All but the last term are simple. For the lazy, an integral table can be used to solve it.

$$ \int_0^y dy = -g \int_0^{t_b} t dt + \frac{v_e}{k} \int_{m_i}^{m_f} ln \left( \frac{m}{m_i} \right) dm $$

$$y = -\frac{g t_b^2}{2} + \frac{v_e}{k} \int_{m_i}^{m_f} ln \left( \frac{m}{m_i} \right) dm $$

Use typical integral substitution technique for the last term, so set

$$ w = \frac{m}{m_i} $$

$$ dw = \frac{dm}{m_i} $$

$$ dm = m_i dw $$

substituting this in leaves us with

$$ \frac{v_e m_i}{k} \int ln(w) dw $$

Integrating the logarithm function is done through integration by parts

\begin{equation} \int u dv = uv - \int v du \end{equation}

setting

$$ u = ln(w)\,\,\,\,\, dv = 1 dw $$

then differentiating and integrating yields

$$ du = \frac{1}{w} dw \,\,\,\,\, v = w $$

plugging these into (6) gives

$$ \int ln(w) dw = w ln(w) - \int \frac{w}{w} dw $$

which is equal to

$$ \int ln(w) dw = w ln(w) - w $$

so going back to the original integral and evaluating the limits we have

$$ y = -\frac{g t_b^2}{2} + \frac{v_e m_i}{k} \left ( \frac{m_f}{m_i} ln \left( \frac{m_f}{m_i} \right) - \frac{m_f}{m_i} - \frac{m_i}{m_i} ln \left( \frac{m_i}{m_i} \right) + \frac{m_i}{m_i} \right) $$

simplfying and using $ln(1) = 0$

$$ y = -\frac{g t_b^2}{2} + \frac{v_e m_i}{k} \left ( \frac{m_f}{m_i} ln \left( \frac{m_f}{m_i} \right) - \frac{m_f}{m_i} + 1 \right) $$

plug in k and expand

$$ y = -\frac{g t_b^2}{2} + \frac{v_e m_f t_b}{m_i - m_f}ln \left( \frac{m_f}{m_i} \right) + \frac{v_e \left( m_f - m_i \right) t_b}{m_i - m_f} $$

finally, except for the last term, we get equation (1)

\begin{equation} y = -\frac{g t_b^2}{2} + \frac{v_e m_f t_b}{m_i - m_f}ln \left( \frac{m_f}{m_i} \right) + v_e t_b \end{equation}

To find how much more altitude is gained after burnout, we note that there is no more thrust and rocket is under constant acceleration, g, so we go back to Newton's 2nd law

$$ F = ma $$

with a constant force due to gravity

$$ F = -mg = ma $$ $$ a = -g $$ $$ \frac{dv}{dt} = -g $$ $$ dv = -g dt $$

integrate to get velocity

$$ \int_{v_b}^0 dv = - \int_0^t g dt $$ $$ 0 - v_b = -gt $$ $$ v_b = gt $$

integrate again to get distance traveled

$$ \frac{dy}{dt} = gt $$ $$ dy = gt \, dt $$ $$ \int_0^y dy = \int_0^t gt \,dt $$

so the distance traveled after burnout is

$$ y = \frac{g t^2}{2} $$

combining this with the distance traveled during the burn (7), yields equation (1), the total distance rocket traveled

\begin{equation} y_{max} = -\frac{g t^2_b}{2} + \frac{v_e m_f t_b}{m_i - m_f} ln \left (\frac{m_f}{m_i} \right ) + v_e t_b + \frac{v^2_b}{2 g} \end{equation}

\subsection{References}

[1] Ellis, R., Gulick, D. "Calculus" Harcourt Brace Jovanovich, Inc., Orlando, FL, 1991.

[2] Etgen, G. "Calculus" John Wiley \& Sons, New York, 1999.

[3] Marion, J., Thornton, S. "Classical \htmladdnormallink{dynamics}{http://planetphysics.us/encyclopedia/NewtonianMechanics.html} of \htmladdnormallink{particles}{http://planetphysics.us/encyclopedia/Particle.html} and \htmladdnormallink{systems}{http://planetphysics.us/encyclopedia/SimilarityAndAnalogousSystemsDynamicAdjointnessAndTopologicalEquivalence.html}" Fourth Edition, Harcourt College Publishers, Fort Worth, 1995.

[4] Ketsdever, A. "Launch Vehicle Analysis", lecture notes, University of Colorado at Colorado Springs, Spring 2006.

\end{document}