Lecturenote - 555581022chapter 2 Compressible Flow
Lecturenote - 555581022chapter 2 Compressible Flow
2
Learning Objectives
After completing this chapter, you should be able to:
▪ Explain speed of sound and Mach number and their practical significance.
3
Introduction to compressible flow
What do you mean by a compressible flow?
❑ Fluids have the capacity to change volume and density, i.e. compressibility.
4
.. cont’d
❑ Therefore, liquid flows will almost always be considered as incompressible,
but gas flows could easily be either incompressible or compressible.
▪ Gas has large compressibility but when its velocity is low compared with the
sonic velocity the change in density is small and it is then treated as an
incompressible fluid.
▪ The compressibility factor has to be considered in the following engineering
problems.
➢ The flight of projectiles and airplanes moving at high altitudes
➢ The flow of gases through nozzles and orifices
➢ In meteorological studies where the gas system involves appreciable
variation in density due to greater heights involved
➢ In quick-return-valves, water hammers, acoustics, rapid & repeated
vibrations and such other fields.
5
.. cont’d
6
.. cont’d
Flows according to Mach number:
▪ Incompressible Ma < 0.3
▪ Subsonic 0.3 < Ma < 0.8
▪ Transonic 0.8 < Ma < 1.2
Compressible flows
▪ Supersonic 1.2 < Ma < 3.0
▪ Hypersonic 3.0 < Ma
❑ Modern aircraft are mainly powered by gas turbine engines that
involve transonic flows.
❑ When a space shuttle reenters the earth’s atmosphere, the flow is
hypersonic. Future aircraft may be expected to operate from subsonic to
hypersonic flow conditions.
7
.. cont’d
▪ Ma < 1 Subsonic ▪ Ma > 1 Supersonic
In general when:
▪ Ma = 1 Sonic ▪ Ma > 5 Hypersonic
➢ Under small Mach number condition, change in fluid density are everywhere
small in the flow field.
➢ If the density change is significant, it follows from the equation of state the
temperature and pressure change are also substantial.
➢ The consequence of compressibility are not limited simply to density changes.
➢ Density change means that, we can have significant compression or expansion
work on a gas, so the thermodynamics states of the fluid will change, meaning
that in general all properties like – temperature, internal energy, entropy,
enthalpy and so on can change.
➢ For this reason, we begin with a review of the thermodynamics need to study
compressible flow.
8
A review of the thermodynamic properties .. cont’d
Compressible flows require:
▪ Continuity equation
▪ Momentum equation
▪ Energy equation
▪ Equation of state
➢ Basic thermodynamics relations:
9
.. cont’d
𝑅𝑈
R=
➢ Characteristic values for air: 𝑀
Where 𝑅𝑈 = 8314J/(kgK) (universal
gas constant)
M = 29.98 ( Molecular weight)
10
.. cont’d
Basic Thermodynamics Process:
▪ Isochoric process (constant volume)
𝑃1 𝑇1 𝑃
= 𝑂𝑟 = 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡 , 𝑃1 𝑉1 = 𝑚𝑅𝑇1 𝑎𝑛𝑑 𝑃2 𝑉2 = 𝑚𝑅𝑇2
𝑃2 𝑇2 𝑇
▪ Isothermal process
𝑃1 𝑉1 = 𝑃2 𝑉2 𝑜𝑟 𝑃𝑉 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
11
.. cont’d
Basic Thermodynamics Process:
▪ Adiabatic process:
𝛾 𝛾
𝐴𝑛 𝑎𝑑𝑖𝑎𝑏𝑎𝑡𝑖𝑐 𝑝𝑟𝑜𝑐𝑒𝑠𝑠 𝑓𝑜𝑙𝑙𝑜𝑤𝑠 𝑡ℎ𝑒 𝑙𝑎𝑤: 𝑃1 𝑉1 = 𝑃2 𝑉2 = 𝑃𝑉 𝛾 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 𝐶
12
.. cont’d
▪ Isentropic Process (reversible adiabatic)
13
Cont…
14
Speed of Sound
15
Speed of Sound … cont’d
▪ Imagine a sound source that creates pressure disturbance
at a certain point (i.e p is changed to P + dP due to the
sound. From that point sound waves will propagate
equally in all directions similar to the wave that occur
when throwing a stone in to water.
▪ Where, C = wave speed
➢ Note that: the wave strength (sound wave) is regarded as a wake
wave and that the change in properties are therefore infinitesimal.
➢ To derive an expression for the speed of sound, we set up the
continuity equation as follows:
17
Speed of Sound … cont’d
𝑃
▪ From equation of state, = 𝑅𝑇
𝜌
▪ From this equation the following are recognized for speed the of sound
➢ The speed of sound is only dependent on gas properties and temperature
➢ The higher the temperature the higher the speed of sound; this is due to
increased activity of the gas molecules
18
… cont’d
Example:
▪ Estimate the speed of sound of carbon monoxide at 200-kPa pressure and 300°C in m/s.
Solution
From Table A.4, for CO, the molecular weight is 28.01 and k =1.40.
Thus RCO = Ru/Mm= 8314/28.01 = 297 m2/(s2K),
For the given temperature T = 300°C + 273 = 573 K.
Thus, from the equation of speed of sound
aCO = (kRT)1/2 = [1.40(297)(573)]1/2 = 488 m/s
Table A.4 Properties of Common Gases at 1 atm and 20°C
(68°F)
▪ The
19
… cont’d
20
Adiabatic and Isentropic Steady Flow
𝑉2
▪ 𝑇𝑂 = 𝑇 + is the stagnation or total head temperature of the flow stream.
2𝐶𝑃
▪ Maximum velocity obtained when enthalpy & temperature drop to(absolute)
zero:
1 1
𝑉𝑚𝑎𝑥 = (2ℎ𝑜 ) ൗ2 = (2𝐶𝑃 𝑇𝑂 ) ൗ2
21
… cont’d
Mach number relations for adiabatic flow (isentropic or not!)
22
… cont’d
• Mach number relations for isentropic flow (must be adiabatic, too!)
23
… cont’d
Critical values (at Ma=1.0)
• Values at critical point, where Mach number is equal to one (sonic
conditions) are of special importance for compressible flow calculations.
For that reason, we mark these values by an asterisk:
24
… cont’d
Some Useful Numbers for Air
▪ For k = 1.4, the following numerical versions of the isentropic and adiabatic
flow formulas are obtained:
▪ Or, if we are given the properties, it is equally easy to solve for the Mach
number (again with k = 1.4)
25
… cont’d
Example:1. Explain the concept of stagnation properties:
▪ A supersonic plane flies at 2000km/hr at an altitude of 9km above sea
level in standard atmospheric. If the pressure and density of air at this
attitude are states to be 30kN/m2 absolute and 0.45kg/m3, make
calculations for the pressure, temperature and density at the stagnation
point on the nose of the plane. Take R = 287 J/kg.K and γ= 1.4
26
… cont’d
27
Isentropic Flow with Area Changes
28
Isentropic Flow with Area Changes …cont’d
How ? proof
30
…cont’d
❑ Flow through converging–diverging nozzle (nozzle with a throat)
(De Laval nozzle)
• It is required to accelerate a flow from subsonic to supersonic flow
conditions.
• Used in supersonic aircrafts
31
…cont’d
Perfect Gas Area Change
❑ Area ratio (relative to the critical area) can be expressed for the
convergent/ divergent nozzle as a function of Mach number only.
After lengthy algebra:
For, k = 1.4
32
Choking …cont’d
▪ Choked flow occurs when the Mach number is 1.0 at the minimum cross-
sectional area.
▪ From the 1D mass conservation, the ratio of mass flow rates per unit area is
inversely proportional to the area ratio:
▪ The maximum flow rate occurs when the narrowest area (the throat) reaches
the sonic conditions (Ma=1). We call such flow choked, since any increase of
incoming velocity will decrease the mass flow rate. The maximum mass
flow rate is thus:
33
…cont’d
▪ The mass flow function starts at 0, for p=p0, and levels off as p/p0
approaches the critical value of 0.5283 (sonic conditions).
34
Isentropic Flow Tables for k = 1.4 …cont’d
35
Isentropic Flow Table (continued)
36
…cont’d
37
…cont’d
38
Normal Shocks
▪ When air undergoes large and rapid compression (e.g. following an
explosion, the release of engine gases into an exhaust pipe, or where
an aircraft or a bullet flies at supersonic velocity) a thin wave of large
pressure change is produced as shown in Fig.1.
40
…cont’d
Normal-Shock Wave Relations:
▪ To compute all property changes rather than just the wave speed, we
use all our basic one-dimensional steady-flow relations.
41
…cont’d
▪ The first successful analyses of these normal-shock relations are credited
to W. J. M. Rankine (1870) and A. Hugoniot (1887), hence the modern
term Rankine-Hugoniot relations.
➢ Assumed upstream conditions (p1, V1, ρ1, h1, T1) and (p2, V2, ρ2, h2,
T2) are unknowns
➢ Due to velocity-squared term, two solutions are found, and from
2nd law of thermodynamic S2 >S1 and again with eliminating V1
and V2 from mass & energy equations the Rankine-Hugoniot
relations obtained as:
42
…cont’d
➢ Also, the actual change in entropy across the shock can be computed
from the perfect gas relation:
43
Mach-Number Relations with Normal shock …cont’d
▪ Property change across the normal shock for a perfect gas are obtained as:
44
…cont’d
Some remarks
▪ The upstream flow is supersonic, and the downstream flow is subsonic.
▪ There is entropy increase across the wave, thus the flow is adiabatic but
not isentropic (because it is irreversible). So, and .
▪ Stagnation temperature remains the same across the shock wave, but the
stagnation pressure and density decrease in the same ratio.
▪ Any stagnation process (e.g. Pitot tube) must induce a shock wave in
front of it, if it is inserted into a supersonic flow
45
Operation of Converging and Diverging Nozzles
❑ Converging Nozzle
• Consider the converging nozzle sketched in Fig. Where, 𝑇0 and 𝑝0
(stagnation properties), pe - exit pressure and pb - back pressure.
(b) Pressure distribution caused by ▪ A nozzle when operating below the design
various back pressures value is said to be under-expanding and
above the design pressure is called over -
expanding
47
…cont’d
48
Fanno line and Rayleigh line
▪ Friction must be included for flow through long ducts, especially if the
cross-sectional area is small.
49
…cont’d
➢ Dividing the first equation by the second one and noting that T01 = T02,
we have:
…………..……….. (a)
…..………... (b)
50
…cont’d
➢ Substituting these into the conservation of mass relation ρ1V1 = ρ2V2 and
noting that, Ma = V/a and , a = √kRT, we have
…..…………………………….…... (c)
➢ Combining Eqs. (a) and (c) gives the pressure ratio across the shock as:
…..………………….…... (d)
51
…cont’d
52
…cont’d
Property relations for Rayleigh Flows:
➢ Recalling the momentum equation as: P1 - P2 = ρ2V22 - ρ1V12
And, ρV2 = PkMa2
➢ After certain manipulation, we obtain
……….……….. (e)
…….…….. (f)
❑ Eq. (f) represents the intersections of the Fanno and Rayleigh lines and
relates the Mach number upstream of the shock to that downstream of the
shock.
53
…cont’d
Fanno & Rayleigh lines plots on h-s diagram
55
Summary …cont’d
56
…cont’d
57
…cont’d
Sound Wave Propagation from a Moving Source
▪ As the particle moves, it continually crashes against fluid particles
and sends out spherical sound waves emanating from every point
along its path.
▪ The semi-angle of the cone is called the Mach angle,𝜇 and it is given
by:
58
…cont’d
Sound wave propagation from a moving source
59
…cont’d
▪ The pressure waves do not reach the region to the left of the normal
front and this region is called the silence zone (out side the cone)
▪ In supersonic flow the pressure disturbance are felt only inside the
cone and this region is called zone of action.
60
…cont’d
.
61
…cont’d
Example 2
▪ An aircraft cruising at 1000-m elevation, above you moves past in a flyby. It
is moving with a Mach number equal to 1.5, speed of sound 343.3 m/s and
the ambient temperature is 20 °C. Find, How many seconds after the plane
passes overhead do you expect to wait before you hear the aircraft?
▪ Solution:
62
…cont’d
Reading Assignment
▪ Oblique shock waves: characteristic features, governing
equations, calculation of properties
▪ Textbook : White and Manson
63
Oblique Shock Waves …cont’d
.
▪ At leading edge, flow is deflected
through an angle called the
turning/deflection angle
64
Oblique Shock Waves …cont’d
Note: the typical supersonic CL is much smaller than the subsonic value CL ≈ 2π𝛼,
but the lift can be very large because of the large value of 1/2ρV2 at supersonic
speeds.
67
…..cont’d
.
68
.
Thank you
End
69