Rocket Propulsion

Variable-Mass Momentum

The mass of a rocket decreases while its engine operates. A calculation must state which material belongs to the system at each instant. The usual rocket system contains the vehicle, its payload, structure, and propellant that has not yet left the vehicle. Exhaust already beyond the nozzle lies outside that system. Its momentum still matters because crossing the boundary transports momentum away.

Let $m (t)$ be the instantaneous mass of the rocket system and let $v (t)$ be its velocity in an inertial frame. During a short interval $d t$ , a positive amount of exhaust mass $d m_{ex}$ leaves the vehicle. The system mass change is therefore

d m = - d m_{ex}, \overset{m}{˙} < 0.

The sign in $\overset{m}{˙}$ carries physical information in every propulsion calculation. Replacing it with a positive burn rate too early is a common source of reversed thrust terms. A separate symbol $R = - \overset{m}{˙} > 0$ is convenient when the engine expels propellant at a positive rate. The two conventions are equivalent only when their signs are carried through every equation.

The rocket system is the vehicle plus unspent propellant at one instant. Once exhaust crosses the dashed boundary it leaves the system, carrying backward momentum as a flux through that boundary.

The boundary is a modelling convention placed around selected material. A larger closed system could include the rocket and all gas expelled during a specified interval. Momentum conservation is then applied to that closed collection. The open-system form used below is shorter because it accounts for the leaving gas through one relative-speed term. Both methods produce the same result when they describe the same interval.

Mass gain uses the same bookkeeping with the opposite sign. Sand falling onto a cart, rain accumulating in a moving wagon, and a spacecraft collecting dust have $\overset{m}{˙} > 0$ . Incoming material may have a laboratory velocity different from the vehicle velocity. It must be accelerated or decelerated after capture, so even a force-free collector can change speed. The rocket is the mass-loss case: expelled material leaves with a speed relative to the vehicle and gives the remaining vehicle a forward momentum change.

An open variable-mass boundary requires more than the substitution $F = d (m v) / d t$ . The momentum of material presently inside the boundary changes as material enters or leaves. Newton's second law applies directly to a closed set of matter; an open-system balance includes both external force and momentum transport.

Momentum balance over a short interval.

Consider one-dimensional motion first. At time $t$ , the rocket has mass $m$ and laboratory speed $v$ . In the interval $d t$ , its mass changes by $d m < 0$ and its speed changes by $d v$ . The expelled mass is $- d m > 0$ . Let $u_{e} > 0$ denote the exhaust speed relative to the rocket, directed backward. The exhaust laboratory speed is therefore $v - u_{e}$ when the rocket moves in the positive direction.

Take the closed system as the rocket plus the small amount of propellant that will leave during the interval. Its initial momentum is $m v$ . At the end, the rocket has momentum $(m + d m) (v + d v)$ and the exhausted material has momentum $(- d m) (v - u_{e})$ . External impulse is $F_{ext} d t$ . Thus,

F_{ext} d t = (m + d m) (v + d v) + (- d m) (v - u_{e}) - m v .

Expanding and discarding the product $d m d v$ , which vanishes faster than either first-order change as $d t \to 0$ , gives

F_{ext} = m \frac{d v}{d t} - u_{e} \frac{d m}{d t} .

The last term is positive for a running engine because $d m / d t$ is negative. With $R = - \overset{m}{˙}$ , the equation becomes

m \overset{v}{˙} = F_{ext} + R u_{e} .

Before gravity, aerodynamic force, or other external forces are specified, this is the one-dimensional rocket equation. The quantity $R u_{e}$ has units of force. It is the momentum per second carried backward by the exhaust relative to the rocket, and its equal-and-opposite contribution to the rocket is the thrust in this ideal model.

The relative speed belongs in the mass-transfer term. In a frame moving at a constant speed $V$ , the rocket and exhaust laboratory speeds both change by $V$ , while their difference remains $u_{e}$ . The thrust prediction therefore has the same value in every inertial frame. Inserting an exhaust speed measured relative to the ground without transforming it to a relative speed produces a frame-dependent answer and signals an inconsistent momentum balance.

A shift to a frame moving at speed V changes both the rocket and exhaust laboratory speeds, but their separation — the relative speed — is unchanged. The momentum-transfer term uses that relative speed, so thrust is frame-independent.

Vector form resolves nozzle-direction changes during flight. If $u_{e}$ is the exhaust velocity relative to the rocket, with its direction opposite the outgoing jet, then

m \frac{d v}{d t} = F_{ext} - u_{e} \frac{d m}{d t} .

A rear-facing nozzle has backward $u_{e}$ and $d m / d t < 0$ , so the thrust term points forward. A gimballed nozzle rotates $u_{e}$ and hence rotates the thrust vector. The equation predicts acceleration from the instantaneous mass, external force, exhaust direction, and mass-flow rate; it does not require the path to be straight.

Thrust and Ideal Velocity Change

A steady engine with constant $u_{e}$ and $R$ produces the ideal thrust

T = R u_{e} .

Thrust separates into two engine characteristics. The burn rate $R$ measures how many kilograms leave per second. The exhaust speed $u_{e}$ measures the momentum change per kilogram. Doubling either quantity doubles ideal thrust, provided the other remains fixed. The two changes have different consequences for mission duration and propellant use. A large burn rate consumes the available propellant quickly; a large exhaust speed raises the velocity change available from a given mass ratio.

At launch, with upward positive and with drag temporarily neglected,

a = \frac{T}{m} - g = \frac{R u _{e}}{m} - g .

The thrust-to-weight ratio is

\frac{T}{m g} .

It must exceed one for upward acceleration from a level launch surface. A ratio of one gives zero initial acceleration in the ideal vertical model. As propellant is expelled, $m$ decreases. If $T$ stays approximately constant, $T / m$ rises and the acceleration rises during the burn. Actual launch vehicles often reduce thrust late in a stage to limit structural load and crew acceleration.

With constant thrust the thrust acceleration T/m rises as the vehicle sheds mass. Subtracting the nearly constant gravitational acceleration gives the net vertical acceleration, which stays positive once thrust exceeds weight.

Thrust is a force, whereas impulse is thrust integrated over time. For a constant thrust interval $Δ t$ , the impulse magnitude is $T Δ t$ . It changes the momentum of the vehicle, but the change in vehicle speed is not generally $T Δ t / m$ with one fixed mass because $m$ changes during the interval. The logarithmic mass dependence derived below is the finite-burn result.

An engine test can determine thrust without measuring the rocket trajectory. A calibrated load cell between an engine mount and a rigid stand records the axial reaction. A separate measurement of propellant mass change over the same time interval gives $R$ . The ratio $T / R$ estimates $u_{e}$ in the ideal momentum model. Stable test conditions matter: a transient ignition trace, changing chamber pressure, or a drifting mass scale should not be averaged without recording the time range used for each quantity.

Nozzle pressure can add a correction when exhaust leaves at a pressure different from the surrounding pressure. The introductory momentum model treats the effective exhaust speed as the quantity that makes $T = R u_{e}$ match the measured thrust. In that form, the pressure contribution and the jet momentum contribution are included together. A detailed nozzle calculation resolves them separately; the variable-mass derivation still uses the total axial momentum transport.

The ideal rocket equation.

In a direction with no external force, substitute $\overset{m}{˙} = - R$ into the momentum balance and separate variables:

m \frac{d v}{d t} = R u_{e}, \frac{d v}{d t} = - u_{e} \frac{1}{m} \frac{d m}{d t} .

Integration from a state with mass $m_{i}$ and speed $v_{i}$ to a state with mass $m_{f}$ and speed $v_{f}$ gives

v_{f} - v_{i} = u_{e} ln (\frac{m _{i}}{m _{f}}) .

The speed increment

Δ v = u_{e} ln (\frac{m _{i}}{m _{f}})

is called the ideal rocket velocity change. It is a scalar only when thrust stays parallel to the velocity-change direction. In a manoeuvre with rotated thrust, the same incremental relation is integrated vectorially along the changing nozzle direction.

Ideal velocity change grows as

u_{e} ln (m_{i} / m_{f})

. Each doubling of the mass ratio adds the same increment

u_{e} ln 2

, so the equal steps from ratio 2 to 4 show why large velocity changes are expensive in propellant.

Expelled reaction mass supplies the momentum partner for a rocket. The combined rocket-plus-exhaust momentum remains conserved in a force-free environment, so a rocket accelerates in vacuum without surrounding air. An aircraft propeller instead accelerates surrounding air and depends on that external working fluid.

The logarithm has an immediate design consequence. A mass ratio of $2$ produces $u_{e} ln 2$ . Raising the ratio from $2$ to $4$ adds another $u_{e} ln 2$ , even though the added propellant is much larger in absolute mass. Large mission velocity changes therefore become expensive in propellant when $u_{e}$ is limited. Structural mass, tanks, engines, and payload all remain in $m_{f}$ and reduce the achievable ratio.

The ideal equation uses a constant effective exhaust speed. A real engine can throttle, change mixture ratio, or operate through an atmosphere whose pressure changes with altitude. For $u_{e} = u_{e} (t)$ , retain the differential form

d v = - u_{e} (t) \frac{d m}{m}

and integrate using the measured or specified function. Segmenting a burn into short intervals with nearly constant $u_{e}$ provides a numerical approximation. Each interval contributes $u_{e, k} ln (m_{k} / m_{k + 1})$ to the ideal axial velocity change. The same segmentation handles a sequence of different engines or stages.

Flight Dynamics

For vertical motion near Earth's surface, take upward as positive. Let the engine operate with constant burn rate $R$ and constant effective exhaust speed $u_{e}$ . With the mass law

m (t) = m_{0} - R t,

and with a constant approximation for $g$ , the equation of motion without drag is

\frac{d v}{d t} = \frac{R u _{e}}{m _{0} - R t} - g .

Integration from $t = 0$ gives

v (t) - v_{0} = u_{e} ln (\frac{m _{0}}{m _{0} - R t}) - g t .

The first contribution is the ideal propulsion gain. The second is the gravity loss: the accumulated external impulse of gravity over the burn duration. A higher-thrust vehicle reaches a given speed in less time and consequently suffers a smaller gravity loss, even when two vehicles have the same ideal $Δ v$ .

Free-body diagram of a rocket in vertical ascent: upward thrust, downward weight, and downward aerodynamic drag along the axis. The exhaust plume lies outside the vehicle because thrust already represents its momentum effect.

The initial launch condition follows directly from the free-body equation. At the pad, a support force can keep the vehicle at rest even when thrust is below weight. After release, a vertical ascent requires

T - D - m g > 0,

where $D$ is the downward drag magnitude. Near zero speed, drag may be small, so $T > m_{0} g$ is the usual first check. It tests positive vertical acceleration after release. An engine can have sub-weight thrust while mounted to a stand, or while accelerating a vehicle horizontally with a supporting track.

The pad reaction changes discontinuously at lift-off. Before release, upward support $N$ satisfies $N + T - m g = 0$ . As $T$ approaches $m g$ from below, $N$ approaches zero. An ideal pad cannot pull downward; once the force balance would require $N < 0$ , contact is lost and the free-flight equation replaces the support equation.

Integrating position adds a second integral. For the constant- $R$ , constant- $u_{e}$ , constant- $g$ model,

y (t) - y_{0} = v_{0} t + u_{e} [t - (\frac{m _{0} - R t}{R}) ln (\frac{m _{0}}{m _{0} - R t})] - \frac{1}{2} g t^{2} .

Compare a numerical time-step model against this expression, but use a changing gravity vector, atmosphere, and steering law for a launch trajectory. It should not be extended past propellant depletion. At $t_{b} = m_{p} / R$ , where $m_{p}$ is the usable propellant mass, the burn ends; the subsequent trajectory is a different dynamical segment with $R = 0$ .

Atmospheric force and trajectory modelling.

At low altitude, aerodynamic force can be comparable with thrust or weight. A simple axial drag model uses

D = \frac{1}{2} ρ C_{D} A v^{2},

where $ρ$ is local air density, $C_{D}$ is a drag coefficient appropriate to the shape and flow regime, and $A$ is a reference area. For upward motion through still air, the drag force is downward. A vector model writes drag opposite the relative velocity through the atmosphere, not necessarily opposite ground-relative velocity. Wind changes that relative velocity.

The vertical differential equation becomes

m \frac{d v}{d t} = T - m g - \frac{1}{2} ρ C_{D} A v ∣ v ∣.

The absolute-value form gives a drag force opposite the signed velocity for either ascent or descent. A constant-density assumption is often adequate for a short laboratory launch but poor for a high ascent. Density decreases strongly with altitude, and $C_{D}$ can change near transonic speed. Record that atmospheric assumption separately from the numerical solution so the result remains traceable to the data used.

Gravity also varies with distance from Earth. The local form $m g$ can be replaced by $GM m / r^{2}$ when the altitude change is large. The mass cancels from the gravity acceleration but does not cancel from the thrust acceleration $T / m$ . A trajectory integrator can update both terms at each time step:

a_{k} = \frac{T _{k} t ^ _{k} + F _{drag, k}}{m _{k}} - \frac{GM}{r _{k}^{3}} r_{k} .

The unit vector $\hat{t}_{k}$ gives nozzle direction. A numerical update must use a time step short enough to resolve the fastest change among thrust, attitude, mass rate, and atmospheric force. Halving the step and comparing the output provides a basic convergence test. A trajectory that shifts materially when the time step is halved has not yet reached a reliable numerical resolution.

The ideal rocket equation is often used as a mission-level accounting relation before a detailed trajectory is constructed. It estimates the propulsion velocity increment available from a mass ratio. A detailed ascent then subtracts gravity and drag losses and includes steering loss, which arises when some thrust is directed away from the desired velocity-change direction. These effects are not corrections to the logarithm; they come from external forces and the changing direction of the thrust term in the vector equation.

During a curved ascent only the thrust component along the desired velocity-change direction adds to it; the perpendicular component turns the velocity vector and shows up as steering loss.

Mass Ratio and Staging

Write the ignition mass as

m_{0} = m_{p} + m_{s} + m_{L},

where $m_{p}$ is usable propellant, $m_{s}$ is stage structure and engine hardware, and $m_{L}$ is payload plus any upper-stage mass carried by the stage. At burnout,

m_{f} = m_{s} + m_{L} .

The ideal velocity change depends on the ratio $m_{0} / m_{f}$ , not on propellant mass alone. Adding propellant increases the numerator, but tanks and feed hardware add to the burnout mass as well. The inert mass fraction is a structural measure,

ϵ = \frac{m _{s}}{m _{s} + m _{p}} .

At fixed payload and exhaust speed, lowering $ϵ$ can raise the available mass ratio. It does not make the structure optional: tanks, engines, thermal protection, guidance hardware, and load paths have minimum mass set by strength, temperature, vibration, and manufacturing constraints.

Ignition mass separates into payload, propellant, and stage hardware. Burnout removes the propellant but keeps payload and hardware, so both count in the burnout mass that sets the mass ratio.

Solving the rocket equation for mass ratio yields

\frac{m _{0}}{m _{f}} = exp (\frac{Δ v}{u _{e}}) .

Payload fraction must also be defined explicitly. A vehicle-level payload fraction can mean $m_{L} / m_{0}$ , whereas an individual stage may carry another stage as part of its payload. Mixing those definitions makes comparisons ambiguous. In a staging calculation, label the mass present immediately before and after each event. The discarded hardware has mass before separation and no mass after separation; its removal changes the starting mass for the next burn but does not itself impart a thrust-derived velocity increment in the ideal separation model.

Stage separation discards the empty lower-stage hardware after its burn. Velocity is continuous across the ideal separation, but the mass carried into the next burn drops sharply, improving the next stage's mass ratio.

Staging and sequences of burns.

Suppose stage 1 starts with mass $m_{0, 1}$ and ends its powered burn at $m_{f, 1}$ . Its ideal velocity change is

Δ v_{1} = u_{e, 1} ln (\frac{m _{0, 1}}{m _{f, 1}}) .

After stage-1 hardware is released, the mass changes discontinuously to $m_{0, 2}$ . Stage 2 then supplies

Δ v_{2} = u_{e, 2} ln (\frac{m _{0, 2}}{m _{f, 2}}) .

Collinear burns have ideal total

Δ v_{total} = k \sum u_{e, k} ln (\frac{m _{0, k}}{m _{f, k}}) .

Each stage may have a different exhaust speed. A lower stage may emphasize high thrust for lift-off, while an upper stage may emphasize a large velocity increment. The sum remains an ideal propulsion accounting total. Gravity, drag, and steering must be integrated over the actual timed sequence before it is compared with a mission requirement.

Stage separation has practical dynamics beyond the ideal mass calculation. The vehicles must avoid collision, exhaust impingement, and recontact. Separation springs or small thrusters can give relative motion. The resulting momentum changes are usually small compared with the main propulsion increments but should be included when predicting close-proximity motion. A detailed stage model also requires the centre of mass and inertia of the retained vehicle because steering torque changes after hardware is discarded.

An orbit manoeuvre is another sequence of burns, often separated by long coasting intervals. Each burn has a velocity-change vector. A burn aligned with the local velocity raises the orbital energy in the simplest two-body picture; a burn opposite the local velocity lowers it. Radial or normal components change the orbit in other ways. The rocket equation gives the propellant relation for each commanded vector magnitude. Orbital mechanics determines the required vectors and timings.

Worked Vertical Burn

Consider an ideal vertical stage with ignition mass $m_{0} = 2.40 \times 1 0^{5} kg$ , burnout mass $m_{f} = 8.00 \times 1 0^{4} kg$ , constant exhaust speed $u_{e} = 3.20 km s^{- 1}$ , and constant burn rate $R = 1.20 \times 1 0^{3} kg s^{- 1}$ . Suppose it starts from rest at a location where $g = 9.81 m s^{- 2}$ . The calculation ignores drag, variation of gravity, and steering. The result is a propulsion and gravity-loss calculation under an idealized vertical flight model.

The two terms refer to different parts of the flight model. The stage supplies $3.52 km s^{- 1}$ of ideal propulsion capability over the mass interval, while gravity lowers the vertical speed during this particular burn.

The worked burn's velocity account: the ideal propulsion gain (

3.52 km/s

) minus the gravity loss (

1.31 km/s

) over the same

133 s

interval gives the simplified burnout speed (

2.21 km/s

A numerical model of the same case should reproduce these values as its time step shrinks. At each step, update mass by $m_{k + 1} = m_{k} - R Δ t$ , compute $a_{k} = T / m_{k} - g$ , update velocity, and then update height. The model should stop its powered update at the prescribed burnout mass. This prevents $m$ from decreasing below dry-plus-payload mass. An automatic stop condition handles intervals for which $t_{b} /Δ t$ is not an integer.

Algorithm:

1
set m <- m0, v <- v0, y <- y0, t <- 0
2
while m > mf do
3
dt_step <- min(dt, (m - mf) / R)
4
a <- T / m - g
5
v <- v + a dt_step
6
y <- y + v dt_step
7
m <- m - R dt_step
8
t <- t + dt_step
9
end while

The update shown is a basic explicit method. A midpoint or higher-order method reduces integration error for a given time step. The physical assumptions dominate the error once the step is sufficiently small; adding numerical precision cannot repair a drag law, mass schedule, or nozzle-direction model that does not match the vehicle.

Burn programs and acceleration limits.

An engine need not operate at constant thrust. If exhaust speed is approximately constant, a chosen thrust program $T (t)$ determines the mass rate through

\overset{m}{˙} = - \frac{T ( t )}{u _{e}} .

During vertical flight with prescribed maximum upward acceleration $a_{max}$ , the ideal thrust is bounded by

T \leq m (a_{max} + g) .

As mass falls, the permitted thrust falls in direct proportion. A throttle schedule can therefore begin at full thrust and reduce later in the burn. This is a load constraint, distinct from the requirement that total propellant and exhaust speed provide enough mission velocity change.

Throttle also changes the burn duration. A lower late-burn thrust uses propellant more slowly, which can increase gravity loss during a vertical ascent. A design comparison should report total ideal $Δ v$ , burn duration, maximum acceleration, and external-force loss over the trajectory together. Thrust or exhaust speed alone cannot show the trade among those quantities.

At fixed total propellant mass and constant $u_{e}$ , varying the thrust schedule does not change the ideal velocity increment in a force-free collinear burn. The integral $\int - u_{e} d m / m$ depends on the mass endpoints. In an atmosphere or a gravitational field, schedule and steering alter the duration and path, so they change the external-force contributions to the delivered velocity.

Measurements and Calibration

Thrust and mass-flow data are time series. A reliable estimate uses a common steady interval for both measurements. Let $T_{j}$ be load-cell samples and let $m_{a}, m_{b}$ be propellant-system masses at the endpoints $t_{a}, t_{b}$ . Over a steady window,

\overset{ˉ}{T} = \frac{1}{N} j = 1 \sum N T_{j}, R = \frac{m _{a} - m _{b}}{t _{b} - t _{a}}, u_{e} = \frac{T ˉ}{R} .

The mass record must represent propellant leaving the engine system. A scale that also includes a changing support load, an unmeasured vent stream, or a different fuel tank from the tested engine does not provide the required $R$ . A thrust trace may contain ignition and shutdown transients. Including them in a steady-engine average changes both the measured thrust and the time interval associated with the mass change.

A thrust test trace has transient rise and fall around a steady interval. Averaging thrust over the dashed steady window and pairing it with the mass change over that same window gives a consistent exhaust-speed estimate.

For independent small uncertainties, a first-order estimate for the inferred exhaust speed is

(\frac{σ _{u_{e}}}{u _{e}})^{2} \approx (\frac{σ _{T}}{T})^{2} + (\frac{σ _{R}}{R})^{2} .

The expression does not replace a calibration record. Load-cell gain, zero offset, temperature drift, vibration, and alignment can introduce systematic error. A misaligned engine transfers a transverse force to the mount; an axial load cell may then report only a component of the thrust. The fixture should define the intended axis, record the sign convention, and use known reference loads before and after a test when drift matters.

Velocity change inferred from masses also has uncertainty. With constant $u_{e}$ ,

Δ v = u_{e} ln (m_{i} / m_{f}) .

Mass errors near burnout matter strongly because the logarithm depends on the ratio. A reported $Δ v$ should identify whether $m_{f}$ includes trapped residual propellant, pressurant, usable reserves, and separation hardware. Each choice can be reasonable for a stated event, but values taken from different event definitions cannot be combined in one mass ratio.

Energy and General Mass Exchange

The rocket equation comes from momentum balance. Kinetic energy is generally not conserved between rocket and exhaust, even in force-free flight, because chemical or stored energy in the propulsion system becomes kinetic and internal energy. The rocket receives forward momentum while the exhaust receives backward momentum. The energy partition depends on the reference frame and on how much propellant has already been expelled.

A small exhaust packet of mass $d m_{ex}$ has laboratory velocity $v - u_{e}$ in the one-dimensional rear-nozzle model. Its kinetic energy contribution is $\frac{1}{2} d m_{ex} (v - u_{e})^{2}$ . The vehicle's kinetic-energy change during the same interval contains both a speed change and a mass change. Equating the rocket's kinetic-energy gain to an exhaust kinetic-energy loss omits the stored energy released by the engine and produces an incorrect propulsion relation.

Momentum transfer gives the rocket and exhaust equal and opposite momentum changes, but their kinetic-energy changes need not match. Stored propulsion energy (chemical or pressure) supplies the difference.

Energy analysis remains essential for engine and vehicle design. It constrains how chemical, electrical, nuclear, or stored-pressure energy becomes exhaust kinetic energy, heat, radiation, and acoustic energy. The momentum balance then converts the effective exhaust speed and mass rate into thrust and velocity change. Separate momentum and energy balances distinguish an engine with adequate energy release but low jet momentum from a high-thrust engine with a short available burn.

The reference-frame dependence of kinetic energy is visible in a simple example. If a rocket and its exhaust are observed from a frame moving forward, both laboratory speeds increase by the same amount. Their momenta and kinetic energies change in that frame, while their relative speed and the ideal thrust prediction remain unchanged. Work or energy claims for an open system therefore require a declared frame and a complete account of energy crossing the boundary.

Incoming mass and the general exchange term.

For comparison with mass loss, suppose a cart of mass $m$ gains material while travelling at speed $v$ . Material with laboratory speed $u$ lands in it at rate $\overset{m}{˙} > 0$ . With external horizontal force $F_{ext}$ , the one-dimensional equation is

m \overset{v}{˙} = F_{ext} + (u - v) \overset{m}{˙} .

For material initially at rest on the ground, $u = 0$ and the transfer term is $- v \overset{m}{˙}$ . The cart must accelerate each incoming kilogram from zero to its own speed, so it slows unless an applied force supplies the missing momentum. For a rocket, $\overset{m}{˙} < 0$ and $u - v$ is the backward relative exhaust speed. The sign reverses naturally when the same general equation is used consistently.

An accumulating rope on a scale offers a vertical example. The scale supports the weight of rope already at rest and also removes the downward momentum of material arriving at the pan. The measured force can therefore exceed the weight of the accumulated portion. A free-body diagram containing only the resting weight misses the momentum flux. Variable-mass systems are often difficult because the material crossing the boundary changes state as it crosses: moving rope becomes resting rope, incoming sand becomes cart mass, and propellant becomes a backward jet.

Model Limits and Checks

Several errors recur in variable-mass calculations.

Fixed-mass impulse for a whole burn. Replacing the changing mass by either the ignition or burnout mass in $Δ v = T Δ t / m$ gives only a rough estimate. The finite-burn result is logarithmic.
Laboratory exhaust speed inserted as relative speed. The exhaust laboratory speed includes rocket speed. The mass-transfer term uses the difference between exhaust and vehicle velocities.
Dry mass omitted from a mass ratio. Tanks, engines, payload, residual fluids, and structures remain after propellant depletion. Removing them from $m_{f}$ creates an impossible performance estimate.
Powered dynamics continued through cutoff. The thrust term exists only while propellant crosses the nozzle boundary. At cutoff, set $R = 0$ and solve the coast equation with the appropriate external forces.
External forces hidden inside the ideal velocity change. Gravity, drag, contact force, and steering effects belong in $F_{ext}$ or in a trajectory integration. The mass-ratio expression describes propulsion alone.

Model limits and richer flight models.

The nonrelativistic rocket equation assumes speeds small compared with the speed of light. At relativistic exhaust or vehicle speeds, momentum and velocity addition require special relativity. The system-boundary analysis still applies, but $p = m v$ and the ordinary logarithmic relation must be replaced by relativistic expressions. The present treatment also uses a point-mass translation model. Attitude dynamics, flexible vibration, sloshing, control delay, and structural load require rotational and deformation equations alongside the translation equation.

An effective exhaust speed summarizes a complicated engine. It can vary with throttle, altitude, nozzle geometry, mixture state, and chamber pressure. A model intended for a particular engine should use measured thrust and mass-flow curves or a validated engine map. The constant- $u_{e}$ analytical result is a limiting check: a simulation should approach it when drag, gravity, steering, and schedule variation are removed from the model.

The force model can be extended in controlled stages. First use mass change and a specified exhaust direction. Add gravity for vertical or orbital motion. Add air force using atmospheric density and relative wind. Add nozzle steering and rotational dynamics once the vehicle attitude affects thrust direction. Each layer adds parameters and measurements. A detailed result is only as credible as its least constrained layer, so the model description should state which forces were included and which were deliberately omitted.

Ground tests have their own boundary conditions. A stationary engine has no vehicle velocity change, yet it expels a jet and transfers force to its mounting structure. A test stand gives an accurate thrust record only after its stiffness, alignment, dynamic response, and thermal environment have been considered. A small rocket flight test introduces wind, guide-rail contact, sensor lag, and finite sampling. The same variable-mass equation applies, while the available measurements and external-force model determine the uncertainty of any inferred exhaust speed or trajectory result.

Calculation checks.

Before accepting a result, record the beginning and ending vehicle masses for every powered segment, the effective exhaust speed or measured thrust and mass rate, the nozzle direction, and every significant external force. Use one declared inertial frame for all velocities in a momentum balance. Then check dimensions:

[R u_{e}] = N, [u_{e} ln (\frac{m _{i}}{m _{f}})] = m s^{- 1}, [g t] = m s^{- 1} .

A negative computed propellant mass, a mass ratio below one during a burn, or a forward-facing exhaust vector paired with forward thrust indicates a sign or event definition error. Check limiting cases as well. As $m_{f} \to m_{i}$ , the ideal velocity change must approach zero. With $R \to 0$ at fixed mass, thrust must approach zero. With $F_{ext} = 0$ and a constant $u_{e}$ , the differential relation must integrate to the standard logarithm.

A measured burn permits independent routes to the same quantities. A load-cell thrust integrated over the recorded interval gives impulse. A mass-flow record and an inferred effective exhaust speed give $\int R u_{e} d t$ . Agreement within the stated calibration and sampling uncertainty provides a consistency test. A trajectory-derived velocity change adds gravity, air force, and steering models; it should be compared with the ideal propulsion value only after those contributions are included explicitly.

A concise report identifies the exhaust-speed data source, the labelled ignition and cutoff masses, and the external-force model used for any reported ascent speed. These details distinguish a propulsion capability from a trajectory outcome and support comparison with a higher-fidelity calculation.

A defined system boundary fixes the material whose mass changes. Momentum transport through that boundary produces thrust, and integration over the changing mass gives the logarithmic velocity relation. External forces, burn scheduling, staging, and measurements determine the observed trajectory from that ideal propulsion result.

Multiple Engines and Thrust Control

A vehicle can carry several engines whose exhaust directions are not identical. For engine $j$ , let $R_{j}$ be its positive mass-flow rate and let $u_{e, j}$ point opposite the jet in the vehicle frame. The total thrust term is the vector sum

T = - j \sum u_{e, j} \overset{m}{˙}_{j} = j \sum R_{j} u_{e, j} .

The total mass rate is $\overset{m}{˙} = - \sum_{j} R_{j}$ . Symmetric engines can provide a large axial thrust while their transverse components cancel. Differential throttling or gimballed nozzles changes the component balance and creates a torque about the vehicle centre of mass. Translational acceleration uses the net force. Attitude control uses the torque of each thrust line about the centre of mass. Treating a multi-engine vehicle as a point mass is adequate only when the attitude is already known or when symmetry makes the torque sum zero.

Two symmetric engines have transverse thrust components that cancel while their axial components add to the net force. Unequal thrust would leave a transverse force and a turning torque about the center of mass.

The efficiency of a commanded manoeuvre can be expressed geometrically. A thrust component parallel to the requested velocity-change vector contributes directly to that request. A perpendicular component changes direction instead. During a launch, the perpendicular component may be necessary for a planned turn. During a precise orbit correction, an unwanted component is an error source. The propulsion mass required for a small vector impulse follows from the magnitude of the commanded change and the applicable mass at that event; direction affects the trajectory, contact constraints, and attitude requirements.

Event bookkeeping.

An engineering mass table should use rows for events. A typical segment has ignition mass, mass after planned burn, mass after residuals or shutdown, mass after separation, and the next-stage ignition mass. The rocket equation for that segment uses the mass immediately before and immediately after propellant expulsion. A separation mass drop belongs between segments. This convention prevents a discarded lower stage from being mistakenly included as upper-stage burnout mass.

A sequence of short burns with the same engine has additive ideal contributions, even when the burns are separated by coast intervals. The coast does not use propellant in this simple model, so the starting mass of the next burn is the ending mass of the previous event. External forces during the coast still alter the orbit or flight path. A velocity budget that sums only propulsion increments is therefore a resource statement; it becomes a trajectory prediction only after the coasts and all external forces are propagated.

Interpreting a rocket result.

The relation

Δ v = u_{e} ln (m_{i} / m_{f})

gives the velocity increment supplied by an ideal exhaust stream between two labelled masses under the stated exhaust-speed assumption. Launch altitude, orbital radius, final ground speed, and structural safety margin require additional equations and data. The system, time interval, and assumptions must remain attached to the mass-ratio calculation before it is used in a trajectory claim.

The ideal increment and a trajectory prediction therefore belong in separate rows of the result. Their input data overlap, but their physical boundaries do not.

Statement	Relation	Additional information needed
Ideal burn increment	$Δ v = u_{e} ln (m_{i} / m_{f})$	exhaust model and labelled masses
Instantaneous acceleration	$a = (T + F_{ext}) / m$ along the chosen axis	external forces and attitude
Trajectory outcome	integrated position and velocity	gravity, drag, guidance, and initial state

The same discipline applies to a thrust value. $T = R u_{e}$ describes the instantaneous axial momentum transfer in the effective-exhaust-speed model. The vehicle response depends on its instantaneous mass and on every other force. A thrust trace with no mass record cannot determine an exhaust speed. A mass record with no force or jet measurement cannot determine thrust. Each measured quantity has a role in the momentum balance, and the relationships among them provide the most diagnostic checks on a propulsion calculation.

Synchronized records are required when $u_{e}$ is inferred from a firing rather than taken from a specification. The thrust channel, propellant-flow measurement, and vehicle-mass record must share a clock and a documented delay correction. A load cell can respond after the chamber pressure changes, while a flow meter can average over a different interval. Dividing unmatched samples can therefore create a false variation in effective exhaust speed. Compute impulse and expelled mass over the same accepted window, propagate the baseline and timing uncertainties, and compare the integrated estimate with the pointwise ratio only where both signals exceed their resolution floors. The pre-ignition load-cell mean supplies a zero check; the post-cutoff tail determines whether late impulse has been omitted. A stated window, sampling rate, filter, and delay correction make the inferred exhaust speed reproducible and separate instrument timing from genuine engine transients.

Event timing resolves a further ambiguity. A reported cutoff mass may be measured when a command is issued, when valves close, when the chamber pressure decays, or when the thrust trace falls below a stated threshold. These events can differ by a measurable propellant mass. The mass-ratio denominator should be tied to the same event used to end the thrust integration. If a report includes residual propellant as unavailable mass, that choice belongs in both the stated mass model and the uncertainty estimate. Clear event labels make comparisons between engine tests, trajectory simulations, and stage budgets meaningful.

For classroom and laboratory work, place the sign convention beside the first diagram, label the two masses with their event times, and retain units in every intermediate line. Those records expose sign and event-definition errors before a numerical result is used and support later changes to the burn rate, payload, or external-force model.

A propulsion result should retain the burn events and the measurement that controls each derived quantity. The logarithmic mass ratio gives the ideal burn increment; trajectory prediction also requires the event-matched records below.

Reported quantity	Relation	Matched event or measurement
Ideal velocity increment	$Δ v = u_{e} ln (m_{i} / m_{f})$	labelled start and cutoff masses
Thrust	$T = R u_{e}$	mass-flow rate and effective exhaust speed
External-force correction	impulse integrated over the burn	same time window as the mass change
Stage performance	post-separation mass and velocity	release event and payload boundary

They also expose which measurement would most improve a later revision of the trajectory model.