Variable-Mass Systems and the Rocket Equation

A fully fueled Falcon 9 stands 70 metres tall and weighs about 550 tonnes — of which roughly 30 tonnes is the rocket and everything else is propellant. That lopsided ratio is not an engineering choice; it is a theorem. The velocity a rocket can gain is chained to the logarithm of its mass ratio, Δv=uln⁡(m0/mf)\Delta v = u\ln(m_0/m_f), and logarithms grow with agonizing slowness: to double your speed gain you must square your mass ratio. Every launch vehicle ever flown — its staging, its slenderness, its cost — is an argument with this one equation, and the equation always wins. The physics underneath is a trap and its resolution. Newton's F=maF = ma silently assumes a body of fixed mass; a rocket sheds mass continuously, and the "obvious" generalization F=d(mv)/dt=mv˙+vm˙F = \mathrm{d}(mv)/\mathrm{d}t = m\dot v + v\dot m is wrong if we consider only the rocket — it predicts different physics for observers in different inertial frames, which no real force can do. The cure is to retreat to what is always true: momentum conservation for a closed system. Track the rocket plus the parcel of exhaust it is about to eject, apply Fext=dp/dtF_{\mathrm{ext}} = \mathrm{d}p/\mathrm{d}t to that fixed collection of matter, and the correct equation of motion falls out — with the thrust u(−m˙)u(-\dot m) appearing automatically, built from the exhaust speed relative to the rocket. Remarkably, the man who first worked this out for spaceflight never launched so much as a firework. Konstantin Tsiolkovsky, a deaf schoolteacher in provincial Russia, published the rocket equation in 1903 — the same year the Wright brothers flew — and drew from it the blueprint of everything that followed: liquid propellants, staging, orbital stations. The same bookkeeping, with signs flipped, governs raindrops swallowing mist and chains piling onto tables. It is all one question: with what velocity does mass cross the system boundary?

Key takeaways
  • F=mv˙F = m\dot v (and its tempting cousin F=mv˙+vm˙F = m\dot v + v\dot m) fail for bodies that shed or absorb mass; Newton's laws govern closed systems of fixed particle content.
  • The fix: apply Fext=dp/dtF_{\mathrm{ext}} = \mathrm{d}p/\mathrm{d}t to the body plus the parcel of mass it exchanges during dt\mathrm{d}t, then take the limit — giving the Meshchersky equation mv˙=Fext+u(−m˙)m\dot v = F_{\mathrm{ext}} + u(-\dot m).
  • The thrust T=u(−m˙)T = u(-\dot m) is the momentum flux across the boundary: exhaust speed relative to the rocket times the burn rate. Only relative velocity can appear — anything else would break Galilean invariance.
  • Integrating in free space gives the Tsiolkovsky rocket equation Δv=uln⁡(m0/mf)\Delta v = u\ln(m_0/m_f): velocity gain depends only on exhaust speed and mass ratio, not on burn schedule.
  • The logarithm is tyranny: m0/mf=eΔv/um_0/m_f = e^{\Delta v/u} grows exponentially with the mission, which is why rockets are mostly propellant, why they stage, and why engineers worship specific impulse.
  • Accretion problems (raindrop in mist, chain on a table) use the same law with the incoming-mass velocity specified — when the swallowed mass is at rest, the equation collapses to d(mv)/dt=Fext\mathrm{d}(mv)/\mathrm{d}t = F_{\mathrm{ext}}.
  • Key idea: never write F=mv˙F = m\dot v for a variable-mass body — conserve momentum for the closed system of body plus exchanged parcel, and the thrust u(−m˙)u(-\dot m) emerges, leading to Δv=uln⁡(m0/mf)\Delta v = u\ln(m_0/m_f).
The core insight

Never write F=mv˙F = m\dot v for a body whose mass changes. Instead apply Fext=dp/dtF_{\mathrm{ext}} = \mathrm{d}p/\mathrm{d}t to a CLOSED system (body + exchanged parcel) over an interval dt\mathrm{d}t, then take the limit. The result is mv˙=Fext+u (−m˙)m\dot v = F_{\mathrm{ext}} + u\,(-\dot m), where uu is the velocity of the exhaust relative to the body and the term u(−m˙)u(-\dot m) is the thrust.

Where this lives today

The rocket equation is the profit-and-loss statement of the modern launch industry. Reusability lives or dies by it: a first stage flies a mission needing only ~4 km/s of the total budget, so reserving propellant to land it costs relatively little payload — while reusing an orbital stage means paying reentry and landing propellant at the top of the exponential, which is why boosters came back years before upper stages. SpaceX's Starship program is the current frontal assault: full reusability plus orbital refueling, which resets m0/mfm_0/m_f in space and effectively splits a deep-space mission into several tame logarithms — tanker flights through 2026 exist purely because equation (4) cannot be repealed, only restarted. The other escape route is uu itself. Electric propulsion abandons chemistry's 4.4 km/s ceiling: ion and Hall-effect thrusters expel xenon at 20–50 km/s (IspI_{\mathrm{sp}} of 2,000–5,000 s), shrinking the required mass ratio dramatically — at the price of thrusts measured in fractions of a newton, hence months-long spirals rather than minutes-long burns. NASA's Psyche is cruising to its metal asteroid on Hall thrusters; hundreds of commercial satellites now raise and hold their orbits the same way. And the equation's accretion twin quietly runs planetary science: protoplanets growing by sweeping up planetesimals, and stars accreting gas, obey the same boundary bookkeeping as the raindrop in the mist — momentum flux across a moving boundary, with the sign of m˙\dot m flipped.

The physical picture

Stand on a frictionless cart with a pile of bricks and start hurling them backward. Each throw sends a brick one way and you the other — Newton's third law, one discrete kick at a time. A rocket is this scene made continuous: the "bricks" are exhaust molecules, thrown backward at kilometres per second, millions of kilograms of them per minute at full throttle. Nothing is pushed against — not air, not launchpad. The rocket brings along everything it will ever push on, and that is precisely its problem.

Figure
Momentum bookkeeping over one instant — the rocket plus its about-to-be-ejected parcel form a closed system
time t — one closed system mass m, velocity v v time t + dt — same matter, two pieces mass m + dm velocity v + dv v + dv −dm v − u total momentum changes only by F_ext dt — equate, cancel, and the thrust u(−dm/dt) appears u is the exhaust speed RELATIVE to the rocket: the only velocity all inertial observers agree the engine controls

Two intuitions carry the whole subject. First: the kick depends on the relative throw speed. Whether your cart is already coasting at 1 m/s or 100 m/s, hurling a brick at 5 m/s relative to yourself delivers exactly the same recoil. That is why the equation of motion will contain the exhaust velocity uu relative to the rocket — never the rocket's own speed vv, which no engine can know or care about. Second: the late propellant pays for the early propellant. At liftoff the engine is accelerating not just the payload but hundreds of tonnes of propellant that will itself be burned later. Every kilogram of fuel you want to burn at the end must first be carried, at full price, through every earlier second of the burn. This compounding is exactly what a logarithm encodes, and it is why velocity gains get exponentially expensive: burn 63% of your mass and gain uu; to gain 2u2u you must burn 86%; for 3u3u, 95%. And the mirror image: a raindrop falling through stationary mist must accelerate every microgram it swallows up to its own speed. Swallowed mass acts as a brake — same bookkeeping, opposite sign of m˙\dot m — which is why a mist-gobbling drop falls with acceleration far below gg even in vacuum-clean air.

The equations of variable-mass motion

Everything descends from the one law that never fails, Newton's second law for a closed system:

Fext=dpsysdt(1)F_{\mathrm{ext}} = \frac{\mathrm{d}p_{\mathrm{sys}}}{\mathrm{d}t} \qquad (1)

where: psysp_{\mathrm{sys}} — total momentum of a fixed collection of matter. The word "closed" is the entire content: momentum bookkeeping is only guaranteed for a system that neither gains nor loses particles. A rocket alone is not such a system; a rocket plus its exhaust is. Applying (1) to the rocket-plus-parcel system (the derivation below) yields the equation of motion for the open body — the Meshchersky equation:

m dvdt=Fext  +  u(−dmdt)(2)m\,\frac{\mathrm{d}v}{\mathrm{d}t} = F_{\mathrm{ext}} \;+\; u\left(-\frac{\mathrm{d}m}{\mathrm{d}t}\right) \qquad (2)

where: m(t)m(t) — the body's instantaneous mass; uu — exhaust speed relative to the body; −m˙>0-\dot m > 0 — the mass-ejection rate. Read it as "mama = external force + thrust": the body responds like an ordinary particle of its current mass, driven by an extra force built entirely from the momentum flowing across its boundary. The second term deserves its own name and display:

T=u(−dmdt)=u m˙exhaust(3)T = u\left(-\frac{\mathrm{d}m}{\mathrm{d}t}\right) = u\,\dot m_{\mathrm{exhaust}} \qquad (3)

where: TT — the thrust. A real engine number: a Falcon 9 Merlin burns about 300 kg of propellant per second at an effective exhaust speed near u≈2.8 km/su \approx 2.8\ \mathrm{km/s}, giving T≈300×2800≈8.4×105 NT \approx 300 \times 2800 \approx 8.4\times10^5\ \mathrm{N} — 840 kN, the published sea-level figure. Engineers quote the same physics as specific impulse Isp=u/g0I_{\mathrm{sp}} = u/g_0, in seconds: chemical engines top out near Isp≈450I_{\mathrm{sp}} \approx 450 s (u≈4.4u \approx 4.4 km/s, hydrogen–oxygen), a hard chemical-energy ceiling that no clever plumbing moves. Integrate (2) in force-free flight and the burn schedule drops out entirely, leaving the most consequential logarithm in engineering — the Tsiolkovsky rocket equation:

Δv=u ln⁡ ⁣m0mf(4)\Delta v = u\,\ln\!\frac{m_0}{m_f} \qquad (4)

where: m0m_0, mfm_f — the masses at the start and end of the burn. Only two things matter: how fast you throw, and what fraction of yourself you throw. Inverted, m0/mf=eΔv/um_0/m_f = e^{\Delta v/u} — the exponential that runs the space industry:

Figure
The tyranny of the rocket equation — required propellant fraction grows exponentially with the mission
mission demand Δv / u propellant fraction 0 1 2 3 4 5 0% 100% Δv = u : burn 63% of yourself Δv = 2.3u : burn 90% Δv = 4.6u : burn 99% propellant fraction = 1 − exp(−Δv/u): doubling the mission squares the mass ratio — the logarithm always wins

For a vertical launch the external force does not vanish — gravity taxes the whole burn. Integrating (2) with Fext=−mgF_{\mathrm{ext}} = -mg (constant gg, constant uu) gives the burnout velocity

v=u ln⁡ ⁣m0mf  −  g tburn(5)v = u\,\ln\!\frac{m_0}{m_f} \;-\; g\,t_{\mathrm{burn}} \qquad (5)

where: tburnt_{\mathrm{burn}} — the duration of the burn. The Tsiolkovsky term depends only on the mass ratio; the gravity loss depends only on time. Burning the same propellant twice as fast halves the tax — which is why launch vehicles light enormous engines and leave the pad at maximum sensible thrust. Real first-stage gravity losses run 1–1.5 km/s (a 160-second burn costs gt≈1.6g t \approx 1.6 km/s). Accretion is the same law with the boundary crossing reversed. If the body sweeps up mass that is at rest in the lab frame (mist, a stationary chain link), the incoming matter arrives with relative velocity −v-v, i.e. u=−vu = -v in (2), and the equation collapses to the compact form

d(mv)dt=Fext(swallowed mass at rest in the lab)(6)\frac{\mathrm{d}(mv)}{\mathrm{d}t} = F_{\mathrm{ext}} \qquad \text{(swallowed mass at rest in the lab)} \qquad (6)

where: the product-rule term vm˙v\dot m now correctly appears — as the drag of accelerating swallowed mass up to speed, not as a fake force. This is the one situation where the naive d(mv)/dt\mathrm{d}(mv)/\mathrm{d}t happens to be right, and the reason is visible: the exchanged mass carries zero lab-frame momentum, so the open body's momentum is the closed system's.

Which formula answers which question

you wantusenotes
equation of motion, any mass exchangeMeshchersky, eq. (2)specify the relative velocity of exchanged mass
engine force from burn rateT=u(−m˙)T = u(-\dot m), eq. (3)u=g0Ispu = g_0 I_{\mathrm{sp}} converts from specific impulse
velocity gain, free spaceTsiolkovsky, eq. (4)burn schedule irrelevant
required mass ratio for a missionm0/mf=eΔv/um_0/m_f = e^{\Delta v/u}the exponential tyranny
vertical launch burnout speedeq. (5)gravity loss =g tburn= g\,t_{\mathrm{burn}}: burn fast
raindrop, chain, conveyor (mass at rest)eq. (6)the only safe home of d(mv)/dt\mathrm{d}(mv)/\mathrm{d}t
Derivation AThe Meshchersky equation from closed-system momentum

We are deriving the variable-mass equation of motion (2) — including the thrust term — from nothing but momentum conservation applied to a closed system. The strategy: at time tt, draw the system boundary around the rocket and the parcel of propellant it will eject during the next dt\mathrm{d}t; that collection of matter is fixed, so equation (1) applies to it without apology. Then compare the total momentum at tt and at t+dtt + \mathrm{d}t, cancel, and keep first-order terms. The step to watch is writing the ejected parcel's velocity: it must be v−uv - u — rocket velocity minus relative exhaust speed — and the discipline of writing all velocities in one inertial frame is exactly what the naive product-rule shortcut violates.

Step 1. Momentum at time tt. The parcel (mass −dm>0-\mathrm{d}m > 0) is still aboard, moving with the rocket:

p(t)=m vp(t) = m\,v

Step 2. Momentum at t+dtt + \mathrm{d}t. The rocket, now of mass m+dmm + \mathrm{d}m, moves at v+dvv + \mathrm{d}v; the ejected parcel moves at v−uv - u in the same inertial frame:

p(t+dt)=(m+dm)(v+dv)  +  (−dm)(v−u)p(t+\mathrm{d}t) = (m + \mathrm{d}m)(v + \mathrm{d}v) \;+\; (-\mathrm{d}m)(v - u)

Why v−uv - u and not −u-u: the exhaust leaves at uu relative to the rocket, and the rocket itself moves at vv; forgetting to add the frame velocity is the classic route to a wrong sign in the thrust.

Step 3. Expand and subtract. The mvmv terms cancel; the second-order piece dm dv\mathrm{d}m\,\mathrm{d}v is discarded in the limit:

dp=m dv+v dm−v dm+u dm=m dv+u dm\mathrm{d}p = m\,\mathrm{d}v + v\,\mathrm{d}m - v\,\mathrm{d}m + u\,\mathrm{d}m = m\,\mathrm{d}v + u\,\mathrm{d}m

Why the v dmv\,\mathrm{d}m pieces cancel: the momentum the rocket "loses" by shrinking is carried away by the parcel at the same base velocity vv — bulk velocity cannot generate force. What survives is only the relative-velocity part u dmu\,\mathrm{d}m: pure Galilean-invariant physics.

Step 4. Impose the closed-system law (1), dp=Fext dt\mathrm{d}p = F_{\mathrm{ext}}\,\mathrm{d}t, and divide by dt\mathrm{d}t:

m dvdt=Fext+u(−dmdt)m\,\frac{\mathrm{d}v}{\mathrm{d}t} = F_{\mathrm{ext}} + u\left(-\frac{\mathrm{d}m}{\mathrm{d}t}\right)

which is equation (2), with the thrust (3) appearing as the momentum flux across the boundary. For accretion, run the same steps with mass arriving at lab velocity ww: the general result replaces uu by the relative velocity of the exchanged mass, and setting w=0w = 0 gives equation (6). Galilean check: boosting all velocities by a constant VV shifts every term of Step 3 identically, leaving the equation unchanged, whereas the naive F=mv˙+vm˙F = m\dot v + v\dot m acquires a spurious Vm˙V\dot m.

Go deeper · Derivation B · 1 minThe Tsiolkovsky rocket equation and the gravity-loss correction
Derivation BThe Tsiolkovsky rocket equation and the gravity-loss correction

The goal is the integrated velocity budget: equation (4) in free space, and its vertical-launch correction (5). The strategy is separation of variables on the Meshchersky equation — velocity on one side, mass on the other — which works because the thrust term is a perfect differential in ln⁡m\ln m. Watch what does not appear in the result: neither the burn rate, nor the thrust level, nor the duration — in free space only the mass ratio survives, a stronger statement than it first looks. Gravity breaks exactly this indifference, which is where the launch-vehicle engineering lives.

Step 1. Free space: set Fext=0F_{\mathrm{ext}} = 0 in (2) and separate:

m dvdt=−u dmdt⟹dv=−u dmmm\,\frac{\mathrm{d}v}{\mathrm{d}t} = -u\,\frac{\mathrm{d}m}{\mathrm{d}t} \quad\Longrightarrow\quad \mathrm{d}v = -u\,\frac{\mathrm{d}m}{m}

Why the separation is legal: uu is constant for a given engine, and m(t)m(t) is monotonic during a burn, so mass can replace time as the integration variable — the schedule drops out here.

Step 2. Integrate from (v0,m0)(v_0, m_0) to (vf,mf)(v_f, m_f):

Δv=vf−v0=−u[ln⁡m]m0mf=u ln⁡ ⁣m0mf\Delta v = v_f - v_0 = -u\big[\ln m\big]_{m_0}^{m_f} = u\,\ln\!\frac{m_0}{m_f}

which is equation (4). Why a logarithm: each kilogram of ejected propellant gives momentum u dmu\,\mathrm{d}m to a rocket of current mass mm, so the incremental gain is u dm/mu\,\mathrm{d}m/m — a fixed fractional cost, and summing fractional costs is what logarithms do.

Step 3. Vertical launch: keep Fext=−mgF_{\mathrm{ext}} = -mg with constant gg:

dv=−u dmm−g dt⟹v=u ln⁡ ⁣m0mf−g tburn\mathrm{d}v = -u\,\frac{\mathrm{d}m}{m} - g\,\mathrm{d}t \quad\Longrightarrow\quad v = u\,\ln\!\frac{m_0}{m_f} - g\,t_{\mathrm{burn}}

which is equation (5). Why the two terms cannot be combined: the Tsiolkovsky term integrates over mass and sees only the ratio; the gravity term integrates over time and sees only the duration. The engineering moral is immediate — you cannot cheat the ratio, but you can starve the tax by burning fast.

Step 4. Sanity checks. Small burns: ln⁡(m0/mf)≈Δm/m\ln(m_0/m_f) \approx \Delta m/m gives Δv≈u Δm/m\Delta v \approx u\,\Delta m/m, the single-brick recoil. And Δv\Delta v exceeds uu itself once m0/mf>e≈2.72m_0/m_f > e \approx 2.72 — a rocket can outrun its own exhaust; late in the burn the exhaust it ejected early is actually moving forward in the ground frame. As a concrete check: a mass ratio of 10 at u=3.0u = 3.0 km/s gives Δv=3.0×ln⁡10=3.0×2.3026=6.91\Delta v = 3.0 \times \ln 10 = 3.0 \times 2.3026 = 6.91 km/s; the small-burn limit reproduces the impulse result; dimensions check throughout.

Limiting cases

CaseResultMeaning
m˙=0\dot m = 0eq. (2) → F=maF = mafixed mass recovers ordinary Newton
small burn, Δm≪m\Delta m \ll mΔv≈u Δm/m\Delta v \approx u\,\Delta m/mthe thrown-brick impulse limit
m0/mf=e≈2.72m_0/m_f = e \approx 2.72Δv=u\Delta v = uburn 63% of yourself to match your exhaust speed
m0/mf→∞m_0/m_f \to \inftyΔv→∞\Delta v \to \infty logarithmicallyno hard ceiling, but exponential cost — staging territory
swallowed mass at resteq. (2) → d(mv)/dt=Fext\mathrm{d}(mv)/\mathrm{d}t = F_{\mathrm{ext}}accretion: the product rule becomes legitimate
instantaneous burn, verticalgravity loss g tburn→0g\,t_{\mathrm{burn}} \to 0the ideal impulsive launch; real stages lose 1–1.5 km/s

Worked examples

1 · Thrust and liftoff of a Falcon-9-class booster (confidence builder)

Nine engines each burn propellant at m˙e=300 kg/s\dot m_e = 300\ \mathrm{kg/s} with effective exhaust speed u=2.8 km/su = 2.8\ \mathrm{km/s}. The fully fueled vehicle has mass m0=550m_0 = 550 tonnes. Find the thrust per engine, the total thrust, and the net acceleration at liftoff.

Step 1. Thrust per engine, equation (3): Te=u(−m˙)=2800×300=8.4×105 N=840 kNT_e = u(-\dot m) = 2800 \times 300 = 8.4\times10^5\ \mathrm{N} = \boxed{840\ \mathrm{kN}} — a good match to the published Merlin sea-level thrust.

Step 2. Total thrust: T=9×8.4×105=7.56×106 N≈7.6 MNT = 9 \times 8.4\times10^5 = 7.56\times10^6\ \mathrm{N} \approx 7.6\ \mathrm{MN}.

Step 3. Liftoff acceleration from equation (2) with gravity: a=T/m0−g=7.56×1065.5×105−9.81=13.7−9.81=3.9 m/s2≈0.4 ga = T/m_0 - g = \dfrac{7.56\times10^6}{5.5\times10^5} - 9.81 = 13.7 - 9.81 = \boxed{3.9\ \mathrm{m/s^2} \approx 0.4\,g}. Rockets leave the pad gently — the thrust barely wins at first.

Step 4. Watch the tyranny work in your favor: ninety seconds later the vehicle has burned 9×300×90=2439 \times 300 \times 90 = 243 tonnes. With m≈307m \approx 307 t the same thrust gives a=24.6−9.81≈14.8 m/s2a = 24.6 - 9.81 \approx 14.8\ \mathrm{m/s^2} — the acceleration nearly quadruples on constant thrust, purely because mm in equation (2) keeps falling.

2 · Why rockets stage: one stage versus two, same hardware (exam level)

A launcher has total mass 100 t, carries a 2 t payload, uses engines with u=3.0u = 3.0 km/s, and its tankage/structure always weighs 10% of the stage's loaded (wet) mass. Compare the total Δv\Delta v of (a) a single stage, and (b) two stages with wet masses 80 t and 18 t.

Step 1. Single stage. Wet vehicle without payload: 98 t, of which structure is 0.10×98=9.80.10 \times 98 = 9.8 t and propellant 88.288.2 t. So m0=100m_0 = 100 t and mf=100−88.2=11.8m_f = 100 - 88.2 = 11.8 t: Δv=3.0 ln⁡(100/11.8)=3.0×ln⁡8.475=3.0×2.137=6.41 km/s\Delta v = 3.0\,\ln(100/11.8) = 3.0 \times \ln 8.475 = 3.0 \times 2.137 = \boxed{6.41\ \mathrm{km/s}}. Short of orbit (need ≈ 9.3–9.5 km/s including losses) — and no amount of extra propellant fixes it fast, because the logarithm has flattened.

Step 2. Two stages. Stage 1: wet 80 t = 8 t structure + 72 t propellant. Stage 2: wet 18 t = 1.8 t structure + 16.2 t propellant. Payload 2 t rides on top; totals check: 80+18+2=10080 + 18 + 2 = 100 t. ✓

Step 3. Burn 1 (everything aboard): m0=100m_0 = 100 t, mf=100−72=28m_f = 100 - 72 = 28 t. Δv1=3.0 ln⁡(100/28)=3.0×ln⁡3.571=3.0×1.273=3.82 km/s\Delta v_1 = 3.0\,\ln(100/28) = 3.0 \times \ln 3.571 = 3.0 \times 1.273 = 3.82\ \mathrm{km/s}.

Step 4. Staging: drop the 8 t of dead stage-1 structure. New start mass 2020 t. Burn 2: mf=20−16.2=3.8m_f = 20 - 16.2 = 3.8 t. Δv2=3.0 ln⁡(20/3.8)=3.0×ln⁡5.263=3.0×1.661=4.98 km/s\Delta v_2 = 3.0\,\ln(20/3.8) = 3.0 \times \ln 5.263 = 3.0 \times 1.661 = 4.98\ \mathrm{km/s}.

Step 5. Total: Δv=3.82+4.98=8.80 km/s\Delta v = 3.82 + 4.98 = \boxed{8.80\ \mathrm{km/s}} — a 2.4 km/s gain over the single stage with identical engines, structural fraction, and payload. The entire profit comes from one move: refusing to keep accelerating 8 tonnes of empty stage-1 tankage. Staging is how engineering answers the logarithm: it cannot change equation (4), so it changes mfm_f mid-flight.

Origins · 1 minKonstantin Tsiolkovsky and how this was discovered
KT
Konstantin Tsiolkovsky
1857–1935 · Kaluga, Russia

Scarlet fever left Tsiolkovsky nearly deaf at ten; barred from schools, he educated himself in Moscow's libraries — living on black bread to buy books — and spent his career as a provincial schoolteacher in Kaluga, doing physics at his kitchen table. In 1903, the year the Wright brothers flew, he published "Exploration of Outer Space by Means of Reaction Devices," containing the rocket equation, the case for liquid hydrogen–oxygen propellants, and the argument that only staged rockets could reach orbital speed — the full blueprint of spaceflight, derived by a man who never launched anything. Dismissed for decades, he was canonized by the Soviet space program; Korolev, chief designer of Sputnik's R-7, called him the movement's grandfather. His epitaph is the sentence every space agency still quotes: "Earth is the cradle of humanity, but one cannot live in a cradle forever."

History
1813
— William Moore, at the Royal Military Academy, publishes the first derivation of the rocket equation, in a treatise on military rockets that is promptly forgotten.
1897
— Ivan Meshchersky's dissertation gives the general equation of motion for bodies of variable mass — equation (2) — as pure mechanics, rockets not required.
1903
— Konstantin Tsiolkovsky publishes the rocket equation with its spaceflight consequences: liquid propellants, mass ratios, staging.
1926
— Robert Goddard flies the first liquid-fueled rocket from a Massachusetts farm: 2.5 seconds, 12.5 metres up — the equation acquires hardware.
1969
— Saturn V reaches the Moon as three stages in series, each discarding dead mass exactly as the logarithm demands; its 2,900 t at liftoff deliver a 45 t payload to translunar injection.

Going further

Texts. Morin, Introduction to Classical Mechanics, ch. 5 (variable-mass problems and the rocket equation) · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., ch. 1 (systems of variable mass) · Landau & Lifshitz, Mechanics, Vol. 1, §7 problems · Taylor, Classical Mechanics, ch. 3 (rockets and the rocket equation) · Tsiolkovsky, K. E. (1903), "Exploration of Outer Space by Means of Reaction Devices". Key numeric anchors: a mass ratio of 10 at u=3u = 3 km/s gives Δv=3000⋅ln⁡10≈6908\Delta v = 3000 \cdot \ln 10 \approx 6908 m/s; standard gravity is g0=9.80665g_0 = 9.80665 m/s²; the worked-example figures (840 kN per engine and 0.4 g0.4\,g liftoff from u=2.8u = 2.8 km/s, m˙=300\dot m = 300 kg/s, m0=550m_0 = 550 t; the staging comparison 6.416.41 vs 3.82+4.98=8.803.82 + 4.98 = 8.80 km/s, from ln⁡8.475=2.137\ln 8.475 = 2.137, ln⁡3.571=1.273\ln 3.571 = 1.273, ln⁡5.263=1.661\ln 5.263 = 1.661) follow directly from the boxed formulas.

Builds on: Newton's Second Law · Newton's Third Law · Relative Motion and Galilean Transformations · Velocity and Acceleration · One-Dimensional Kinematics. Leads to: Conservation of Linear Momentum · Center of Mass and Many-Body Systems · Orbital Mechanics and Delta-v Budgets · The Tsiolkovsky Equation and Staging.

What does a rocket actually push against in empty space?

Its own propellant — and that is enough. The engine walls and nozzle push the combustion gas backward, accelerating it to speed uu relative to the ship; by Newton's third law the gas pushes the walls forward with equal force, and that contact force is the thrust, equation (3). No air, road, or launchpad is involved: the closed system of rocket-plus-exhaust conserves momentum exactly, the exhaust carrying backward exactly what the rocket gains forward. The intuition "there's nothing to push against in space" has it backwards — the rocket brings along the thing it pushes against, which is precisely why it must be mostly propellant and why the vacuum, far from being a problem, improves engine performance (no back-pressure on the nozzle).

Why exactly is F = d(mv)/dt wrong for a rocket? It looks like the most general form of Newton's law

Because it answers a different question. d(mv)/dt=F\mathrm{d}(mv)/\mathrm{d}t = F is correct for the momentum of a closed system — but a rocket alone is an open system whose boundary leaks momentum-carrying matter, and "the rocket's momentum" changes both because forces act and because matter leaves. The failed form charges the departing bulk momentum v dmv\,\mathrm{d}m to the force account. The clean diagnostic is Galilean invariance: boost to a frame moving at VV and the candidate force changes by Vm˙V\dot m, which no physical force may do — two observers cannot disagree about the reading of a force gauge. Done properly (Derivation A), the bulk terms cancel and only the invariant relative velocity uu survives. The product-rule form is redeemed only when the exchanged mass is at rest in your frame, equation (6) — accretion, never a rocket.

Where does the logarithm come from, intuitively?

From compounding. When the rocket's mass is mm, ejecting dm\mathrm{d}m of propellant at relative speed uu buys a velocity increment dv=u dm/m\mathrm{d}v = u\,\mathrm{d}m/m — a fractional cost: the same kilogram of propellant is worth more later, when the rocket is lighter, and less earlier, when it must also accelerate all the propellant still aboard. Summing fixed fractional gains is exactly what a logarithm does (the same reason compound interest is exponential). The brutal converse: the propellant needed for the last km/s must be carried through every earlier km/s, so mission demands multiply rather than add — m0/mf=eΔv/um_0/m_f = e^{\Delta v/u}, and a doubled mission squares the ratio. That one integral explains why rockets are balloons of fuel with a payload thimble on top.

What is specific impulse, and why is it quoted in seconds?

Specific impulse is exhaust velocity in disguise: Isp=u/g0I_{\mathrm{sp}} = u/g_0, with g0=9.80665 m/s2g_0 = 9.80665\ \mathrm{m/s^2} the standard gravity. The seconds arise from an old engineering convention — impulse delivered per unit weight (not mass) of propellant — which made the number identical in imperial and metric units and let 1950s engineers compare engines without unit fights. Physically, an IspI_{\mathrm{sp}} of 300 s means each newton of thrust costs a newton-weight of propellant every 300 seconds; multiply by g0g_0 and you recover u=2.94u = 2.94 km/s for use in equations (2)–(5). Rough map of the territory: solid boosters ~250 s, kerosene–oxygen ~300–350 s, hydrogen–oxygen ~450 s (the chemical ceiling), electric thrusters 2,000–5,000 s.

Why does staging help? The propellant is the same either way

Because the rocket equation charges you for every kilogram you carry, not just every kilogram you burn. A single-stage vehicle drags its entire tankage — built to hold all the propellant — all the way to burnout, so the final mass mfm_f in equation (4) is bloated by structure that long ago finished its job. Staging amputates dead mass mid-flight: Worked Example 2 shows identical hardware fractions and payload yielding 6.41 km/s in one stage but 8.80 km/s in two, a 2.4 km/s profit purely from dropping 8 tonnes of empty tank before the second burn. Each stage also gets an engine sized for its own weight rather than the liftoff weight. The cost is complexity and separation risk — but the logarithm is merciless enough that every orbital launcher in history has paid it.

Does the thrust depend on how fast the rocket is already moving?

No — and yes, delightfully. The force is u(−m˙)u(-\dot m), built from frame-invariant quantities; the engine cannot know its speed, so thrust is identical at liftoff and at 7 km/s. But the power delivered to the vehicle is P=TvP = Tv, which grows linearly with speed: the same newtons do more work per second when the ship moves faster. That is the seed of the Oberth effect — burns executed deep in a gravity well, where the ship moves fastest, extract more useful orbital energy from the same propellant, which is why interplanetary missions light their engines at closest approach. Energy books still balance: at high vv the exhaust is left moving slowly in the ground frame, so less of the chemical energy is wasted on exhaust kinetic energy and more lands on the ship.

Can the rocket end up moving faster than its own exhaust?

Yes — nothing in equation (4) caps Δv\Delta v at uu. Once the mass ratio exceeds e≈2.72e \approx 2.72 (burning about 63% of initial mass), Δv>u\Delta v > u. It sounds paradoxical only if you imagine the exhaust must always recede backward in the ground frame: late in a long burn, the rocket moves so fast that gas ejected at uu relative to the ship still drifts forward in the ground frame — just slower than the rocket. Momentum bookkeeping stays exact: the rocket's gain is balanced against exhaust laid down along the whole trail, the early (slow-frame) exhaust carrying large backward momentum. Multi-stage vehicles do this routinely: 8.8 km/s in Worked Example 2 from engines whose exhaust never exceeds 3 km/s.

In the raindrop problem, why does swallowing mist act like a drag force even with no air resistance?

Because every microgram of mist is at rest and must be yanked up to the drop's speed vv the instant it merges. That momentum, v dmv\,\mathrm{d}m, has to come from the drop, exactly as if a small backward force acted on it. Equation (6) makes this automatic: d(mv)/dt=mg\mathrm{d}(mv)/\mathrm{d}t = mg expands to mv˙=mg−vm˙m\dot v = mg - v\dot m, and the term −vm˙-v\dot m is the accretion drag — the mirror image of thrust. The effect is dramatic: for a drop growing in proportion to its size (a classic Morin problem), the acceleration settles to g/7g/7, not gg. Same law as the rocket, opposite sign of m˙\dot m, and the same golden question decides everything: with what velocity does mass cross the system boundary?

Before this

Where this leads