Nozzels
Nozzels
Nozzels
13
Nozzles
13.1 INTRODUCTION
A nozzle is a flow passage of varying cross sectional area in which the velocity of fluid increases and
pressure drops in the direction of flow. Thus in nozzle the fluid enters the variable cross section area
duct with small velocity and high pressure and leaves it with high velocity and small pressure. During
flow through nozzle the enthalpy drops and heat drop in expansion is spent in increasing the velocity of
fluid. Similar to nozzle a duct with variable cross-section area will be called diffuser if the fluid gets
decelerated, causing a rise in pressure along the direction of flow. Nozzles are generally used in turbines,
jet engines, rockets, injectors, ejectors etc.
Fig. 13.2
The one dimensional steady flow through a duct is shown above. For control volume shown the
principle of conservation of mass, energy and Newton’s law of motion may be applied.
By Newton’s law of motion, F = m × a where F is the resultant force acting on system of mass
‘m’ and acceleration ‘a’.
Pressure waves and Sound waves: Let us consider a cylindrical duct having piston on one end for
generating the pressure wave. Figure 13.3 shows the arrangement for producing a pressure wave
moving to right with velocity ‘a’. Sound wave is also a small pressure disturbance that propagates
through a gas, liquid or solid with velocity ‘a’ that depends on the properties of medium.
Figure 13.3 shows how the generation of pressure wave causes displacement of fluid thereby
causing rise in pressure, density and temperature by Dp, Dr and DT in respect to the region on the right
of wave (undisturbed region). In the undisturbed region say pressure, density, temperature and fluid
velocity be p, r, T and C = 0 respectively. Due to piston movement fluid velocity increases by DC and
other properties also change by elemental values as shown. For analysing there are two approaches
available as shown in Figs. 13.3 (a) and (b). One approach considers observer to be stationary and gas
moving and second approach considers observer to be moving along with wave i.e. relative velocity of
observer with respect to wave is zero.
566 ________________________________________________________ Applied Thermodynamics
Respective values of fluid velocity, wave propagation velocity, pressure, density and temperature
are labelled on figure. For an observer at rest relative to wave (observer moving with wave) it seems as
if the fluid is moving towards the stationary wave from right with velocity a, pressure p, density r and
temperature T and moves away on left with velocity ‘a – DC’, pressure ‘p + Dp’, density ‘r + Dr’ and
temperature ‘T + DT’.
From conservation of mass, applying continuity equation upon control volume we get
m¢1 = m¢2 = m¢
r ×A×a = (r + Dr)×A(a – DC)
where A is constant cross section area of duct.
r×A×a = (r×A×a) – (r×A×DC) + (Dr×A×a) – (Dr×A×DC)
Upon neglecting higher order terms and rearranging we get,
(a×Dr) – (r×DC) = 0
a×Dr
or, DC =
r
Applying momentum equation to the control volume;
(p×A) – ((p + Dp)×A) = (m¢ (a – DC)) – (m¢×a)
– Dp × A = m¢×(– DC)
for mass flow rate m¢ we can write, m¢ = r×A×a
so, Dp×A = r×A×a×DC
Dp
or, DC =
r×a
Equating two values obtained for ‘DC’ we get
a×Dr Dp
=
r r×a
Dp
a=
Dr
Thus, velocity of wave propagation comes out as the square root of the ratio of change in
pressure and change in density.
In case of sound waves the magnitude of changes in pressure, density and temperature are
infinitesimal and so these may also be called as infinitesimal pressure wave. It is also seen that thermo-
dynamic process occurring across an infinitesimal pressure wave may be considered nearly isentropic.
Therefore the velocity of sound can be given as square root of derivative of pressure with respect to
density across the wave under isentropic conditions.
∂ p
a=
∂ ρ s =constt.
− dv
in terms of specific volume values; dr =
v2
∂ p
so, a= − v2
∂ v s =constt
Nozzles _______________________________________________________________________ 567
Let us consider fluid to be a perfect gas following isentropic process given by pvk = constt. Taking
log of both sides and then partially differentiating we get,
FG ¶ p IJ -k× p
H ¶v K s
=
v
Substituting in expression for sound velocity
a= k pv
For ideal gas,
a= k R T . In case of air, a = g RT
Using the velocity of sound and fluid velocity a non dimensional parameter called Mach number
can be defined. Mach number is given by the ratio of velocity of fluid (object) to the velocity of sound.
It is generally denoted by M.
C
M=
a
Based upon Mach no. value flow can be classified as given below.
For
M < 1 flow is called subsonic flow.
M = 1 flow is called sonic flow.
M > 1 flow is called supersonic flow.
Nozzle flow analysis: Let us consider one dimensional steady flow in nozzles. Let us take a varying
cross-section area duct such that velocity increases and pressure decreases from inlet to exit.
Fig. 13.4
From conservation of mass, upon applying continuity equation, it can be given that,
r×A×C = constant
Taking log of both the sides,
ln r + ln A + ln C = ln constant
Differentiating partially we get,
dr d A dC
+ + =0
r A C
Let us now apply steady flow energy equation assuming no change in potential energy, no heat
interaction and no work interaction during the flow across control volume.
S.F.E.E. yields, dq = dh + dw + d(KE) + d(PE)
Applying assumptions,
568 ________________________________________________________ Applied Thermodynamics
dh + d(KE) = 0
C2 RS UV
It can be rewritten for section 1 and 2 as, KE = 2
T W
C12 C2
or h1 + = h2 + 2
2 2
ho1 = ho2
Stagnation enthalpy at section 1 = Stagnation enthalpy at section 2.
C2F I
From differential form, dh + d
2H K =0
C dC
dh + 2 =0
2
or, dh = – CdC
From first and second law combined we know,
dh = Tds + vdp
Using the adiabatic flow considerations, ds = 0, so
dp
dh = vdp =
r
Above shows that with increase or decrease in pressure along the direction of flow the specific
enthalpy also change in same way.
From thermodynamic property relations pressure can be given as function of density and entropy
i.e. p = p(r, s).
FG ¶ p IJ × dr + FG ¶ p IJ
or, dp =
H ¶r K s
H ¶s K r
× ds
dp = a2×dr
Combining two expressions for dh we get
dp
– C × dC = . This shows that as pressure increases in direction of flow then velocity must
r
decrease. Substituting from dp as obtained above, it yields,
a2 × d r
– C × dC =
r
Nozzles _______________________________________________________________________ 569
or,
FH dC IK = a FH d r IK
2
F dr I = - C e dC j
2
-
C C 2
r
Þ
H rK a C 2
Fig. 13.5
dC
Diffuser: For negative velocity gradient i.e., is − ve the velocity decreases along the direc-
C
dA
tion of flow as in case of diffuser, it yields as + ve. Positive area gradient means duct has diverging
A
cross section area along the direction of flow.
570 ________________________________________________________ Applied Thermodynamics
Case 2 For supersonic flow i.e. M > 1
dC dA
Nozzle: For positive velocity gradient i.e. being + ve, it yields as + ve. It means that in
C A
supersonic flow the nozzle duct shall have diverging cross-sectional area along the direction of flow.
Fig. 13.6
dC dA
Diffuser: For negative velocity gradient i.e. being – ve it yields as – ve. It means in
C A
supersonic flow the diffuser duct shall have converging cross-sectional area along the direction of flow.
From above discussion it can be concluded that
(i) Nozzle must be of convergent duct type in subsonic flow region and such nozzles are called
subsonic nozzles or convergent nozzles.
(ii) Nozzle must be of divergent duct type in supersonic flow region and such nozzles are called
supersonic nozzles or divergent nozzles.
(iii) For acceleration of fluid flow from subsonic to supersonic velocity the nozzle must be first
of converging type till flow becomes sonic and subsequently nozzle should be of diverging
type in supersonic flow. The portion of duct at which flow becomes sonic (M = 1) and dA
is zero i.e. duct is constant cross-section area duct, is called throat. Thus in this type of
flow from subsonic to supersonic the duct is of converging type followed by throat and a
diverging duct. Such nozzles are also called convergent-divergent nozzles. Throat gives the
minimum cross-section area in convergent-divergent nozzles.
Let us consider the expansion through a nozzle between sections 1 and 2. In nozzle the velocity
of fluid is so high that there is no time available for heat exchange with the surroundings and the
expansion may be considered adiabatic. Also the change in potential energy may be negligible if the
elevation does not change from inlet to exit. Work done during flow is absent.
Application of steady flow energy equation yields,
C12 C2
h1 + = h2 + 2
2 2
Velocity at exit from nozzle:
C2 = a f
2 h1 - h2 + C12 , m/s
For negligible velocity of fluid at inlet to nozzle, C1 » 0
C2 = a f
2 h1 - h2 , m/s, where h1 and h2 are enthalpy in J/kg at sections 1 and 2 respectively.
Nozzles _______________________________________________________________________ 571
Expansion of fluid on p-v diagram is shown below.
Fig. 13.9 T-s and h-s representation for steam flow through nozzle
In above representations the isentropic heat drop shown by 1–2 and 3–4 is also called ‘Rankine
heat drop’.
Mass flow through a nozzle can be obtained from continuity equation between sections 1 and 2.
A1C1 AC
m¢ = = 2 2
v1 v2
m¢ C2
Mass flow per unit area; =
A2 v2
From different from of S.F.E.E.
dq = dh + dw + d(K.E.) + d(P.E.)
or, dh + d(K.E.) = 0
572 ________________________________________________________ Applied Thermodynamics
du + pdv + vdp + d(K.E.) = 0
also as dq = du + pdv = 0, so d(K.E.) = – vdp
z
p2
C22 - C12
or = - vdp
2
p1
For the expansion through a nozzle being governed by process pvn = constt.,
FG n IJ p v FG1 - p v IJ
H n -1K H p v K
2 2
C22 – C12 = 2 1 1
1 1
or,
FH n IK p v FG1 - p v IJ
H pv K
2 2
C2 = 2 1 1
n -1 1 1
If the working fluid is perfect gas then n = g and for air g = 1.4. However, if working fluid is steam
a good approximation for n can be obtained from some polytropic considerations. For steam being dry
saturated initially and process of expansion occurring in wet region the index n can be approximated as
1.135. For steam being initially superheated and expanded in superheated region the index n can be
approximated as 1.3.
Looking at mathematical expression for exit velocity it could be concluded that maximum exit
velocity is possible only when fluid is expanded upto zero pressure. The maximum velocity is,
FG n IJ p v
Cmax = 2
H n -1K 1 1
n p2 v 2
2 p1v1 1 −
m¢ n −1 p1v1
=
A2 v2
F p IJ
=G
1/ n
Hp K × v1
1
or, v2
2
Nozzles _______________________________________________________________________ 573
or,
FG n IJ p v FG 1 - p v IJ
H n -1K H p v K
2 2
2 1 1
m¢
FG p IJ × v
1 1
= 1/ n
A2
Hp K
1
1
2
or,
LM F Ip R|F p I F p I FH
( n + 1) I U O1/ 2
K |P
S|H p K V|P
2/n
m¢
MN H Kv HpK
n n
= 2 1 2
- 2
A2 n -1 1
T 1 1
WQ
This expression for mass flow rate through nozzle depends upon inlet and exit pressures, initial
specific volume and index of expansion. It has been seen earlier that the mass flow per unit area is
maximum at throat and nozzle should be designed for maximum discharge per unit area. Thus there will
be some value of throat pressure (p2) which offers maximum discharge per unit area. Mathematically
this pressure value can be obtained by differentiating expression for mass flow per unit area and equat-
ing it to zero. This pressure at throat for maximum discharge per unit area is also called ‘critical
pressure’ and pressure ratio with inlet pressure is called ‘critical pressure ratio’.
p
Let pressure ratio 2 = r, then mass flow per unit area can be re-written as;
p1
RS F I p FG r ( n + 1)
I UV 1/ 2
m¢
A2
= 2
n
TH
n -1 Kv H
1
1
2/n
-r n
JK W
F m¢ I
d
HA K RS F I p FG r ( n + 1)
I UV 1/ 2
dr
2
=
d
dr
2
n
n -1 TH Kv H
1
1
2/n
-r n
JK W
Here p1, v1 are inlet conditions and remain constant. Also n being index of expansion remains
constant so differentiating and putting equal to zero.
2
×r
(2 - n)
n - FH n + 1 IK r 1/ n
=0
n n
or,
2
×r
(2 - n)
n =FH n + 1 IK r 1/ n
n n
=F
H IK
1- n
n n +1
r
2
or, Critical pressure ratio,
n
2 (n − 1)
r=
n + 1
Let critical pressure at throat be given by pc or pt then,
n n
pc 2 (n − 1) pt 2 (n − 1)
= Þ =
p1 n + 1 p1 n + 1
574 ________________________________________________________ Applied Thermodynamics
Here subscript ‘c’ and ‘t’ refer to critical and throat respectively.
While designing a nozzle the critical pressure ratio at throat is equal to the one obtained above.
Critical pressure ratio value depends only upon expansion index and so shall have constant value. Value
of adiabatic expansion index and critical pressure ratio are tabulated ahead;
Table 13.1: Adiabatic expansion index and critical pressure ratio for selected fluids
F 2 I
n
pc
H n + 1K
( n - 1)
=
p1
Wet steam 1.135 (n = 1.035 + 0.1x, where
x is dryness fraction of wet steam) 0.577
Superheated 1.3 0.545
steam
Air 1.4 0.528
The maximum discharge per unit area can be obtained by substituting critical pressure ratio in
expression for mass flow per unit area at throat section.
LM F n I p R|F 2 I U|OP
-F
2 I
2 ( n + 1)
- 1V P
m¢ ( n - 1) ( n - 1)
L 2n I p F 2 I OP
= MF
( n + 1) 1/ 2
MNH n - 1 K × v H n + 1 K
1 ( n - 1)
× {
n+1
-1 }P
1 2
Q
L p F 2 I OP ( n + 1) 1/ 2
= Mn ×
m¢ ( n - 1)
MN v H n + 1 K PQ
1
At 1
F 2 I
( n + 1)
p
H n +1K
( n - 1)
Maximum discharge per unit area = n× 1
v1
For this maximum discharge per unit area at throat the velocity at throat can be obtained for
critical pressure ratio. This velocity may also be termed as ‘critical velocity’.
n
C2 = 2 ( p1 v1 − p2 v2 )
n − 1
At throat
n p1v1
Ct = 2 pt vt − 1
n − 1 pt vt
Nozzles _______________________________________________________________________ 575
F n I
R|F p I U| (1 - n )
H n - 1 K ST|H p K V| n
C = 2 × p v -t
1
W
t t t
1
FpI
Substituting critical pressure ratio G J
HpK
t
C = 2G
F n IJ × p v {F n + 1 I - 1}
t
H n -1K H 2 K t t
FG g IJ × a p v f
C2 = 2
H g -1K 1 1 - p2 v 2
or,
FG g IJ R aT - T f
C2 = 2
H g -1K 1 2
or,
C2 = a
2 c p T1 - T2 f as cp =
gR
g -1
or,
C2 = a
2 h1 - h2 f
Critical velocity at throat, Ct = g R Tt
Mass flow rate per unit area,
LM F g I p R|F p I 2/g
Fp I
( g + 1)
U|OP 1/ 2
m¢
=
MN2 H g - 1 K × v
1
S|H p K 2
-
HpK
2 g
V|P
A2 1
T 1 1
WQ
Maximum discharge per unit area at throat for critical conditions,
F 2 I
( g + 1)
m¢ p ( g - 1)
At
= g× 1
v1 H g + 1K
576 ________________________________________________________ Applied Thermodynamics
Critical pressure ratio,
F 2 I
g
pc ( g - 1)
p1
=
H g +1K
13.3 CHOKED FLOW
Let us consider a converging nozzle as shown in Fig. 13.10 with arrangement for varying back pres-
sure. A valve is provided at exit of nozzle for regulating the back pressure at section 2-2. Let us denote
back pressure by pb. Expansion occurs in nozzle from pressure p1 to pb.
Initially when back pressure pb is equal to p1 there shall be no flow through the nozzle but as back
pressure pb is reduced the mass flow through nozzle increases. With the reduction in back pressure a
situation comes when pressure ratio equals to critical pressure ratio (back pressure attains critical
pressure value) then mass flow through nozzle is found maximum. Further reduction in back pressure
beyond critical pressure value does not affect the mass flow i.e. mass flow rate does not increase
beyond its’ limiting value at critical pressure ratio. Thus under these situations flow is said to be choked
flow or critical flow.
aT - T f
aT - T f
2¢
hNozzle = 1
1 2
In case of working fluid being steam the friction causes heating of steam flowing through nozzle
thereby increasing dryness fraction. The volume of steam at exit also increases due to increase in
dryness fraction. T–s and h-s representation of steam expanding through nozzle is given in Fig. 13.14.
Fig. 13.14 T-s and h-s representation for steam expanding through nozzle
Nozzles _______________________________________________________________________ 581
Due to friction the velocity at exit from nozzle gets modified by nozzle efficiency as given below.
Velocity at exit, C2 = a f
2 h1 - h2 + C12 , for no friction
In case of nozzle with friction the enthalpy drop, (h1 – h2¢) gives velocity at exit as,
C2¢ = a f
2 h1 - h2 ¢ + C12
-
C22 C12
or, (h1 – h2) =
2
C - C12
2
and (h1 – h2¢) = 2 ¢
2
bC 2
2¢ - C12 g
Substituting in nozzle efficiency, hNozzle =
bC 2
2 - C12 g
For negligible inlet velocity i.e. C1 » 0
C22¢
Nozzle efficiency, hNozzle =
C22
Thus it could be seen that friction loss will be high with higher velocity of fluid. Generally frictional
losses are found to be more in the downstream after throat in convergent-divergent nozzle because of
simple fact that velocity in converging section upto throat is smaller as compared to after throat. Expansion
upto throat may be considered isentropic due to small frictional losses. Apart from velocity considerations
the significantly high frictional loss in diverging portion of nozzle compared to converging portion can
be attributed to the contact surface area. Length of converging section upto throat is quite small compared
to length of diverging portion after throat as it has subsonic acceleration which can be completed in
short length. Diverging section of nozzle is designed comparatively longer than converging section so as
to avoid flow separation due to adverse duct geometry (diverging type). Turbulence losses are also
significant in diverging portion compared to converging portion. Due to the different factors discussed
above the frictional losses are found to be more in diverging portion compared to converging portion.
Normally angle of divergence in divergent portion is kept between 10° and 25° so as to avoid flow
separation. But small divergence angle causes increase in length of diverging portion therefore increas-
ing frictional losses. Thus a compromise should be striked in selecting angle of divergence as very small
angle is desirable from flow separation point of view but undesirable due to long length and larger
frictional losses point of view. Length of diverging portion of nozzle can be empirically obtained as
below
L= 15 ⋅ At
where At is cross-sectional area at throat.
While designing the nozzle parameters due care should be taken for smoothness of nozzle profile,
surface finish for minimum friction and ease of manufacturing etc. Thus finally, it can be concluded that
nozzle efficiency depends upon nozzle material, size and shape of nozzle, angle of divergence, nature of
fluid flowing and its properties etc.
Coefficient of velocity: The ‘coefficient of velocity’ or the ‘velocity coefficient’ can be given by
the ratio of actual velocity at exit and the isentropic velocity at exit. Thus it measures the amount of
deviation from ideal flow conditions. Mathematically,
Cactual at exit
Coefficient of velocity =
Cisentropic at exit
582 ________________________________________________________ Applied Thermodynamics
Coefficient of discharge: The ‘coefficient of discharge’ or ‘discharge coefficient’ is given by the
ratio of actual discharge and the discharge during isentropic flow through nozzle. Mathematically,
mactual
Cofficient of discharge =
misentropic
Here m refers to discharge rate.
Fig. 13.15 Expansion of steam on T-s and h-s diagram under equilibrium
Superheated steam undergoes continuous change in state and becomes dry saturated steam at
state 3 and subsequently wet steam leaving steam turbine at state 2. Some times expansion of steam
occurs in metastable equilibrium or in equilibrium in which change of steam state could not maintain its
pace with expanding steam. This phenomenon in which change of steam state could not occur simulta-
neously with expanding steam in nozzle is called phenomenon of supersaturation and flow is called
supersaturated flow or metastable flow.
EXAMPLES
1. Dry steam at 10 bar and 100 m/s enters a nozzle and leaves it with velocity of 300 m/s at 5 bar.
For 16 kg/s of steam mass flow rate determine heat drop in nozzle and final state of steam leaving nozzle
assuming heat loss to surroundings as 10 kJ/kg.
Solution:
Given,C1 = 100 m/s, C2 = 300 m/s
p1 = 10 bar = 1 MPa, p2 = 0.5 MPa
By steady flow energy equation between sections 1 and 2.
C12 C2
h1 + + q = h2 + 2 + w
2 2
Nozzles _______________________________________________________________________ 585
H 2 K
1 2
(h2 – h1) = – 104 +
+ S
R (100) - (300) UV
2 2
= – 104
T 2 W
h2 – h1 = – 30,000
Heat drop in nozzle = 30 kJ/kg Total heat drop = 480 kJ/s Ans.
or, h2 = 2748.1 kJ/kg
At 5 bar, pressure,
hf = 640.23 kJ/kg, hfg = 2108.5 kJ/kg
Let dryness fraction at exit be x2,
2748.1 = 640.23 + x2 × 2108.5
x2 = 0.99
h2 = h400°C, 6 bar +
ch 500 °C,6 bar - h400 °C , 6 bar h (418.45 – 400)
(500 - 400)
(3482.8 - 3270.3)
= 3270.3 + ´ 18.45
100
h2 = 3309.51 kJ/kg×K
Neglecting inlet velocity, C1 » 0, assuming no heat loss,
a
C2 = 2 h1 - h2 f
C2 = 581.36 m/s
A2 ´ C2
Mass flow rate, =
v2
Specific volume at exit (2)
v2 = v400°C, 6 bar +
cv
500 °C , 6 bar - v 400 °C , 6 bar h ´ (418.45 – 400)
( 500 - 400 )
20 ´ 10 -4 ´ 581.36
Mass flow rate = = 2.2017
0.5281
Mass flow rate = 2.202 kg/s Ans.
3. In a nozzle steam expands from 12 bar and 300°C to 6 bar with flow rate of 5 kg/s. Determine
throat and exit area if exit velocity is 500 m/s and velocity at inlet to nozzle is negligible. Also find
coefficient of velocity at exit. Coefficient of velocity is the ratio of actual velocity of fluid at nozzle exit
to the velocity at exit considering isentropic flow through nozzle.
Solution:
FG p IJ = F 2 I n
( n - 1)
H p K H n + 1K
2
1
For superheated steam, n = 1.3
FH p IK = FG 2 IJ 1. 3 / 0 .3
H 1 + 1.3 K
2
= 0.545
12
Throat pressure, p2 = 6.54 bar
From steam table;
At 6.54 bar, h2 = 2900.05 kJ/kg
T 2 = 224.48°C, v2 = 0.3466 m3/kg
Velocity at throat C2 = a
2 h1 - h2 f
= 2 ´ 10 3 (3045.8 - 2900.05)
C2 = 539.9 m/s
From continuity equation, m1 = m2 = 5 kg/s
A2 C2
m2 =
v2
588 ________________________________________________________ Applied Thermodynamics
A2 ´ 539.9
5=
0.3466
A2 = 3.209 ´ 10–3 m2
a
Velocity at exit of nozzle, C3 = 2 h1 - h3 f
From steam tables, Temperature at exit of nozzle, T3 = 215.15°C
Enthalpy at exit of nozzle, h3 = 2882.55 kJ/kg
Specific volume at exit of nozzle, v3 = 0.3647 m3/kg
Ideal velocity at exit, C3 = 2 ´ 10 3 (3045.8 - 2882.55)
C3 = 571.40 m/s
Actual velocity at exit, C3, actual = 500 m/s
A3 ´ C3, actual
Area at exit, m1 = m2 = m3 = 5 kg/s =
v3
5 ´ 0.3647
A3 = = 3.647 ´ 10–3 m2
500
Cross-sectional area at exit = 3.647 ´ 10–3 m2 Ans.
C3, actual 500
Coefficient of velocity = = = 0.875
C3 571.40
Enthalpy at throat and exit of nozzle can be looked from steam table as,
h2 = 2891.39 kJ/kg
h3 = 2777.0 kJ/kg
v2 = 0.2559 m3/kg
v3 = 0.3882 m3/kg
Nozzles _______________________________________________________________________ 589
Fig. 13.22
Heat drop up to throat section = h1 – h2
Dq12 = 143.41 kJ/kg
Velocity at throat, a f
C2 = 2 h1 - h2 = 535.56 m/s
Heat drop from throat to exit = h2 – h3
Dq23 = 114.39 kJ/kg
Velocity at exit, C3 = a
2 h2 - h3 + C22f
C3 = 718.06 m/s
m2 ´ v 2 1 ´ 0.2559
Throat area, A2 = = = 4.78 ´ 10–4 m2
C2 535.56
m3 ´ v3 1 ´ 0.3882
Exit area, A3 = = = 5.41 ´ 10–4 m2
C3 718.06
For frictionless expansion, Ans.
Throat area = 4.78 cm2
Exit area = 5.41 cm2
(ii) Considering expansion to have 10% friction loss:
Due to frictional heat loss actual heat drop shall get reduced.
Actual heat drop upto throat, Dq¢12 = 0.9 ´ q12
Dq¢12 = 129.07 kJ/kg
Thus, enthalpy at throat, h2, actual = 2905.73 kJ/kg
from steam table, v2, actual = 0.2598 m3/kg corresponding to the modified state as indicated by
enthalpy, h2, actual. This new state at throat shall have temperature of 231.35°C, which can be calculated
from known enthalpy h2, actual at respective pressure.
Actual velocity at throat, C2, actual = 2 ´ 129.07 ´ 10 3
C2, actual = 508.1 m/s
1 ´ 0.2598
Actual throat area, A2, actual =
508.1
Fig. 13.23
Considering negligible inlet velocity, the velocity at throat,
C2 = 2 ´ 10 3 ´ 142
C2 = 532.9 m/s
Nozzles _______________________________________________________________________ 591
A2 × C2
Mass flow rate through nozzle =
v2
=
b g
p ´ 10 -2
2
´ 532.9
4 ´ 0.20
= 0.209 kg/s
The number of nozzles can be determined by using the capacity of turbine, steam consumption
and mass flow rate through a nozzle.
1 ´ 10 3 ´ 8
Number of nozzles = = 10.63 » 11 nozzles
3600 ´ 0.209
Fig. 13.24
Considering velocity at inlet to be negligible, velocity at throat,
C2 = 2 ´ 138 ´ 10 3
C2 = 525.35 m/s
p
4
b g
´ 0.5 ´ 10 - 2 ´ 525.35
2
Fig. 13.25
At inlet to nozzle, from steam table, h1 = hg at 16 bar = 2794 kJ/kg
s1 = sg at 16 bar = 6.4218 kJ/kg×K
At exit of nozzle,
s1 = s2 = 6.4218 kJ/kg×K
Since s2 < sg at 8 bar so the state at the end of expansion lies in wet region. Let dryness fraction at
state 2 be x2, then,
s2 = sf at 8 bar + x2 × sfg at 8 bar
6.4218 = 2.0462 + x2 × 4.6166
or, x2 = 0.9478
h2 = hf at 8 bar + x2 × hfg at 8 bar
= 721.11 + (0.9478 ´ 2048)
h2 = 2662.2 kJ/kg
Theoretically, change in enthalpy, Dh12 = 2794 – 2662.2
Dh12 = 131.8 kJ/kg
D h12 actual
= 0.90 or Dh12 actual = 0.9 ´ 131.8 = 118.62 kJ/kg
D h 12
h2 actual = 2794 – 118.62 = 2675.38 kJ/kg
h2, actual = 2675.38 = hf at 8 bar + x2 actual × hfg at 8 bar
594 ________________________________________________________ Applied Thermodynamics
2675.38 = 721.11 + (x2 actual ´ 2048)
x2 actual = 0.9542
Specific volume, v2 actual = vf at 8 bar + x2 actual ´ vfg at 8 bar
v2 actual = 0.001115 + (0.9542 ´ (0.2404 – 0.001115))
v2 actual = 0.2294 m3/kg
Neglecting the velocity at inlet to nozzle,
C2 = 2 ´ D h12 actual
C2 = 487.1 m/s
By continuity equation,
A2 C2
m=
v2 actual
5000 ´ 6 A ´ 487.1
= 2
3600 0.2294
A2 = 3.925 ´ 10–3 m2
Area at exit of nozzle » 39.25 cm2
60 ´ p
Approximate length of nozzle arc, = = 62.83 cm
3
Length of nozzle arc
Number of nozzles =
Pitch
62.83
= » 13 nozzles.
5
Correct length of nozzle arc = 13 ´ 5 = 65 cm
Flow area at exit of each nozzle = (5 sin 12 – 0.3) ´ h
39.25 = (5 sin 12 – 0.3) ´ h ´ 13
or, h= 4.08 cm
Length of nozzle = 65 cm
Ans.
Radial height of nozzle = 4.08 cm
8. Air is expanded reversibly and adiabatically in a nozzle from 13 bar and 150°C to a pressure of
6 bar. The inlet velocity of the nozzle is very small and the process occurs under steady flow conditions.
Calculate the exit velocity of the nozzle. [U.P.S.C., 1992]
Solution:
Let inlet and exit states in nozzle be given as 1 and 2 respectively. Assuming no heat interaction and
work interaction, no change in potential energy during flow across the nozzle. Steady flow energy
equation gets modified as under;
C12 C2
h1 + = h2 + 2
2 2
Neglecting velocity at inlet to nozzle C1.
C2 = a
2 h1 - h2 f
Nozzles _______________________________________________________________________ 595
C2 = a
2 c p × T1 - T2 f
Assuming adiabatic expansion in nozzle, between states 1 and 2.
FPI
( g - 1)
T2 g
HPK
2
=
T1 1
FPI
( g - 1)
g
=T ×
HPK
2
or, T2 1
1
(1. 4 - 1)
= (273 + 150)
6
13 e j 1. 4
T2 = 339.16 K
Substituting in expression for exit velocity, C2,
C2 = 2 ´ 1.005 ´ ( 423 - 339.16 )
C2 = 12.98 m/s
Exit Velocity = 12.98 m/s Ans.
9. During a test on steam nozzle steam impinges a stationary flat plate which is perpendicular to
the direction of flow and the force on the plate is measured. The force is found to be 350 N when dry
saturated steam at 8 bar is expanded to 1 bar. Throat cross-section area is 5 cm2 and exit area is such that
the complete expansion is achieved under these conditions.
Determine
(i) the discharge at throat.
(ii) the efficiency of nozzle assuming that all the losses occur after throat and that n = 1.13
for isentropic expansion.
Solution:
For isentropic expansion up to throat.
F 2 I
n
P2 ( n - 1)
P1
=
H n + 1K
8F
2 I
1.13
H 113
. +1K
0 .13
P2 =
P 2 = 4.62 bar
From steam table,
Enthalpy at 1, h1 = hg at 8 bar = 2769.1 kJ/kg
s1 = sg at 8 bar = 6.6628 kJ/kg.K
Fig. 13.26
596 ________________________________________________________ Applied Thermodynamics
At throat section let the dryness fraction be x2, then
s1 = s2 = 6.6628 = sf at 4.62 bar + sfg at 4.62 bar × x2
6.6628 = 1.7872 + x2 × 5.0178
x2 = 0.9717 or 97.17%
h2 = hf at 4.62 bar + x2 × hfg at 4.62 bar
= 627.33 + (0.9717 ´ 2117.77)
h2 = 2685.17 kJ/kg
Specific volume at throat, v2 = vf at 4.62 bar + x2 × vfg at 4.62 bar
= 0.001089 + (0.9717 ´ 0.4035)
v2 = 0.3932 m3/kg
Enthalpy change up to throat Dh12 = h1 – h2
Dh12 = 83.93 kJ/kg
Velocity at throat, C2 = 2 ´ D h12 = 2 ´ 83.93 ´ 10 3 = 409.71 m/s
A2 × C2
Discharge at throat, m2 =
v2
5 ´ 10 -4 ´ 409.71
=
0.3932
m2 = 0.521 kg/s
Fig. 13.27
Enthalpy drop up to throat,
F n IJ P v
( n − 1)
– h¢ = G
1 − 2P n
h1 2 H n -1K 1 1 P1
FH 1.3 IK (10 ´ 10 LM F 5.45 I 0 .3 / 1.3
OP
N H 10 K
´ 0.2060) 1 -
Q
= 2
0 .3
h1 – h2¢ = 116.67 kJ/kg
h2¢ = 2711.23 kJ/kg
Let dryness fraction be x2 at throat, so using steam tables it can be obtained as under,
h2¢ = hf at 5.45 bar + (x2 ´ hfg at 5.45 bar)
598 ________________________________________________________ Applied Thermodynamics
2711.23 = 663.36 + (x2 ´ 2098.15)
x2 = 0.9760
s2¢ = sf at 5.45 bar + (x2 × sfg at 5.45 bar)
= (1.8936 + (0.9760 ´ 4.8989) = 6.6749 kJ/kg×K
Let dryness fraction at exit be x3, so
s3 = s2¢ = 6.6749 kJ/kg×K
s3 = sf at 1 bar + x3 × sfg at 1 bar
6.6749 = 1.3026 + x3 ´ 6.0568
x3 = 0.8869
Enthalpy at exit, h3 = hf at 1 bar + x3 × hfg at 1 bar
= 417.46 + (0.8869 ´ 2258)
h3 = 2420.08 kJ/kg
Specific volume at exit, v 3 = vf at 1 bar + x3 × vfg at 1 ba
= (0.001043 + 0.8869 ´ 1.6929)
v3 = 1.5025 m3/kg
Velocity at exit, C3 = 2 × D h13 = 903.13 m/s
A3 × C3
Mass flow rate, m3 =
v3
5 A ´ 903.13
= 3
60 1.5025
A3 = 1.386 ´ 10–4 m2 or 1.386 cm2
Diameter of nozzle at exit d 3 = 0.0133 m or 1.33 cm
FP I
( n - 1)
T2
HPK
n
2
=
T1 1
(1. 3 - 1)
T2 = 473 ( 0.545) 1.3 = 411.18 K = 138.18°C
Saturation pressure corresponding to 411.18 K,
Psat, 138.18°C = 3.44 bar
P2
Degree of supersaturation =
Psat ,138.18°C
5.45
= = 1.58
3.44
Amount of undercooling = (Tsat at 5.45 bar – T2)
= (155.12 – 138.18)
= 16.94°C
Nozzles _______________________________________________________________________ 599
Fig. 13.28
Using governing equation given,
P1 v11.3 = P2 v21.3
4 ´ 102 ´ (0.5088)1.3 = 1.5 ´ 102 ´ v21.3
v2 = 1.082 m3/kg
Enthalpy at 2,
P2 v2 ´ 102 = (h2 – 2614)
(1.5 ´ 1.082 ´ 102) = (h2 – 2614)
or, h2 = 2776.3 kJ/kg
Actual heat drop = hN ´ (h1 – h2)
= 0.95 ´ (2817.52 – 2776.3)
= 39.16 kJ/kg
F 2 I
n
P2 ¢ ( n - 1)
P1
=
H n + 1K
F IK
F 2 IH
1. 3
P ¢ = 14 ´ G
H 1.3 + 1 JK
0 .3
2
= 2 ´ 10 3 ´ 88.92
= 421.71 m/s
m × v2
Cross-sectional area at exit of nozzle A 2 =
C2 × N
5 ´ 0.2827
= ´ 104
421.71 ´ 16
= 2.13 cm2
= b 2 ´ 10
´ 88.92 + (100 ) 2
3
g
C2¢ = 433.41 m/s
A2 ´ C ¢2 ´ 16
Discharge with modified velocity =
v2
16 ´ 2.13 ´ 10 -4 ´ 433.41
=
0.2827
= 5.22 kg/s
602 ________________________________________________________ Applied Thermodynamics
Fig. 13.29
But saturation temperature corresponding to throat pressure,
Tsat 11.6 bar = 186.43°C
Saturation pressure corresponding to 428.14K (155.14°C)
Psat at 155.14°C = 0.5452 MPa
Psat at 155.14°C = 5.452 bar
11.6
Degree of supersaturation =
5.452
= 2.13
Degree of undercooling = Tsat at 11.6 bar – T2
= 186.43 – 155.14 = 31.29°C
Nozzles _______________________________________________________________________ 603
Degree of supersaturation = 2.13
Ans.
Degree of undercooling = 31.29°C
Enthalpy at state 1, h1 = hf at 20 bar = 2799.5 kJ/kg
v 1 = vg at 20 bar = 0.009963 m3/kg
s1 = sg at 20 bar = 6.3409 kJ/kg×K
Isentropic enthalpy drop,
(n − 1)
n P2 n
h 1 – h2 = ⋅ Pv
1 1 1 −
n −1 1
P
1.3 LM F IK FH
1.3 - 1 IK O
P
MN H
11.6
´ 20 ´ 10 ´ 0.009963 1 -
PQ
2 1. 3
h 1 – h2 =
1.3 − 1 20
88.55
= ´ 100
(2799.5 - 2544.21)
= 34.68%
Fig. 13.30
Pressure at inlet of nozzle, P1 = 20 bar,
Temperature T1 = Tsat at 20 bar = 212.42°C
Pressure at throat, P2 = 0.7 ´ 20 = 14 bar
Enthalpy at state 1, h1 = hf at 20 bar + x1 × hfg at 20 bar
= (908.79) + (0.95 ´ 1890.7)
h1 = 2704.95 kJ/kg
s1 = sf at 20 bar + x1 × sfg at 20 bar
= 2.4474 + (0.95 ´ 3.8935)
s1 = 6.1462 kJ/kg×K
At throat s2= s1 = 6.1462 kJ/kg×K
Let dryness fraction at throat be x2
s2 = sf at 14 bar + x2 × sfg at 14 bar
6.1462 = 2.2842 + (x2 ´ 4.1850)
x2 = 0.923
h2 = hf at 14 bar + x2 × hfg at 14 bar
= 830.30 + (0.923 ´ 1959.7)
h2 = 2639.10 kJ/kg
Specific volume v2 = vf + x2 × vfg at 14 bar = 0.13 m3/kg
a f
at 14 bar
Velocity of steam at throat, C2 = 2 h1 - h2
= 2 ´ (2704.95 - 2639.10 ) ´ 10 3
606 ________________________________________________________ Applied Thermodynamics
C2 = 362.9 m/s
Applying Bernoulli’s equation between states 3 and 4,
P3 C32 P C2
gZ3 + + = gZ4 + 4 + 4
r 2 r 2
C32 1.013 ´ 10 5 20 ´ 10 5 ( 20 ) 2
+ = + (0.8 ´ 9.81) +
2 10 3 10 3 2
C3 = 64.91 m/s
Let mass of water lifted by per kg of steam be m kg, i.e. m kg of water per kg of steam.
C2 – m 2 gh2 w = (m + 1) ´ C3
C2 – m ´ 2 ´ 9.81 ´ 5 = (m + 1) ´ C3
362.9 – 9.90 m = (m + 1) ´ 64.91
m = 3.98 kg of water per kg steam
150 1
Mass of steam required for given water flow rate = ´
60 3.98
msteam = 0.628 kg/s
Asteam × Csteam
msteam =
vsteam
A2 × C2
msteam =
v2
A2 ´ 362.9
0.628 =
0.13
A2 = 2.25 cm2
d2 = 1.69 cm
F 2 I
n
( n - 1)
Pressure at throat, P2 = P1
H n + 1K , here n = 1.13, P2 = 0.87 bar
Fig. 13.31
Temperature at state 1, T1 = Tsat at 1.5 bar = 111.37°C = 384.37 K
Enthalpy at state 1, h1 = hf at 1.5 bar + x1 × hfg at 1.5 bar
= 467.11 + (0.9 ´ 2226.5)
h1 = 2470.96 kJ/kg
s1 = 1.4336 + (0.9 ´ 5.7897)
s1 = 6.6443 kJ/kg×K
At state 2, s2 = s1 = 6.6443 = sf at 0.87 bar + x2 × sfg at 0.87 bar
6.6443 = 1.2560 + x2 × 6.1538
x2 = 0.88
h2 = hf at 0.87 bar + x2 ´ hfg at 0.87 bar
= 400.26 + (0.88 ´ 2268.71)
h2 = 2396.72 kJ/kg
Specific volume, v 2 = vf at 0.87 bar + x2 × vfg at 0.87 bar
608 ________________________________________________________ Applied Thermodynamics
= 0.001039 + (0.88 ´ 1.9649), v2 = 1.7302 m3/kg
Steam velocity, C2 = a
2 h1 - h2 f
C2 = 385.33 m/s
By Bernoulli’s equation,
C32 P
= 4
2 r
1.2 ´ 20 ´ 2 ´ 10 5
C3 =
1000
C3 = 69.28 m/s
Let mass of water entrained per kg of steam be m kg water per kg of steam.
Neglecting momentum of entering water the momentum balance can be given as under.
1 ´ C2 = (m + 1) ´ C3
385.33 = (m + 1) ´ 69.28
m = 4.56 kg water per kg of steam
5000
Mass of steam supplied per second = = 0.3046 kg/s
3600 ´ 4.56
m v2
Area of steam nozzle, A2 =
C2
0.3046 ´ 1.7302
=
385.33
A2 = 13.68 cm2
-:-4+15-5