Tsiolkovsky rocket equation

Tsiolkovsky's rocket equation, named after Konstantin Tsiolkovsky who first derived it, considers the principle of a rocket: a device that can apply an acceleration to itself (a thrust) by expelling part of its mass with high speed in the opposite direction, due to the conservation of momentum.

It says that for any maneuver or any journey involving a number of maneuvers:

<math>\Delta v = v_e \ln \frac {m_0} {m_1}<math>

or equivalently

<math>m_1=m_0 e^{-\Delta v\ / v_e}<math> or <math>m_0=m_1 e^{\Delta v\ / v_e}<math>

where <math>m_0<math> is the initial total mass, and <math>m_1<math> the final total mass and <math>v_e<math> the velocity of the rocket exhaust with respect to the rocket (the specific impulse).

<math>\Delta v<math> (delta v) is the integration over time of the magnitude of the acceleration produced by using the rocket engine (not the acceleration due to other sources such as gravity or drag). For the typical case of an acceleration in the direction of the velocity, this is the increase of the speed. In the case of an acceleration in opposite direction (deceleration) it is the decrease of the speed. Note that gravity or drag also change velocity, but they are not part of the quantity delta-v. Hence delta-v is not simply the change in speed or velocity. However, thrust is often applied in short bursts, and during these short periods the other sources of acceleration may be negligable, and the delta-v of one burst may be simply approximated by the speed change. The total delta-v can simply be found by addition, even though between bursts the magnitude and direction of the velocity changes due to gravity, e.g. in an elliptic orbit.

Note that, as mentioned, at any time the magnitude of the acceleration contributes to the delta-v, hence always a non-negative value, regardless of whether the rocket is used for acceleration or deceleration. This again demonstrates that delta-v is not simply the change in speed or velocity: the latter may be zero if we first accelerate and than decelerate, but the delta-v accumulates.

The equation is obtained by integrating the conservation of momentum equation

for a simple rocket that emits mass at a constant velocity (<math>dm<math> is here the reaction mass; if it is the change of the rocket mass then there is a minus sign in the latter equation).

The equation does not take into account the discarding of rocket stages; these reduce the mass even more, but unfortunately without providing the corresponding momentum (they can not be ejected with the exhaust speed); however, for the resulting payload this is more advantageous than having to subtract the mass of used fuel tanks and engines from the final <math>m_1<math>.

Although an extreme simplification, the rocket equation captures the essentials of rocket flight physics in a single short equation. It happens that delta-v is one of the most important quantities in orbital mechanics, that quantifies how difficult it is to get from one trajectory to another.

Clearly, to achieve a large delta-v, either <math>m_0<math> must be huge (growing exponentially as delta-v rises), or <math>m_1<math> must be tiny, or <math>v<math> must be very high, or some combination of all of these.

In practice, this has been achieved by using very large rockets (increasing <math>m_0<math>), with multiple stages (decreasing <math>m_1<math>), and rockets with very high exhaust velocities. The Saturn V rockets used in the Apollo space program and the ion thrusters used in long-distance unmanned probes are good examples of this.

The rocket equation shows a kind of "exponential decay" of mass, but not as a function of time, but as a function of delta-v produced. The delta-v that is the corresponding "half-life" is <math>v_e \ln 2 \approx 0.693 v_e<math>

Energy

In the ideal case <math>m_1<math> is useful payload and <math>m_0-m_1<math> is reaction mass (this corresponds to empty tanks having no mass, etc.). The energy required can simply be computed as

Seemingly this is just the kinetic energy of the reaction mass and not the kinetic energy required for the payload, but if e.g <math>v_e<math>=10 km/s and the speed of the rocket is 3 km/s, then the speed of the reaction mass changes only from 3 to 7 km/s; the energy thus "saved" corresponds to the increase of the specific kinetic energy (kinetic energy per kg) for the rocket. In general:

Thus <math>\Delta \epsilon = \int v\, d (\Delta v)<math>

where <math>\epsilon<math> is the specific energy of the rocket (potential plus kinetic energy) and <math>\Delta v<math> is a separate variable, not just the change in <math>v<math>. In the case of using the rocket for deceleration, i.e. expelling reaction mass in the direction of the velocity, <math>v<math> should be taken negative.

The formula is for the ideal case again, with no energy lost on heat, etc. The latter causes a reduction of thrust, so it is a disadvantage even when the objective is to lose energy (deceleration).

If the energy is produced by the mass itself, as in a chemical rocket, the fuel value has to be :<math>\frac{1}{2}v_e^2<math>, where for the fuel value also the mass of the oxidizer has to be taken into account. A typical value is <math>v_e=4.5 km/s<math>, corresponding to a fuel value of 10.1 MJ/kg. The actual fuel value is higher, but part of the energy is lost on heat that flows off as radiation.

The required energy is

<math>E = \frac{1}{2}m_1(e^{\Delta v\ / v_e}-1)v_e^2<math>

Conclusions:

for <math>\Delta v <
for a given <math>\Delta v<math>, the minimum energy is needed if <math>v_e=0.6275 \Delta v<math>, requiring an energy of

<math>E = 0.772 m_1(\Delta v)^2<math>.

Starting from zero speed this is 54.4 % more than just the kinetic energy of the payload. Starting from a nonzero speed the required energy may be less than the increase in energy in the payload. This can be the case when the reaction mass has a lower speed after being expelled than before. For example, from a LEO of 300 km altitude to an escape orbit is an increase of 29.8 MJ/kg, which, using a specific impulse of 4.5 km/s, has a net cost of 20.6 MJ/kg (<math>\Delta v<math> = 3.20 km/s; the energies are per kg payload).

This optimization does not take into account the masses of various kinds of engines.

Also, for a given objective such as moving from one orbit to another, the required <math>\Delta v<math> may depend greatly on the rate at which the engine can produce <math>\Delta v<math> and maneuvers may even be impossible if that rate is too low. For example, a launch to LEO normally requires a <math>\Delta v<math> of ca. 9.5 km/s (mostly for the speed to be acquired), but if the engine could produce <math>\Delta v<math> at a rate of only slightly more than g, it would be a slow launch requiring altogether a very large <math>\Delta v<math> (think of hovering without making any progress in speed or altitude, it would cost a <math>\Delta v<math> of 9.8 m/s each second). If the possible rate is only <math>g<math> or less, the maneuver can not be carried out at all with this engine.

The power is given by

where <math>F<math> is the thrust and <math>a<math> the acceleration due to it. Thus the theoretically possible thrust per unit power is 2 divided by the specific impulse in m/s. The thrust efficiency is the actual thrust as percentage of this.

If e.g. solar power is used this restricts <math>a<math>; in the case of a large <math>v_e<math> the possible acceleration is inversely proportional to it, hence the time to reach a required delta-v is proportional to <math>v_e<math>; with 100% efficiency:

for <math>\Delta v <

Examples:

power 1000 W, mass 100 kg, <math>\Delta v<math>= 5 km/s, <math>v_e<math>= 16 km/s, takes 1.5 months.
power 1000 W, mass 100 kg, <math>\Delta v<math>= 5 km/s, <math>v_e<math>= 50 km/s, takes 5 months.

Thus <math>v_e<math> should not be too large.

Examples

Assume a specific impulse of 4.5 km/s and a <math>\Delta v<math> of 9.7 km/s (Earth to LEO).

SSTO rocket: <math>1-e^{-9.7/4.5}<math> = 0.884, therefore 88.4 % of the initial total mass has to be propellant. The remaining 11.6 % is for the engines, the tank, and the payload. In the case of a space shuttle, it would also include the orbiter.

Two stage to orbit: suppose that the first stage should provide a <math>\Delta v<math> of 5.0 km/s; <math>1-e^{-5.0/4.5}<math> = 0.671, therefore 67.1% of the initial total mass has to be propellant. The remaining mass is 32.9 %. After deposing of the first stage, a mass remains equal to this 32.9 %, minus the mass of the tank and engines of the first stage. Assume that this is 8 % of the initial total mass, then 24.9 % remains. The second stage should provide a <math>\Delta v<math> of 4.7 km/s; <math>1-e^{-4.7/4.5}<math> = 0.648, therefore 64.8% of the remaining mass has to be propellant, which is 16.2 %, and 8.7 % remains for the tank and engines of the second stage, the payload, and in the case of a space shuttle, also the orbiter. Thus together 16.7 % is available for all engines, the tanks, the payload, and the possible orbiter.

Tsiolkovsky rocket equation

Energy

Examples

See also