09.01.08 · classical-mech / newtonian

Variable-Mass Systems: The Rocket Equation (Tsiolkovsky) and Momentum Transport

shipped3 tiersLean: none

Anchor (Master): Goldstein-Poole-Safko, Classical Mechanics 3e, Ch. 1; Ackeret, 'Zur Theorie der Raketen,' Helv. Phys. Acta 19 (1946)

Intuition Beginner

Stand on a frozen pond in ice skates. Hold a heavy backpack. Throw it forward as hard as you can. You slide backward. This is the rocket principle in its simplest form: to move in one direction, throw mass in the opposite direction.

A rocket does exactly this, continuously. It burns fuel and expels the hot exhaust gas backward at high speed. Each gram of gas pushed backward pushes the rocket forward. The rocket does not push against the air behind it -- it pushes against its own exhaust. This is why rockets work in the vacuum of space, where there is nothing to push against ^{[source pending]}.

Consider the momentum. Before the engine fires, the rocket plus its fuel sit still: total momentum is zero. After some fuel is expelled backward at high speed, the exhaust carries momentum in one direction and the rocket carries an equal and opposite momentum in the other. The total momentum remains zero. This is momentum conservation -- the same principle that governs billiard ball collisions and ice skaters pushing off each other ^{[source pending]}.

A heavier rocket needs more fuel. A faster exhaust needs less fuel. The relationship between these quantities is the Tsiolkovsky rocket equation, derived independently by Konstantin Tsiolkovsky in 1903. It says that the change in velocity a rocket can achieve depends on two things: how fast the exhaust leaves the engine (the exhaust velocity), and the ratio of the rocket's starting mass to its final mass (the mass ratio). Specifically:

Δ v = v_{e} ln \frac{m _{0}}{m _{f}},

where $v_{e}$ is the exhaust velocity, $m_{0}$ is the initial mass (rocket plus all fuel), and $m_{f}$ is the final mass (rocket after fuel is spent). The natural logarithm appears because each kilogram of fuel has to accelerate not just the rocket but all the remaining fuel still on board. The first gram of fuel expelled accelerates the entire full rocket by a tiny amount. The last gram accelerates a much lighter rocket by a larger amount ^{[source pending]}.

The logarithm is punishing. Doubling the fuel does not double the velocity change -- it only adds $v_{e} ln 2$ more. This is why single-stage rockets have such limited performance and why staging is essential for reaching orbit. A staged rocket discards empty tanks and engines once their fuel is spent, shedding dead weight. The upper stages then accelerate a much lighter vehicle. The total velocity change is the sum of each stage's contribution, and staging can achieve far more than any single stage could ^{[source pending]}.

Visual Beginner

Figure: Rocket momentum conservation at two instants. At time $t$ , a rocket of mass $m$ moves to the right at velocity $v$ . At time $t + d t$ , the rocket has mass $m - d m$ and velocity $v + d v$ (moving faster, lighter), while a small element of exhaust gas of mass $d m$ moves to the left at velocity $v_{e}$ relative to the rocket. The total momentum at each time is shown as horizontal arrows: at time $t$ , a single arrow $m v$ to the right; at time $t + d t$ , two arrows -- $(m - d m) (v + d v)$ to the right and $d m (v - v_{e})$ to the left. Momentum conservation requires the sum of the two arrows at $t + d t$ to equal the single arrow at $t$ .

Worked example Beginner

A small rocket has dry mass (structure plus payload) of 1000 kg and carries 4000 kg of fuel. The exhaust velocity is 3000 m/s. What velocity change can it achieve?

Step 1: Identify the masses.

Initial mass $m_{0} = 1000 + 4000 = 5000$ kg. Final mass $m_{f} = 1000$ kg (dry mass after all fuel is burned). Mass ratio $m_{0} / m_{f} = 5$ .

Step 2: Apply the Tsiolkovsky equation.

Δ v = v_{e} ln \frac{m _{0}}{m _{f}} = 3000 \times ln 5 = 3000 \times 1.609 = 4828 m/s .

Step 3: Check staging.

Suppose instead the rocket has two stages. Each stage carries half the fuel (2000 kg each) and has a dry mass of 500 kg (structure + engine). The payload is 500 kg.

Stage 1: initial mass $= 500 + 2000 + 500 + 2000 + 500 = 5500$ kg (everything). Final mass after burning $= 5500 - 2000 = 3500$ kg. Delta-v from stage 1: $3000 \times ln (5500/3500) = 3000 \times 0.452 = 1355$ m/s. Stage 1 structure (500 kg) is then discarded.

Stage 2: initial mass $= 500 + 2000 + 500 = 3000$ kg (stage 2 + payload). Final mass $= 3000 - 2000 = 1000$ kg. Delta-v from stage 2: $3000 \times ln (3000/1000) = 3000 \times 1.099 = 3296$ m/s.

Total delta-v: $1355 + 3296 = 4651$ m/s. This is less than the single-stage case because of the additional structural mass of the second-stage engine. But real staged rockets carry much more fuel relative to structure, and staging then provides large gains. The key insight is that staging discards structure, not just empty tanks ^{[source pending]}.

Check your understanding Beginner

Exercise (easy, multiple choice).

A rocket fires its engine in the vacuum of space. What provides the force that accelerates the rocket?

A. The exhaust pushing against the vacuum behind the rocket

B. The rocket pushing against its own expelled exhaust gas

C. The engine nozzle pushing against the spacecraft structure

D. Gravity from a nearby planet

Hint

Think about Newton's third law: the force on the rocket is the reaction to the force that accelerates the exhaust gas.

Answer

The rocket pushing against its own expelled exhaust gas -- option B, more precisely: the engine accelerates exhaust gas backward, and by Newton's third law the exhaust exerts an equal and opposite forward force on the engine. This reaction force is called thrust. No external medium is required, which is why rockets work in vacuum.

Formal definition Intermediate+

Momentum conservation for an open system

A rocket is an open system: mass flows out through the engine nozzle. Newton's second law in the form $F = m a$ applies to closed systems (constant mass). For an open system, the correct statement is

F_{ext} = \frac{d p}{d t},

where $p$ is the total momentum of the material currently inside the system boundary, and $F_{ext}$ includes the momentum flux carried by mass entering or leaving. For a rocket ejecting mass at rate $\overset{m}{˙} < 0$ (mass decreasing) with exhaust velocity $v_{e}$ relative to the rocket:

F_{ext} + \overset{m}{˙} v_{e} = m \frac{d v}{d t} .

The term $\overset{m}{˙} v_{e}$ is the momentum flux carried away by the exhaust. Rearranging:

m \frac{d v}{d t} = F_{ext} + ∣ \overset{m}{˙} ∣ v_{e} \hat{e}_{exhaust},

where $\hat{e}_{exhaust}$ points opposite to the exhaust direction. The quantity $∣ \overset{m}{˙} ∣ v_{e}$ is the thrust ^{[source pending]}.

The thrust equation

For a rocket engine, the thrust is

T = ∣ \overset{m}{˙} ∣ v_{e} = \frac{d m _{prop}}{d t} v_{e},

where $d m_{prop} / d t$ is the propellant mass flow rate and $v_{e}$ is the effective exhaust velocity. In engineering practice, thrust is often written in terms of the specific impulse $I_{sp}$ :

T = \frac{d m _{prop}}{d t} g_{0} I_{sp},

where $g_{0} = 9.80665$ m/s $^{2}$ is standard gravity and $I_{sp}$ has units of seconds. The specific impulse is a measure of engine efficiency: how much thrust is produced per unit weight flow of propellant. The relation is $v_{e} = g_{0} I_{sp}$ ^{[source pending]}.

The Tsiolkovsky rocket equation

For a rocket in free space (no external forces), the equation of motion is

m \frac{d v}{d t} = - v_{e} \frac{d m}{d t},

where $d m / d t < 0$ because mass is decreasing, and the minus sign gives thrust in the forward direction. Separating variables and integrating from initial state $(m_{0}, v_{0})$ to final state $(m_{f}, v_{f})$ :

\int_{v_{0}}^{v_{f}} d v = - v_{e} \int_{m_{0}}^{m_{f}} \frac{d m}{m},

Δ v \equiv v_{f} - v_{0} = v_{e} ln \frac{m _{0}}{m _{f}} .

This is the Tsiolkovsky rocket equation (1903). It relates the achievable velocity change to the exhaust velocity and the mass ratio $R = m_{0} / m_{f}$ . The mass ratio is the single most important performance parameter of a rocket vehicle ^{[source pending]}.

Key features:

The velocity change scales linearly with exhaust velocity. Better engines (higher $v_{e}$ ) directly translate to more delta-v.
The velocity change scales logarithmically with mass ratio. Diminishing returns: to double delta-v, the mass ratio must be squared.
The result is independent of the burn rate. Whether the fuel is expelled quickly or slowly, the total delta-v is the same (in the absence of external forces). Only the mass ratio and exhaust velocity matter.
For a mass ratio $R = e \approx 2.718$ , the delta-v equals the exhaust velocity. For $R = e^{2} \approx 7.389$ , the delta-v is twice the exhaust velocity.

Staging

A single-stage rocket's mass ratio is limited by structural constraints: the tanks, engines, and interstage structure all contribute to the dry mass $m_{f}$ . Typical single-stage mass ratios are 8--15 for well-designed vehicles ^{[source pending]}.

Staging overcomes this limitation by dividing the vehicle into stages, each with its own engines and propellant. When a stage is expended, its structure is discarded. An $n$ -stage rocket achieves total delta-v

Δ v_{total} = i = 1 \sum n v_{e, i} ln \frac{m _{0, i}}{m _{f, i}},

where $m_{0, i}$ is the mass of the entire vehicle at ignition of stage $i$ , and $m_{f, i}$ is the mass after stage $i$ has burned its propellant (just before stage $i$ is jettisoned).

For identical stages (same $v_{e}$ and same structural fraction), the optimal distribution of propellant among stages is obtained when the mass ratios of each stage are equal. This follows from the convexity of the logarithm: equal mass ratios minimise the total propellant required for a given $Δ v$ , or equivalently maximise $Δ v$ for a given total propellant mass ^{[source pending]}.

Gravity losses

When a rocket launches from a planetary surface, gravity opposes the motion. The equation of motion becomes

m \frac{d v}{d t} = T - m g,

and the gravity loss is the velocity penalty incurred by thrusting against gravity:

Δ v_{gravity loss} = \int_{0}^{t_{burn}} g d t .

For a vertical launch from Earth's surface with a burn time of 150 seconds, the gravity loss is approximately $150 \times 9.8 \approx 1470$ m/s. This is a significant fraction of the total delta-v budget for reaching orbit (about 7800 m/s for low Earth orbit). Gravity losses motivate launching with high thrust-to-weight ratio and pitching over quickly to gain horizontal velocity, since only the vertical component of thrust fights gravity ^{[source pending]}.

Key derivation Intermediate+

Key result (Tsiolkovsky equation from momentum conservation). Consider a rocket in free space at time $t$ with mass $m$ and velocity $v$ . In time $d t$ , it expels propellant of mass $d m > 0$ at exhaust velocity $v_{e}$ relative to the rocket. By momentum conservation:

Δ v = v_{e} ln \frac{m _{0}}{m _{f}} .

Derivation (momentum-flux method). At time $t$ , the total momentum is $p (t) = m v$ . At time $t + d t$ , the rocket has mass $m - d m$ and velocity $v + d v$ , while the expelled propellant has mass $d m$ and velocity $v - v_{e}$ (in the inertial frame). Momentum conservation requires

m v = (m - d m) (v + d v) + d m (v - v_{e}) .

Expanding: $m v = m v + m d v - d m \cdot v - d m \cdot d v + d m \cdot v - d m \cdot v_{e}$ . Cancel $m v$ on both sides, and note $d m \cdot d v$ is second-order small:

0 = m d v - d m \cdot v_{e} .

Since the rocket loses mass, $d m > 0$ and $d m = - d m_{rocket}$ where $d m_{rocket}$ is the change in rocket mass. Thus $m d v = v_{e} d m$ , and since $d m_{rocket} = - d m < 0$ :

d v = - v_{e} \frac{d m}{m},

where now $m$ is the (decreasing) rocket mass. Integrating from $m_{0}$ to $m_{f}$ :

\int_{v_{0}}^{v_{f}} d v = - v_{e} \int_{m_{0}}^{m_{f}} \frac{d m}{m} = v_{e} ln \frac{m _{0}}{m _{f}},

which gives $Δ v = v_{e} ln (m_{0} / m_{f})$ ^{[source pending]}.

Alternative derivation: thrust from Newton's second law

Write Newton's second law for the rocket as a variable-mass body. The rate of change of the rocket's momentum is

\frac{d p}{d t} = \frac{d ( m v )}{d t} = m \frac{d v}{d t} + v \frac{d m}{d t} .

But the rocket is losing mass at rate $\overset{m}{˙} = d m / d t < 0$ , and the expelled mass carries away momentum at rate $- \overset{m}{˙} (v - v_{e})$ . The net force (zero in free space) equals the rate of change of the total momentum of rocket plus exhaust:

0 = \frac{d p _{rocket}}{d t} + \frac{d p _{exhaust}}{d t} = (m \frac{d v}{d t} + v \frac{d m}{d t}) + (- \overset{m}{˙}) (v - v_{e}) .

Simplifying: $m d v / d t + v \overset{m}{˙} - \overset{m}{˙} v + \overset{m}{˙} v_{e} = 0$ , giving $m d v / d t = - \overset{m}{˙} v_{e} = ∣ \overset{m}{˙} ∣ v_{e}$ . This identifies the thrust as $T = ∣ \overset{m}{˙} ∣ v_{e}$ and reproduces the Tsiolkovsky equation upon integration ^{[source pending]}.

Worked example: two-stage rocket optimization

A two-stage rocket carries a 500 kg payload. Both stages use the same propellant with $v_{e} = 3500$ m/s. Each stage has a structural mass equal to 10% of its propellant mass. The total propellant budget is 20,000 kg. How should the propellant be divided between stages to maximise delta-v?

Let $m_{p 1}$ and $m_{p 2}$ be the propellant masses of stages 1 and 2, with $m_{p 1} + m_{p 2} = 20, 000$ kg. Structural masses are $0.1 m_{p 1}$ and $0.1 m_{p 2}$ .

Stage 2 (upper): initial mass $= m_{p 2} + 0.1 m_{p 2} + 500 = 1.1 m_{p 2} + 500$ . Final mass $= 0.1 m_{p 2} + 500$ . Mass ratio: $R_{2} = (1.1 m_{p 2} + 500) / (0.1 m_{p 2} + 500)$ .

Stage 1 (lower): initial mass $= (1.1 m_{p 2} + 500) + m_{p 1} + 0.1 m_{p 1} = 1.1 m_{p 2} + 500 + 1.1 m_{p 1}$ . Final mass $= 1.1 m_{p 2} + 500 + 0.1 m_{p 1}$ . Mass ratio: $R_{1} = (1.1 m_{p 2} + 500 + 1.1 m_{p 1}) / (1.1 m_{p 2} + 500 + 0.1 m_{p 1})$ .

Total delta-v: $Δ v = 3500 (ln R_{1} + ln R_{2})$ .

For the equal mass-ratio optimum (which applies when stages are identical), numerical evaluation with $m_{p 1} = m_{p 2} = 10, 000$ kg gives $R_{1} = (11000 + 500 + 11000) / (11000 + 500 + 1000) = 22500/12500 = 1.8$ and $R_{2} = 11500/1500 = 7.667$ . The mass ratios are unequal, so the even split is not optimal. The optimum shifts more propellant to stage 1 to balance the mass ratios ^{[source pending]}.

Worked example: vertical launch with gravity

A sounding rocket launches vertically from rest with initial mass 500 kg, propellant mass 400 kg, exhaust velocity 2500 m/s, and constant burn rate of 10 kg/s. What is the velocity at burnout, including gravity losses?

Burn time: $t_{b} = 400/10 = 40$ s. Thrust: $T = 10 \times 2500 = 25, 000$ N. Initial weight: $W_{0} = 500 \times 9.8 = 4900$ N. Thrust-to-weight ratio: $T / W_{0} = 5.1$ .

The equation of motion is $m d v / d t = T - m g$ where $m (t) = 500 - 10 t$ . Substituting $d m = - 10 d t$ :

d v = \frac{T}{m} d t - g d t = - \frac{T}{10} \frac{d m}{m} - g d t .

The velocity change from thrust alone is $Δ v_{thrust} = - (T /10) \int_{500}^{100} d m / m = - (2500) ln (100/500) = 2500 \times 1.609 = 4024$ m/s. The gravity loss is $Δ v_{grav} = - g t_{b} = - 9.8 \times 40 = - 392$ m/s. Net burnout velocity: $4024 - 392 = 3632$ m/s. The gravity loss is about 10% of the ideal delta-v ^{[source pending]}.

Bridge. The Tsiolkovsky equation governs the velocity budgets of every space mission. The delta-v it computes determines whether a trajectory is achievable, and the trajectory itself is governed by the central-force problem 09.03.01. The staging optimisation is a constrained optimisation that connects to the calculus of variations 09.02.01. Gravity losses introduce the interaction between thrust and gravitational fields, which is the domain of orbit transfer theory and the Oberth effect discussed in the Master section.

Exercises Intermediate+

Exercise 4 (medium, symbolic).

A two-stage rocket has payload mass $m_{L}$ , structural coefficient $ϵ$ (structural mass equals $ϵ$ times propellant mass for each stage), exhaust velocity $v_{e}$ for both stages, and total propellant mass $m_{P}$ divided equally between the stages. Find the total delta-v in terms of these parameters.

Hint

Write the mass ratio for each stage. The second stage's initial mass is the payload plus the second stage's structure and propellant.

Answer

Each stage has propellant mass $m_{P} /2$ . Structural mass per stage: $ϵ m_{P} /2$ .

Stage 2 initial mass: $m_{0, 2} = m_{P} /2 + ϵ m_{P} /2 + m_{L} = (1 + ϵ) m_{P} /2 + m_{L}$ . Final mass: $m_{f, 2} = ϵ m_{P} /2 + m_{L}$ . Mass ratio: $R_{2} = [(1 + ϵ) m_{P} /2 + m_{L}] / [ϵ m_{P} /2 + m_{L}]$ .

Stage 1 initial mass: $m_{0, 1} = m_{P} + ϵ m_{P} + m_{L} = (1 + ϵ) m_{P} + m_{L}$ . Final mass (stage 1 structure remains at burnout): $m_{f, 1} = ϵ m_{P} /2 + m_{P} /2 + m_{0, 2} = (1 + ϵ) m_{P} + m_{L} - m_{P} /2 = m_{P} /2 + ϵ m_{P} + m_{L}$ . Wait -- let me redo this cleanly.

Stage 1 initial: everything $= m_{P} + ϵ m_{P} + m_{L}$ . Stage 1 final (stage 1 propellant burned, structure jettisoned): $= m_{0, 2} = (1 + ϵ) m_{P} /2 + m_{L}$ . So $R_{1} = [(1 + ϵ) m_{P} + m_{L}] / [(1 + ϵ) m_{P} /2 + m_{L}]$ .

$Δ v = v_{e} [ln R_{1} + ln R_{2}] = v_{e} ln (R_{1} R_{2})$ .

Exercise 5 (medium, symbolic).

Derive the Tsiolkovsky equation by considering the total momentum of rocket plus all exhaust ever expelled, showing that it equals the initial momentum (zero for a rocket starting from rest).

Hint

At any time $t$ , the rocket has mass $m (t)$ and velocity $v (t)$ . The exhaust expelled before time $t$ has total mass $m_{0} - m (t)$ and a range of velocities (earlier exhaust moves slower). Integrate the exhaust momentum.

Answer

At time $t$ , rocket momentum is $m (t) v (t)$ . The exhaust element expelled at time $t^{'}$ (when rocket velocity was $v (t^{'})$ ) has mass $- d m$ and velocity $v (t^{'}) - v_{e}$ . Total exhaust momentum is $\int_{0}^{t} [v (t^{'}) - v_{e}] (- d m) = \int_{m_{0}}^{m (t)} [v (t^{'}) - v_{e}] d m^{'}$ where $d m^{'}$ runs over the rocket mass at expulsion time.

Total momentum conservation: $m (t) v (t) + \int_{m (t)}^{m_{0}} [v (t^{'}) - v_{e}] d m^{'} = m_{0} \cdot 0 = 0$ (starting from rest). Differentiating with respect to $t$ :

$m \overset{v}{˙} + \overset{m}{˙} v + [v - v_{e}] (- \overset{m}{˙}) = 0$ , giving $m \overset{v}{˙} = - \overset{m}{˙} v_{e}$ , which integrates to the Tsiolkovsky equation.

Exercise 6 (medium, symbolic).

A rocket launches vertically in a uniform gravitational field $g$ with constant thrust $T$ and constant mass flow rate. Show that the burnout velocity is $v_{f} = v_{e} ln (m_{0} / m_{f}) - g t_{b}$ where $t_{b}$ is the burn time. Verify that a higher thrust (shorter burn time) reduces gravity losses.

Hint

Integrate the equation of motion $m d v / d t = T - m g$ with $m (t) = m_{0} - (T / v_{e}) t$ .

Answer

With $m (t) = m_{0} - (T / v_{e}) t$ and $d m = - (T / v_{e}) d t$ :

$d v = (T / m) d t - g d t = - v_{e} d m / m - g d t$ .

Integrating: $Δ v = v_{e} ln (m_{0} / m_{f}) - g t_{b}$ .

Since $t_{b} = (m_{0} - m_{f}) v_{e} / T$ , the gravity loss is $g (m_{0} - m_{f}) v_{e} / T$ . Higher thrust $T$ reduces $t_{b}$ and hence the gravity loss. This is why launch vehicles use high-thrust engines even though the Tsiolkovsky delta-v is independent of thrust in free space.

Exercise 7 (hard, symbolic).

A conveyor belt moves horizontally at speed $v$ . Sand is dumped onto it at rate $σ$ (kg/s) from a stationary hopper. Find the force required to keep the belt moving at constant speed, and the power input. Where does the energy go?

Hint

This is a variable-mass problem where the system gains mass. The incoming sand has zero horizontal momentum and must be accelerated to $v$ .

Answer

The belt-plus-sand system gains mass at rate $σ$ . The incoming sand has zero horizontal velocity. The momentum balance is $F = d (m v) / d t = v d m / d t = σ v$ . This force accelerates the new sand from rest to speed $v$ .

The power input is $P = F v = σ v^{2}$ . But the rate of increase of kinetic energy of the sand is $d (\frac{1}{2} m v^{2}) / d t = \frac{1}{2} σ v^{2}$ . Half the input power becomes kinetic energy; the other half is dissipated as heat in the inelastic collision between the sand and the belt. This is a characteristic feature of inelastic mass accretion: the energy efficiency is exactly 50%.

Exercise 8 (hard, symbolic).

An $n$ -stage rocket has identical stages, each with structural mass fraction $ϵ$ (structural mass = $ϵ$ times propellant mass) and exhaust velocity $v_{e}$ . The payload mass is $m_{L}$ . Show that for a fixed total propellant mass $m_{P}$ , the total delta-v is maximised when each stage has the same mass ratio, and find the optimal delta-v.

Hint

Let $R_{i}$ be the mass ratio of stage $i$ . Express the total delta-v as $v_{e} \sum_{i} ln R_{i}$ subject to the constraint that the product $R_{1} R_{2} \dots R_{n}$ is fixed by the total propellant and structural mass budget. Use the AM-GM inequality or Lagrange multipliers.

Answer

For identical stages, the mass ratio of each stage is $R = (1 + ϵ) / ϵ$ in the limit of zero payload. More precisely, with payload $m_{L}$ , the overall mass ratio is $R_{total} = [(1 + ϵ)^{n} / ϵ^{n}] \cdot f (m_{L} / m_{P})$ where $f$ encodes the payload fraction.

For fixed $m_{P}$ and $n$ , the constraint is that $R_{1} R_{2} \dots R_{n} = C$ (a constant determined by $m_{P}$ , $ϵ$ , and $m_{L}$ ). By the AM-GM inequality, $\sum ln R_{i}$ is maximised when all $R_{i}$ are equal: $R_{i} = C^{1/ n}$ for all $i$ .

Then $Δ v = n v_{e} ln (C^{1/ n}) = v_{e} ln C$ , which is independent of $n$ for fixed $C$ . But increasing $n$ increases $C$ because more staging allows more propellant to be used effectively. The practical optimum is limited by the increasing complexity and structural overhead of additional stages. Real rockets rarely exceed three stages.

Lean formalization Intermediate+

The Tsiolkovsky equation is a clean target for formalisation. It requires: a time-dependent mass function $m (t)$ satisfying $d m / d t = - α$ for constant burn rate $α > 0$ ; a velocity function $v (t)$ satisfying $m d v / d t = α v_{e}$ ; and the theorem that $v_{f} - v_{0} = v_{e} ln (m_{0} / m_{f})$ . Mathlib's ODE existence theorems and logarithm properties suffice for the final statement, but the coupling of variable mass with Newton's second law is not represented. The momentum-flux formulation of $F = d p / d t$ for open systems would require modelling system boundaries with mass flow, which is a general framework not currently in Mathlib.

Advanced results Master

The relativistic rocket equation

When the rocket's velocity becomes a significant fraction of the speed of light, the Tsiolkovsky equation must be modified. The relativistic equation of motion for a rocket in free space is

\frac{d}{d t} (γ m v) = γ_{e} v_{e} \frac{d m}{d t},

where $γ = (1 - v^{2} / c^{2})^{- 1/2}$ is the Lorentz factor of the rocket and $γ_{e} = (1 - v_{e}^{2} / c^{2})^{- 1/2}$ is the Lorentz factor corresponding to the exhaust velocity in the rocket's rest frame. Integration yields the relativistic Tsiolkovsky equation:

\frac{m _{0}}{m _{f}} = (\frac{1 + v _{f} / c}{1 - v _{f} / c})^{c / (2 v_{e})} .

For $v_{f} ≪ c$ , this reduces to the nonrelativistic form $Δ v = v_{e} ln (m_{0} / m_{f})$ . The relativistic form shows that arbitrarily high final velocities are in principle achievable with sufficient mass ratio, but the mass ratio required grows exponentially with $v_{f} / c$ ^{[source pending]}.

For a photon rocket (exhaust velocity $v_{e} = c$ ), the equation simplifies to

\frac{m _{0}}{m _{f}} = \frac{1 + v _{f} / c}{1 - v _{f} / c},

or equivalently $v_{f} / c = (m_{0} - m_{f}) / (m_{0} + m_{f})$ . Even converting the entire mass to photons ( $m_{f} = 0$ ) gives $v_{f} = c$ , as expected. A photon rocket is the most efficient possible reaction drive in terms of exhaust velocity, but the thrust per unit power is extremely low: $T = P / c$ where $P$ is the radiated power ^{[source pending]}.

The Oberth effect

The Oberth effect, described by Hermann Oberth in 1927, is the counterintuitive result that a rocket manoeuvre produces more kinetic energy change when performed at higher velocity. Specifically, a burn of delta-v at velocity $v$ changes the kinetic energy by

Δ T = \frac{1}{2} m (v + Δ v)^{2} - \frac{1}{2} m v^{2} = m v Δ v + \frac{1}{2} m (Δ v)^{2} .

The first term $m v Δ v$ dominates for $v ≫ Δ v$ . At high speed, a given propellant expenditure converts more propellant chemical energy into vehicle kinetic energy because the exhaust is left at lower velocity (and hence lower kinetic energy) in the inertial frame ^{[source pending]}.

The practical consequence: orbital manoeuvres are most efficient at periapsis (closest approach), where the spacecraft moves fastest. A retrograde burn at apoapsis (slowest point) is least efficient. This is why transfer orbits like the Hohmann transfer 09.03.01 use brief burns at periapsis and apoapsis rather than continuous thrust.

The Oberth effect is not a violation of energy conservation. The total energy change (vehicle plus exhaust) is the same regardless of when the burn occurs. What changes is the partition of energy between vehicle and exhaust. At high vehicle velocity, the exhaust is expelled at lower velocity in the inertial frame and carries away less kinetic energy, leaving more for the vehicle ^{[source pending]}.

Effective exhaust velocity and nozzle physics

The ideal exhaust velocity for a chemical rocket is given by

v_{e} = \frac{2 γ}{γ - 1} \frac{R T _{c}}{M} [1 - (\frac{p _{e}}{p _{c}})^{(γ - 1) / γ}],

where $γ$ is the ratio of specific heats, $R$ is the gas constant, $T_{c}$ is the combustion chamber temperature, $M$ is the molecular weight of the exhaust, $p_{e}$ is the nozzle exit pressure, and $p_{c}$ is the chamber pressure. Lower molecular weight (hydrogen is ideal) and higher chamber temperature increase $v_{e}$ . The square-root dependence means that doubling exhaust velocity requires quadrupling the chamber temperature ^{[source pending]}.

The effective exhaust velocity also includes a pressure-thrust term. The total thrust is

T = \overset{m}{˙} v_{e} + (p_{e} - p_{a}) A_{e},

where $p_{a}$ is the ambient pressure and $A_{e}$ is the nozzle exit area. The effective exhaust velocity is $v_{e, eff} = T / \overset{m}{˙}$ , which exceeds the actual gas exit velocity when $p_{e} > p_{a}$ (underexpanded nozzle) and is less than it when $p_{e} < p_{a}$ (overexpanded nozzle). Optimum expansion ( $p_{e} = p_{a}$ ) maximises thrust for a given mass flow rate ^{[source pending]}.

Variable-mass systems beyond rockets

The momentum-flux formulation applies to any system with mass flow. Several classical examples:

Leaky freight car. A freight car of initial mass $m_{0}$ rolls without friction at initial velocity $v_{0}$ . It leaks sand at rate $σ$ . Since the sand leaves with the same horizontal velocity as the car (it falls through the floor), it carries away momentum but exerts no horizontal force on the car. The car maintains constant velocity $v_{0}$ while its mass decreases. The kinetic energy of the car decreases, but the total kinetic energy of car plus leaked sand is conserved (each grain retains the velocity it had when it left) ^{[source pending]}.

Raindrop falling through cloud. A raindrop falls through a cloud and accumulates mass by condensation at rate $d m / d t = k m^{2/3}$ (proportional to surface area). The equation of motion is

\frac{d}{d t} (m v) = m g,

since the accreted water droplets have negligible initial velocity. This gives $m d v / d t + v d m / d t = m g$ , a different structure from the rocket equation because the accreted mass brings zero momentum.

Chain falling onto a scale. A chain of linear density $λ$ is held with its lower end just touching a scale, then released. As the chain falls, links pile up on the scale. The force on the scale has two contributions: the weight of the accumulated chain ( $λ g y$ where $y$ is the fallen length) and the impulse per unit time from the arriving links ( $λ v^{2}$ where $v = 2 g y$ is the impact velocity). The total force is $3 λ g y$ -- three times the weight of the chain on the scale ^{[source pending]}.

Historical context

Konstantin Eduardovich Tsiolkovsky (1857--1935) derived the rocket equation in his 1903 paper "Exploration of the World Space with Reactive Devices" (Issledovanie mirovykh prostranstv reaktivnymi priborami), published in the Russian journal Nauchnoe Obozrenie (Scientific Review), No. 5. Tsiolkovsky was a self-taught schoolteacher in Kaluga, Russia, who had never attended university. His paper passed almost unnoticed in the scientific community due to its publication in Russian and in an obscure journal ^{[source pending]}.

Tsiolkovsky's 1903 paper also proposed liquid hydrogen and liquid oxygen as rocket propellants, described the benefits of staging, and calculated the velocities needed to reach orbit and escape Earth's gravity. He continued publishing on rocketry throughout his life, including a 1911 revision of his 1903 paper and a 1929 work on multistage rockets. His famous epitaph, composed by himself, reads: "Mankind will not remain tied to Earth forever."

Robert Hutchings Goddard (1882--1945) independently derived similar results in his 1919 Smithsonian publication "A Method of Reaching Extreme Altitudes." Goddard went further than Tsiolkovsky in one crucial respect: he built and flew the world's first liquid-fueled rocket on March 16, 1926, in Auburn, Massachusetts. The rocket burned for 2.5 seconds, reaching an altitude of 12.5 meters and a speed of about 97 km/h. Goddard's subsequent rockets achieved altitudes of over 2 km by 1935 ^{[source pending]}.

Hermann Julius Oberth (1894--1989) published Die Rakete zu den Planetenraumen ("The Rocket into Planetary Space") in 1923, independently deriving the rocket equation and proposing liquid-fueled rockets, staging, and space stations. Oberth's work, unlike Tsiolkovsky's, reached a wide audience and inspired the Verein fur Raumschiffahrt (German Rocket Society) and a generation of rocket engineers. His 1929 follow-up Wege zur Raumschiffahrt ("Ways to Spaceflight") described the Oberth effect and proposed electric and ion propulsion ^{[source pending]}.

The three men -- Tsiolkovsky, Goddard, and Oberth -- are regarded as the founding fathers of rocketry. Each worked independently, each derived the rocket equation from momentum conservation, and each recognised that liquid propellants and staging were essential for reaching space. Their theoretical work laid the foundation for Wernher von Braun's V-2 rocket (first ballistic missile, 1944), Sergei Korolev's R-7 launcher (first orbital launch, Sputnik, 1957), and the Saturn V (first human lunar landing, Apollo 11, 1969) ^{[source pending]}.

Modern rocketry extends the Tsiolkovsky framework to ion engines ( $v_{e} \sim 30, 000$ m/s but very low thrust), nuclear thermal rockets ( $v_{e} \sim 9000$ m/s), and proposed nuclear pulse propulsion (Project Orion, theoretical $v_{e} \sim 100, 000$ m/s). The Tsiolkovsky equation remains the fundamental performance equation for all reaction drives: it sets the mass-ratio requirements for any mission profile, from satellite station-keeping to interstellar probes ^{[source pending]}.

Synthesis. The Tsiolkovsky rocket equation is the bridge between Newtonian momentum conservation and practical spaceflight. Its logarithmic dependence on mass ratio is the fundamental constraint that drives the design of every launch vehicle and spacecraft. Staging partially circumvents this constraint by discarding structure, and the Oberth effect recovers propellant efficiency by exploiting the nonlinear relationship between velocity and kinetic energy. The relativistic extension shows how the equation deforms near the speed of light, maintaining the same physical structure (exhaust velocity and mass ratio) but with relativistic kinematics. Variable-mass dynamics appears beyond rockets in conveyor belts, raindrops, and chain problems, all unified by the momentum-flux formulation of Newton's second law. The central insight -- that a system can propel itself by ejecting part of its own mass -- is a direct consequence of momentum conservation and is independent of any external medium.

Connections Master

09.01.01 Kinematics provides the velocity and acceleration descriptions that the thrust equation and gravity-loss calculations require; the rocket's position and velocity as functions of time are the kinematic outputs of the variable-mass equations of motion.
09.01.02 Newton's second law $F = d p / d t$ is the foundation from which the thrust equation and Tsiolkovsky equation are derived; this unit extends Newton's law to open systems with mass flux.
09.02.01 The action principle provides an alternative formulation of mechanics; the Oberth effect, derived here from kinetic energy considerations, has a natural variational interpretation in terms of the work done by thrust along a trajectory.
09.03.01 Central-force motion and orbital mechanics determine the gravitational environment in which rockets operate; the Tsiolkovsky delta-v budget computed here determines whether orbit transfer manoeuvres are achievable, and the Oberth effect optimises where along the orbit to burn.
09.06.01 Special relativity modifies the Tsiolkovsky equation for velocities approaching $c$ ; the relativistic rocket equation derived in the Master section is the correct high-velocity extension that reduces to the classical form in the $v ≪ c$ limit.
09.07.01 Continuum mechanics treats propellant as a fluid with its own equations of motion; the effective exhaust velocity and nozzle physics discussed here are coarse-grained descriptions of the compressible gas dynamics that govern real engine performance.

Historical & philosophical context Master

The Tsiolkovsky equation is one of the rare instances where a foundational physics result was derived not by a professional scientist but by a self-taught outsider. Konstantin Tsiolkovsky lost most of his hearing to scarlet fever at age ten and never received formal higher education. He taught himself mathematics, physics, and chemistry from books in his father's library. His 1903 derivation of the rocket equation predates practical rocketry by decades -- the Wright brothers' first powered flight occurred the same year -- and his vision of liquid-fueled, staged rockets was not realised until Goddard's 1926 flight and the German V-2's 1944 combat deployment ^{[source pending]}.

The equation itself is a direct consequence of momentum conservation, which was understood since Newton's Principia (1687). The question of why it took more than two centuries to derive the rocket equation is partly a question of motivation: before the twentieth century, there was no practical reason to consider rocket propulsion as a means of reaching space. Tsiolkovsky, Goddard, and Oberth were driven not by military or commercial applications but by the vision of interplanetary travel.

Robert Goddard's contribution was engineering rather than theoretical. His 1919 paper derived the rocket equation independently and proposed the liquid-fueled rocket, but his greater achievement was building one. Goddard pioneered the use of gyroscopic stabilisation, regenerative cooling of the combustion chamber, and the de Laval nozzle for exhaust acceleration. His secrecy and isolation -- he worked alone with a small team, publishing minimal details -- meant that his innovations had limited influence on the broader rocketry community. The New York Times famously editorialised in 1920 that Goddard lacked "the knowledge ladled out daily in high schools" because he believed rockets could operate in vacuum. The newspaper issued a retraction the day after Apollo 11 launched in 1969 ^{[source pending]}.

Hermann Oberth's 1923 doctoral dissertation, rejected by the University of Heidelberg as "utopian," became the foundational text of practical astronautics. Unlike Tsiolkovsky and Goddard, Oberth engaged with a broad community: his students and collaborators included Wernher von Braun, who would lead the development of the V-2 and later the Saturn V. The V-2, the first ballistic missile and the first vehicle to reach the edge of space (84 km altitude), demonstrated that the physics Tsiolkovsky had calculated on paper could be engineered into working hardware. The V-2 also demonstrated the moral ambiguity of technological progress: it killed an estimated 9,000 civilians in London and Antwerp, while more slave labourers died producing it than were killed by its attacks ^{[source pending]}.

A philosophical point: the rocket equation illustrates a deep asymmetry in momentum conservation. A rocket can accelerate in any direction by expelling mass, but it can never stop expelling mass without stopping acceleration. There is no "momentum battery" -- the only way to change momentum is to transfer it to expelled propellant. This is why every spacecraft manoeuvre consumes irreplaceable propellant, and why mission design centres on delta-v budgets. The logarithmic dependence means that the propellant cost of a manoeuvre grows exponentially with its velocity requirement, making high-energy missions (like interstellar probes) extraordinarily expensive in mass-ratio terms. This constraint drives the search for propellantless propulsion concepts (solar sails, magnetic sails, gravitational assists) that circumvent the Tsiolkovsky equation by transferring momentum to or from an external body rather than carrying all reaction mass onboard ^{[source pending]}.

Bibliography Master

Tsiolkovsky, K. E., "Issledovanie mirovykh prostranstv reaktivnymi priborami" [Exploration of the World Space with Reactive Devices], Nauchnoe Obozrenie [Scientific Review], No. 5 (1903). English translation in NASA Technical Translation F-15 (1966).
Goddard, R. H., "A Method of Reaching Extreme Altitudes," Smithsonian Miscellaneous Collections 71(2) (1919).
Oberth, H., Die Rakete zu den Planetenraumen [The Rocket into Planetary Space] (R. Oldenbourg, Munich, 1923).
Oberth, H., Wege zur Raumschiffahrt [Ways to Spaceflight] (R. Oldenbourg, Munich, 1929).
Ackeret, J., "Zur Theorie der Raketen," Helvetia Physica Acta 19, 103--112 (1946). [Relativistic rocket equation.]
Taylor, J. R., Classical Mechanics (University Science Books, 2005), Ch. 3.
Goldstein, H., Poole, C. P. & Safko, J. L., Classical Mechanics, 3rd ed. (Pearson, 2002), Ch. 1.
Sutton, G. P. & Biblarz, O., Rocket Propulsion Elements, 7th ed. (Wiley, 2001).
Turner, M. J. L., Rocket and Spacecraft Propulsion: Principles, Practice and New Developments (Springer, 2005).
Clark, J. D., Ignition! An Informal History of Liquid Rocket Propellants (Rutgers University Press, 1972).

Prerequisites

09.01.01
09.01.02

Used in

09.03.01
09.01.09

Tier anchors

beginner: Taylor, Classical Mechanics (2005), Ch. 3 Momentum and CM motion
intermediate: Taylor, Classical Mechanics (2005), Ch. 3; Goldstein-Poole-Safko, Classical Mechanics 3e, Ch. 1
master: Goldstein-Poole-Safko, Classical Mechanics 3e, Ch. 1; Ackeret, 'Zur Theorie der Raketen,' Helv. Phys. Acta 19 (1946)

References

Taylor — Classical Mechanics (University Science Books, 2005) · Ch. 3 Momentum (rocket as variable-mass system, momentum conservation derivation of Tsiolkovsky equation, thrust).
Goldstein, Poole & Safko — Classical Mechanics, 3rd ed. (Pearson, 2002) · Ch. 1 Survey of the Elementary Principles. Variable-mass systems, momentum flux, thrust equation.
Tsiolkovsky — 'Exploration of the World Space with Reactive Devices' (1903, in Russian) · Original derivation of the rocket equation relating exhaust velocity, mass ratio, and final velocity. Published in Nauchnoe Obozrenie (Scientific Review), No. 5.
Sutton — Rocket Propulsion Elements, 7th ed. (Wiley, 2001) · Ch. 2 Fundamentals of Rocket Propulsion. Thrust equation, effective exhaust velocity, specific impulse, staging, gravity losses.

Estimated time

beginner: 15m
intermediate: 45m
master: 80m