Isentropic
Isentropic
Isentropic
C HAPTER 4: ISENTROPIC
Version: 0.4.8.4
October 8, 2008
Fundamentals of Compressible
Flow Mechanics
T HE
Project
Name
Progress
Compressible Flow
beta
Die Casting
Remarks
P OTTO P ROJECT
Version
Availability
for
Public
Download
0.4.8.3
alpha
0.0.3
Dynamics
NSY
0.0.0
Fluid Mechanics
alpha
0.1.8
Heat Transfer
NSY
0.0.0
Mechanics
NSY
0.0.0
NSY
0.0.0
0.0.1
Open
Flow
Channel
Based
on
Eckert
Statics
early
alpha
Strength of Material
NSY
0.0.0
Thermodynamics
early
alpha
0.0.01
Two/Multi
flow
NSY
0.0.0
phases
first
chapter
TelAvivs
notes
CHAPTER 5
Isentropic Flow
In this chapter a discussion on a steady
state flow through a smooth and continuous area flow rate is presented.
A discussion about the flow through
a convergingdiverging nozzle is also
part of this chapter. The isentropic flow
models are important because of two
main reasons: One, it provides the information about the trends and important
parameters. Two, the correction factors
can be introduced later to account for
deviations from the ideal state.
PB = P0
P
P0
Subsonic
M <1
Supersonic
M >1
distance, x
49
50
U2
= h0
2
(5.1)
Perfect gas is an ideal gas with a constant heat capacity, Cp . For perfect gas
equation (5.1) is simplified into
Cp T +
U2
= Cp T0
2
(5.2)
(5.3)
kR
k1
(5.4)
T0
U2
=
2Cp T
T
(5.5)
T0
kRU 2
=
2
2Cp c
T
(5.6)
By utilizing the definition of k by equation (2.24) and inserting it into equation (5.6)
yields
1+
k 1 U2
T0
=
2
2
c
T
(5.7)
It very useful to convert equation (5.6) into a dimensionless form and denote Mach number as the ratio of velocity to speed of sound as
M
U
c
(5.8)
Inserting the definition of Mach number (5.8) into equation (5.7) reads
T0
k1 2
=1+
M
T
2
(5.9)
51
Q
B
A
The usefulness of Mach number and equation (5.9) can be demonT0
T0
P0
P0
strated by this following simple example.
velocity
0
0
In this example a gas flows through a
tube (see Figure 5.2) of any shape can
be expressed as a function of only the
Fig. -5.2: Perfect gas flows through a tube
stagnation temperature as opposed to
the function of the temperatures and velocities.
The definition of the stagnation state provides the advantage of compact
writing. For example, writing the energy equation for the tube shown in Figure (5.2)
can be reduced to
Q = Cp (T0B T0A )m
(5.10)
T0
T
k
k1
k
k 1 2 k1
1+
M
2
(5.11)
T0
T
1
k1
k 1 2 k1
M
= 1+
2
(5.12)
A new useful definition is introduced for the case when M = 1 and denoted by
superscript . The special case of ratio of the star values to stagnation values are
dependent only on the heat ratio as the following:
T
c 2
2
= 2 =
T0
c0
k+1
P
=
P0
=
0
k
k1
(5.14)
1
k1
(5.15)
2
k+1
2
k+1
(5.13)
52
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
1
2
3
Mon Jun 5 17:39:34 2006
4
5
Mach number
Fig. -5.3: The stagnation properties as a function of the Mach number, k = 1.4
n(n 1)x2
n(n 1)(n 2)x3
+
+
2!
3!
(5.16)
P
4
8
48
(5.17)
0
(k 1)M 2
kM 4
2(2 k)M 6
=1+
+
+
4
8
48
(5.18)
}|
{
z
M2
(2 k)M 4
P0 P
=1+
+
+
1
2
4
24
2 U
(5.19)
From the above equation, it can be observed that the correction factor approaches
zero when M 0 and then equation (5.19) approaches the standard equation
for incompressible flow.
The definition of the star Mach is ratio of the velocity and star speed of
sound at M = 1.
U
M = =
c
k+1
M
2
k1 2
1
M +
4
(5.20)
M2
+
1+
4
(5.21)
kM 2
+
1
4
(5.22)
P0 P
kM 2
=
P
2
0
M2
=
kP0 2 M 2
k1 2
m
1+
=
M +
A
RT0
4
(5.23)
2
k+1
k+1
2(k1)
k+1
(3 k)(k + 1) 3
1
+
M+
M +
M
4
32
(5.24)
54
T+dT
+d
P
U
P+dP
U+dU
(5.25)
A
U
(5.26)
dP
(5.27)
For isentropic process ds 0 and combining equations (5.25) with (5.27) yields
dP
+ U dU = 0
(5.28)
(5.29)
}|
{
dA
d
dP
U2
+
=0
(5.30)
Rearranging equation (5.30) so that the density, , can be replaced by the static
pressure, dP/ yields
1
c2
z}|{
dA
d
dP
dA
dP
d dP
2
2
=U
+
+
=U
A
dP
dP
A
(5.31)
Recalling that dP/d = c2 and substitute the speed of sound into equation (5.31)
to obtain
"
2 #
dA
U
dP
(5.32)
= U2
1
c
A
Or in a dimensionless form
dA
dP
1 M2 = U2
(5.33)
Equation (5.33) is a differential equation for the pressure as a function of the cross
section area. It is convenient to rearrange equation (5.33) to obtain a variables
separation form of
dP =
dA
U 2
A 1 M2
(5.34)
56
one dP and dA have the same sign. For the subsonic branch M < 1 the term
1/(1 M 2 ) is positive hence
dA > 0 = dP > 0
dA < 0 = dP < 0
From these observations the trends are similar to those in incompressible fluid.
An increase in area results in an increase of the static pressure (converting the
dynamic pressure to a static pressure). Conversely, if the area decreases (as a
function of x) the pressure decreases. Note that the pressure decrease is larger in
compressible flow compared to incompressible flow.
For the supersonic branch M > 1, the phenomenon is different. For M > 1
the term 1/1 M 2 is negative and change the character of the equation.
dA > 0 dP < 0
dA < 0 dP > 0
This behavior is opposite to incompressible flow behavior.
For the special case of M = 1 (sonic flow) the value of the term 1 M 2 = 0
thus mathematically dP or dA = 0. Since physically dP can increase only
in a finite amount it must that dA = 0.It must also be noted that when M = 1
occurs only when dA = 0. However, the opposite, not necessarily means that
when dA = 0 that M = 1. In that case, it is possible that dM = 0 thus the diverging
side is in the subsonic branch and the flow isnt choked.
The relationship between the velocity and the pressure can be observed
from equation (5.28) by solving it for dU .
dU =
dP
PU
(5.35)
From equation (5.35) it is obvious that dU has an opposite sign to dP (since the
term P U is positive). Hence the pressure increases when the velocity decreases
and vice versa.
From the speed of sound, one can observe that the density, , increases
with pressure and vice versa (see equation 5.36).
d =
1
dP
c2
(5.36)
It can be noted that in the derivations of the above equations (5.35 - 5.36), the
equation of state was not used. Thus, the equations are applicable for any gas
(perfect or imperfect gas).
The second law (isentropic relationship) dictates that ds = 0 and from
thermodynamics
dP
dT
R
ds = 0 = Cp
T
P
k1 2
dT0 = 0 = dT 1 +
+ T (k 1)M dM
(5.38)
M
2
and simplifying equation (5.38) yields
dT
(k 1)M dM
=
2
T
1 + k1
2 M
(5.39)
(5.40)
U 2
dA
A 1M 2
M dM
2
1 + k1
2 M
(5.41)
z}|{ P
P
U 2 z }| {
U 2 z }| {
2
kM P = k 2 RT = k
RT = U 2
c
kRT
Using the identity in equation (5.42) changes equation (5.41) into
1 The
dA
M2 1
dM
=
2
A
M 1 + k1
2 M
(5.42)
(5.43)
58
croon
k1
2
,
A cti
M 1 + 2 M and A are positive rese
gardless of the values of M or A. TheredA = 0
fore, the only factor that affects relationdx
ship between the cross area and the
dM = 0
2
Mach number is M 1. For M < 1
dx
the Mach number is varied opposite to
M, Much nubmer
the cross section area. In the case of
M > 1 the Mach number increases with
x
the cross section area and vice versa.
The special case is when M = 1 which
Fig. -5.5: The relationship between the cross
requires that dA = 0. This condition
section and the Mach number on
imposes that internal2 flow has to pass
the subsonic branch
a convertingdiverging device to obtain
supersonic velocity.
This minimum
area is referred to as throat.
Again, the opposite conclusion that when dA = 0 implies that M = 1 is
not correct because possibility of dM = 0. In subsonic flow branch, from the
mathematical point of view: on one hand, a decrease of the cross section increases
the velocity and the Mach number, on the other hand, an increase of the cross
section decreases the velocity and Mach number (see Figure (5.5)).
T
T0
A
A
P
P0
AP
A P0
F
F
With these values the static temperature and the density can be calculated.
T = 0.86420338 (273 + 21) = 254.076K
0
z }| {
P0
5 106 [P a]
h
i
=
= 0.69428839
J
0 RT0
287.0 kgK
294[K]
kg
= 41.1416
m3
The velocity at that point is
c
z }| {
60
S OLUTION
With the known Mach number at point A all the ratios of the static properties to
total (stagnation) properties can be calculated. Therefore, the stagnation pressure
at point A is known and stagnation temperature can be calculated.
At M = 2 (supersonic flow) the ratios are
M
2.0000
T
T0
A
A
P
P0
AP
A P0
F
F
= 0.12780453
= 0.17040604
P0
P0
PB
1.5
The corresponding Mach number for this pressure ratio is 1.8137788 and TB =
0.60315132 PPB0 = 0.17040879. The stagnation temperature can be bypassed to
calculate the temperature at point B
M =2
M =1.81..
z}|{
z}|{
TB
T0
TB = TA
TA
T0
= 250[K]
1
0.60315132 271.42[K]
0.55555556
Example 5.3:
Gas flows through a convergingdiverging duct. At point A the cross section area
is 50 [cm2 ] and the Mach number was measured to be 0.4. At point B in the duct
the cross section area is 40 [cm2 ]. Find the Mach number at point B. Assume that
the flow is isentropic and the gas specific heat ratio is 1.4.
S OLUTION
To obtain the Mach number at point B by finding the ratio of the area to the critical
area. This relationship can be obtained by
AB
AB
AA
40
=
=
A
AA
A
50
B
With the value of A
A from the Table (5.2) or from Potto-GDC two solutions can
be obtained. The two possible solutions: the first supersonic M = 1.6265306 and
second subsonic M = 0.53884934. Both solution are possible and acceptable. The
supersonic branch solution is possible only if there where a transition at throat
where M=1.
T
T0
A
A
P
P0
AP
A P0
0.22617 0.28772
0.82071 1.0440
Example 5.4:
Engineer needs to redesign a syringe for medical applications. The complain in
the syringe is that the syringe is hard to push. The engineer analyzes the flow
and conclude that the flow is choke. Upon this fact, what engineer should do with
syringe increase the pushing diameter or decrease the diameter? Explain.
S OLUTION
This problem is a typical to compressible flow in the sense the solution is opposite
the regular intuition. The diameter should be decreased. The pressure in the choke
flow in the syringe is past the critical pressure ratio. Hence, the force is a function
of the cross area of the syringe. So, to decrease the force one should decrease
the area.
End Solution
P
UA
RT
(5.44)
This parameter is studied here, to examine the maximum flow rate and to see what
is the effect of the compressibility on the flow rate. The area ratio as a function of
the Mach number needed to be established, specifically and explicitly the relationship for the chocked flow. The area ratio is defined as the ratio of the cross section
at any point to the throat area (the narrow area). It is convenient to rearrange the
equation (5.44) to be expressed in terms of the stagnation properties as
f (M,k)
P P0 U
m
=
A
P0 kRT
k
R
P
T0 1
= 0 M
T
T0
T0
}| {
z r
k P
T0
R P0 T
(5.45)
=
A
kM P0
kRT0
k1 2
1+
M
2
k+1
2(k1)
(5.46)
62
as by A . It can be noticed that at the throat when the flow is chocked or in other
words M = 1 and that the stagnation conditions (i.e. temperature, pressure) do
not change. Hence equation (5.46) obtained the form
!
k+1
2(k1)
kP0
m
k1
(5.47)
=
1+
A
2
RT0
Since the mass flow rate is constant in the duct, dividing equations (5.47) by
equation (5.46) yields
A
1
=
A
M
1+
k1
2
2 M
k+1
2
k+1
! 2(k1)
(5.48)
Equation (5.48) relates the Mach number at any point to the cross section area
ratio.
The maximum flow rate can be expressed either by taking the derivative of
equation (5.47) in with respect to M and equating to zero. Carrying this calculation
results at M = 1.
r
k+1
2(k1)
k k+1
P0
m
=
(5.49)
A max T0
R
2
For specific heat ratio, k = 1.4
m
P
0.68473
0
A max T0
R
The maximum flow rate for air (R = 287j/kgK) becomes,
m
T0
= 0.040418
A P 0
(5.50)
(5.51)
Equation (5.51) is known as Fliegners Formula on the name of one of the first engineers who observed experimentally the choking phenomenon. It can be noticed
that Fliengners equation can lead to definition of the Fliengners Number.
c0
z }| { z F}|n {
p
mc
0
m
T0
m
kRT0
1
=
=
A P0
RA P0 k
kRA P0
(5.52)
mc
0
RA P0
(5.53)
k1 2
M
F n = kM 1 +
2
(5.54)
k+1
2
k+1
2(k1)
(5.55)
results in
energy equation (5.9) and substituting for Mach number M = m/Ac
k1
T0
=1+
T
2
Ac
(5.56)
1/kR
z }| {
2
z}|{
T k1 m
2
T0 = T +
c2
2
A
(5.57)
P
k1
+
T0 R 2kRT0
(5.58)
64
Equation (5.58) is quadratic equation for density, when all other variables are
known. It is convenient to change it into
2
P
k1 m
2
=0
(5.59)
T0 R 2kRT0 A
The only physical solution is when the density is positive and thus the only solution
is
2
2
u
u P
1
k
1
m
=
(5.60)
+2
+u
u
2 RT0 t RT0
kRT0 A
{z
}
|
(M 0)0
For almost incompressible flow the density is reduced and the familiar form of
perfect gas model is seen since stagnation temperature is approaching the static
P
temperature for very small Mach number ( = RT
). In other words, the terms
0
for the group over the underbrace approaches zero when the flow rate (Mach
number) is very small.
It is convenient to denote a new dimensionless density as
=
p
RT0
RT0
1
=
P
T
(5.61)
s
2
(k 1)RT0 m
1
1+ 1+2
=
2
kP 2
A
(5.62)
The dimensionless density now is related to a dimensionless group that is a function of Fn number and Mach number only! Thus, this dimensionless group is function of Mach number only.
A P 0
AP
F n2
=f (M )
}|
{z
}|
{
z
2 2 2
2
2
A
RT0 m
m
1 c0
P0
=
P2
A
k P0 2 A
A
P
(5.63)
Thus,
RT0
P2
F n2
k
A P0
AP
2
(k 1)F n2 A P0
1
1+ 1+2
=
2
k2
AP
(5.64)
(5.65)
F n2
=
k
A P0
AP
F n2
=
k
T
T0
P0
P
(5.66)
P0
P
(5.67)
A P 0
AP
(5.68)
The right hand side is tabulated in the regular isentropic Table such (5.2). This
example shows how a dimensional analysis is used to solve a problems without
actually solving any equations. The actual solution of the equation is left as exercise (this example under construction). What is the legitimacy of this method?
The explanation simply based the previous experience in which for a given ratio of
area or pressure ratio (etcetera) determines the Mach number. Based on the same
arguments, if it was shown that a group of parameters depends only Mach number
than the Mach is determined by this group.
A
The method of solution for given these parameters is by calculating the PP0 A
and then using the table to find the corresponding Mach number.
The case of 0 and T or P
The last case sometimes referred to as the naughty professors question case
dealt here is when the stagnation density given with the static temperature/pressure.
First, the dimensionless approach is used later analytical method is discussed (under construction).
c0 2
1
R0 P
z }| { 2
2
m
c0 2
F n2 P0
kRT0
=
=
=
k
P
kRP0 P0 PP0 A
kRP0 2 PP0 A
(5.69)
66
The last case dealt here is of the stagnation density with static pressure and the
following is dimensionless group
c0 2
1
R0 2 T
z }| {
2
2
kRT0 T0 m
c0 2 T0
F n 2 T0
=
=
=
k
T
kRP0 2 T A
kRP0 2 T A
(5.70)
It was hidden in the derivations/explanations of the above analysis didnt explicitly state under what conditions these analysis is correct. Unfortunately, not all the
analysis valid for the same conditions and is as the regular isentropic Table, (5.2).
The heat/temperature part is valid for enough adiabatic condition while the pressure condition requires also isentropic process. All the above conditions/situations
require to have the perfect gas model as the equation of state. For example the
first naughty professor question is sufficient that process is adiabatic only (T0 , P ,
mass flow rate per area.).
Table -5.1: Fliegners number and other parameters as a function of Mach number
Fn
0.00E+001.400E06
1.000
0.050001 0.070106 1.000
0.10000 0.14084 1.000
0.20000 0.28677 1.001
0.21000 0.30185 1.001
0.22000 0.31703 1.001
0.23000 0.33233 1.002
0.24000 0.34775 1.002
0.25000 0.36329 1.003
0.26000 0.37896 1.003
0.27000 0.39478 1.003
0.28000 0.41073 1.004
0.29000 0.42683 1.005
0.30000 0.44309 1.005
0.31000 0.45951 1.006
0.32000 0.47609 1.007
0.33000 0.49285 1.008
0.34000 0.50978 1.009
0.35000 0.52690 1.011
0.36000 0.54422 1.012
0.37000 0.56172 1.013
0.38000 0.57944 1.015
0.39000 0.59736 1.017
0.40000 0.61550 1.019
P0 A
AP
0.0
0.00747
0.029920
0.12039
0.13284
0.14592
0.15963
0.17397
0.18896
0.20458
0.22085
0.23777
0.25535
0.27358
0.29247
0.31203
0.33226
0.35316
0.37474
0.39701
0.41997
0.44363
0.46798
0.49305
RT0
P2
2
A
0.0
2.62E05
0.000424
0.00707
0.00865
0.010476
0.012593
0.015027
0.017813
0.020986
0.024585
0.028651
0.033229
0.038365
0.044110
0.050518
0.057647
0.065557
0.074314
0.083989
0.094654
0.10639
0.11928
0.13342
1
R0 P
2
A
0.0
0.00352
0.014268
0.060404
0.067111
0.074254
0.081847
0.089910
0.098460
0.10752
0.11710
0.12724
0.13796
0.14927
0.16121
0.17381
0.18709
0.20109
0.21584
0.23137
0.24773
0.26495
0.28307
0.30214
1
R0 2 T
2
A
0.0
0.00351
0.014197
0.059212
0.065654
0.072487
0.079722
0.087372
0.095449
0.10397
0.11294
0.12239
0.13232
0.14276
0.15372
0.16522
0.17728
0.18992
0.20316
0.21703
0.23155
0.24674
0.26264
0.27926
Fn
0.41000
0.42000
0.43000
0.44000
0.45000
0.46000
0.47000
0.48000
0.49000
0.50000
0.51000
0.52000
0.53000
0.54000
0.55000
0.56000
0.57000
0.58000
0.59000
0.60000
0.61000
0.62000
0.63000
0.64000
0.65000
0.66000
0.67000
0.68000
0.69000
0.70000
0.71000
0.72000
0.73000
0.74000
0.75000
0.76000
0.77000
0.78000
0.79000
0.63386
0.65246
0.67129
0.69036
0.70969
0.72927
0.74912
0.76924
0.78965
0.81034
0.83132
0.85261
0.87421
0.89613
0.91838
0.94096
0.96389
0.98717
1.011
1.035
1.059
1.084
1.109
1.135
1.161
1.187
1.214
1.241
1.269
1.297
1.326
1.355
1.385
1.415
1.446
1.477
1.509
1.541
1.574
1.021
1.023
1.026
1.028
1.031
1.035
1.038
1.042
1.046
1.050
1.055
1.060
1.065
1.071
1.077
1.083
1.090
1.097
1.105
1.113
1.122
1.131
1.141
1.151
1.162
1.173
1.185
1.198
1.211
1.225
1.240
1.255
1.271
1.288
1.305
1.324
1.343
1.362
1.383
P0 A
AP
0.51882
0.54531
0.57253
0.60047
0.62915
0.65857
0.68875
0.71967
0.75136
0.78382
0.81706
0.85107
0.88588
0.92149
0.95791
0.99514
1.033
1.072
1.112
1.152
1.194
1.236
1.279
1.323
1.368
1.414
1.461
1.508
1.557
1.607
1.657
1.708
1.761
1.814
1.869
1.924
1.980
2.038
2.096
RT0
P2
2
A
0.14889
0.16581
0.18428
0.20442
0.22634
0.25018
0.27608
0.30418
0.33465
0.36764
0.40333
0.44192
0.48360
0.52858
0.57709
0.62936
0.68565
0.74624
0.81139
0.88142
0.95665
1.037
1.124
1.217
1.317
1.423
1.538
1.660
1.791
1.931
2.081
2.241
2.412
2.595
2.790
2.998
3.220
3.457
3.709
1
R0 P
2
A
0.32220
0.34330
0.36550
0.38884
0.41338
0.43919
0.46633
0.49485
0.52485
0.55637
0.58952
0.62436
0.66098
0.69948
0.73995
0.78250
0.82722
0.87424
0.92366
0.97562
1.030
1.088
1.148
1.212
1.278
1.349
1.422
1.500
1.582
1.667
1.758
1.853
1.953
2.058
2.168
2.284
2.407
2.536
2.671
1
R0 2 T
2
A
0.29663
0.31480
0.33378
0.35361
0.37432
0.39596
0.41855
0.44215
0.46677
0.49249
0.51932
0.54733
0.57656
0.60706
0.63889
0.67210
0.70675
0.74290
0.78062
0.81996
0.86101
0.90382
0.94848
0.99507
1.044
1.094
1.147
1.202
1.260
1.320
1.382
1.448
1.516
1.587
1.661
1.738
1.819
1.903
1.991
68
Table -5.1: Fliegners number and other parameters as function of Mach number (continue)
Fn
0.80000
0.81000
0.82000
0.83000
0.84000
0.85000
0.86000
0.87000
0.88000
0.89000
0.90000
0.91000
0.92000
0.93000
0.94000
0.95000
0.96000
0.97000
0.98000
0.99000
1.000
1.607
1.642
1.676
1.712
1.747
1.784
1.821
1.859
1.898
1.937
1.977
2.018
2.059
2.101
2.144
2.188
2.233
2.278
2.324
2.371
2.419
1.405
1.427
1.450
1.474
1.500
1.526
1.553
1.581
1.610
1.640
1.671
1.703
1.736
1.771
1.806
1.843
1.881
1.920
1.961
2.003
2.046
P0 A
AP
2.156
2.216
2.278
2.340
2.404
2.469
2.535
2.602
2.670
2.740
2.810
2.882
2.955
3.029
3.105
3.181
3.259
3.338
3.419
3.500
3.583
RT0
P2
2
A
3.979
4.266
4.571
4.897
5.244
5.613
6.006
6.424
6.869
7.342
7.846
8.381
8.949
9.554
10.20
10.88
11.60
12.37
13.19
14.06
14.98
1
R0 P
2
A
1
R0 2 T
2.813
2.963
3.121
3.287
3.462
3.646
3.840
4.043
4.258
4.484
4.721
4.972
5.235
5.513
5.805
6.112
6.436
6.777
7.136
7.515
7.913
2
A
2.082
2.177
2.277
2.381
2.489
2.602
2.720
2.842
2.971
3.104
3.244
3.389
3.541
3.699
3.865
4.037
4.217
4.404
4.600
4.804
5.016
Example 5.5:
A gas flows in the tube with mass flow rate of 0.1 [kg/sec] and tube cross section
is 0.001[m2 ]. The temperature at Chamber supplying the pressure to tube is 27 C.
At some point the static pressure was measured to be 1.5[Bar]. Calculate for that
point the Mach number, the velocity, and the stagnation pressure. Assume that the
process is isentropic, k = 1.3, R = 287[j/kgK].
S OLUTION
The first thing that need to be done is to find the mass flow per area and it is
m
= 0.1/0.001 = 100.0[kg/sec/m2 ]
A
It can be noticed that the total temperature is 300K and the static pressure is
1.5[Bar]. The solution is based on section equations (5.60) through (5.65). It is
fortunate that Potto-GDC exist and it can be just plug into it and it provide that
M
T
T0
A
A
P
P0
AP
A P0
0.98116 3.4102
F
F
1.5392
P
= 1.5/0.98116 = 1.5288[Bar]
P/P0
2(k1)
k1 2
kP0 A
1+
M
= constant
(5.71)
RT0
2
From equation (5.71), it is clear that the function f (P0 , T0 , A ) = constant. There
are two possible models that can be used to simplify the calculations. The first
model for neglected heat transfer (adiabatic) flow and in which the total temperature
remained constant (Fanno flow like). The second model which there is significant
heat transfer but insignificant pressure loss (Rayleigh flow like).
If the mass flow rate is constant at any point on the tube (no mass loss occur) then
s
k+1
k1
k
2
m
=A
P0
(5.72)
RT0 k + 1
For adiabatic flow, comparison of mass flow rate at point A and point B leads to
P0 A |A = P0 A |B
;
And utilizing the equality of A =
A |
P0 |A
= A
P0 |B
A |B
(5.73)
A
A A
P0 |A
P0 |B
leads to
A
A MA A|A
= A
A|B
A M
B
(5.74)
70
For a flow with a constant stagnation pressure (frictionless flow) and non adiabatic
flow reads
" B
#2
T 0 |A
A MB A|B
(5.75)
=
A
T 0 |B
A|A
A M
A
Example 5.6:
At point A of the tube the pressure is 3[Bar], Mach number is 2.5, and the duct
section area is 0.01[m2 ]. Downstream at exit of tube, point B, the cross section
area is 0.015[m2 ] and Mach number is 1.5. Assume no mass lost and adiabatic
steady state flow, calculated the total pressure lost.
S OLUTION
Both Mach numbers are known, thus the area ratios can be calculated. The total pressure can be calculated because the Mach number and static pressure are
known. With these information, and utilizing equation (5.74) the stagnation pressure at point B can be obtained.
M
1.5000
2.5000
T
T0
A
A
P
P0
AP
A P0
F
F
3
P
= 51.25781291[Bar]
P0 |A = =
P
0.058527663
P0
| {z }
M =2.5 A
1.1761671
0.01
15.243[Bar]
2.6367187 0.015
Hence
P0 |A P0 |B = 51.257 15.243 = 36.013[Bar]
Note that the large total pressure loss is much larger than the static pressure loss
(Pressure point B the pressure is 0.27240307 15.243 = 4.146[Bar]).
71
M
0.000
0.050
0.100
0.200
0.300
0.400
0.500
0.600
0.700
0.800
0.900
1.00
1.100
1.200
1.300
1.400
1.500
1.600
1.700
1.800
1.900
2.000
2.500
3.000
3.500
4.000
4.500
5.000
5.500
6.000
6.500
7.000
7.500
8.000
8.500
9.000
9.500
10.00
T
T0
1.00000
0.99950
0.99800
0.99206
0.98232
0.96899
0.95238
0.93284
0.91075
0.88652
0.86059
0.83333
0.80515
0.77640
0.74738
0.71839
0.68966
0.66138
0.63371
0.60680
0.58072
0.55556
0.44444
0.35714
0.28986
0.23810
0.19802
0.16667
0.14184
0.12195
0.10582
0.092593
0.081633
0.072464
0.064725
0.058140
0.052493
0.047619
1.00000
0.99875
0.99502
0.98028
0.95638
0.92427
0.88517
0.84045
0.79158
0.73999
0.68704
0.63394
0.58170
0.53114
0.48290
0.43742
0.39498
0.35573
0.31969
0.28682
0.25699
0.23005
0.13169
0.076226
0.045233
0.027662
0.017449
0.011340
0.00758
0.00519
0.00364
0.00261
0.00190
0.00141
0.00107
0.000815
0.000631
0.000495
A
A
5.8E+5
11.59
5.822
2.964
2.035
1.590
1.340
1.188
1.094
1.038
1.009
1.000
1.008
1.030
1.066
1.115
1.176
1.250
1.338
1.439
1.555
1.688
2.637
4.235
6.790
10.72
16.56
25.00
36.87
53.18
75.13
1.0E+2
1.4E+2
1.9E+2
2.5E+2
3.3E+2
4.2E+2
5.4E+2
P
P0
AP
A P0
1.0000
5.8E + 5
0.99825 11.57
0.99303
5.781
0.97250
2.882
0.93947
1.912
0.89561
1.424
0.84302
1.130
0.78400
0.93155
0.72093
0.78896
0.65602
0.68110
0.59126
0.59650
0.52828
0.52828
0.46835
0.47207
0.41238
0.42493
0.36091
0.38484
0.31424
0.35036
0.27240
0.32039
0.23527
0.29414
0.20259
0.27099
0.17404
0.25044
0.14924
0.23211
0.12780
0.21567
0.058528 0.15432
0.027224 0.11528
0.013111 0.089018
0.00659
0.070595
0.00346
0.057227
0.00189
0.047251
0.00107
0.039628
0.000633 0.033682
0.000385 0.028962
0.000242 0.025156
0.000155 0.022046
0.000102 0.019473
6.90E5 0.017321
4.74E5 0.015504
3.31E5 0.013957
2.36E5 0.012628
F
F
2.4E+5
4.838
2.443
1.268
0.89699
0.72632
0.63535
0.58377
0.55425
0.53807
0.53039
0.52828
0.52989
0.53399
0.53974
0.54655
0.55401
0.56182
0.56976
0.57768
0.58549
0.59309
0.62693
0.65326
0.67320
0.68830
0.69983
0.70876
0.71578
0.72136
0.72586
0.72953
0.73257
0.73510
0.73723
0.73903
0.74058
0.74192
72
5.3.1
(5.76)
RT dP
=0
P
(5.77)
Integration of equation (5.77) yields the Bernoullis equation for ideal gas in isothermal process which reads
;
P2
U2 2 U1 2
+ RT ln
=0
2
P1
2RT ln
P2
U1 2
P1
(5.78)
(5.79)
(5.80)
(5.81)
kR
T
P2
ln
;
T
= R
2
P1
(5.82)
6 The one dimensional momentum equation for steady state is U dU/dx = dP/dx+0(other effects)
which are neglected here.
73
k(M2 2 M1 2 )
P2
= ln
2
P1
(5.83)
eM1
2
eM2
! k2
(5.84)
As oppose to the adiabatic case (T0 = constant) in the isothermal flow the stagnation temperature ratio can be expressed
1
T01
T1 1 +
=
T02
T2 1 +
2
k1
2 M1
2
k1
2 M2
1+
=
1+
2
k1
2 M1
2
k1
2 M2
(5.85)
(5.86)
eM2
2
eM1
! k2
(5.87)
1+
1+
2
k1
2 M2
2
k1
2 M1
k
! k1
"
eM1
=
2
eM1
# k2
(5.88)
The critical point, at this stage, is unknown (at what Mach number the nozzle is
choked is unknown) so there are two possibilities: the choking point or M = 1 to
normalize the equation. Here the critical point defined as the point where M = 1
so results can be compared to the adiabatic case and denoted by star. Again it
has to emphasis that this critical point is not really related to physical critical point
but it is arbitrary definition. The true critical point is when flow is choked and the
relationship between two will be presented.
The critical pressure ratio can be obtained from (5.84) to read
(1M 2 )k
P
= =e 2
(5.89)
74
A
M
(5.90)
k
2 k1
2 1 + k1
T0
2 M1
=
T0
k+1
(5.91)
k
! k1
2
2 1 + k1
M
1
2
k+1
(5.92)
k1
k
P0
2
= e2
(5.93)
P0 M =0
k+1
For specific heat ratio of k = 1.4, this maximum value is about two. It can be noted
that the stagnation pressure is monotonically reduced during this process.
Of course in isothermal process T = T . All these equations are plotted in Figure
(5.6). From the Figure 5.3 it can be observed that minimum of the curve A/A isnt
on M = 1. The minimum of the curve is when area is minimum and at the point
where the flow is choked. It should be noted that the stagnation temperature is not
constant as in the adiabatic case and the critical point is the only one constant.
The mathematical procedure to find the minimum is simply taking the derivative
and equating to zero as following
d
A
A
dM
kM 2 e
k(M 2 1)
2
k(M 2 1)
2
M2
=0
(5.94)
(5.95)
Uthroatmax = RT
(5.96)
75
Isothermal Nozzle
k=14
4
*
P/P
*
A/A
*
P0 / P0
3.5
3
T 0 / T0
T/T
2.5
2
1.5
1
0.5
0
0.5
2
M
1.5
2.5
3.5
On the other hand, the pressure loss in adiabatic flow is milder as can be seen in
Figure (5.7(a)).
It should be emphasized that the stagnation pressure decrees. It is convenient to
find expression for the ratio of the initial stagnation pressure (the stagnation pressure before entering the nozzle) to the pressure at the throat. Utilizing equation
(5.89) the following relationship can be obtained
P Pthroat
Pthroat
=
=
P0initial
P0initial P
1
(102 )k
2
k
1 1
2
k
e 2 = 0.60653
(5.97)
Notice that the critical pressure is independent of the specific heat ratio, k, as
opposed to the adiabatic case. It also has to be emphasized that the stagnation
values of the isothermal model are not constant. Again, the heat transfer is ex-
76
Isothermal Nozzle
k=14
k=14
5
4.5
3.5
M isoT
M isentropic
Uisntropic/UisoT
4
3
3.5
2.5
2.5
*
A / A iso
*
A / A adiabatic
*
P / P iso
*
P / P adiabatic
1.5
1
1.5
1
0.5
0
0.5
0
0.5
1.5
2
M
2.5
3.5
1
0.5
1.5
Distance (normalized distance two scales)
pressed as
Q = Cp (T02 T02 )
(5.98)
77
For
comparison
between
Comparison between the two models
the adiabatic model and the
k=14
isothermal a simple profile of
nozzle area as a function of
1
the distance is assumed. This
P / P0 isentropic
profile isnt an ideal profile
0.8
T / T0 isentropic
but rather a simple sample
P / P0 isothermal
T/T0 isothermal
just to examine the difference
0.6
between the two models so in
an actual situation it can be
0.4
bounded. To make sense and
eliminate unnecessary details
0.2
the distance from the entrance
to the throat is normalized (to
0
0
1
2
0.5
1.5
one (1)). In the same fashion
Distance (normalized distance two scales)
the distance from the throat to
Fri Apr 8 15:11:44 2005
the exit is normalized (to one
(1)) (it doesnt mean that these
distances are the same). In this Fig. -5.8: Comparison of the pressure and temperature
drop as a function of the normalized length
comparison the entrance area
(two scales)
ratio and the exit area ratio are
the same and equal to 20. The Mach number was computed for the two models
and plotted in Figure (5.7(b)). In this comparison it has to be remembered that
critical area for the two models are different by about 3% (for k = 1.4). As can be
observed from Figure (5.7(b)). The Mach number for the isentropic is larger for
the supersonic branch but the velocity is lower. The ratio of the velocities can be
expressed as
Ms kRTs
Us
(5.99)
=
UT
MT kRTs
It can be noticed that temperature in the isothermal model is constant while temperature in the adiabatic model can be expressed as a function of the stagnation
temperature. The initial stagnation temperatures are almost the same and can be
canceled out to obtain
Us
Ms
q
UT
2
MT 1 + k1
2 Ms
(5.100)
By utilizing equation (5.100) the velocity ratio was obtained and is plotted in Figure
(5.7(b)).
Thus, using the isentropic model results in under prediction of the actual results for
the velocity in the supersonic branch. While, the isentropic for the subsonic branch
will be over prediction. The prediction of the Mach number are similarly shown in
Figure (5.7(b)).
78
Two other ratios need to be examined: temperature and pressure. The initial stagnation temperature is denoted as T0int . The temperature ratio of T /T0int can be
obtained via the isentropic model as
T
T0 int
1
1+
k1
2
2 M
(5.101)
While the temperature ratio of the isothermal model is constant and equal to one
(1). The pressure ratio for the isentropic model is
P
P0int
1
1+
k1
2
2 M
k1
k
(5.102)
and for the isothermal process the stagnation pressure varies and has to be taken
into account as the following:
isentropic
Pz
P0 P0z
=
P0int
P0int P0
z}|{
Pz
P0z
(5.103)
where z is an arbitrary point on the nozzle. Using equations (5.88) and the isentropic relationship, the sought ratio is provided.
Figure (5.8) shows that the range between the predicted temperatures of the two
models is very large, while the range between the predicted pressure by the two
models is relatively small. The meaning of this analysis is that transferred heat
affects the temperature to a larger degree but the effect on the pressure is much
less significant.
To demonstrate the relativity of the approach advocated in this book consider the
following example.
Example 5.7:
Consider a divergingconverging nozzle made out of wood (low conductive material) with exit area equal entrance area. The throat area ratio to entrance area is
1:4 respectively. The stagnation pressure is 5[Bar] and the stagnation temperature
is 27 C. Assume that the back pressure is low enough to have supersonic flow
without shock and k = 1.4. Calculate the velocity at the exit using the adiabatic
model. If the nozzle was made from copper (a good heat conductor) a larger heat
transfer occurs, should the velocity increase or decrease? What is the maximum
possible increase?
S OLUTION
The first part of the question deals with the adiabatic model i.e. the conservation
of the stagnation properties. Thus, with known area ratio and known stagnation
PottoGDC provides the following table:
T
T0
0.99572
0.36644
79
0.98934
0.08129
A
A
4.0000
4.0000
P
P0
0.98511
0.02979
AP
A P0
3.9405
0.11915
With the known Mach number and temperature at the exit, the velocity can be
calculated. The exit temperature is 0.36644 300 = 109.9K. The exit velocity, then,
is
T
T0
1.9910
1.4940
0.51183
A
A
4.0000
P
P0
0.12556
AP
A P0
0.50225
The exit Mach number is known and the initial temperature to the throat temperature ratio can be calculated as the following:
1
1
T0ini
=
= 0.777777778
=
k1
1
T0
1+ 2 k
1 + k1
k
Thus the stagnation temperature at the exit is
T0ini
= 1.4940/0.777777778 = 1.921
T0exit
The exit stagnation temperature is 1.92 300 = 576.2K. The exit velocity can be
determined by utilizing the following equation
(5.104)
80
The net force is denoted here as Fnet . The mass conservation also can be applied
to our control volume
m
= 1 A1 U1 = 2 A2 U2
(5.105)
Combining equation (5.104) with equation (5.105) and by utilizing the identity in
equation (5.42) results in
Fnet = kP2 A2 M2 2 kP1 A1 M1 2 + P2 A2 P1 A1
(5.106)
f (M1 )
f (M1 )
2)
z }| { z f (M
}|
{ z }| { z
}|
{
P2 A2
P
A
Fnet
1
1
2
2
=
1
+
kM
1
+
kM
2
1
P 0 A
P0 A
P0 A
(5.107)
F = P A 1 + kM2 2
(5.108)
In the Impulse function when F (M = 1) is denoted as F
F = P A (1 + k)
(5.109)
F
P1 A1 1 + kM1 2
=
=
F
P A
(1 + k)
1
P
P0
|{z}
P1 A1
1
2
1 + kM1
P0 A
(1 + k)
(5.110)
2
) k1
( k+1
This ratio is different only in a coefficient from the ratio defined in equation (5.107)
which makes the ratio a function of k and the Mach number. Hence, the net force
is
F1
k + 1 k1 F2
(5.111)
Fnet = P0 A (1 + k)
2
F
F
81
S OLUTION
The solution is obtained by getting the data for the Mach number. To obtained the
Mach number, the ratio of P1 A1 /A P0 is needed to be calculated. To obtain this
ratio the denominator is needed to be obtained. Utilizing Fliegners equation (5.51),
provides the following
m
RT
1.0 400 287
A P0 =
=
70061.76[N ]
0.058
0.058
and
A2 P2
500000 0.003
=
2.1
A P0
70061.76
M
T
T0
A
A
A
A
P
P0
AP
A P0
0.94934 2.1000
F
F
0.96666
A2 A1
0.009
A1
=
= 2.2121
= 5.2227
A
A A2
0.003
And utilizing again Potto-GDC provides
M
T
T0
A
A
P
P0
AP
A P0
0.99132 5.1774
F
F
2.1949
P0 P1
= 5.0times0.94934/0.99380 4.776[Bar]
P2 P0
82
P 0 A
F1
k + 1 k1 F2
= P2 A2
(1 + k)
P2 A2
2
F
F
1
3.5
= 500000
2.4 1.2 (2.1949 0.96666) 614[kN ]
2.1
(5.112)
Utilizing the ideal gas model for density and some rearrangement results in
P2 A2 1 +
F2
=
F1
P1 A1 1 +
U2 2
RT
U1 2
RT
(5.113)
Since U 2 /RT = kM 2 and the ratio of equation (5.86) transformed equation into
(5.113)
F2
M1 1 + kM2 2
=
F1
M2 1 + kM1 2
(5.114)
F
M2 1 + k
(5.115)
M
0.00
T0
T0
0.52828
P0
P0
1.064
A
A
5.0E + 5
P
P
2.014
AP
A P0
F
F
1.0E+6 4.2E+5
83
M
0.05
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.00
1.10
1.20
1.30
1.40
1.50
1.60
1.70
1.80
1.90
2.00
2.50
3.000
3.500
4.000
4.500
5.000
5.500
6.000
6.500
7.000
7.500
8.000
8.500
9.000
9.500
10.00
T0
T0
0.52921
0.53199
0.54322
0.56232
0.58985
0.62665
0.67383
0.73278
0.80528
0.89348
1.000
1.128
1.281
1.464
1.681
1.939
2.245
2.608
3.035
3.540
4.134
9.026
19.41
40.29
80.21
1.5E + 2
2.8E + 2
4.9E + 2
8.3E + 2
1.4E + 3
2.2E + 3
3.4E + 3
5.2E + 3
7.7E + 3
1.1E + 4
1.6E + 4
2.2E + 4
P0
P0
A
A
1.064
9.949
1.064
5.001
1.064
2.553
1.063
1.763
1.062
1.389
1.059
1.183
1.055
1.065
1.047
0.99967
1.036
0.97156
1.021
0.97274
1.000
1.000
0.97376
1.053
0.94147
1.134
0.90302
1.247
0.85853
1.399
0.80844
1.599
0.75344
1.863
0.69449
2.209
0.63276
2.665
0.56954
3.271
0.50618
4.083
0.22881 15.78
0.071758 90.14
0.015317 7.5E + 2
0.00221
9.1E + 3
0.000215 1.6E + 5
1.41E5 4.0E + 6
0.0
1.4E + 8
0.0
7.3E + 9
0.0
5.3E+11
0.0
5.6E+13
0.0
8.3E+15
0.0
1.8E+18
0.0
5.4E+20
0.0
2.3E+23
0.0
1.4E+26
0.0
1.2E+29
P
P
AP
A P0
2.010
20.00
2.000
10.00
1.958
5.000
1.891
3.333
1.800
2.500
1.690
2.000
1.565
1.667
1.429
1.429
1.287
1.250
1.142
1.111
1.000
1.000
0.86329
0.90909
0.73492
0.83333
0.61693
0.76923
0.51069
0.71429
0.41686
0.66667
0.33554
0.62500
0.26634
0.58824
0.20846
0.55556
0.16090
0.52632
0.12246
0.50000
0.025349 0.40000
0.00370
0.33333
0.000380 0.28571
2.75E5 0.25000
1.41E6 0.22222
0.0
0.20000
0.0
0.18182
0.0
0.16667
0.0
0.15385
0.0
0.14286
0.0
0.13333
0.0
0.12500
0.0
0.11765
0.0
0.11111
0.0
0.10526
0.0
0.100000
F
F
8.362
4.225
2.200
1.564
1.275
1.125
1.044
1.004
0.98750
0.98796
1.000
1.020
1.047
1.079
1.114
1.153
1.194
1.237
1.281
1.328
1.375
1.625
1.889
2.161
2.438
2.718
3.000
3.284
3.569
3.856
4.143
4.431
4.719
5.007
5.296
5.586
5.875
84
(5.116)
Again, Gibbs equation (5.27) is reused to related the entropy change to the
change in thermodynamics properties and applied on non-ideal gas. Since ds = 0
and utilizing the equation of the state dh = dP/. The enthalpy is a function of the
temperature and pressure thus, h = h(T, P ) and full differential is
h
h
dT +
dP
(5.117)
dh =
T P
P T
h
The definition of pressure specific heat is Cp T
and second derivative is
Maxwell relation hence,
h
s
=vT
(5.118)
P T
T P
First, the differential of enthalpy is calculated for real gas equation of state as
dP
z
T
(5.119)
dh = Cp dT
Z
T P
Equations (5.27) and (4.19) are combined to form
ds
z
dP
Cp dT
T
=
z 1+
R
R T
Z
T P P
The mechanical energy equation can be expressed as
Z
Z 2
dP
U
=
d
2
(5.120)
(5.121)
At the stagnation the definition requires that the velocity is zero. To carry the
integration of the right hand side the relationship between the pressure and the
density has to be defined. The following power relationship is assumed
=
0
P
P0
n1
(5.122)
85
Notice, that for perfect gas the n is substituted by k. With integration of equation
(5.121) when using relationship which is defined in equation (5.122) results
U2
=
2
P1
dP
=
P0
P
P0
1
0
P0
P
n1
dP
(5.123)
z0 RT0
P0
P0
P0
P
n1
dP
(5.124)
2z0 RT0 ln
U=
P0
P
(5.125)
(5.126)
It must be noted that n is a function of the critical temperature and critical pressure.
The mass flow rate is regardless to equation of state as following
m
= A U
(5.127)
Where is the density at the throat (assuming the chocking condition) and A is
the cross area of the throat. Thus, the mass flow rate in our properties
U
}|
{
z
z
}|
{v
"
n1 #
n1 u
(
)
u
n
P
P0
tz RT 2n 1 P
m
= A
0
0
z0 RT0 P0
n1
P0
(5.128)
P0
m
= A
z0 RT0
}|
P
P0
z
{s
1
n
}|
2z0 RT0 ln
P0
P
(5.129)
The Mach number can be obtained by utilizing equation (4.34) to defined the Mach
number as
M=
U
znRT
(5.130)
86
T2
T1
Cp dT
=
R T
P2
P1
dP
z
T
z 1+
Z
T P P
(5.131)
P/P0
0
d
P
P0
(5.133)
It was shown in Chapter (4) that (4.33) is applicable for some ranges of relative
temperature and pressure (relative to critical temperature and pressure and not the
stagnation conditions).
v
#
u
"
n1
u
n
2n
P
U = tz0 RT0
(5.134)
1
n1
P0
When n = 1 or when n 1
U=
2z0 RT0 ln
P0
P
(5.135)
2n
n1
A P 0
m
=
z0 RT0
2n
n1
P
P0
n1
(5.136)
P0
P
(5.137)
P
P0
(5.138)
P
P0
And for n = 1
2z0 RT0 ln
2n
n1
P
P0
n1
87
mc
0
=2
A P0
P
P0
ln
P
P0
(5.139)
2
n+1
n
n1
(5.140)
P
= e
P0
(5.141)
To obtain the relationship between the temperature and pressure, equation (5.131)
can be integrated
T0
=
T
z
CR [z+T ( T
)P ]
p
P0
P
(5.142)
1n #
u
2
z 0 T0
P 0 n
t
M=
1
n1 z T
P
(5.144)
M=
z0 T0 P0
ln
z T
P
(5.145)
The pressure ratio at any point can be expressed as a function of the Mach number
as
z
n1 z+T ( T
)P ]
n 1 2 ( n )[
T0
(5.146)
= 1+
M
T
2
for n = 1
T0
=
T
eM [z+T (
z
T
)P ]
(5.147)
88
)P ]
T
1 + n ( 1n )[
=
T0
2
(5.148)
and for n = 1
T
=
T0
e[z+T (
z
T
)P ]
(5.149)
n+1
n 1 2 n1
P0 n
1+
M
M
m
=
c0
2
(5.150)
M2
n 1 2 n1
P 0 A n
M
1+
m
=
c0
2
(5.151)
Example 5.9:
A design is required that at a specific point the Mach number should be M = 2.61,
the pressure 2[Bar], and temperature 300K.
i. Calculate the area ratio between the point and the throat.
ii. Calculate the stagnation pressure and the stagnation temperature.
iii. Are the stagnation pressure and temperature at the entrance different from the
point? You can assume that k = 1.405.
S OLUTION
1. The solution is simplified by using Potto-GDC for M = 2.61 the results are
M
2.6100
T
T0
0.42027
0.11761
A
A
2.9066
P
P0
0.04943
AP
A P0
0.14366
T0 =
T0
300
T =
713.82K
T
0.42027