AE2230 II Rocket Propulsion - 2019
AE2230 II Rocket Propulsion - 2019
Rocket Propulsion
February 2017
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Preface
The "rocket propulsion" part of TU-Delft’s Propulsion and Power course (AE2230-II) is
about the theory, technology and performance of rocket and spacecraft propulsion systems.
Some earlier courses can be of interest for this part, in particular:
Aerospace Design & Systems Engineering Elements I, part on spacecraft design and sizing
Physics I, part on thermodynamics
The level of detail of this course allows for a more in depth sizing of rocket propulsion
systems, taking into account physical and technological aspects not taught by previous
courses. Some typical high level learning objectives are given in the block below.
Specific Learning Objectives (At the end of this course, the student is able to…)
Understand the basic principles of thrust producing mechanisms for rockets and space
vehicles.
Perform basic sizing of thrust generation systems suitable for rockets and space
vehicles.
List/describe/explain:
the main thrust generation options available;
the components that make up the propulsion system and their function;
the current limits to thrust generation.
Apply control volume analysis and integral momentum equations to estimate the
produced thrust.
Develop a system model from a schematic system description.
Select the appropriate propulsion system from basic types depending on system
requirements.
Asses effect of changes in design/operating parameters on system performance.
Thanks to Raphael Klein for editing, modifying and correcting the old reader, and to Maria
Fernandez Jimenez for her transcripts of the video lectures.
ii
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Nomenclature
Abbreviation Description
ESA European Space Agency
EU European Union
IRT Ideal Rocket Theory
LEO Low Earth Orbit
LH2 Liquid Hydrogen
LOX Liquid Oxygen
LRE Liquid Rocket Engine
MMH Mono/Methyl Hydrazine
NTO Nitrogen Tetroxide
RP-1 Rocket Propellant-1
rpm Revolutions per Minute
SL Sea Level
SRM Solid Rocket Motor
SSME Space Shuttle Main Engine
UDMH Unsymmetrical Di/Methyl Hydrazine
iii
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
pC Chamber pressure Pa
pe Nozzle exhaust pressure Pa
P Power W
Pjet Jet power W
PT Thrust power W
r Regression rate m/s
t Time s
tb Burn time s
T Temperature K
TC Chamber temperature K
T0 Initial propellant temperature K
v Velocity m/s
ve Jet velocity m/s
veq Equivalent (or effective) jet velocity m/s
v0 Flight velocity m/s
(Greek)
p Pressure increase Pa
v Velocity change m/s
Nozzle expansion ratio -
Specific heat ratio -
Vandenkerckhove function -
Efficiency -
C Cycle efficiency -
P Propulsive efficiency -
Propellant mass density kg/m3
Important Constants
Definition Symbol Value
Earth’s gravitational acceleration at sea level g0 9.81 m/s2
Universal gas constant RA 8314 J/(K∙kmol)
iv
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Contents
PREFACE .............................................................................................................................................. II
NOMENCLATURE ................................................................................................................................ III
CONTENTS ...........................................................................................................................................V
1 – INTRODUCTION ........................................................................................................................ 1
1.1 How does a rocket work? ........................................................................................................ 1
1.2 Rocket propulsion features and applications .......................................................................... 1
1.3 Functions and components of rocket systems .......................................................................... 3
1.4 Rocket engines classification ................................................................................................... 4
2 - FUNDAMENTALS....................................................................................................................... 7
2.1 The working principle of rockets: exchange of momentum ..................................................... 7
2.2 The rocket equation ................................................................................................................. 9
2.3 The thrust equation................................................................................................................ 12
2.4 Total impulse and specific impulse ........................................................................................ 15
2.5 Thrust power and propulsive efficiency................................................................................. 18
3 – THE IDEAL ROCKET THEORY .............................................................................................. 20
3.1 Ideal Rocket Theory: assumptions......................................................................................... 20
3.2 Ideal Rocket Theory: building blocks .................................................................................... 22
3.3 Why a convergent-divergent nozzle? ..................................................................................... 25
3.4 Nozzle temperature and velocity equations ........................................................................... 26
3.5 Mach number relationships ................................................................................................... 28
3.6 Mass flow rate equation ........................................................................................................ 29
3.7 Expansion ratio equation ...................................................................................................... 30
3.8 Nozzle thrust .......................................................................................................................... 32
3.9 Characteristic velocity and thrust coefficient ........................................................................ 33
3.10 Ideal Rocket Theory: summary ............................................................................................ 36
4 – COLD GAS ROCKETS ............................................................................................................. 38
4.1 General characteristics ......................................................................................................... 38
4.2 Cold gas propulsion systems ................................................................................................. 38
4.3 Cold gas propellants ............................................................................................................. 39
4.4 Blow-down systems ............................................................................................................... 40
5 – LIQUID PROPELLANT ROCKETS ......................................................................................... 42
5.1 Mono-propellant rockets ....................................................................................................... 42
5.2 Components of a mono-propellant system............................................................................. 42
5.3 Mono-propellant fluids .......................................................................................................... 43
5.4 Mono-propellant chamber temperature ................................................................................ 44
5.5 Bi-propellant rockets ............................................................................................................. 44
5.6 Pressure-fed and pump-fed bi-propellant systems ................................................................ 45
5.7 Bi-propellant fluids ............................................................................................................... 48
5.8 Bi-propellant chamber temperature ...................................................................................... 49
5.9 Molar-averaged fluid properties ........................................................................................... 50
6 – SOLID PROPELLANT ROCKETS ........................................................................................... 52
6.1 General characteristics ......................................................................................................... 52
6.2 Components of a solid propellant rocket ............................................................................... 52
6.3 Regression rate and chamber pressure ................................................................................. 53
6.4 Influence of the grain geometry ............................................................................................. 54
6.5 Hybrid propellant rockets ..................................................................................................... 55
7 – SUMMARY PROBLEMS .......................................................................................................... 57
7.1 Rocketdyne F1 Engine ........................................................................................................... 57
7.2 Space Shuttle propulsion system ........................................................................................... 59
REFERENCES ................................................................................................................................. 62
v
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
1 – Introduction
Aircraft engines and rocket engines are based on two different types of propulsion, with some
similarities but also several big differences. In a nutshell, the most important difference is
that, even if they are both based on accelerating a working fluid, in an aircraft engine the fluid
to be accelerated comes from outside the vehicle, while in a rocket the fluid to be accelerated
is carried on board of the rocket. Actually, in larger rockets and in most of the space
launchers, the mass of the accelerated fluid (usually called propellant), can be up to 80% of
the total rocket mass at launch.
1
Note that “nearly independent” does not mean completely independent: a weak dependence on
altitude and flight velocity might still be present, as it will be shown in the following chapters.
1
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
This happens mainly because both the oxidizer and the fuel have to be carried on board, while
in aircraft engines only the fuel has to be carried on board.
Table 1: Comparison between rocket engines and turbojet engines (adapted from [1]).
Feature Rocket engine or rocket Turbojet engine
motor
T/W, typical 75:1 5:1
Specific fuel/propellant 0,8-1,4 kg/(Nhr) 0,05-0,15 kg/(Nhr)
consumption
Thrust density 375000 N/m2 125000 N/m2
Thrust versus altitude Nearly constant Decreases with
increasing altitude
Thrust versus flight velocity Nearly constant Decreases with
increasing flight velocity
Thrust versus air temperature Constant Decreases with
increasing air
temperature
Flight velocity versus exhaust Unrelated; flight velocity Flight velocity always
velocity can be greater less than exhaust
velocity
Altitude limitation None; suited to space 14-17 km
travel
Adapted from: Rocket Propulsion Elements
Although rocket engines are sometimes used in military applications to propel ballistic
missiles and anti-tank weapons, their main applications are related to space vehicles.
The most typical example are space launchers, which require high flight velocities (in the
order of 7 km/s or higher), but also operate at high altitudes with high thrust levels to
overcome the Earth’s gravity. As an example, the European Ariane 5 launcher is capable of
bringing a payload of about 40 tons into LEO (Low Earth Orbit). The total launcher mass at
lift-off is about 746 tons, of which 642 tons is propellant. The main stage is powered by a
single rocket engine (the Vulcain engine), which produces a vacuum thrust of 1145 kN and
has a nominal burn time of 590 s, and two large booster rockets providing a lift-off thrust of
5500 kN for a burn time of 130 s. Similar to launchers are the applications in ascent and
lander vehicles, from/to the surface of the Moon or other celestial bodies.
In spacecraft applications, rocket engines are used for tasks such as orbit transfer, acquisition,
repositioning, de-orbiting or plane changes, station-keeping, wheel momentum dumping, but
also (in pairs of two) for attitude control. Figure 2 shows some features of the rocket
propulsion system on a specific spacecraft, including engines to provide spin, axial and radial
control.
2
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Other less common but still important applications of rocket engines are sounding rockets,
amateur rockets, ejection seats, rocket-assisted take off on some heavy military transport
aircraft, special race cars, propulsion belts. Table 2 provides an overview of a number of
primary and secondary rocket propulsion applications, with some of their typical
characteristics and requirements.
3
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Various types of components are used to accomplish these different functions. Starting from
the propellant storage function, Figure 4 shows some specific types of tanks, from smaller
ones (left hand side) to the very large liquid propellant tank used for the Space Shuttle, shown
on the right hand side. The central figure also shows that, inside the tank, additional features
can sometimes be present to handle the propellant, for example when it is required to
correctly direct the propellant towards the exhaust and avoid the formation of large bubbles of
gas during flight, especially when the vehicle operates in micro-gravity or in space.
Other components of the system include the propellants and the pressurant gas (when
needed), as well as all propellant handling devices, such as:
pressure regulators, to make sure that the pressure at which the propellant is used in the
thrusters is at the right level (often different to the pressure in the tanks);
filters, to avoid that potentially dangerous solid particles, that could be present in the
propellant, reach the thrusters and thus procure big damage to the propulsion system;
pressure transducers, to measure the pressure at different points of the feed system;
flow regulators, to ensure the right level of propellant mass flow rate, either actively (for
example by means of proportional valves) or passively (for example by means of a
calibrated nozzle);
gas generators, sometimes used to generate the gas which pressurizes the propellant in the
tank, a function that in larger engines can also be accomplished by turbopumps;
a wide range of valves, used for different purposes: we might have fill/drain valves, to fill
the tanks with propellant or to drain it if necessary, check valves, a particular type of valve
letting the propellant flow only in one direction but not in the opposite one, pyrotechnic
valves, opened only once by means of a pyrotechnic device and then always staying open,
thrust valves, the last ones to be opened to make the propellant flow into the thruster where
it is accelerated, usually built with particularly good accuracy and very fast response time.
4
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Rocket engines can also be classified based on the type of energy source used, which can
simply be propellant pressurization (in the so-called cold gas rockets) but also, more
commonly, chemical, nuclear, solar or electric energy.
A typical way to compare and classify different types of propulsion is to put them in a chart
like the one shown in Figure 5, where their performance range is identified in terms of thrust
(on the horizontal axis) and specific impulse2 (on the vertical axis).
The propulsion systems that will be discussed more in detail in this course (coloured boxes in
the figure) are cold gas rockets, characterized by very low specific impulse and low-moderate
thrust levels, and several different types of rockets using chemical energy, such as liquid
mono-propellant rockets (characterized by slightly higher thrust and specific impulse levels
than cold gas ones), solid propellant rockets (comparable specific impulse to mono-
propellant rockets, but significantly higher thrust), and liquid bi-propellant rockets (higher
specific impulse than solid propellant rockets, while still staying at a very high thrust level).
Other propulsion options, not studied in detail by this course, belong to the field of electric
propulsion (either electrostatic, electromagnetic or electro-thermal), typically characterized
by very good specific impulse but low thrust levels, rarely in excess of a few Newton. More
advanced concepts, such as nuclear, solar or even antimatter and laser propulsion, usually
offer an intermediate specific impulse between chemical and electric propulsion, but higher
thrust levels than those achievable with electric propulsion.
Let’s now try to associate the most common types of propulsion to some of the applications
that have been previously discussed, recalling that higher thrust means higher acceleration,
while higher specific impulse means lower propellant mass consumption and thus, as a
consequence, longer possible thrust durations for a given amount of available propellant.
Starting from launchers and missiles, where (see Figure 6) typical requirements ask for high
acceleration and relatively short operational time, we need in this case a propulsion system
providing high thrust levels, but not necessarily very high specific impulse. Therefore, for
2
These performance parameters will be discussed and explained more in detail in the next chapter. For
now, we just mention (in a simplified way) that higher thrust means higher acceleration, while higher
specific impulse means lower propellant mass consumption.
5
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
these applications, solid or liquid bi-propellant rockets are chosen, since they are the options
providing the highest amount of thrust.
In amateur rockets or ascent and lander vehicles, requirements are typically based on
moderate-high acceleration and, again, relatively short operational time. Therefore, in this
case, together with solid and liquid bi-propellant rockets, other options characterized by
slightly lower thrust levels (such as liquid mono-propellant systems) can also be considered.
For orbit insertion and transfer, the main driving requirement is typically asking for moderate
to long, or in some cases even very long, operational time. Therefore, systems with medium
to high specific impulse are chosen, ranging from solid and liquid bi-propellant rockets to, in
many cases, even electric propulsion.
In orbit maintenance and attitude control of satellites, the acceleration provided by the
propulsion system shall not be too high, to allow for precise positioning and pointing of the
spacecraft. This means that any kind of propulsion system allowing for moderate and low
thrust levels can be used; this can be, depending on the specific requirements for the given
spacecraft, electric propulsion, cold gas systems, liquid mono-propellant rockets but even
smaller types of liquid bi-propellant rockets.
Finally, for a more exotic application such as propulsion belts, the acceleration shall not be
too high (to avoid physical issues to the person using the belt) and, usually, short operational
times are required. Therefore, in this case, a propulsion system with intermediate performance
characteristics in terms of both specific impulse and thrust is typically chosen; looking at the
chart in Figure 5, this can usually be a liquid mono-propellant rocket.
6
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
2 - Fundamentals
In this chapter, we will look more in detail at the physics behind the working principles of
rocket propulsion. Two fundamental equations will be derived: the rocket equation and the
thrust equation. Finally, some other important propulsion performance parameters will be
defined and discussed, including the specific impulse.
Also in the case of a system of more bodies, when there are no external forces acting on the
system (meaning that there is only exchange of internal forces between the system bodies),
there cannot be any change in the total momentum and it will remain constant.
Let’s see how this works with a very simple example, see Figure 7. The figure shows a
system made of two bodies that, in this case, are two skaters.
Assume that, initially, the system is at rest, so the skaters are not moving (left hand side of
Figure 7). This means the total momentum of the system is initially zero, because the velocity
of the two bodies (skaters) is zero. Additionally, if we assume that both friction between
skaters and ground and air drag forces are negligible, no external forces act on this system:
the only forces are internal to the system, exchanged between the two skaters. This means that
the total momentum of the system shall always remain equal to zero. Assume now that the
first skater pushes the second one (right hand side of Figure 7): they will start moving in
opposite directions, the first one going left with velocity v1 , the second one going right with
velocity v2.
7
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Since the total momentum of the system shall always be zero and the velocities of the two
skaters are in opposite directions, assuming that their masses are respectively M1 and M2, we
can write:
M2
I system M1 v1 M 2 v2 0 v1 v2
M1
This means, among other things, that the skater with higher mass will have lower velocity.
Why is this example relevant to rocket propulsion? Because a rocket works with exactly the
same principle. Assume that the first skater is the rocket and the second skater is the
propellant: the more propellant mass M2 is expelled and the higher the expulsion velocity v2
of this propellant mass, the higher the final velocity v1 of the rocket will be.
Two dwarfs are on a skateboard, initially at rest. The mass of each dwarf is M1 = 5 kg and
the skateboard mass is M2 = 2 kg. The dwarfs jump off the skateboard with horizontal
velocity (relative to the skateboard) v1 = 1 m/s. We can neglect friction and air drag forces.
(1) What is the final skateboard velocity if the dwarfs jump off together, at the same time?
(2) How does the final skateboard velocity change if the dwarfs jump off separately, one
after each other?
-----------------------------------------------------------------------------------------------------------
If the dwarfs jump off together, we can apply exactly the same equation as in the case of
the two skaters (the first body of the system is the skateboard, the second body are the two
dwarfs jumping off together). The final skateboard velocity v2 is:
2M1 10 kg
M 2 v2 2M1 v1 v2 v1 1 m/s 5 m/s
M2 2 kg
To answer question (2), we need to consider the jumps of the two dwarfs separately. When
the first dwarf jumps off, the two bodies of the system are the first dwarf (mass M1) and the
skateboard with the second dwarf on it (mass M2+M1). The skateboard velocity v2_a is:
M1 5
I in I fin 0 M 2 M1 v2 _ a M1 v1 v2 _ a v1 m/s
M 2 M1 7
When the second dwarf jumps off, we have now to perform the same calculations with a
new two-body system, in which the first body is the skateboard alone, the second body is
the dwarf jumping off. The initial momentum of the system is not zero anymore, since the
skateboard is initially moving at velocity v2_a . Furthermore, when the dwarf jumps off, its
absolute velocity needs to be considered (equal to v1-v2_a ), not the relative velocity v1 with
which it jumps off the skateboard. The final skateboard velocity v2_fin is then:
25
I in I fin M 2 M1 v2 _ a M 2 v2 _ fin M1 v1 v2 _ a v2 _ fin 5 m/s 5 m/s
14
This result, which could in principle be extended (with slightly longer calculations) to the
generic case of N dwarfs, shows that the maximum skateboard velocity is obtained when
the dwarfs jumps off all together, and we always get a lower velocity when they jump off
separately.
When we look at the result of the previous exercise, we can once again find an analogy to
rocket propulsion. If we assume that the skateboard is our “rocket” and the dwarfs are its
“propellant”, we can conclude that, in order to achieve higher rocket (= skateboard) velocity,
we need to expel all the propellant (= dwarfs) as fast as possible, ideally (for maximum
velocity) instantaneously at a time. Thus, it clearly appears that not only the amount of
expelled propellant and its exhaust velocity are important, but also the rate at which the
propellant is expelled. We will better analyse these aspects in the following sections, where
the most important equations for characterizing the performance of a rocket will be derived.
8
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Another assumption we are making here is that the jet velocity is parallel to the flight
velocity (horizontal direction in the case of Figure 8), meaning that the propellant is expelled
in a direction that is parallel and opposite to the flight direction.
Under these assumptions, we can now write the total momentum I+dI of the system (rocket +
expelled propellant) in the direction of flight, at the time t+dt, in the following way:
Note that, when writing the momentum of the expelled propellant, we need to consider its
absolute velocity, not the relative velocity with respect to the rocket (recall the previous
exercise with the two dwarfs on a skateboard). This is the reason why, in the above equation,
we are not simply using the relative jet velocity ve , but the absolute jet velocity ve-v .
With a few simple calculations, it is now possible to write in a more compact way the
momentum change dI:
I dI M dM P v dv dM P ve v
Mv M dv dM P v dM P dv dM P v dM P ve
Which, recalling that the initial total momentum at the time t was simply M∙v, leads to the
following expression:
The term dMP∙dv, in the above expression, is a lower-order term given by the product of two
small quantities and thus, in a first order approximation, can be neglected compared to the
other ones. Now, since we have assumed that no external forces act on this system, the total
momentum of the system cannot change, therefore the variation of total momentum dI shall
be zero:
dI 0 M dv dM P ve 0 M dv dM P ve
9
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
The next step is variables separation, by putting all the velocity variables on the left side of
the equation and all the mass variables on the right side. Also note, looking at Figure 8, that
the propellant mass variation dMP is equal and opposite to the rocket mass variation dM:
dv dM P dM
M dv dM P ve
ve M M
Finally we can integrate the above expression, making the additional assumption that the jet
velocity ve is constant with time:
1 dM 1 M M in
ve dv
M
v ln in
ve M v ve ln
M in M P
fin
Where Min is the initial mass of the rocket and Mfin is its final mass, which is simply equal to
the initial mass Min minus the expelled propellant mass MP.
We have now derived what is commonly called the rocket equation or, alternatively,
Tsiolkovsky equation (from the name of the scientist who first applied it to the study of the
motion of a rocket, in the first years of XX century):
M in
v ve ln (3)
M in M P
This equation gives the total velocity change of a rocket of given initial mass Min when a
propellant mass MP is expelled by the rocket with jet velocity ve , provided that all the
simplifying assumptions that have been made to derive the equation are valid. In summary,
these assumptions are:
1. No gravity forces;
2. No atmospheric drag;
3. No other external forces acting on the rocket;
4. The jet velocity ve is constant with time;
5. The jet velocity ve is parallel and opposite to the flight velocity.
If at least one of the above assumptions is not true, the value calculated by the rocket equation
is not the actual velocity change of the rocket. However the v calculated by the rocket
equation can still be used, even in this case, as a good indicator of the total energy provided
by the propulsion system.
As a sample case when the assumptions made to derive the rocket equation are not valid
anymore, we can mention a launcher that takes off from the surface of the Earth or another
celestial body; in this case, the first two assumptions in the above list are not valid anymore,
because there is certainly a gravity force acting on the rocket and, in most cases, atmospheric
drag too. The fourth assumption in the list (constant jet velocity with time) is not valid when
the operational conditions of the rocket vary during the period considered, which (as we will
see in the following chapters) actually happens in some particular types of rockets. The fifth
assumption, finally, is not valid when the rocket axis is not parallel to the flight direction, as it
usually happens in most launcher trajectories.
One particular example, in which it is still possible to derive a relatively simple equation for
the actual velocity change of the rocket, is the case when a gravity force of constant
gravitational acceleration g is present, but all the other assumptions made to derive the rocket
equation are still valid. In this case, it is possible to show that the actual velocity change of
the rocket is given by:
M in
v ve ln g tb (4)
M in M P
10
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Where tb is the burn time (total amount of time during which the propellant is expelled). The
difference between the velocity change calculated by equation (4) and the one calculated by
the rocket equation is given by the term g∙tb, usually called gravity loss because it is a
measure of how much rocket velocity is lost, with respect to the ideal case, due to the
presence of gravity.
A small rocket ascends from Moon’s surface: thus, no atmosphere is present, and the only
external force is Moon’s gravity (g = 1.62 m/s2, assumed constant).
The rocket flights in vertical direction and the propellant is expelled vertically at constant
jet velocity ve = 1000 m/s. The initial rocket mass is Min = 100 kg, of which 70 kg are
propellant. After a burn time tb = 100 s, all the propellant has been expelled.
What is the v calculated by the rocket equation, and what is the actual velocity change of
the rocket?
-----------------------------------------------------------------------------------------------------------
From the rocket equation:
M in m 100 kg m 100 m
v ve ln 1000 ln 1000 ln 1204
M in M P s 100 kg - 70 kg s 30 s
While the actual velocity change of the rocket is, from equation (4):
M in m m m
vactual ve ln g tb 1204 1.62 2 100 s 1042
M in M P s s s
The difference between these two values is 162 m/s, which is therefore the gravity loss:
this v is needed to overcome the gravity force, and is therefore not used to accelerate the
rocket.
The rocket equation shows us, once again, that to achieve higher flight velocity using rocket
propulsion it is preferable to expel the propellant at the highest possible jet velocity, as clearly
illustrated by the sample case shown in Figure 9.
2000
1800
Launch mass [kg]
1600
1400
1200
1000
800
600
400
0 5 10 15 20
Exhaust velocity [km/s]
Figure 9: Initial launch mass Min calculated by the rocket equation, as a function of the
exhaust velocity ve, for a rocket of required final mass Mfin = 400 kg and v = 5000 m/s.
The jet velocity ve is, however, strongly dependant on the type of propulsion used by the
rocket and its specific acceleration mechanism, as it has briefly been mentioned in the
previous chapter and will be better explained in the following ones.
11
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Figure 10: Pressure forces acting on a rocket flying in the horizontal direction.
On the area Ae where the propellant is expelled, the pressure acting on the rocket is exactly
the pressure pe at which the propellant is expelled, which we assume to be uniform over the
whole surface. On the remaining rocket surface, we have to consider the pressure po of the
fluid surrounding the rocket, for instance air in case the rocket is flying inside Earth’s
atmosphere. This pressure is not necessarily uniform everywhere around the rocket, so the
total force it exerts is calculated by integrating it over the whole remaining external area of
the rocket (indicated as Ao in this case). The total pressure force is therefore:
FP p
Ao
o dAo pe Ae (7)
To write our equations in a more useful way, we recall from fundamental physics that, when
we have a uniform pressure acting on a closed volume, its integral over the surface enclosing
that volume shall be zero. If we now assume that this uniform pressure is the ambient pressure
pa (or, in other words, the pressure of unperturbed atmosphere far enough from the rocket,
where the perturbation given by the rocket is not present anymore), we can simply write that
the integral of this pressure over the surface of the same control volume used for the previous
calculation is zero:
12
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
p
Ao
a dAo pa Ae 0 (8)
Note that equation (8) has no real physical meaning and does not reflect the situation of the
surface around the rocket, where the pressure is not uniform and severely perturbed by the
rocket motion; it is therefore just a mathematical ideal expression, that will be used here in
order to simplify the final equation. This aspect becomes more clear if we combine equations
(7) and (8) in the following way:
FP 0 p
Ao
o dAo pe Ae pa dAo pa Ae
Ao
FP p
Ao
o pa dAo pe pa Ae
The first term on the right hand side of this expression, the integral over Ao of po-pa, can be
easily recognized as the general definition of aerodynamic drag Fdrag, since it gives an
indication of the force exerted over the rocket surface, in the direction of flight, due to the
presence of a perturbed pressure field with respect to the unperturbed ambient pressure. Thus,
we can finally write the expression for the pressure forces in the following way:
p
Ao
o pa dAo Fdrag FP Fdrag pe pa Ae
13
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
at which it is expelled, the more force will be exerted by the propellant to the rocket. The
second term on the right hand side of equation (10) is usually called pressure term, generated
by the difference between the pressure at which the propellant is expelled and the pressure of
the external environment. Note that the pressure term can be either equal to zero, negative, or
positive, depending on whether the propellant exhaust pressure pe is the same, higher or lower
than the ambient pressure pa.
It should also be noticed that the thrust calculated by equation (10) is an external force for the
rocket itself (since it is exerted by the expelled propellant on the rocket), but it is an internal
force for the system made by the rocket and the expelled propellant.
Looking at equation (10), we can also see that there is no explicit dependence of the thrust on
the flight velocity. Conversely, there is a clear dependence on the altitude at which the rocket
is flying, through the term pa present in the equation. As a consequence, the thrust increases
with the altitude, since an increase in altitude is directly associated to a decrease in ambient
pressure (see Figure 11).
Figure 11: Ambient pressure variation with altitude (standard atmosphere values).
This is the reason why, when we indicate the thrust of a given propulsion system, we also
need to specify at which altitude it is referred. Most typical cases are the sea level thrust,
corresponding to an ambient pressure pa = 101325 Pa, and the vacuum thrust, corresponding
to an ambient pressure pa = 0. An example is shown in Figure 12 for the SSME (Space Shuttle
Main Engine), where it is evident that, indeed, there is a significant change in rocket thrust at
altitudes from 0 to about 30 km, above which the thrust stays approximately constant.
14
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
A rocket engine delivers a thrust equal to 5 kN when the ambient pressure is 105 Pa, and a
thrust of 5.5 kN when the ambient pressure is 5∙104 Pa. How much is the vacuum thrust of
this engine (i.e., the thrust it delivers when the ambient pressure is zero)?
-----------------------------------------------------------------------------------------------------------
Assuming that the propellant mass flow rate and the jet velocity do not vary with altitude,
as it usually happens in most propulsion systems of practical interest, the difference
between the thrust levels FT1 and FT2 at the two given altitudes is:
FT m ve pe pa Ae FT 2 FT 1 pa1 pa 2 Ae
In order to write the rocket thrust in a more compact way, just as the product of a mass flow
rate times a velocity, it is possible to define a fictitious performance parameter usually called
equivalent jet velocity veq (or, sometimes, effective jet velocity), which includes both the
momentum and the pressure term present in the thrust equation:
veq
FT p pa Ae F mv
ve e (11)
T eq
m m
The equivalent jet velocity is a pure mathematical entity, which does not have any real
physical meaning. It is exactly equal to the jet velocity when the pressure term in the thrust
equation is identically zero or, in other terms, the propellant exit pressure is exactly equal to
the ambient pressure. It is possible to show that, taking into account the effects of the pressure
term in the thrust equation, the rocket equation should be written using the equivalent jet
velocity instead of the “pure” jet velocity. For this reason, equation (3) should be replaced by
the following, more correct one:
M in
v veq ln (12)
M in M P
I tot FT dt (13)
0
In particular, when thrust is constant over time, the total impulse is simply equal to:
Itot FT tb (14)
We can also consider the case when the equivalent jet velocity is constant over time. In this
case, recalling from equation (11) that thrust is simply mass flow rate times equivalent jet
velocity, it is possible to write:
15
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
tb tb
dM P
I tot mveq dt veq dt veq M P (15)
0 0
dt
A rocket engine delivers during a burn a total impulse Itot = 104 Ns. Assuming that the
engine is characterized by an equivalent jet velocity (constant during the burn) of 2500
m/s, how much propellant has been expelled during this burn?
-----------------------------------------------------------------------------------------------------------
From equation (15), since the equivalent jet velocity is assumed constant:
I 104
I tot veq M P M P tot 4 kg
veq 2500
Based on these considerations, we can now define the gravimetric specific impulse, or simply
specific impulse Isp of a rocket, as a quantity proportional to the total impulse divided by the
total mass of propellant used to generate that total impulse:
tb
1 0
FT dt
I sp (16)
g 0 tb
m dt
0
Note that the term g0 in equation (16) is always the gravitational acceleration on Earth at sea
level (equal to 9.81 m/s2), and is therefore a constant, whatever is the environment and the
gravity level into which the rocket or spacecraft is operating.
In case the equivalent jet velocity is constant over time, equation (16) can be significantly
simplified and becomes:
tb tb
Which gives the following, very commonly used simplified expression for the specific
impulse:
veq
I sp (17)
g0
valid in all cases when we can consider the equivalent jet velocity as a constant over the
whole burn time tb.
The specific impulse is probably the most important figure of merit for rocket performance.
From the definition shown in equation (16), we can note that higher specific impulse means
that a larger momentum change (or total impulse) can be generated with the same amount of
propellant or, in other words, the same momentum change can be obtained with a smaller
amount of propellant. The specific impulse is therefore directly associated to the propellant
mass consumption of the rocket.
Since, from equation (17), the specific impulse is proportional to the equivalent jet velocity
and the same applies to the v and the rocket thrust, see equations (11) and (12), higher
specific impulse means, in principle, higher thrust and higher v. All these considerations,
taken together, show very clearly why the specific impulse is such an important performance
parameter for rocket propulsion.
As we have already seen for the thrust and the equivalent jet velocity, the specific impulse is a
function of altitude too, increasing when altitude increases. This means that, as for the thrust,
16
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
when we indicate the specific impulse of a given propulsion system we always need to
specify the altitude to which it is referred.
But for which reason is the constant g0 included in the definition of specific impulse, equation
(16)? This is mainly due to historical reasons and, more specifically, to the fact that without
this constant the unit for the specific impulse would be a velocity and, thus, it might be
expressed in several different units (m/s, inches/s, feet/s, etc.) depending on the unit
convention used. By including g0, the specific impulse has the unit of a time and, therefore, is
more universally expressed in seconds.
Two rocket engines A and B are characterized by the following values of the mass flow
rate and jet velocity:
Rocket A: mass flow rate = 100 kg/s, jet velocity = 3000 m/s
Rocket B: mass flow rate = 200 kg/s, jet velocity = 1500 m/s
Assuming that the two rockets are working at a nozzle exit pressure equal to the ambient
pressure and their jet velocity is constant over time, which rocket provides more thrust?
Which one has higher specific impulse?
-----------------------------------------------------------------------------------------------------------
Under the given assumptions, the rocket thrust and specific impulse are simply given by:
FT m ve
ve
I sp
g0
This leads to the same thrust for the two rockets, equal in both cases to 300 kN,
Conversely, the specific impulse is equal to 306 s for the first rocket and 153 s for the
second one. Rocket A has therefore higher specific impulse.
It is also possible to define a volumetric specific impulse I, exactly in the same way as the
gravimetric specific impulse but considering, in this case, the propellant volume instead of its
mass:
tb
1 0
FT dt
I (18)
g 0 tb m
dt
0
where is the mass density of the propellant. Higher volumetric specific impulse means that
a given momentum change can be obtained using a smaller volume of propellant and, thus,
the volumetric specific impulse is an important figure of merit to define the size of a rocket.
A rocket is characterized by a specific impulse of 200 s and equivalent jet velocity of 1962
m/s on Earth at sea level (g = 9.81 m/s2). What is the specific impulse of the same rocket:
(1) On the surface of a planet with gravity = 50% of Earth’s gravity at sea level?
(2) In an ideal case where no gravity is present (g = 0)?
Assume, for simplicity, that the equivalent jet velocity of the rocket (constant over time)
remains the same in all the cases considered in this problem.
-----------------------------------------------------------------------------------------------------------
Under the given assumptions, the specific impulse is given by equivalent jet velocity
divided by g0, which is a constant and is always equal to 9.81 m/s2. Thus, the specific
impulse of the rocket will always be the same (200 s), under any gravity conditions.
17
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Finally we can consider another power parameter, the absolute jet power Pjet_abs , which is a
measure of the amount of jet power not effectively used for generating thrust and, thus, lost in
the expelled propellant. The absolute jet power is basically the same as the jet power, using
however the absolute jet velocity instead of the jet velocity relative to the rocket:
By using the definitions of thrust power and absolute jet power given in equations (21) and
(22), it is possible to write in a more interesting form the propulsive efficiency expression:
v
2 0
PT m veq v0 2veq v0 veq
P 2 (24)
PT Pjet _ abs m veq v0 m veq v0
1 2 veq v0
2
v
2
1 0
2 v
eq
Equation (24) shows very clearly that the propulsive efficiency is a function of the ratio of
rocket flight velocity to equivalent jet velocity. If we plot this function and compare it to what
typically happens in air-breathing engines (see Figure 13), we note that there is a maximum
for the propulsive efficiency, equal to exactly 100%, when the equivalent jet velocity is equal
to the flight velocity. Furthermore, differently to what happens in air-breathing engines,
rocket propulsion can also work when the flight velocity is higher that the equivalent jet
velocity.
18
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Figure 13: Propulsive efficiency of rocket and air-breathing jet engines [2].
Finally, we can combine the two efficiencies defined so far in order to obtain the overall
energy efficiency of the rocket:
P C (25)
19
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
FT m ve pe pa Ae veq
FT
ve e
p pa Ae
m m
M in veq
v veq ln I sp
M in M P g0
Let’s now take a closer look at the variables in these equations, and the exact role played by
the physical quantities they represent:
The pressure pa depends on the ambient conditions (atmospheric pressure of the
environment in which the rocket is operating);
The nozzle exit area Ae and the masses Min and MP depend on the design of the rocket (or
spacecraft) and on the nozzle geometry;
The jet velocity ve , the nozzle exit pressure pe and the mass flow rate m depend on the
way how the propellant is heated and accelerated in the rocket.
The objective of this chapter is to evaluate the three quantities listed in the last bullet point
above (jet velocity, nozzle exit pressure, mass flow rate) by means of a set of simplified
equations describing the flow in a rocket nozzle. In order to do it we will use the so-called
Ideal Rocket Theory (IRT), starting from a simplified rocket geometry and a set of physical
assumptions, illustrated in the next section.
20
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
An important parameter to characterize the nozzle geometry, which will be widely used in the
following, is the expansion ratio , simply defined as the ratio of nozzle exit area to nozzle
throat area:
Ae
(26)
A*
As previously mentioned, the Ideal Rocket Theory is also based on a set of physical
assumptions, necessary to simplify the physics of the nozzle flow at a sufficient level to
derive meaningful equations for the parameters which characterize it. These assumptions can
be summarized as follows:
1. The propellant in the chamber and in the nozzle is a perfect gas.
2. The propellant in the chamber and in the nozzle is also a calorically ideal gas, meaning
that its specific heats are not dependent on temperature.
3. The propellant is homogeneous and has uniform chemical properties in every point of the
chamber and the nozzle.
4. The flow in the nozzle is steady, meaning that all the quantities that characterize it are not
dependent on time.
5. The flow in the nozzle is isentropic: no external energy is provided and, in the same way,
no energy is lost by the propellant.
6. The flow in the nozzle is one-dimensional, meaning that all quantities vary only along the
nozzle axis and do not vary in the other two directions.
7. The flow velocity in the nozzle is purely axial: everywhere in the nozzle, the velocity is
parallel to the nozzle axis.
8. No external forces act on the propellant flowing in the nozzle; in particular, there is no
friction between the propellant and the nozzle walls.
9. The propellant in the chamber has negligible velocity, allowing in all equations to assume
the chamber velocity vC equal to zero.
21
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
This long list of assumptions clearly indicates that the final equations we are going to derive
in this chapter describe in an extremely simplified way the actual flow in a rocket nozzle.
They can certainly be used for an initial characterization of the rocket performance, but in
defining its real performance an important role will also be played by all the aspects neglected
as a consequence of the above assumptions.
To derive the mass conservation equation, we can observe that the propellant mass flow rate
in any section of the nozzle can be written as the product of gas density, times its velocity,
times the frontal area of the nozzle section through which it is flowing. Since no additional
mass of propellant is injected in the nozzle or extracted from it, no mass flow rate variations
are possible and, thus, the mass flow rate variation through our small portion of nozzle shall
also be zero:
d v A 0 (27)
This equation can also be integrated, in order to write it in a form that applies to the entire
nozzle and not only to the small portion we have considered here:
v A m constant (28)
For what concerns the momentum conservation equation, we start by recalling that the
propellant momentum per unit time is the product of its mass flow rate times its velocity.
Now, since one of our initial assumptions is that no external forces act on the propellant at
any place in the nozzle, the only remaining forces that can cause momentum variations are the
pressure forces. In other terms, we can write that the variation of momentum in our small
22
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
portion of nozzle is balanced by the difference between the pressure forces acting on its inlet
and outlet sections:
d m v p A p dp A m dv A dp
where we have neglected the term dp∙dA since it is a lower-order term given by the product of
two small quantities and, thus, negligible compared to the other ones.
Finally, for what concerns the energy conservation equation, we recall from our initial
assumptions that the propellant does not exchange any energy with the external environment.
This means that there cannot be any variations in the total enthalpy, defined as the enthalpy
plus the kinetic energy per unit mass:
1
d h v2 0
2
By slightly elaborating this expression, we obtain the energy conservation equation:
dh v dv 0 (31)
or, in its integral form:
23
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
1
h v 2 constant (32)
2
Another set of Ideal Rocket Theory building blocks is given by the ideal gas equations.
Since we have assumed that the gas in the nozzle can be considered a perfect gas, it is
possible to characterize it by means of the following equation of state:
RA
p T (33)
MW
where MW is the molecular mass of the gas and RA is the universal gas constant, equal to 8314
J/(K∙kmol). We will also extensively use the specific heat relationships for an ideal gas, which
relate the constant pressure specific heat cp and the constant volume specific heat cv ,through
the universal gas constant and the molecular mass. An important role in these equations is
played by the specific heat ratio , which is simply defined as the ratio of constant pressure
to constant volume specific heat.
RA
c p cv (34)
MW
cp
(35)
cv
RA
cp (36)
1 MW
In addition, the enthalpy for an ideal gas is simply the product of its constant pressure specific
heat times its temperature:
h cp T (37)
which, since we have assumed that the gas is calorically ideal and, thus, its constant pressure
specific heat is constant with temperature, can also be written as:
dh c p dT (38)
Another one of our initial assumptions was that the flow in the nozzle is isentropic. For an
ideal gas, this leads to the following relationships between its pressure, density and
temperature:
p
constant (39)
1
pT constant (40)
We can also recall the definition of Mach number M, that is simply the ratio of the gas
velocity to its speed of sound at the given flow conditions:
v
M (41)
a
where the speed of sound a is defined as:
dp
a2 (42)
d Constant entropy
It is possible to show that, for an ideal gas, this leads to the following equation for the speed
of sound:
24
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
RA p
a2 T (43)
MW
Figure 17 shows, in summary, all the equations that we have introduced as building blocks for
the Ideal Rocket Theory.
dp d v 2 dv dv dA d dv dA dv dv
a2 v dv 2 M 2 M2 0
d Const. entropy dp a v v A v A v v
M 2 1
dA dv
(44)
A v
25
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
This equation gives a relationship between the variation of nozzle frontal area, the Mach
number and the variation of velocity, and is valid everywhere in the nozzle. Applying it to all
possible situations in the convergent or the divergent part of the nozzle, we obtain the
overview shown in Figure 18.
Figure 18: Velocity variations in the convergent and divergent sections of a nozzle.
In the convergent part of the nozzle the frontal area decreases and, thus, the term dA is
negative. This means that we can have a velocity increase only if the Mach number is lower
than one or, which is the same, only if the flow in the convergent is subsonic. We have an
exactly opposite situation in the divergent, where the frontal area increases and, thus, dA is
positive: in this case, we can have a velocity increase only if the Mach number is higher than
one and, thus, only if the flow in the divergent is supersonic.
In order to accelerate the flow in the entire nozzle, we need to achieve a situation in which the
term dv is always positive. As a consequence, recalling that our assumed chamber conditions
correspond to a velocity (and, thus, Mach number) equal to zero, the flow at the convergent
inlet is certainly subsonic; therefore, it is accelerated along the convergent part of the nozzle.
If this acceleration leads to a Mach number at the throat exactly equal to one, the flow will
continue to be accelerated also in the divergent and a Mach number higher than one will be
achieved at the nozzle exit. On the other hand, if the Mach number at the throat is not yet
equal to one, the flow at the divergent inlet will still be subsonic and, thus, the velocity will
decrease along the divergent and a Mach number higher than one will never be achieved.
This means that, in order to have a properly operational rocket in which the flow continues to
be accelerated from the combustion chamber to the nozzle exit, the nozzle needs to be
designed in such a way that a Mach number exactly equal to one is achieved at the throat (or,
in other words, the nozzle is characterized by a sonic throat).
26
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Figure 19: Nozzle geometry and chamber conditions used for the derivation of the
nozzle flow equations.
For our derivations we start from the energy conservation, equation (32), combined with the
enthalpy expression for an ideal gas, equation (37):
1 1
h v 2 constant c pT v 2 constant
2 2
Since the total enthalpy shall be constant in every section of the nozzle, we can equate its
value at the nozzle inlet (combustion chamber) with the value at a generic nozzle section
where the flow is characterized by temperature T and velocity v. Recalling that the flow
velocity at the nozzle inlet (combustion chamber) is assumed to be zero, this leads to:
1
c pTC 0 c pT v 2 v 2c p TC T (45)
2
This expression, together with equation (36) for the constant pressure specific heat and
equation (40) for the isentropic flow, allows to find a useful expression for the flow velocity v
in a generic nozzle section:
RA 2 RA 2 RA T
cp v TC T TC 1
1 MW 1 MW 1 MW TC
1
1 T p
pT constant (46)
TC pC
1
2 RA p
v TC 1
pC
(47)
1 MW
In particular, if we write equation (47) at the nozzle exit section, we obtain an equation for the
jet velocity:
1
2 RA p
ve TC 1
pC
e
(48)
1 MW
The jet velocity is therefore a function of the specific heat ratio, the molecular mass, the
chamber temperature and the ratio of chamber pressure to nozzle exit pressure, often also
called “nozzle pressure ratio”. In particular, equation (48) shows that to obtain higher jet
velocity (which, recall, plays a fundamental role in the thrust, specific impulse and v
equations), we can either:
27
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Using again equation (36) for the constant pressure specific heat and equation (43) for the
speed of sound, it is possible to write this equation as:
RA 1 R a2 T 1 v2 1 2
cp c pT A T C 1 2 1 M
1 MW 1 MW 1 T 2 a 2
From which, using equation (40) for the isentropic flow:
1
p T 1 1 1
pT constant C C 1 M2
p T 2
28
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
In this way, we have derived two useful equations that allow for calculating the pressure and
temperature at any generic section of the nozzle, as functions of the chamber pressure and
temperature, the specific heat ratio and the Mach number M at that nozzle section:
TC 1 2
1 M (50)
T 2
pC 1 1
1 M2 (51)
p 2
In particular, when written at the nozzle throat (where the flow is sonic and the Mach number
is exactly equal to one), these equations become:
TC 1 1
*
1 (52)
T 2 2
pC 1 1
(53)
p* 2
Equations (52) and (53) show that, for given chamber pressure and chamber temperature, the
pressure and temperature at the nozzle throat depend only on the specific heat ratio of the gas
flowing in the nozzle.
A C 1 M W MW 1 C pC
Using the isentropic flow equation (39) and the equation of state of an ideal gas, equation
(33), we obtain:
1
p p
constant
C pC
RA pC
pC C TC C
MW R A
TC
MW
1
2
m pC 2 p p
1
pC
(54)
A RA 1 pC
TC
MW
Remember, now, that the mass flow rate is constant along the whole nozzle. This means that
the mass flow rate calculated by means of this equation shall be the same at any nozzle
section and, in particular, at the nozzle throat with p=p* and A=A*.
From equation (54) for a sonic throat:
29
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
pC 1 1 p* 1 1
p* 2 pC 2
And, using this expression in the previous equation for the mass flow rate, written at the
nozzle throat:
1
2
2
2 1 1 1
1
m pC 2 p* p *
pC
1
pC 1 2
1
A* RA 1 pC RA
2
TC TC
MW MW
1
1 1
(56)
2
where is a function of the specific heat ratio, often called Vandenkerckhove function,
very widely used in the theory of rocket propulsion.
Equation (55) shows that, for given nozzle inlet conditions (specific heat ratio, molecular
mass, chamber temperature and chamber pressure) and given nozzle throat area, there is only
one possible value of mass flow rate that makes sonic throat conditions possible. When this
mass flow rate value is achieved we say that we are in presence of a chocked flow, because
the mass flow rate is “fixed” and cannot be increased anymore.
30
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Ae
(58)
A*
2 1
2 pe pe
1
1 pC pC
Equation (58) provides, for a given value of specific heat ratio , a direct relationship between
the nozzle expansion ratio and the nozzle pressure ratio. However, it is unfortunately an
implicit equation for the nozzle pressure ratio, meaning that it is not possible to find a closed-
form solution for it as a function of the nozzle expansion ratio; in this case, the equation needs
to be solved either graphically or numerically.
Figure 20 shows the results obtained from equation (58), in terms of nozzle expansion ratio as
a function of pressure ratio for some given values of specific heat ratio. It is possible to note,
in particular, a generally weak dependence of the results on the specific heat ratio, with the
three curves in the plot very close to each other. We also notice that higher expansion ratio
leads, as it could be expected, to higher pressure ratio which, in turn, means lower nozzle exit
pressure for a given value of chamber pressure. Thus, for given initial chamber pressure, the
higher the nozzle expansion ratio, the lower the pressure that will be achieved at the nozzle
exit.
Figure 20: Nozzle expansion ratio as a function of pressure ratio, for different values of
specific heat ratio.
Depending on how the nozzle exit pressure compares to the ambient pressure, three different
situations are possible for the flow expansion in a supersonic nozzle, see Figure 21:
When the nozzle exit pressure is lower than the ambient pressure (case 1 in Figure 21), the
flow in the nozzle has expanded more than necessary and we thus call it an over-
expanded nozzle.
31
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
When the nozzle exit pressure is exactly the same as the ambient pressure (case 2) we have
an adapted nozzle, because its shape is perfectly adapted to the altitude and, thus, the
ambient pressure at which the rocket is flying.
When the nozzle exit pressure is higher than the ambient pressure (case 3), the flow in the
nozzle has not expanded enough and we call it an under-expanded nozzle.
In case the nozzle is over-expanded or under-expanded, the flow pressure usually adapts to
the external pressure through a set of shock waves outside the nozzle, as schematically shown
in Figure 21.
A rocket nozzle is characterized by chamber pressure pC = 20 bar and specific heat ratio
= 1.3. The nozzle expansion ratio is = 10. What are the expansion conditions of this
nozzle at sea level (pa = 1 bar) and in vacuum (pa = 0)?
-----------------------------------------------------------------------------------------------------------
For a nozzle expansion ratio of 10 and a specific heat ratio of 1.3, we can graphically find
from Figure 20 that the nozzle pressure ratio is equal to 100. Thus:
pC p 20 bar
100 pe C 0.2 bar
pe 100 100
This means that the nozzle is over-expanded at sea level (0.2 bar < 1 bar), and under-
expanded in vacuum (0.2 bar > 0 bar).
1
2 p
p p A*
FT pC A
*
1 e
(59)
1 pC e a
2
1
2 pe pe
1
1 pC pC
32
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Using equation (59) it is possible to show that, for given values of chamber pressure, ambient
pressure, specific heat ratio and nozzle throat area, the thrust is maximum when the nozzle is
adapted, so when the exit pressure is exactly equal to the ambient pressure.
Thus, in summary, putting together all the information we have found on the thrust produced
by a rocket, we can conclude that:
For fixed nozzle geometry and variable ambient pressure, thrust is maximum in vacuum
(zero ambient pressure), see equation (10).
At fixed ambient pressure (or, in other words, fixed altitude), thrust is maximum for the
particular nozzle geometry that leads to an adapted nozzle (nozzle exit pressure equal to
the ambient pressure), see equation (59).
A rocket engine is characterized by vacuum thrust FTvac = 180 kN , nozzle throat area A* =
0.1 m2, nozzle expansion ratio Ae /A* = 50, chamber pressure pc = 106 Pa.
Assuming that the nozzle expansion ratio and the specific heat ratio stay constant, how
much vacuum thrust will be obtained if the throat area is increased to A* = 0.2 m2 and the
chamber pressure becomes pc = 2∙106 Pa?
-----------------------------------------------------------------------------------------------------------
From the vacuum thrust equation:
1
pC A* 2 RA p
FTvac m ve pe Ae TC 1 e
pe Ae
RA 1 M C
p
TC W
MW
1
FTvac 2 p
pe Ae
1 e
pC A *
1 pC pC A*
Note that the term on the right hand side of the above expression stays constant under the
given conditions, because the nozzle pressure ratio remains the same when the nozzle
expansion ratio and the specific heat ratio don’t change, see equation (58). Thus:
33
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
For the characteristic velocity, using the mass flow rate equation (55), it is also possible to
write it as:
pC A* pC A* 1 RA
m c* TC (63)
RA m MW
TC
MW
Equation (63) shows that the characteristic velocity is only a function of the combustion
chamber conditions (chamber temperature, specific heat ratio, molecular mass) and does not
depend on the nozzle geometry. Therefore, the characteristic velocity is a very important
parameter to measure the performance of the propellant used by the rocket: in particular,
higher characteristic velocity generally means higher specific impulse and, thus, a more
performant propellant.
For the thrust coefficient, using the rocket thrust equation (59), it is possible to write:
1
FT 2 pe pe pa
CF
1 (64)
pC A* 1 pC pC pC 2
1
2 pe pe
1
1 pC pC
Equation (64) shows that, although the thrust coefficient is also a function of the specific heat
ratio and the ratio of ambient pressure to chamber pressure, it is mainly dependent on the
nozzle geometry through the nozzle pressure ratio. Thus, the thrust coefficient is a very good
measure of the nozzle geometry effects on the thrust.
The plot in Figure 22 gives a very good idea of the way how the thrust coefficient varies as a
function of several parameters. Note that the plot is obtained for a fixed specific heat ratio
value, equal to 1.3. It clearly shows that, for a given value of the ambient pressure to chamber
pressure ratio, the thrust coefficient has a maximum exactly corresponding to the case of an
adapted nozzle (red curve in the plot). It can also be observed that typical values of the thrust
coefficient are in the range of 1.2 to 1.9.
Figure 22: Thrust coefficient as a function of the nozzle expansion ratio, for fixed
specific heat ratio and different values of the ambient/chamber pressure ratio.
34
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
From the above considerations, it should now be clear why the characteristic velocity and the
thrust coefficient are important parameters in defining the performance of a rocket engine.
Looking at the way how we have expressed the rocket thrust in equation (62), it is evident
that it allows to separate the contributions of the propellant (taken into account by the
characteristic velocity) and the nozzle geometry (taken into account by the thrust coefficient).
The LR91 rocket motor (2nd stage of the Titan 23G launcher) is characterized by a
chamber pressure pC = 5.6 MPa, mass flow rate = 132.9 kg/s, throat area A* = 0.041 m2,
nozzle expansion ratio = 49, adapted nozzle thrust (when pa = pe ) of 401491 N.
Calculate the characteristic velocity and the thrust coefficient at adapted nozzle conditions.
-----------------------------------------------------------------------------------------------------------
Characteristic velocity:
FT _ AdaptedNozzle 401491 N
CF _ AdaptedNozzle 1.75
pC A *
5.6 106 Pa 0.041 m2
For a rocket engine it is known that the vacuum thrust is equal to 106 N and the vacuum
thrust coefficient is CF_vac = 1.8. It is also known that the molecular mass of the gas in the
nozzle is Mw = 20.8 g/mol, the specific heat ratio of the gas in the nozzle is = 1.3, and the
mass flow rate is 310 kg/s.
What is the combustion chamber temperature of this engine?
-----------------------------------------------------------------------------------------------------------
For the Vandenkerckhove function, from equation (56):
1 11.3
1 1 1 1.3 11.3
1.3 0.66726
2 2
Then, from the definition of thrust coefficient, equation (61):
FT FT _ vac 106
CF p A*
5.55 105 N
pC A* C
CF _ vac 1.8
And from the expression for the mass flow rate, equation (55):
2
pC A* M p A*
m TC W C
RA
TC
RA m
MW
2
20.8 5.55 105
0.66726 3577 K
8314 310
35
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
1
pC A*
2 RA
p m
ve TC 1 e RA
1 MW pC TC
MW
Ae
A* 2
1
2 pe p
1
pC
e
1 pC
Looking more in detail at the equations in Table 4, it can be observed that they all depend on
the same parameters. Two of these parameters, the nozzle throat area A* and the nozzle exit
area Ae , are functions of the nozzle geometry. The other parameters are the specific heat ratio
, the molecular mass MW, the chamber temperature TC and the chamber pressure pC , which
depend on:
The type of propellant used by the rocket (liquid mono-propellant, liquid bi-propellant,
solid, hybrid, etc.);
The type of rocket considered (cold gas, chemical rocket, electro-thermal etc.).
In the next chapters some specific examples of rockets and propellants will be discussed,
looking in particular at how chamber conditions and propellant characteristics can be
calculated for these different types of rockets.
Recalling the classification of rocket engines already introduced in chapter 1, it is now
possible to discuss this classification more in detail starting from the additional information
we have found on the rocket thrust and specific impulse. In particular, we have found that the
specific impulse is proportional to the equivalent jet velocity, while the thrust is proportional
to the mass flow rate and the equivalent jet velocity. This means, see Figure 23, that even if
the equivalent jet velocity is high, the thrust is not necessarily high too, because it is still
possible that the mass flow rate is very small. A very good example is electrical propulsion
(blue area in Figure 23), typically characterized by high specific impulse but very low thrust
levels. The opposite applies to chemical propulsion (green area in Figure 23), characterized
by lower specific impulse (and, thus, lower equivalent jet velocity) but the possibility to
achieve very high thrust levels, since the mass flow rate can be very high in these rockets.
Two rockets (rocket A and rocket B) have the same specific impulse.
The mass flow rate of rocket A is 10 times higher than mass flow rate of rocket B.
What is the relationship between the thrust FTA of rocket A and FTB of rocket B?
-----------------------------------------------------------------------------------------------------------
Since the two rockets have the same specific impulse, they also have the same equivalent
jet velocity. Thus, for the thrust:
FT m veq
36
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
37
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Figure 24: Example (left) and cut-off drawing (right) of typical cold gas thrusters.
Cold gas rockets are typically used for attitude control tasks, especially in small and very
small satellites. Their specific impulse is typically not higher than 100 s, with a thrust varying
from the N-mN level, for miniaturized thrusters, to the 1-100 N range, for larger ones.
38
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
A pressure regulator, needed because the required pressure at which the propellant is
used in the thruster is usually much lower than the very high pressure at which the gas is
stored in the tank;
Finally a low pressure isolation valve, opened to allow the propellant flow into the
thruster where it is accelerated and the thrust force is generated.
A cold gas propulsion system is usually very compact, especially when it is intended for use
in miniaturized satellites where volume constraints are even more important than mass ones.
39
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
The gravimetric and volumetric specific impulse of the same propellants, calculated under the
same conditions as in Table 5, are shown in a comparative way in Figure 26. It can be
observed that the propellants with lower molecular mass are also the ones with higher specific
impulse, not surprising if we recall that, from the Ideal Rocket Theory, the jet velocity is
inversely proportional to the square root of the molecular mass, see equation (48). However, it
can also be noticed that propellants characterized by higher molecular mass are typically
leading to a higher value of the volumetric specific impulse, as a direct consequence of the
fact that higher molecular mass normally means higher density and, thus, the same mass of
propellant can be stored in a much smaller volume.
Figure 26: Gravimetric and volumetric specific impulse of typical cold gas propellants.
40
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Since the tank volume is constant in the case of non-flexible walls, it is also possible to write
equation (65) in the following way:
pT
constant (66)
MT k
where MT is the mass of gas in the tank.
The tank of a cold gas rocket is initially at a pressure of 40 bar and temperature of 400 K.
Assuming an isothermal blow-down process, what are the final tank pressure and
temperature when 3/4 of the initial propellant mass has been expelled?
-----------------------------------------------------------------------------------------------------------
Final tank pressure:
pT p pTf M Tf 1
constant T0 pTf pT0 pT0 10 bar
MT M T0 M Tf M T0 4
Since the process is isothermal, the tank temperature remains constant and the final tank
temperature is therefore still equal to 400 K.
41
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Figure 28: Example (left) and schematic (right) of a liquid mono-propellant thruster.
42
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
A filter, used to avoid that any solid impurities present in the propellant eventually reach
the thruster and the nozzle;
In some cases a pressure regulator might be used (not shown in the figure), when the
pressure in the propellant tank is very high compared to the required pressure in the
thruster.
43
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Assuming now that no thermal losses occur and the entire decomposition power is directly
used to heat the decomposition products, we can equate equations (67) and (68) and obtain an
expression for the chamber temperature:
HD
Pheat Pdec m H D m c p TC T0 TC T0 (69)
cp
1 RA 1 8314 J K kmol m
c* TC 1544 K 1613
MW 0.679 10.7 g mol s
44
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
enter in contact with each other, without any additional input energy provided, for example,
by means of an igniter to accelerate the combustion process.
An important parameter for the characterization of these engines is the oxidizer to fuel ratio
O/F , usually also called mixture ratio, defined as the ratio of the mass flow rate of oxidizer
to the mass flow rate of fuel:
moxidizer
O /F (70)
m fuel
Note that, since the oxidizer and the fuel flow together in the combustion chamber at the same
time, the mixture ratio can also be written as the ratio of the total mass of oxidizer to the total
mass of fuel burned by the engine in a given amount of time.
Bi-propellant rocket engines are used in a wide range of applications ranging from launchers
to missiles, amateur rockets, orbit insertion and orbital maneuvering. They are typically
characterized by the highest possible specific impulse values among chemical rockets,
normally in the range from 300 to 450 s, and are usually intended for high or very high thrust
levels, in the kN or even the MN range. The left hand side of Figure 30 shows a picture of the
Space Shuttle Main Engine, where the divergent part of the nozzle and the feeding system are
particularly visible. The right hand side of Figure 30 shows a schematic of the engine parts
close to the combustion chamber in another bi-propellant engine, including the convergent,
the throat and the first part of the divergent section of the nozzle. The fuel (shown in red) and
oxidizer (shown in yellow) are injected in the combustion chamber through an injector,
generating a fine spray of fluid and ensuring proper distribution and mixing of the propellants
in the combustion chamber. In this particular case, before being injected in the combustion
chamber, the fuel flows through some coolant channels around the nozzle and combustion
chamber walls.
Figure 30: The Space Shuttle Main Engine (left, courtesy NASA); schematic of the
Aestus liquid rocket thrust chamber (right, courtesy EADS).
45
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
As usual, pressure sensors are used to measure the pressure in the tank or in other points of
the system, fill/drain valves are attached to the two propellant tanks and to the pressurizing
gas tank, and filters are included to avoid that solid particles reach the thruster and eventually
block the nozzle. Finally, the propellant valves (also called thrust valves) are opened to make
the propellants flow in the combustion chamber and start producing thrust.
In the case of a pump-fed system (Figure 32), the propellants are stored in their tanks at a
relatively low pressure and, although sometimes a small quantity of pressurizing gas (helium
in the example shown in the figure) is used to initiate the propellant flow, the main pressure
increase is provided to the propellant by means of a turbo-pump system. The hot gas turbine
can be driven in several different ways; in the case shown in the figure this is done by a gas
generator, in which combustion of a small quantity of propellants is used to generate the hot
gases to drive the turbine. The turbine drives the fuel and oxidizer pumps, either directly on
the same shaft or through a gearbox as in the case shown in the figure. Finally, also in this
case, thrust valves are opened to allow the pressurized propellants flow in the combustion
chamber.
46
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
As already mentioned, pump-fed systems are particularly convenient for use in large rockets,
because in this case a pressurizing gas tank would be very large and heavy, given the size of
the engine and the amount of propellants used.
Figure 33 shows, as an example, the turbo-pump assembly of the Viking rocket engine used
in the Ariane 4 launcher. It is a single-shaft turbo pump in which, similarly to what shown in
the schematic of Figure 32, both the fuel and the oxidizer pump are driven by the same
turbine. In this particular case the pumps run at a rotational speed of 10000 rpm and require a
power of 2.5 MW. The figure also shows the gas generator, which in this case provides hot
gas to the turbine at a mass flow rate of 1.2 kg/s.
In most practical cases, turbo-pumps are built in the same body of the combustion chamber
and the nozzle, to form a single assembly. This can be clearly seen, as an example, in the
picture of the Space Shuttle Main Engine shown on the left hand side of Figure 30.
In pump-fed rockets, the power Ppump required by each pump can be estimated by:
mp
Ppump (71)
where p is the pressure increase provided by the pump to the fluid, is the density of the
liquid propellant and is the pump efficiency, which is usually in the order of 50% for the
pumps typically used in space rocket engines.
47
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
In the Space Shuttle Main Engine, the total mass flow rate of liquid oxygen is 1264 kg/s.
The liquid oxygen pumps need to increase the propellant pressure from 1 MPa to the value
required in the combustion chamber, which is 20.6 MPa. The density of liquid oxygen is
1140 kg/m3.
Assuming a pump efficiency of 50%, calculate the total power required by these pumps.
-----------------------------------------------------------------------------------------------------------
The required pressure increase is:
p 20.6 MPa 1 MPa 19.6 MPa 19.6 106 Pa
And the required pump power, from equation (71), is:
kg
1264 19.6 106 Pa
mp s
Ppump 43.46 106 W 43.46 MW
0.5 1140 3
kg
m
48
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
propellants, it is not always practical to use them and, to avoid the design complexity
introduced by the need of maintaining the very low temperature required by these propellants,
it is preferred in some cases to accept a lower specific impulse and opt for a different, non-
cryogenic or only partially cryogenic, combination of propellants.
49
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
temperature of the combustion products from its initial value T0 to a final value TC is
calculated by taking into account the total mass flow rate of propellants (fuel and oxidizer):
Pheat m fuel moxidizer c p TC T0 m fuel 1 O /F c p TC T0 (73)
where cp is, in this case, the constant pressure specific heat of the combustion products.
Also in this case, we assume that no thermal losses occur and the entire combustion power is
directly used to heat the combustion products. We can therefore equate equations (71) and
(72) and obtain an expression for the chamber temperature:
HV
Pheat Pdec m fuel HV m fuel 1 O /F c p TC T0 TC T0 (74)
c p 1 O /F
Note, however, that this equation has been obtained under the assumption that the entire
amount of fuel flowing in the combustion chamber is burned by the oxidizer. This does not
happen in the so-called “fuel-rich” mixtures, in which the amount of fuel is in excess to what
can be combined in the chemical reaction with the oxidizer, and part of it is not effectively
burned. In this case, the value given by equation (74) is an over-estimation of the actual
combustion chamber temperature achieved in the rocket.
A rocket uses kerosene and liquid oxygen as propellants, injected in the combustion
chamber at an initial temperature of 300 K with mixture ratio O/F = 2.3. The heating value
of kerosene is 46200 kJ/kg. The combustion products have constant pressure specific heat
cp = 4.48 kJ/(kg∙K), specific heat ratio = 1.12 and molecular mass MW = 21.3 g/mol.
Calculate the chamber temperature and the characteristic velocity of this rocket.
-----------------------------------------------------------------------------------------------------------
For the chamber temperature, from equation (74):
kJ
46200
HV kg
TC T0 300 K 3425 K
c p 1 O /F 4.48
kJ
1 2.3
K kg
And for the characteristic velocity, from equation (63):
1 11.12
1 1 1 1.12 11.12
1.12 0.633
2 2
1 RA 1 8314 J K kmol m
c* TC 3425 K 1827
MW 0.633 21.3 g mol s
50
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
37
C12 H26 O2 12CO2 13H 2 O
2
The combustion products of this reaction are water vapour (H2O) and carbon dioxide (CO2);
in particular, nCO2 = 12 moles of CO2 and nH2O = 13 moles of H2O are formed for each mole of
kerosene. These are the weights used to evaluate the molar averaged values of constant
pressure specific heat and molecular mass. In particular, if cpCO2 , cpH2O , MWCO2 and MWH2O
are, respectively, the constant pressure specific heats and molecular masses of CO2 and H2O:
nCO2 c pCO2 nH2O c pH 2O 12 c pCO2 13 c pH 2O
cp
nCO2 nH2O 12 13
A similar method can be applied to calculate the mixture ratio corresponding to a given
chemical reaction, but in this case (recalling the definition of mixture ratio) the molecular
masses of the combustion reactants need to be considered, using their respective number of
moles as weights. As an example, for the above reaction of oxygen (oxidizer) and kerosene
(fuel):
nO2 M WO2 37 2 MWO
O/F 2
A liquid bi-propellant engine operates with liquid oxygen and liquid methane, according to
the following chemical reaction in the combustion chamber:
CH4 +2O2 CO2 +2H2O
The molecular mass of oxygen (O2) is 32 g/mol and the molecular mass of methane (CH4)
is 16 g/mol. The constant pressure specific heat cp is 45 J/(mol∙K) for carbon dioxide
(CO2) and 39 J/(mol∙K) for water vapour (H2O). Assuming that the initial propellant
temperature is 300 K and the heating value of methane is HV = 7.5∙105 J/mol, how much is
the chamber temperature?
-----------------------------------------------------------------------------------------------------------
The mixture ratio O/F is calculated based on the combustion reactants:
2 M WO2 64
O/F 4
1 M WCH 4 16
While the constant pressure specific heat is calculated using the combustion products:
1 c pCO2 2 c pH 2O 1 45 2 39 J
cp 41
1 2 1 2 mol K
And the chamber temperature, from equation (74), is:
HV 7.5 105
TC T0 300 3958.5 K
c p 1 O /F 41 1 4
51
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
52
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
53
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
To understand the importance of equation (78) note that, assuming that combustion takes
place at constant chamber temperature, all the variables on the right hand side of the equation
are constant with time, except the propellant grain burning surface Ab. The equation thus
shows a clear dependence of chamber pressure on the burning surface: if the burning surface
changes with time, the chamber pressure changes with time accordingly, as it will be
explained more in detail in the next section.
In the same way, however, it is also possible to design grain shape geometries in which the
burning surface is constant with time or decreases with time. Some examples of possible grain
burning behaviours are shown in Figure 37.
Figure 37: Different types of solid propellant grain geometries, in terms of burning
surface as a function of time.
Looking back at equation (78), and recalling that in practical cases the exponent n is always
lower than one to avoid unstable combustion, the burning surface behaviour as a function of
time shown in Figure 37 also reflects on the chamber pressure behaviour as a function of
time: increasing, constant or decreasing burning surface also means, in turn, increasing,
constant or decreasing chamber pressure.
54
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
A solid propellant rocket has cylindrical grain shape similar to the one shown in Figure 36.
The initial grain diameter db0 is 0.5 m, its final diameter (at the end of burning) dbf is 1 m,
its length L (assumed constant during the burn) is 5 m. The solid propellant has a density
of 1600 kg/m3, a characteristic velocity of 1600 m/s, and its regression rate coefficients
are n = 0.5 and a = 10-5 m/s*Pa0.5. The throat area A* is 0.2 m2.
Calculate the initial and final chamber pressure, and the initial and final regression rate.
-----------------------------------------------------------------------------------------------------------
The initial and final grain burning surface is:
Ab0 db0 L 0.5 m 5 m 7.85 m2
Abf dbf L 1 m 5 m 15.7 m2
Thus, from equation (78), the initial and final chamber pressure is:
1
5 m kg m 10.5
10 s Pa 0.5 1600 m3 1600 s
1
a s c* 1n
pC0 *
Ab0 2
7.85 m 2 10.1 105 Pa 10.1 bar
A 0.2 m
1
5 m kg m 10.5
10 s Pa 0.5 1600 m3 1600 s
1
a s c* 1n
pCf *
Abf 2
15.7 m 2 40.4 105 Pa 40.4 bar
A 0.2 m
And finally, from equation (76), the initial and final regression rate is:
0.5
10.1 105 Pa 0.01 1
m 0.5 m cm
r0 a pC0 n 105
s Pa s s
0.5
40.4 105 Pa 0.02
m 0.5 m cm
rf a pCf n 105 2
s Pa s s
55
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
56
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
7 – Summary problems
7.1 Rocketdyne F1 Engine
The Rocketdyne F1 liquid bi-propellant
engine was used for the first stage of the
Saturn V rocket (the launcher for the Apollo
missions to the Moon). The propellants are
RP-1 kerosene (C12H26) and oxygen (O2), and
the following chemical reaction is assumed in
the combustion chamber:
C12 H26 18.5O2 12CO2 13H2O
The following data are given on this bi-propellant engine:
• Combustion chamber pressure: pC = 7∙106 Pa
• Initial propellant temperature: T0 = 300 K
• Circular nozzle, exit diameter: de = 3.75 m
• Specific heat ratio of the combustion products: = 1.13
• Heating value of RP-1: HV = 46200 kJ/kg
• Constant pressure specific heat of the combustion products: cp =2.34 kJ/(K*kg)
• Molecular masses: 44 g/mol (CO2) , 18 g/mol (H2O), 32 g/mol (O2) , 170 g/mol (C12H26)
• Nozzle pressure ratio: pC / pe = 119
• Atmospheric pressure at sea level: pa = 101325 Pa
Assuming that the equivalent jet velocity is constant with time, calculate the following
quantities for this engine:
1. Mixture ratio
2. Nozzle exit pressure
3. Nozzle throat area
4. Thrust at sea level
5. Specific impulse at sea level
6. Number of F1 engines required by the Saturn V first stage, knowing that the required v is
3.75 km/s, the burn time is 150 s, the rocket take-off mass is 2.8106 kg, and assuming for
simplicity that the specific impulse is constant and equal to its sea level value.
Solution
Mixture ratio:
18.5 M W (O 2) 18.5 32
O/F 3.482
1 M W (C12 H 26) 1 170
57
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
in
M M P
P in
MP 2.03 106 kg
M P N m tB N 5 engines
m tB 2707 150
58
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Solution
Combustion chamber temperature of the Space Shuttle Main Engine:
kJ
141900
HV kg
TC T0 75 K 3600 K
c p 1 O /F 5.75
kJ
1 6
K kg
59
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Total mass flow rate of one single Space Shuttle Main Engine:
pC A* 20.6 106 Pa 0.053 m 2 kg
m 0.631 471.2
RA J s
TC 8314
MW K kmol 3600 K
g
14
mol
Fuel mass flow rate in one single Space Shuttle Main Engine:
moxidizer m m fuel m m 471.2 kg
O/F 1 m fuel 67.3
m fuel m fuel m fuel O / F 1 6 1 s
Oxidizer mass flow rate in one single Space Shuttle Main Engine:
kg
moxidizer m fuel O / F 67.3 6 403.9
s
Nozzle expansion ratio in the Space Shuttle Main Engine:
Ae 0.631
77.45
A* 2 1
2
1.115 1
2 pe pe
2 1.115 1 1.115
1 1.115
1 1
1 pC pC 1.115 1 753 753
Nozzle exit pressure and exit area in the Space Shuttle Main Engine:
pC p 20.6 106 Pa
753 pe C 27357 Pa
pe 753 753
60
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
Sea level specific impulse of one single Space Shuttle Main Engine:
veq ( sea level ) FT ( sea level ) 1.83 106 N
I sp ( sea level ) 395.9 s
g0 m g0 471.2
kg m
9.81 2
s s
Take-off thrust of each Solid Rocket Booster:
kg m
FT (TakeOff ) mTakeOff veq (TakeOff ) 5265 2374 12.5 106 N
s s
Take-off regression rate in the Solid Rocket Boosters:
mTakeOff mTakeOff 5265 m cm
m s r Ab rTakeOff 18.2 103 1.82
s Ab (TakeOff ) s Db i L 1920 1.2 40 s s
61
Faculty of Aerospace Engineering, Delft University of Technology
Lecture Notes AE2230-II
Rocket Propulsion
References
[1] Sutton, G.P., 1978, Rocket Propulsion Elements, 7th edition, John Wiley & Sons Inc.
[2] Ward, T.A., 2010, Aerospace Propulsion Systems, John Wiley & Sons Inc.
[3] Hill, P.G., Peterson C.R., 1992, Mechanics and Thermodynamics of Propulsion, Addison-
Wesley Publishing Co.
[4] Fortescue, P., Swinerd, G., Stark, J., 2011, Spacecraft Systems Engineering, 3rd Ed., John
Wiley & Sons
[5] Larson, W. J., Wertz, J.R., 1992, Space Mission Analysis and Design, 3rd edition,
Microcosm, Inc., Torrance, CA
62