# 9.11: Rocket Propulsion

- Page ID
- 4675

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)- Describe the application of conservation of momentum when the mass changes with time, as well as the velocity
- Calculate the speed of a rocket in empty space, at some time, given initial conditions
- Calculate the speed of a rocket in Earth’s gravity field, at some time, given initial conditions

Now we deal with the case where the mass of an object is changing. We analyze the motion of a rocket, which changes its velocity (and hence its momentum) by ejecting burned fuel gases, thus causing it to accelerate in the opposite direction of the velocity of the ejected fuel (Figure \(\PageIndex{1}\)). Specifically: A fully fueled rocket ship in deep space has a total mass m_{0} (this mass includes the initial mass of the fuel). At some moment in time, the rocket has a velocity \(\vec{v}\) and mass m; this mass is a combination of the mass of the empty rocket and the mass of the remaining unburned fuel it contains. (We refer to m as the “instantaneous mass” and \(\vec{v}\) as the “instantaneous velocity.”) The rocket accelerates by burning the fuel it carries and ejecting the burned exhaust gases. If the burn rate of the fuel is constant, and the velocity at which the exhaust is ejected is also constant, what is the change of velocity of the rocket as a result of burning all of its fuel?

## Physical Analysis

Here’s a description of what happens, so that you get a feel for the physics involved.

- As the rocket engines operate, they are continuously ejecting burned fuel gases, which have both mass and velocity, and therefore some momentum. By conservation of momentum, the rocket’s momentum changes by this same amount (with the opposite sign). We will assume the burned fuel is being ejected at a constant rate, which means the rate of change of the rocket’s momentum is also constant. By Equation 9.4.17, this represents a constant force on the rocket.
- However, as time goes on, the mass of the rocket (which includes the mass of the remaining fuel) continuously decreases. Thus, even though the force on the rocket is constant, the resulting acceleration is not; it is continuously increasing.
- So, the total change of the rocket’s velocity will depend on the amount of mass of fuel that is burned, and that dependence is not linear.

The problem has the mass and velocity of the rocket changing; also, the total mass of ejected gases is changing. If we define our system to be the rocket + fuel, then this is a closed system (since the rocket is in deep space, there are no external forces acting on this system); as a result, momentum is conserved for this system. Thus, we can apply conservation of momentum to answer the question (Figure \(\PageIndex{2}\)).

At the same moment that the total instantaneous rocket mass is m (i.e., m is the mass of the rocket body plus the mass of the fuel at that point in time), we define the rocket’s instantaneous velocity to be \(\vec{v}\) = v \(\hat{i}\) (in the +x-direction); this velocity is measured relative to an inertial reference system (the Earth, for example). Thus, the initial momentum of the system is \(\vec{p}_{i}\) = mv \(\hat{i}\).

The rocket’s engines are burning fuel at a constant rate and ejecting the exhaust gases in the −x-direction. During an infinitesimal time interval dt, the engines eject a (positive) infinitesimal mass of gas dm_{g} at velocity \(\vec{u}\) = −u \(\hat{i}\); note that although the rocket velocity v \(\hat{i}\) is measured with respect to Earth, the exhaust gas velocity is measured with respect to the (moving) rocket. Measured with respect to the Earth, therefore, the exhaust gas has velocity (v − u) \(\hat{i}\).

As a consequence of the ejection of the fuel gas, the rocket’s mass decreases by dm_{g}, and its velocity increases by dv \(\hat{i}\). Therefore, including both the change for the rocket and the change for the exhaust gas, the final momentum of the system is

\[\begin{split} \vec{p}_{f} & = \vec{p}_{rocket} + \vec{p}_{gas} \\ & = (m - dm_{g})(v + dv) \hat{i} + dm_{g} (v - u) \hat{i} \ldotp \end{split}\]

Since all vectors are in the x-direction, we drop the vector notation. Applying conservation of momentum, we obtain

\[\begin{split} p_{i} & = p_{f} \\ mv & = (m - dm_{g})(v + dv) + dm_{g} (v - u) \\ mv & = mv + mdv - dm_{g} v - dm_{g} dv + dm_{g} v - dm_{g} u \\ mdv & = dm_{g} dv + dm_{g} u \ldotp \end{split}\]

Now, dm_{g }and dv are each very small; thus, their product dm_{g}dv is very, very small, much smaller than the other two terms in this expression. We neglect this term, therefore, and obtain:

\[mdv = dm_{g} u \ldotp\]

Our next step is to remember that, since dm_{g}represents an increase in the mass of ejected gases, it must also represent a decrease of mass of the rocket:

\[dm_{g} = - dm \ldotp\]

Replacing this, we have

\[mdv = -dmu\]

or

\[dv = -u \frac{dm}{m} \ldotp\]

Integrating from the initial mass m_{0} to the final mass m of the rocket gives us the result we are after:

\[\begin{split} \int_{v_{i}}^{v} dv & = -u \int_{m_{0}}^{m} \frac{1}{m} dm \\ v - v_{i} & = u \ln \left(\dfrac{m_{0}}{m}\right) \end{split}\]

and thus our final answer is

\[\Delta v = u \ln \left(\dfrac{m_{0}}{m}\right) \ldotp \label{9.38}\]

This result is called the **rocket equation**. It was originally derived by the Soviet physicist Konstantin Tsiolkovsky in 1897. It gives us the change of velocity that the rocket obtains from burning a mass of fuel that decreases the total rocket mass from m_{0} down to m. As expected, the relationship between \(\Delta\)v and the change of mass of the rocket is nonlinear.

In rocket problems, the most common questions are finding the change of velocity due to burning some amount of fuel for some amount of time; or to determine the acceleration that results from burning fuel.

- To determine the change of velocity, use the rocket equation Equation \ref{9.38}.
- To determine the acceleration, determine the force by using the impulse-momentum theorem, using the rocket equation to determine the change of velocity

A spacecraft is moving in gravity-free space along a straight path when its pilot decides to accelerate forward. He turns on the thrusters, and burned fuel is ejected at a constant rate of \(2.0 \times 10^2\, kg/s\), at a speed (relative to the rocket) of \(2.5 \times 10^2 \,m/s\). The initial mass of the spacecraft and its unburned fuel is \(2.0 \times 10^4\, kg\), and the thrusters are on for 30 s.

- What is the thrust (the force applied to the rocket by the ejected fuel) on the spacecraft?
- What is the spacecraft’s acceleration as a function of time?
- What are the spacecraft’s accelerations at t = 0, 15, 30, and 35 s?

**Strategy**

- The force on the spacecraft is equal to the rate of change of the momentum of the fuel.
- Knowing the force from part (a), we can use Newton’s second law to calculate the consequent acceleration. The key here is that, although the force applied to the spacecraft is constant (the fuel is being ejected at a constant rate), the mass of the spacecraft isn’t; thus, the acceleration caused by the force won’t be constant. We expect to get a function \(a(t)\), therefore.
- We’ll use the function we obtain in part (b), and just substitute the numbers given. Important: We expect that the acceleration will get larger as time goes on, since the mass being accelerated is continuously decreasing (fuel is being ejected from the rocket).

###### Solution

- The momentum of the ejected fuel gas is $$p = m_{g} v \ldotp$$The ejection velocity v = 2.5 x 10
^{2 }m/s is constant, and therefore the force is $$F = \frac{dp}{dt} = v \frac{dm_{g}}{dt} = -v \frac{dm}{dt} \ldotp$$Now, \(\frac{dm_{g}}{dt}\) is the rate of change of the mass of the fuel; the problem states that this is 2.0 x 10^{2}kg/s. Substituting, we get $$\begin{split} F & = v \frac{dm_{g}}{dt} \\ & = (2.5 \times 10^{2}\; m/s)(2.0 \times 10^{2}\; kg/s) \\ & = 5 \times 10^{4}\; N \ldotp \end{split}$$ - Above, we defined m to be the combined mass of the empty rocket plus however much unburned fuel it contained: m = m
_{R}+ m_{g}. From Newton’s second law, $$a = \frac{F}{m} = \frac{F}{m_{R} + m_{g}} \ldotp$$The force is constant and the empty rocket mass m_{R}is constant, but the fuel mass m_{g}is decreasing at a uniform rate; specifically: $$m_{g} = m_{g} (t) - m_{g_{0}} - \left(\dfrac{dm_{g}}{dt}\right)t \ldotp$$This gives us $$a(t) = \frac{F}{m_{g_{1}} - \left(\dfrac{dm_{g}}{dt}\right)t} = \frac{F}{M - \left(\dfrac{dm_{g}}{dt}\right)t} \ldotp$$Notice that, as expected, the acceleration is a function of time. Substituting the given numbers: $$a(t) = \frac{5 \times 10^{4}\; N}{(2.0 \times 10^{4}\; kg) - (2.0 \times 10^{2}\; kg/s)t} \ldotp$$ - At t = 0 s: $$a(0\; s) = \frac{5 \times 10^{4}\; N}{(2.0 \times 10^{4}\; kg) - (2.0 \times 10^{2}\; kg/s)(0\; s)} = 2.5\; m/s^{2} \ldotp$$

At t = 15 s, a(15 s) = 2.9 m/s^{2}.

At t = 30 s, a(30 s) = 3.6 m/s^{2}.

Acceleration is increasing, as we expected.

**Significance**

Notice that the acceleration is not constant; as a result, any dynamical quantities must be calculated either using integrals, or (more easily) conservation of total energy

What is the physical difference (or relationship) between \(\frac{dm}{dt}\) and \(\frac{dm_{g}}{dt}\) in this example?

## Rocket in a Gravitational Field

Let’s now analyze the velocity change of the rocket during the launch phase, from the surface of Earth. To keep the math manageable, we’ll restrict our attention to distances for which the acceleration caused by gravity can be treated as a constant g.

The analysis is similar, except that now there is an external force of \(\vec{F}\) = −mg \(\hat{j}\) acting on our system. This force applies an impulse d\(\vec{J}\) = \(\vec{F}\)dt = −mgdt \(\hat{j}\), which is equal to the change of momentum. This gives us

\[\begin{split} d \vec{p} & = d \vec{J} \\ \vec{p}_{f} - \vec{p}_{i} & = -mgdt\; \hat{j} \\ \big[ (m - dm_{g})(v + dv) + dm_{g}(v - u) - mv \big] \hat{j} & = -mgdt\; \hat{j} \end{split}\]

and so

\[mdv - dm_{g} u = -mgdt\]

where we have again neglected the term dm_{g}dv and dropped the vector notation. Next we replace dm_{g} with −dm:

\[\begin{split} mdv + dmu & = -mgdt \\ mdv & = -dmu - mgdt \ldotp \end{split}\]

Dividing through by \(m\) gives

\[dv = -u \frac{dm}{m} - gdt\]

and integrating, we have

\[\Delta v = u \ln \left(\dfrac{m_{0}}{m}\right) - g \Delta t \ldotp \label{9.39}\]

Unsurprisingly, the rocket’s velocity is affected by the (constant) acceleration of gravity.

Remember that \(\Delta\)t is the burn time of the fuel. Now, in the absence of gravity, Equation \ref{9.38} implies that it makes no difference how much time it takes to burn the entire mass of fuel; the change of velocity does not depend on \(\Delta\)t. However, in the presence of gravity, it matters a lot. The −g\(\Delta\)t term in Equation \ref{9.39} tells us that the **longer **the burn time is, the **smaller **the rocket’s change of velocity will be. This is the reason that the launch of a rocket is so spectacular at the first moment of liftoff: It’s essential to burn the fuel as quickly as possible, to get as large a \(\Delta\)v as possible.