Tsiolkovsky rocket equation

Tsiolkovsky rocket equation

Tsiolkovsky's rocket equation, or "ideal rocket equation" is named after Konstantin Tsiolkovsky, who independently derived it and published in his 1903 work [К. Э. Циолковский, Исследование мировых пространств реактивными приборами, 1903. It is available online [http://epizodsspace.testpilot.ru/bibl/dorev-knigi/ciolkovskiy/issl-03st.html here] ] , 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.

The equation relates the delta-v with the effective exhaust velocity and the initial mass and the end mass of a rocket.

For any such maneuver (or journey involving a number of such maneuvers):

:Delta v = v_ ext{e} ln frac {m_0} {m_1}

where::m_0 is the initial total mass, including propellant, in kg (or lb):m_1 is the final total mass in kg (or lb):v_ ext{e} is the effective exhaust velocity in m/s or (ft/s) or v_ ext{e} = I_ ext{sp} cdot g_0:Delta v is the delta-v in m/s (or ft/s)

History

This equation was independently derived by Konstantin Tsiolkovsky towards the end of the 19th century and is widely known under this name and ideal rocket equation. However a recently discovered pamphlet "A Treatise on the Motion of Rockets" by William Moore [Johnson W., "Contents and commentary on William Moore's a treatise on the motion of rockets and an essay on naval gunnery", International Journal of Impact Engineering, Volume 16, Number 3, June 1995, pp. 499-521] shows that the earliest known derivation of this kind of equation was in fact in Royal Military Academy at Woolwich in England in 1813, and was used for weapons research.

tages

In the case of sequentially thrusting rocket stages, the equation applies for each stage, where for each stage the initial mass in the equation is the total mass of the rocket after discarding the previous stage, and the final mass in the equation is the total mass of the rocket just before discarding the stage concerned. For each stage the specific impulse may be different.

For example, if 80% of the mass of a rocket is the fuel of the first stage, and 10% is the dry mass of the first stage, and 10% is the remaining rocket, then

:egin{align}Delta v & = v_ ext{e} ln { 100 over 100 - 80 }\ & = v_ ext{e} ln 5 \ & = 1.61 v_ ext{e} \end{align}

With three similar, subsequently smaller stages with the same v_e for each stage, we have

:Delta v = 3 v_ ext{e} ln 5 = 4.83 v_ ext{e}

and the payload is 10%*10%*10% = 0.1% of the initial mass.

A comparable SSTO rocket, also with a 0.1% payload, could have a mass of 11% for fuel tanks and engines, and 88.9% for fuel. This would give

:Delta v = v_ ext{e} ln(100/11.1) = 2.20 v_ ext{e}

If the motor of a new stage is ignited before the previous stage has been discarded and the simultaneously working motors have a different specific impulse (as is often the case with solid rocket boosters and a liquid-fuel stage), the situation is more complicated.

Energy

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

:frac{1}{2}(m_0-m_1)v_ ext{e}^2

This is the kinetic energy of accelerating the reaction mass.

However, this is not the kinetic energy gained by the rocket and payload. To see this, if, for example, v_e=10 km/s and the speed of the rocket is 3 km/s, then the speed of a small amount of expended reaction mass changes from 3 forwards to 7 km/s rearwards. From conservation of energy, the energy difference must corresponds to the increase of the specific kinetic energy (kinetic energy per kg) for the rocket and payload.

In general:

:dleft(frac{1}{2}v^2 ight)=vdv=vv_ ext{e}dm/m=frac{1}{2}left(v_ ext{e}^2-(v-v_ ext{e})^2+v^2 ight)dm/m

Thus the specific energy gain of the rocket in any small time interval is the energy gain of the rocket including the remaining fuel, divided by its mass, where the energy gain is equal to the energy produced by the fuel minus the energy gain of the reaction mass. The larger the speed of the rocket, the smaller the energy gain of the reaction mass; if the rocket speed is more than half of the exhaust speed the reaction mass even loses energy on being expelled, to the benefit of the energy gain of the rocket; the larger the speed of the rocket, the larger the energy loss of the reaction mass.

We have:Delta epsilon = int v, d (Delta v) where epsilon is the specific energy of the rocket (potential plus kinetic energy) and Delta v is a separate variable, not just the change in v. In the case of using the rocket for deceleration, i.e. expelling reaction mass in the direction of the velocity, v 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 scriptstyle{v_ ext{e}^2/2}, where for the fuel value also the mass of the oxidizer has to be taken into account. A typical value is v_ ext{e} = 4.5 km/s, corresponding to a fuel value of 10.1 MJ/kg. The actual fuel value is higher, but much of the energy is lost as waste heat in the exhaust that the nozzle was unable to extract.

The required energy E is

:E = frac{1}{2}m_1left(e^{Delta v / v_ ext{e-1 ight)v_ ext{e}^2

Conclusions:

*for Delta v ll v_e we have Eapprox frac{1}{2}m_1 v_ ext{e} Delta v
*for a given Delta v, the minimum energy is needed if v_ ext{e}=0.6275 Delta v, requiring an energy of:E = 0.772 m_1(Delta v)^2.: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 (Delta v = 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 Delta v may depend greatly on the rate at which the engine can produce Delta v and maneuvers may even be impossible if that rate is too low. For example, a launch to LEO normally requires a Delta v of ca. 9.5 km/s (mostly for the speed to be acquired), but if the engine could produce Delta v at a rate of only slightly more than "g", it would be a slow launch requiring altogether a very large Delta v (think of hovering without making any progress in speed or altitude, it would cost a Delta v of 9.8 m/s each second). If the possible rate is only g or less, the maneuver can not be carried out at all with this engine.

The power is given by:P= frac{1}{2} m a v_ ext{e} = frac{1}{2}F v_ ext{e} where F is the thrust and a 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 a; in the case of a large v_ ext{e} the possible acceleration is inversely proportional to it, hence the time to reach a required delta-v is proportional to v_ ext{e}; with 100% efficiency:
*for Delta v ll v_ ext{e} we have tapprox frac{m v_ ext{e} Delta v}{2P}

Examples:
*power 1000 W, mass 100 kg, Delta v= 5 km/s, v_ ext{e}= 16 km/s, takes 1.5 months.
*power 1000 W, mass 100 kg, Delta v= 5 km/s, v_ ext{e}= 50 km/s, takes 5 months.

Thus v_ ext{e} should not be too large.

Derivation

Consider the following system:

Newton's second law of motion relates external forces and change in linear momentum as follows:

:egin{array}{l} sum F = lim_{Delta t o 0} frac{Delta P}{Delta t} \\ Delta P = P_2 - P_1 \\ P_1 = left( {m + Delta m} ight)V_ ext{b/i} \\ end{array}

:V_ ext{b/i}: Velocity of body with respect to inertial frame.

:egin{array}{l} P_2 = mleft( {V_ ext{b/i} + Delta V_ ext{b/i} } ight) + Delta mV_ ext{ex/i} \\ Delta P = mDelta V_ ext{b/i} + Delta mleft( {V_ ext{ex/i} - V_ ext{b/i} } ight) \\ V_ ext{ex/i} - V_ ext{b/i} = V_ ext{ex/b} Rightarrow Delta P = mDelta V_ ext{b/i} + Delta mV_ ext{ex/b} \\ Rightarrow sum F = mfracdV_ ext{b/i} dt + V_ ext{ex/b} fracdmdt end{array}

:V_ ext{ex/i}: Velocity of exhaust with respect to inertial frame

:V_ ext{ex/b}: Velocity of exhaust with respect to body

If there are no external forces:

:egin{array}{l} sum F = 0 Rightarrow m,dV_ ext{b/i} = - V_ ext{ex/b} dm Rightarrow dV_ ext{b/i} = - V_ ext{ex/b} fracdm{m} \\ Rightarrow Delta V_ ext{b/i} = V_ ext{ex/b} ln frac {m_0} {m_1} \\ end{array}

Substitute the following:

:egin{array}{l}Delta v = Delta V_ ext{b/i} \ v_ ext{e} = V_ ext{ex/b} end{array}

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

:Delta v = v_ ext{e} ln frac {m_0} {m_1}

or equivalently

:m_1=m_0 e^{-Delta v / v_ ext{e or m_0=m_1 e^{Delta v / v_ ext{e

where m_0 is the initial total mass including propellant, and m_1 the final total mass and v_ ext{e} the velocity of the rocket exhaust with respect to the rocket (the specific impulse, or, if measured in time, that multiplied by gravity-on-Earth acceleration).

:1-frac {m_1} {m_0}=1-e^{-Delta v / v_ ext{e is the mass fraction (the part of the initial total mass that is spent as reaction mass).

Delta v (delta v) is the integration over time of the magnitude of the acceleration produced by using the rocket engine (what would be the actual acceleration if external forces were absent). In free space, for the case of 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.

Of course gravity and drag also accelerate the vehicle, and they can add or subtract to the change in velocity experienced by the vehicle. Hence delta-v is not usually the actual change in speed or velocity of the vehicle.

The equation is obtained by integrating the conservation of momentum equation

: frac{}{}m dv = v_ ext{e} dm

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

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

Clearly, to achieve a large delta-v, either m_0 must be huge (growing exponentially as delta-v rises), or m_1 must be tiny, or v_ ext{e} must be very high, or some combination of all of these.

In practice, this has been achieved by using very large rockets (increasing m_0), with multiple stages (decreasing m_1), 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 m_1, not as a function of time, but as a function of delta-v produced. The delta-v that is the corresponding "half-life" is v_ ext{e} ln 2 approx 0.693 v_ ext{e}

Examples

Assume an exhaust velocity of 4.5 km/s and a Delta v of 9.7 km/s (Earth to LEO).

*SSTO rocket: 1-e^{-9.7/4.5} = 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 Delta v of 5.0 km/s; 1-e^{-5.0/4.5} = 0.671, therefore 67.1% of the initial total mass has to be propellant. The remaining mass is 32.9 %. After disposing 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 Delta v of 4.7 km/s; 1-e^{-4.7/4.5} = 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.

Applicability

The rocket equation is based on physical laws, and therefore must be obeyed by all reaction engines where all the propulsion is from propellant taken from the vehicle's mass. However, it does not apply to other technologies such as gun launches, space elevators, launch loops, tether propulsion and air-breathing engines.

See also

* Delta-v
* Delta-v budget
* Specific impulse
* Spacecraft propulsion
* Mass ratio
* Working mass
* Relativistic rocket

References

* [http://ed-thelen.org/rocket-eq.html How to derive the rocket equation]
* [http://www.relativitycalculator.com/rocket_equations.shtml Relativity Calculator - Learn Tsiolkovsky's rocket equations]


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Rocket — This article is about vehicles powered by rocket engines. For other uses, see Rocket (disambiguation). A Soyuz U, at Baikonur Site 1/5 A rocket is a missile, spacecraft, aircraft or other vehicle which obtains thrust from a rocket engi …   Wikipedia

  • Konstantin Tsiolkovsky — Born 17 September [O.S. 5 September] …   Wikipedia

  • Relativistic rocket — A relativistic rocket is any spacecraft that is travelling at a velocity close enough to light speed for relativistic effects to become significant. What significant means is a matter of context, but generally speaking a velocity of at least 50%… …   Wikipedia

  • Vis-viva equation — Astrodynamics Orbital mechanics Equations …   Wikipedia

  • Rocket engine — RS 68 being tested at NASA s Stennis Space Center. The nearly transparent exhaust is due to this engine s exhaust being mostly superheated steam (water vapor from its propellants, hydrogen and oxygen) …   Wikipedia

  • Konstantin Eduardovich Tsiolkovskii — Infobox Scientist name = Konstantin Tsiolkovsky box width = image width =150px caption = Konstantin Tsiolkovsky birth date = OldStyleDate|September 17|1857|September 5 birth place = Izhevskoye death date = September 19, 1935 death place = Kaluga… …   Wikipedia

  • Spacecraft propulsion — A remote camera captures a close up view of a Space Shuttle Main Engine during a test firing at the John C. Stennis Space Center in Hancock County, Mississippi Spacecraft propulsion is any method used to accelerate spacecraft and artificial… …   Wikipedia

  • Mass ratio — Astrodynamics Orbital mechanics Equations …   Wikipedia

  • Specific impulse — (usually abbreviated I sp) is a way to describe the efficiency of rocket and jet engines. It represents the impulse (change in momentum) per unit of propellant. The higher the specific impulse, the less propellant is needed to gain a given amount …   Wikipedia

  • Orbital mechanics — A satellite orbiting the earth has a tangential velocity and an inward acceleration. Orbital mechanics or astrodynamics is the application of ballistics and celestial mechanics to the practical problems concerning the motion of rockets and other… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”