Lecture 2 - Compressible Flow (2022 Notes)
Lecture 2 - Compressible Flow (2022 Notes)
COMPRESSIBLE FLOW
B. Huyssen
15 March 2022
OUTLINE
Converging-diverging nozzle
Chapter 10
1
15/03/2022
dA dP
A u 2
1 M 2
This eq governs the shape of a nozzle or
dA
A
du
u
1 M 2 a diffuser in subsonic or supersonic
isentropic flow
QUASI ONE-DIMENSIONAL
ONE-DIMENSIONAL ISENTROPIC
ISENTROPIC FLOW
FLOW
diverging
2
15/03/2022
3
15/03/2022
QUASI ONE-DIMENSIONAL
ONE-DIMENSIONAL ISENTROPIC
ISENTROPIC FLOW
FLOW
Example illustrates
M = 1 at the location of the smallest flow area, called the throat
Velocity continues to increase past the throat, and is due to decrease in
density
Area decreases, and then increases. Known as a converging - diverging
nozzle. Used to accelerate gases to supersonic speeds.
To take a gas to supersonic speed we must accelerate the gas subsonically through a
convergent duct reach sonic condition in the throat and than expand the gas to
supersonic speed through a divergent duct.
To slow down isentropically a supersonic flow to subsonic flow we must decelerate the
gas through a convergent duct reach sonic condition in the throat and than decelerate to
subsonic speed.
4
15/03/2022
Here, we will study the effects of area change on flow properties, as the
exit velocity, mass flow rate, and pressure distribution along the
nozzle.
RT
P
m AV A M RT PAM M
V
RT RT a
a RT
Substituting T and P from the expressions
of total conditions T0 1 2 p0 1 2
1
1 M1 1 M1
T1 2 p1 2
5
15/03/2022
4
P0 AM P0 A* 2
RT0 RT0
m 1 m 1
2 1
1 2 2 1
1 M 1
2 1
2
=Pb=P*
The maximum mass flow rate through a
nozzle (point 3) with a given throat area is fixed by the m AV
stagnation pressure and temperature of the inlet flow!
The flow at the throat is frozen, once the flow As Pe decreases V increases
becomes sonic the disturbances cannot work their way and ρ decreases slightly
upstream of the throat. No communication between
convergent section and divergent section.
This principal is important for chemical processes,
medical devices, flow meters, and anywhere the mass
flux of a gas must be known and controlled.
6
15/03/2022
EXAMPLE
GIVEN the converging nozzle below with throat area of 50cm2 and
back pressure of 0.7MPa
m AV PAM
RT
cp
1.289
cv
P1=1MPa P0 AM
T1=600C c p cv R m
RT0
V1=150m/s 1
2 1
1 2
1 M
2
R=287J/kgK=287m2/Ks2
A * a*
m Au * A*u *
A* u * * 0
1
* 2 1 0
0 1 2 2 2 2 2
A * 0 a* a* 1
* *2
1
A 0 u u M
0 1 2 1
1 M1
1 2
M *2
1M 2
2 1M 2
Area Mach number relation 1
This eq is tabulated in appendix A, A
2
1 2
* 2 1
1 M 2
1
1
isentropic flow. It yields to two
A M 1 2
solutions of M
7
15/03/2022
A M 1 2
A M 1 2
3
A 1 1 0.2M 2
A* M 2 1.728
For γ=1.4
P0 A* * For isentropic flow through a duct, the
RT0 P0 A 0.6847
m maximum mass flow rate possible is
1
1
2 1 RT0 1/ 2 proportional to the throat area and stagnation
pressure and inversely proportional to the
2 square root of stagnation temperature
8
15/03/2022
A M 1 2
When we get the variation of M we can get any other fluid
ratio parameter as p/p0, T/T0. To be able to get Me>1, pe
must have a specific value for which
3
the eq Ae/A* is satisfied.
For all converging-diverging nozzles,
the maximum value of Pb/P0 is For a given shape of nozzle there
always 1 and the minimum is always 0. is only one supersonic solution
This leaves three other pressure ratios. and depends from A/A*
9
15/03/2022
EXAMPLE
1
m AV PAM A
2
1 2
1 M 2
1
A* M 2 1 1
RT
2
p0 1 2 1
1 Me P0 AM
pe 2 RT0
m 1
2 1
1 2
1 M
2
Here, we will study the conditions under which shock waves develop
and how they affect the flow.
10
15/03/2022
DISTURBANCES
IN SUBSONIC & SUPERSONIC FLOW
MACH CONE
Beeper is stationary
V 0
M 0
a a
At time t the beeper, in 0, emits a sound wave that propagate a distance at. The beeper
has zero velocity. The wave fronts spread out spherically and can be drawn as
concentric circles with position corresponding to times t0, t1, t2, t3.
11
15/03/2022
MACH CONE
Beeper is moving slow, Subsonic compressible Flow
V
M 1
a
At time t the beeper emits a sound wave that propagate a distance at. In the meantime A
moved to B with v<a while emitting several sound waves. The beeper is slower than the
sound wave. It is inside the family of circular sound waves. The circles can not overlap
MACH CONE
Beeper is moving at sound speed
V
M 1
a
B A
At time t the beeper emits a sound wave that propagate to distance at. In the meantime
A moved to B while emitting several sound waves. The rate of spread of the disturbance
upstream is exactly balanced by the beeper velocity.
12
15/03/2022
MACH CONE
Beeper is moving faster than sound speed, SUPERSONIC Flow
Zone of V
silence M 1
a
Zone of action
Mach angle
at a 1
sin
Vt V M
At time t the beeper emits a sound wave that propagate a distance at. In the meantime A moved to
B while emitting several sound waves. The beeper is faster than the sound wave, disturbances are
swept downstream at a greater speed than they can propagate upstream. The beeper is outside the
family of circular sound waves tangent to BC. BC line of disturbances is called Mach wave. The
disturbances(messages) emitted by the beeper do not have the opportunity to warn the fluid ahead
of its arrival. The fluid that knows about the presence of the beeper is restricted to the Mach cone.
The Mach wave is the envelope of disturbances in the supersonic flow and it is oblique to the
direction of motion. The effect of the disturbances does not go beyond the Mach lines into the
surrounding fluid.
EXAMPLE
GIVEN a slender missile flying at Mach 1.5 at low altitude. Assume the wave
generated by the nose of the missile is a mach wave. This wave intersects the
ground 170m behind the nose.
1
sin
M
13
15/03/2022
Wave 2
Gas at rest
Wave 1
14
15/03/2022
SHOCK WAVES
NORMAL SHOCKS
15
15/03/2022
16
15/03/2022
There is no external heat transfer across the control volume so the shock is adiabatic.
cpT0,2 =cp T0,1
By applying the thermodynamics eq we arrive to the result of T0,2 = T0,1 s2 > s1
17
15/03/2022
Typical shock wave is thin region, thickness 1/1,000 mm, highly viscous flow. Normal
shock wave is one dimensional supersonic flow.
Conservation of energy
Conservation of momentum
Increase in entropy
18
15/03/2022
a p /
a1
2
1 1
u12 a *2
1 2 2 1 a* stays constant because the flow is adiabatic
a2
2
1 1
u 22 a *2 u1 u2
1 2 2 1 a *2 u1u 2 1
a* a*
M *1M * 2
u1 u 2 1
1 M *1M * 2 M *2 1
a* a* M *1
1
1 M 2
1M 2 M 22 2
1
1
*2
M
2 1M 2 M 12
2
Mach number after the shock is function only of Mach number ahead the
shock. The upstream Mach number dictates the shock wave parameters.
If M1=1 than M2=1
19
15/03/2022
1
1 M 2 M behind a normal shock
Normal shock is adiabatic
wave is a function only of
but non isentropic
1
M 22 2 M ahead the wave!
M 12
1 M1>1 than M2<1 more M1>1 Equations are functions of
2 more stronger becomes M1, only
2 u1
1M 12 the normal shock wave and
Mach number behind a normal
less is M2
1 u2 2 1M 12 shock wave is always subsonic
(M2 < 1)
2
p2
p1
1
1
M 12 1 Density, static pressure, and
temperature increase across
T2 h2 2 2 1M 12
1 M 12 1 a normal shock wave
T1 h1 1 1M 1 Velocity and total pressure
2
20
15/03/2022
0.5 10
0.4 8
0.3 6
0.2 4
0.1 2
0 0
1 2 3 4 5 6 7 8 9 10
Upstream Mach Number, M1
21
15/03/2022
EXAMPLE
GIVEN the pressure upstream of a normal shock wave is 1 atm. The pressure and
the temperature downstream of the wave are 10.33 atm and 772.2 K˚ respectively.
FIND Mach number and the temperature upstream of the wave and the
total temperature and total pressure downstream of the wave.
22