Fluid 11
Fluid 11
Fluid 11
Jyh-Cherng Shieh
Department of Bio-Industrial Mechatronics Engineering National Taiwan University 1
MAIN TOPICS
Ideal Gas Relationships Mach Number and Speed of Sound Categories of Compressible Flow Isentropic Flow of an Ideal Gas Nonisentropic Flow of an Ideal Gas Two-Dimensional Compressible Flow
Introduction
Fluid compressibility is a very important consideration in numerous engineering applications of fluid mechanics. For example, The measurement of high-speed flow velocities requires compressible flow theory. The flows in gas turbine engine components are generally compressible. Many aircraft fly fast enough to involve a compressible flow field. In this study of compressibility effects, we mainly consider the steady, one-dimensional, constant (including zero) viscosity, compressible flow of an ideal gas.
4
p = RT
(1)
R is gas constant
The gas constant, R, represents a constant for each distinct ideal gas or mixture of ideal gases
R= M gas
(2) is the universal gas constant. Mgas is the molecular weight of the ideal gas or gas mixture.
( ( du u cv = = T v dT
(3)
Where the subscript v on the partial derivative refers to differentiation at constant specific volume, v=1/.
( ( ( du = c v dT u 2 u1 = c v (T2 T1 )
T2 ( ( u 2 u1 = cvdT T1
(5)
(4)
6
Enthalpy 1/2
The fluid property enthalpy
( ( p h=u+
(6)
( ( h = u + RT
p = RT ( ( h = h()
(9)
Enthalpy 2/2
( ( (10) h = u + RT ( ( ( dh du ( dh = du + RdT = + R (11) dT dT
( ( h dh cP = = T p dT
( ( du u cv = = T v dT
cp cv = R
(12)
cP k= cv
(12) + (13)
(13)
k cP = R k 1
(14)
1 cV = R k 1
(15)
8
Example 11.1 Internal Energy, Enthalpy, and Density for an Ideal Gas
Air flows steadily between two sections in a long straight portion of 4-in.-diameter pipe as is indicated in Figure E11.1. The uniformly distributed temperature and pressure at each section are T1=540R, p1=100 psia, and T2=453R, p2=18.4 psia. Calculate the (a) change in internal energy between sections (1) and (2), (b) change in enthalpy between sections (1) and (2), and (c) change in density between sections (1) and (2).
( ( u 2 u1 = c v (T2 T1 )
cV = R 1 (15) k 1
k = 1.4 R = 1716(ft lb ) /(slug R ) = 53.3(ft lb ) /(lbm R )
cP k= cv
(13)
10
p = RT
(1)
11
Entropy 1/3
For any pure substance including ideal gases, the first T ds equation is
1 ( (16) Tds = du + pd T is absolute temperature ( 1 1 ( dh = du + pd + dp (17)
(16)+(17)
( 1 Tds = dh dp
(18)
second T ds equation
dT R 1 (1)+(3)+(16) ds = c v + d T y
(1)+(7)+(18) ds = c P
(19) (20)
12
dp dT R P T
Entropy 2/3
If cp and cv are assumed to be constant for a given gas
(19)
1 T2 s 2 s1 = c v ln + R ln T1 2
(21)
(20)
p2 T2 s 2 s1 = c P ln R ln T1 p1
(22)
13
Entropy 3/3
For the adiabatic and frictionless flow
ds = 0, c v ln s1 = s 2
Constant entropy flow T2 T p + R ln 1 = c P ln 2 R ln 2 = 0 (23) or isentropic flow T1 2 T1 p1
k / k 1 k
T2 T 1
p2 2 = = p 1 1
(25)
(24)
p = cons tan t k
14
15
T2 p2 s 2 s1 = c P ln R ln T1 p1
(22)
17
Mach Number
The Mach number, Ma, was a dimensionless measure of compressibility in a fluid flow. The Mach number is defined as the ratio of the value of the local flow velocity, V, to the local speed of sound, c.
V Ma = c
What we perceive as sound generally consists of weak pressure pulse that move through air. When our ear drums respond to a succession of moving pulse, we hear sounds.
18
Leaving CV
20
Ac = ( + ) A(c V )
(26)
c = c V ()(V)
small
(27)
V = c
(28)
21
ccA + (c V )( + )(c V ) A = A ( p + p) A
ccA + (c V )Ac = pA
(26)+(29)
(29)
p V = c
(30)
(28)+(30)
c2 =
c=
(31)
This expression for the speed of sound results from application of the conservation of mass and conservation of linear momentum principles to the flow through the control volume.
22
(32)
(V )2 << cV
p (c V ) 2 c 2 + =0 2 2
2
p (33) V = c
p c =
p c=
(31)
This expression for the speed of sound results from application of the conservation of mass and conservation of energy principles to the flow through the control volume.
23
p p 0
(25)
p c= s
(34)
(35)
c = RTk
(36)
c = RTk =
Ev =
Ev
(38) (37)
24
p dp = d / s
c = RTk
k = 1.4 R = 286.9J /( kg K )
26
27
28
29
30
31
33
The pressure information can still travel unrestricted throughout the flow field, but not symmetrically or instantaneously.
34
Mach wave
The pressure waves are all tangent to a plane that is perpendicular to the flow and that passes through the point source. This plane is called a Mach wave. The communication of pressure information is restricted to the region of flow downstream of the Mach wave.
35
c 1 sin = = V Ma
(39)
36
The schlieren visualization of flow (supersonic to subsonic) through a row of compressor airfoils.
37
38
39
1 sin = Ma
The speed of the aircraft V = Ma c
Ma =
1 sin[tan 1 (1000 / Vt )]
41
42
(40)
In chapter 3, Newtons second law was applied to the inviscid (frictionless) and steady flow of a fluid particle. For either compressible or incompressible flow along the streamwise direction,
1 dp + d( V 2 ) + dz = 0 2
(41)
Can be dropped because of its small size in comparison to the other terms
43
d dA dV + + =0 A V
dV d dA = + V A
(45)
(44)
(42)+(44)
dp V 2 dA 1 = 2 V dp / d A
44
V c
(46)
p c= s
(47)
(34)
dp dA 2 (1 Ma ) = 2 V A
dV dA 1 = V A (1 Ma 2 )
(48)
(42)+(47)
(44)+(48)
d dA Ma 2 = A (1 Ma 2 )
(49)
45
When the flow is subsonic (Ma<1), velocity and section area changes are in opposite directions. When the flow is supersonic (Ma>1), velocity and area changes in the same direction.
(a) A diverging duct. (b) A converging duct.
46
47
For subsonic flow (Ma<1), density and area changes are in the same direction. For supersonic flow (Ma>1), density and area changes are in the opposite direction.
48
For Ma =1, this equation requires that dA/dV=0. This means that the area associated with Ma=1 is either a minimum or a maximum amount. If the flow enter converging-diverging duct were subsonic, the fluid velocity would increase in the converging portion of the duct, and achievement of a sonic condition (Ma=1) at the minimum area location appears possible.
49
51
Conclusion 10/10
1. 2. 3. 4. 5. 6. Sonic condition (Ma=1) can be attained in a converging-diverging duct at the minimum area location. This minimum area location is often called the throat of the converging-diverging duct. To achieve supersonic flow from a subsonic state in a duct, a converging-diverging area variations is necessary. A converging-diverging duct is referred as a converging-diverging nozzle. A converging-diverging duct can also decrease a supersonic flow to subsonic conditions. A converging-diverging duct can be a nozzle or a diffuser depending on whether the flow in the converging portion of the duct is subsonic or supersonic.
52
53
p p0 = cons tan t = k k 0
(25)
It is convenient to use the stagnation state of the fluid as a reference state for compressible flow calculations.
The streamwise equation of motion for steady and frictionless flow can be expressed for an ideal gas
(41)
V2 dp =0 + d 2
(51)
Note: The stagnation state is can also be achieved by isentropically decelerating a flow to zero velocity.
V2 p1 / k dp 0 + d 1/ k 2 =0 o p
(52)
54
k p0 p V 2 2 = 0 k 1 0
(53)
kR V2 (To T ) = 0 k 1 2
V2 c P (To T ) =0 2
(14)
( ( h 2 h1 = c P (T2 T1 )
(54)
k cP = R k 1
kR V2 (To T ) = 0 k 1 2
T 1 = 2 To 1 + k 1 M a 2
(56)
Pressure variation ?
With this equation we can calculate the temperature of an ideal gas anywhere in the converging-diverging duct if the flow is steady, one-dimensional, and isentropic, provided we know the value of the local Mach number and the stagnation temperature.
56
p o T = po To
T p = T po o
(57)
k / k 1
(58)+(56)
p = po 1 +
1
k 1 2
Ma
k / k 1
(59)
1 / k 1
(56)+(57)+(59)
= o 1 +
1
k 1 2
Ma
(60)
57
58
(56)
Ma
p = po 1 +
1 2 k 1 Ma 2
k / k 1
(59)
Ma
Lower fluid temperature and pressures are associated with higher Mach number in an isentropic converging-diverging duct.
59
p = p0 1 +
1
k 1 2
Ma
k / k 1
= 0 1 +
1
k 1 2
Ma
1 / k 1
Used to describe steady + T-s diagram flow through the converging-diverging duct.
60
Ma
(59)
k / k 1
Ma=1
p* 2 = po k + 1
k / k 1
(61)
62
Ma=1
T* 2 = To k + 1
(65)
(63)
* p* To 2 = * = Po T k + 1 o
63
64
65
Stagnation density
0 = 1.23kg / m 3
k = 1 .4
th = o 1 +
1
k 1 2
Ma
2 th
(11.5.2)
(59)
p th = po 1 +
1
k 1 2
Ma
2 th
k / k 1
(11.5.3)
(62)
66
Vth = Ma th c th = Ma th RTth k
(11.5.4)
Tth 1 = To 1 + k2 1 Ma 2 th
(11.5.5)
67
Ma
2 th
(11.5.2)
th = 1.04kg / m3
Tth 1 = To 1 + k2 1 Ma 2 th
(11.5.5)
Tth = 269 K
68
& m = th A th Vth
(11.5.1)
& m = 0.0201kg / s
69
Ma
2 th
1 / k 1
(11.5.2)
th = 0.780kg / m3
Tth 1 = To 1 + k2 1 Ma 2 th
Vth = 310( J / kg)1 / 2 = 310m / s
(11.5.5)
Tth = 240K
& m = 0.0242kg / s
70
71
72
& m = th A th Vth
(11.6.1)
Since the receiver pressure, pre=80 kPa (abs) > 53.3 kPa (abs) = p*, we have pth= pre.
p th 80 kPa = = 0 . 792 po 101 kPa
Ma th = 0.59
th = 0.85 (11.6.3)
73
th = (0.85)(1.23kg / m 3 ) = 1.04 kg / m 3
74
th = (0.64)(1.23kg / m 3 ) = 0.79 kg / m 3
75
76
V = Ma RTk
k = 1.4
77
Ma
k / k 1
(59)
Ma = 0.499
T = To 1 +
1
k 1 Ma 2 2
(56)
AV = A V
* *
* V * = A * V A
(68)
(67)
(36)+(46)
V* = RT*k
A
V = Ma RTk
T * / T0 T/T 0
(69)
(67)+(68)+(69)
1 * = A* Ma o
(70)
79
A 1 1 + [( k 1) / 2]Ma = * A Ma 1 + ( k 1) / 2
2
( k +1) / 2 ( k 1)
(71)
The variation of area ratio with Mach number for isentropic flow of an ideal gas (k = 1.4, linear coordinate scales).
80
Figure D.1
Value of p/p0, T/T0, /0, and A/A* are graphed in Figure D.1 as a function of Mach number for air (k=1.4).
T = To 1 + 1
k 1 Ma 2 2
(56)
k / k 1
= o 1 +
(59)
( k +1) / 2 ( k 1)
1
k 1 2
Ma
1 / k 1
(60)
p = po 1 +
k 1 Ma 2 2 1
A 1 1 + [( k 1) / 2]Ma = * A Ma 1 + ( k 1) / 2
2
(71)
81
Example 11.8 Isentropic Choked Flow in a Converging-Diverging Duct with Subsonic Entry
Air enters subsonically from standard atmosphere and flows isentropically through a choked converging-diverging duct having a circular cross-sectional area, A, that varies with axial distance from the throat, x, according to the formula A=0.1+x2 where A is in square meters and x is in meters. The duct extends from x = -0.5 m to x = +0.5 m. For this flow situation, sketch the side view of the duct and graph the variation of Mach number, static temperature to stagnation temperature ratio, T/T0, and static pressure to stagnation pressure ratio, p/p0, through the duct from x = -0.5 m to x = +0.5 m. Also show the possible fluid states at x = -0.5 m, 0 m, and +0.5 m using temperature entropy coordinates.
82
0.1 + x r=
(11.8.3)
83
Values of A/A* can be used in Eq. (71) to calculate corresponding values of Mach number, Ma.
A 1 1 + [( k 1) / 2]Ma = * A Ma 1 + ( k 1) / 2
2
( k +1) / 2 ( k 1)
(71)
84
(56)
p = po 1 +
2 k 1 Ma 2 1
k / k 1
T/T0 p/p0
(59)
85
A 0.1 + x 2 = A* 0.1
Figure D.1
86
87
89
91
92
93
A A* = = 0.07 m 2 ( A / A*)
With known values of duct area at different axial locations, we can calculate corresponding area ratios, A/A*, knowing A*=0.07m2. With values of area ratio A/A*, we can use Figure D.1 and obtain related values of Ma, T/T0, and p/p0.
95
96
(a) Subsonic to subsonic isentropic flow (not choked). (b) Subsonic to subsonic isentropic flow (Choked).
97
(c) Subsonic to supersonic isentropic flow (choked), (d) Supersonic to supersonic isentropic flow (choked).
98
(e) Supersonic to subsonic isentropic flow (choked). (f) Supersonic to subsonic isentropic flow (not choked).
99
Solution 1/3
For a given stagnation state (i.e., T0 and p0 fixed), ideal gas (k=constant), and converging-diverging duct geometry, an infinite number of isentropic subsonic to subsonic (not choked) and supersonic to supersonic (not choked) flow solutions exist. In contrast, the isentropic subsonic to supersonic (choked), subsonic to subsonic (choked), supersonic to subsonic (choked), and supersonic to supersonic (choked) from solutions are each unique. Given stagnation state (i.e., T0 and p0 fixed), ideal gas (k=constant), and converging-diverging duct geometry.
Solution?
100
Solution 2/3
When the pressure at x=+0.5 (exit) is greater than or equal to pI, isentropic flows are possible. When the pressure at x=+0.5 (exit) is less than or equal to pII, isentropic flows are possible. When the exit pressure is less than pI and greater than pII,
101
Solution 3/3
(a) The variation of duct radius with axial distance. (b) The variation of Mach number with axial distance. (c) The variation of temperature with axial distance. (d) The variation of pressure with axial distance.
102
103
104
Nonisentropic Flow
105
Fanno flow: adiabatic flow with friction. Rayleigh flow: Flows with heat transfer (diabatic flows) without friction.
106
Fanno Flow
Adiabatic flow with friction
107
Small
( ( h h 0 = c P (T T0 )
(73) (74)
For particular Fanno flow, the stagnation temperature T0 is fixed Used to calculate values of fluid (V) 2 T 2 T+ = T0 = cons tan t (75) temperature corresponding to 2c P ( p 2 / R 2 ) value of pressure in the Fanno flow. 109
Determined until later
p2 T2 s 2 s1 = c P ln R ln T1 p1
(22)
If the temperature, T1, and entropy s1, at the entrance of the Fanno flow duct as reference values
(22)
T p s s1 = c P ln R ln T1 p1
(76)
110
(76)
= T0 = cons tan t
(75)
This curve involves a given gas (cp and R) with fixed values of stagnation temperature, densityvelocity product, and inlet temperature, pressure, and entropy.
T-s diagram for Fanno flow
111
112
= T0 = cons tan t
(11.11.1) (11.11.2)
T p s s1 = c P ln R ln T1 p1
k = 1.4
(11.11.3) (11.11.4)
113
T = To 1 +
1
k 1 Ma 2 2
(56)
114
115
( dh = c P dT (7)
p = RT
(1)
dp d dT = + T
(78) (79)
(77)
d dT Tds = c PdT RT + T
dV dT Tds = c P dT RT + T V ds c P 1 dV 1 = R + dT T V dT T
(80)
116
dV c = P dT V
(81) (82)
cp 1 ds c P = R 2 + V dT T T
ds =0 dT cp = Rk k 1
At state a
Va = RTa k
(83)
Va = RTa k
ds =0 dT
vs.
c = RTk
(a) Subsonic Fanno flow. (b) Supersonic Fanno flow. (c) Normal shock occurrence in Fanno flow (an abruot change from supersonic to subsonic flow in the Fanno duct).
119
120
121
& = m(V2 V1 )
Rx is the frictional force exerted by the inner pipe wall on the fluid
RX = V( V2 V1 ) A
122
f=
8 w V 2
(86)
123
V 2 dx dp f = VdV 2 D dp f V 2 dx d ( V 2 ) + + =0 p p 2 D p 2
(1)+(36)+(46)+(88)
(87) (88)
dp fk Ma 2 d ( V 2 ) 2 dx + Ma +k = 0 (89) 2 p 2 D 2 V
V = Ma RTk
2 2
V = Ma c = Ma RTk
d( V 2 ) V2
d ( Ma 2 ) Ma 2
dT (90) T
124
dT d ( V 2 ) + =0 T 2c P T
2 dT k 1 2 d( V ) + Ma =0 2 T 2 V
(91)
(92)
(92)+(90)
d( V 2 ) V
2
1 + [( k 1) / 2]Ma
d ( Ma 2 ) / Ma 2
(93)
125
dp 1 d( V 2 ) d ( Ma 2 ) = 2 p 2 V Ma 2
(94)
(94)+(89)
2 dx 1 d( Ma 2 ) fk 2 d( V ) (1 + kMa ) + Ma 2 =0 2 2 2 V Ma 2 D
(95)
(1 Ma 2 )d ( Ma 2 ) {1 + [( k 1) / 2]Ma 2 }kMa 4
=f
dx D
(96)
(96) can be integrated from one section to another in Fanno flow duct. We elect to use the critical (*) state as a reference and to integrate (96) from an upstream state to the critical state.
126
M* a
(1 Ma 2 ) d ( Ma 2 ) {1 + [( k 1) / 2 ] Ma }kMa
2 4
Ma
l*
dx D
(97)
l is length measured from an arbitrary but fixed upstream reference location to a section on the Fanno flow. Friction factor is constant at an average value over the integration length l*- l. K is also constant.
1 (1 Ma 2 ) k + 1 [( k + 1) / 2]Ma 2 f (l l *) + = ln 2 2 D k Ma 2k 1 + [( k 1) / 2]Ma
(98)
127
NEXT PAGE
128
130
(k 1) dT d ( Ma 2 ) = T 2{1 + [( k 1) / 2]Ma 2 }
(100)
Integrating Eqs. (100) from any state upstream in a Fanno flow to the critical (*) state
T ( k + 1) / 2 = T * 1 + [( k 1) / 2]Ma 2
(101)
= Ma T T*
(68)+(69)
V V*
Ma RTk RT k
*
(102)
1/ 2
V [( k + 1) / 2]Ma 2 = V * 1 + [( k 1) / 2]Ma 2
(103)
131
V = * V *
2
(104)
V * 1 + [( k 1) / 2]Ma = = * V [( k + 1) / 2]Ma 2
1/ 2
(105)
p T = p* * T*
(106)
1/ 2
(106)+(105)+(101)
p 1 ( k + 1) / 2 = p * Ma 1 + [( k 1) / 2]Ma 2
(107)
132
1 2 k 1 2 = Ma 1 + * 2 p o Ma k + 1 po
[( k +1) / 2( k 1)]
(109)
133
Figure D.2
Value of f(l*- l)/D, T/T*, V/V*, p/p*, and p0/p0* for Fanno flow of air (k=1.4) are graphed as a function of Mach number in Figure D.2.
f (l * l 2 ) f (l * l1 ) f = (l1 l 2 ) D D D
T ( k + 1) / 2 = T * 1 + [( k 1) / 2]Ma 2
(99)
1/ 2
(101)
1/ 2
V [( k + 1) / 2]Ma 2 = V * 1 + [( k 1) / 2]Ma 2
(103)
p 1 ( k + 1) / 2 = p * Ma 1 + [( k 1) / 2]Ma 2
(107)
[( k +1) / 2( k 1)]
1 2 k 1 2 Ma = 1 + * 2 p o Ma k + 1 po
(109)
134
135
136
137
p 0,1 p0 *
= 1.16
138
Since T0=288K
T1 = 0.93T0 = 267.84 K
T* = T1 / 1.1 = 243.49 K
The stagnation temperature, T0, remain constant through this adiabatic flow at a value of
T0,1 = T0,2 = T0 = 288K
140
The duct exit stagnation pressure p * 1 p 0,2 = p 0 * = 0 p 0,1 = 101kPa (abs) = 84 kPa (abs) p 0,1 1.16
141
142
Entering Figure D.2 with f ( l * l 1 ) / D = 0.2 We read p1 T1 V1 Ma1 = 0.7 = 1 .5 = 1.05 = 0.73 p* T* V*
143
T1 = 0.92 T0
p1 = 0.72 p 0,1
1 = 0.79 0,1
Since T0=288K
T1 = 0.92T0 = 264.96 K
T* = T1 / 1.05 = 253.34 K
146
147
Rayleigh flow
Flows with heat transfer without friction
149
(110)
Frictionless
150
For a given Rayleigh flow, the constant, V, and R are all fixed. Eqs. (111) can be used to determine values of fluid temperature corresponding to the local pressure in a Rayleigh flow.
151
(111) (76)
p T s s1 = c P ln R ln p1 T1
These two equations can be solved simultaneously to obtain the Rayleigh lines.
Rayleigh line.
152
dp = VdV
(113)
dp = VdV
(112)
( Tds = dh + VdV
( ( h dh cP = = (7) T dT p
(114)
(114)+(112)+(1)+(77)+(79)
ds c P V dV c P V 1 = + = + dT T T dT T T [( T / V ) ( V / R ) ]
(115)
153
Va = RTa k
(116)
dT 1 = 0 Ma b = ds k
(118)
154
& + Wshaft
net in
or in differential form for Rayleigh flow through the seminfinitesimal control volume (
dh + VdV = q
(119)
The heat transfer per unit mass of fluid in the semi-infinitesimal control volume
155
dV q V dT V ( k 1) = + V c P T T dV kRT
2
(120)
(36)+(46)+(1)+(77)+(79)+(112)+(120)
dV q 1 = V c P T (1 Ma 2 )
(121)
From (121), we see clearly that when the Rayleigh flow is subsonic (Ma<1), fluid heating (q>0) increases fluid velocity while fluid cooling decreases fluid velocity. When the Rayleigh flow is supersonic (Ma>1), fluid heating decrease fluid velocity and fluid cooling increases fluid velocity.
156
(a) Subsonic Rayleigh flow. (b) Supersonic Rayleigh flow. (c) Normal shock in a Rayleigh flow.
157
a 2 p V 2 + =1+ Va pa pa pa
(122)
159
T p a = Ta p a
p 1+ k = p a 1 + kMa 2
(124)
(123)
a V = Va
(125) (126)
a T = Ma Ta
T p = Ma Ta p a
2
(124)+(126)
(127)
160
(128)
(125)+(126)+(128)
a V (1 + k ) Ma = = Ma Va 1 + kMa 2
(129)
The energy equation tells us that because of the heat transfer involved in Rayleigh flows, the stagnation temperature varies.
To To T = To, a T Ta
Ta T o ,a
(130)
161
To T o, a
k 1 2( k + 1) Ma 1 + Ma 2 2 = (1 + kMa 2 ) 2
2
(131)
p o p po = p p p 0 ,a a
(132)+(123)+(59)
p a p 0 ,a
(132)
p0 1+ k = p 0,a 1 + kMa 2
2 k 1 2 Ma 1 + 2 k + 1
k / k 1
(133)
162
Figure D.3
Value of p/pa, T/Ta, /a or V/Va, T0/T0,a, and p0/p0,a are graphed in Figure D.3 as a function of Mach number for Rayleigh flow of air (k=1.4).
p 1+ k = p a 1 + kMa 2
(123)
T (1 + k ) Ma = (128) Ta 1 + kMa 2
2
k 1 Ma 2 2( k + 1) Ma 1 + 2 = (131) 2 2 (1 + kMa )
p0 1+ k = p 0,a 1 + kMa 2
2 k 1 1+ Ma 2 2 k + 1
k / k 1
(133)
163
164
k 1
2
V = 16.7lbm /(ft s)
165
167
Example 11.16 Effect of Mach Number and Heating/Cooling for Rayleigh Flow
The information in Table 11.2 shows us that subsonic Rayleigh flow accelerates when heated and decelerates when cooled. Supersonic Rayleigh flow behaves just opposite to subsonic Rayleigh flow; it decelerates when heated and accelerates when cooled. Using Figure D.3 for air (k=1.4), state whether velocity, Mach number, static temperature, static pressure, and stagnation pressure increase or decrease as subsonic and supersonic Rayleigh flow is (a) heated, (b) cooled.
168
169
p = RT
(V )2 T 2 T+ = To = cons tan t 2 2 2c P ( p / R )
172
174
We should be able to use equations developed earlier for Fanno and Rayleigh flows to quantify normal shock wave.
175
p y pa = p x pa p x
py
(137)
px 1+ k = p a 1 + kMa 2 x
py px
1 + kMa 2 x 1 + kMa 2 y
(140)
py
(107)
1 + kMa 2 x = p x 1 + kMa 2 y py
(140)
Ty T * = Tx T * Tx
Ty
(141)
T ( k + 1) / 2 = From T * 1 + [( k 1) / 2]Ma 2
(101)
Ty T*
( k + 1) / 2 1 + [( k 1) / 2]Ma 2 y
(142)
Tx ( k + 1) / 2 = T * 1 + [( k 1) / 2]Ma 2 x
Ty Tx
1 + [( k 1) / 2]Ma 2 x 1 + [( k 1) / 2]Ma 2 y
(144)
179
py px
Ty y Tx x
py
x Vx = y Vy
(145)
Ty Vx = p x Tx Vy
(147)
py
(146)
(146)+(36)+(46)
Ty = p x Tx
1/ 2
Ma x Ma y
1/ 2
(147)+(144)
p y 1 + [( k 1) / 2]Ma 2 x = p x 1 + [( k 1) / 2]Ma 2 y
Ma x Ma y
(148)
180
Ma 2 = y
Ma 2 + [2 /( k 1)] x [2 k /( k 1)]Ma 2 x 1
(149)
Equation (149) can be used to calculate values of Mach number downstream of a normal shock from a known Mach number of the shock. Max>1..May<1
181
k 1 2k 2 = Ma x px k + 1 k +1
Ty Tx = {( k + 1)
(152)
2
py
(150)
(149)+(144)
(151)
(40)
Vx = x Vy
(152)+(1)
p y Tx = x p x Ty
(153)
182
Vx ( k + 1) Ma 2 x = = x Vy ( k 1) Ma 2 + 2 x
(154)
(59)+(149)+(150)
p 0,y p 0,x
+ ( k 2 1 Ma
2 k /( k 1 ) ) (1 + k 2 1 Ma 2 ) k /( 1 k ) x x k k ( k2+ 1 Ma 2 k + 1 ) 1 /( k 1 ) x 1
(156)
183
Figure D.4
Value of downstream Mach numbers, May, pressure ration, py/px, temperature ratio, Ty/Tx, density ratio, y/x or velocity ratio Vy/Vx, and stagnation pressure ratio, p0,y/p0,x are graphed in Figure D.4 as a function of upstream Mach number, Max, for the steady flow across a normal shock wave of an ideal gas having a specific heat ratio k=1.4.
Ma 2 y
Ty Tx =
Ma 2 + [2 /( k 1)] x [2 k /( k 1)]Ma 2 x
{( k + 1)
2
(149)
py px
2k k 1 Ma 2 x k +1 k +1
y
(150)
Vx ( k + 1) Ma 2 x = = (151) x Vy ( k 1) Ma 2 + 2 x
(154)
p 0,y p 0,x
2 k /( k 1 ) (1 + k 2 1 Ma 2 ) k /( 1 k ) x) x k k ( k2+ 1 Ma 2 k + 1 ) 1 /( k 1 ) x 1
(156)
184
185
186
k=1.4
p 0,y p 0,x
2k k 1 2 = Ma x px k + 1 k +1
py
187
188
Where p0,y is the stagnation pressure measured by the probe, and px is the static pressure measured by the wall trap. The stagnation pressure upstream of the shock, p0,x, is not measured.
(11.18-1)+(156)+(59)
p 0,y px
+ ( k 2 1 Ma k ( k2+ 1 Ma
2 k /( k 1 ) x) 2 k k + 1 ) 1 /( k 1 ) x 1
189
Max=1.9
To determine the flow velocity we need to know the static temperature upstream of the shock.
V x = Ma x c x = Ma
RT x k
The stagnation temperature downstream of the shock was measured and found to be
T o , y = 1000 R T o , x = T o , y = 1000 R
190
or Figure D.1
V x = Ma
RT x k = 1 . 9
191
192
M a x = 2 .8
px = 0 . 04 p 0 ,x
py px
= 9 .0
py p 0 ,x p y p x = p p x 0,x = ( 9 . 0 )( 0 . 04 ) = 0 . 36 = p III p 0 ,x
193
px = 0 . 10 p 0,x
194
For x = +0.3m, the ratio of duct exit area to local area (A2/Ay) is
Ay A2 A2 0 .1 + ( 0 .5 ) 2 A2 = = 1 . 842 = = 2 . 28 2 Ay 0 .1 + ( 0 .3 ) A * A * Ay
195
196
198
gy
(157) (159)
Voc Fr = c oc
Voc is the velocity of the channel flow coc is the velocity of a small amplitude wave on the surface of an open-channel flow
199
For open-channel flows, Continuity equation ybVoc = cons tan t (161) y and b are the depth and width of the open-channel flow. Comparing Eqs. (160) and (161) we note that if flow velocities are considered similar and flow area, A, and channel width, b, are considered similar, the compressible flow density, , is analogous to open-channel depth, y.
200
c = (cons tan t ) k k 1
(162)
201
should be equal to 2. Typically k=1.4 or 1.67, not 2. This limitation to exactness is, however, usually not serious enough to compromise the benefits of the analogy between compressible and open-channel flows.
202
208
Supersonic flow over a wedge: (a) Smaller wedge angle results in attached oblique shock.
209
Supersonic flow over a wedge: (b) Large wedge angle results in detached curve shock.
210