Hydraulic Machines Textbook
Hydraulic Machines Textbook
Hydraulic Machines Textbook
Historical Review . . . . . . . . . . . . . 2
1.1 General Introduction . . . . . . . . . . . . 4
1.2 Velocity Diagram . . . . . . . . . . . . . . 5
1.3 Momentum Transfer Principles . . . . . . . . 6
1.4 Energy Equation . . . . . . . . . . . . . . 9
1.5 Theories of Turbomachines . . . . . . . . . . 11
Euler Theory (Elementary) 11
Modern Theory 14
Necessity for flow unsteadiness 18
1.5.4 Approximate calculation of deviation 19
after Stodola
1.6 Some Practical Considerations (Actual
Machine Design) . . . . . . . . . . . . . . 20
Friction 20
Disk Friction 20
Leakage 21
Pre-rotation of the fluid 21
1.7 Coefficients and Efficiencies . . . . . . . . . 22
Circulatory Flow Coefficient 22
Manometric Efficiency 22
Mechanical Efficiency 23
Volumetric Efficiency 23
Hydraulic Efficiency (Turbine) 23
2.1 Introduction . . . . . . . . . . . . . . . . 24
2.2 Dimensional Analysis . . . . . . . . . . . . 24
2.3 Hydraulic Similarity . . . . . . . . . . . . . 25
2.4 Application of Dimensional Analysis on
A
Turbomachines . . . . . . . . . . . . . . . 26
2.4.1 Discussion 27
2.4.2 Performance Curves 28
2.5 Scale Effect . . . . . . . . . . . . . . . . 29
Reynolds Number effect 29
Scale effects in Hydraulic Machines 30
Scale effects in compressible machines 35
2.6 Affinity Laws . . . . . . . . . . . . . . . . 39
2.7 Specific Speed . . . . . . . . . . . . . . . 40
Pumps 41
Compressors and Blowers 43
Hydraulic Turbines 43
2.8 Pressure and Flow Coefficients . . . . . . . . 44
2.9 Specific Diameter . . . . . . . . . . . . . . 44
3.1 Introduction . . . . . . . . . . . . . . . . 47
3.2 Cascade Nomenclature . . . . . . . . . . . . 47
3.3 Analysis of Cascade Forces . . . . . . . . . . 49
3.4 Lift and Drag . . . . . . . . . . . . . . . . 51
3.5 Cascades in Motion . . . . . . . . . . . . . 54
3.6 Cascade Performance . . . . . . . . . . . . 55
General Approach 55
Fluid Deviation 57
Off-Design Performance 59
Turbine Cascade Performance 59
3.7 Mach Number Effect . . . . . . . . . . . . 60
3.8 Ideal Characteristics . . . . . . . . . . . . 62
Zero Lift Angle 62
Impulse Flow Angle 63
3.9 The Head-Capacity Curve of Straight Cascade . 64
3.10 Radial Cascade . . . . . . . . . . . . . . 65
3.11 Cascade Characteristics Analysis . . . . . . . 67
3.12 Singularity Method . . . . . . . . . . . . 67
Method of Solution for Single Airfoil 69
Conformal Transformation Method 71
B
CHAPTER (IV) INCOMPRESSIBLE FLOW
TURBOMACHINES ( PUMPS )
4.1 Introduction . . . . . . . . . . . . . . . . . 78
4.2 Centrifugal Pumps (Radial) . . . . . . . . . . 78
4.2.1 General Considerations 78
a. Volute type pump 79
b. Diffuser type pump 79
4.2.2 Effect of Impeller Exit Angle β2 80
4.2.3 Efficiencies and Coefficients of Centrifugal Pumps 82
i. Efficiencies 82
ii. Coefficients 83
iii. Affinity Laws 83
iv. Specific Speed 83
4.2.4 Centrifugal Pump Actual Performance 83
4.2.4.1 Actual Head Capacity Curve 83
4.2.4.2 Brake Horsepower and Efficiency Curves 85
4.2.4.3 Analysis of Characteristic Curves 86
4.2.4.4 Influence of Physical
Properties on 87
Performance 87
i. Viscosity Effect 88
ii. Density 88
4.2.5 Some Design Features of Centrifugal Pumps 88
4.2.5.1 Leakage Calculation 89
4.2.5.2 Disk Friction 89
4.2.5.3 Diffuser Losses 90
4.2.5.4 Mechanical Seals 90
a. Single Seals 93
b. Tandem Seals 95
c. Double Seals 96
4.2.5.5 Bearing Losses 96
4.2.5.6 Axial Thrust 99
4.2.5.7 Impeller Design 100
a. Impeller Inlet Dimensions and 101
Angles 105
b. Impeller Exit Dimensions and 105
Angles 106
C
4.2.6 Centrifugal Pump Types 106
4.2.6.1 Fire Pump 106
4.2.6.2 Dredge Pumps 107
4.2.6.3 Slurry Pumps 107
4.2.6.4 Deep Well Pumps 108
4.2.6.5 Circulating Pumps
4.2.6.6 Boiler Feed Pumps
4.2.6.7 Pumping Liquid/Gas Mixtures
4.3 Axial Pumps (Propeller Pumps) . . . . . . . . 109
Degree of Reaction 110
Pressure and Flow Coefficients 110
Study of Flow Inside the Rotor (Radial
Equilibrium) 112
Performance of Axial Flow Propeller Pumps 113
p Selection and Applications . . . . . . . 116
Pumps in Parallel 117
Pumps in Series 117
Economic Considerations 118
Design of the Intake Chamber of Vertical Pumps 118
4.4.4.1 General 118
4.4.4.2 Open Intake Chambers 119
4.4.4.3 Covered Intake Chambers 122
4.4.4.4 Inlet Elbows 123
Pressure Surges (Water Hammer) in Piping
Systems 124
Pump Installation 126
Centrifugal Pump Trouble Shooting 134
D
5.2.2.2 Power,
Efficiency and 153
Coefficients 156
5.2.2.3 Head 157
Delivered by 157
Turbine and 160
Draft Tube 161
5.2.2.4 Types of 162
Draft Tube 163
5.2.2.5 Net Head 163
5.2.2.6 Cavitation in
Turbines 164
5.2.2.7 Power and
Speed
Regulation
5.2.2.8 Francis
Turbine
Performance
5.2.3 Axial Flow Reaction Turbines
a. Propeller Turbine
b. Kaplan Turbine
5.2.4 Some Design Characteristics for Hydraulic
Turbines
5.3 Some Turbines Installations . . . . . . . . . 165
a. Impulse Turbine 165
b. Francis Turbine 166
c. Axial Turbine 173
5.4 Fluid Coupling and Torque Converters . . . . 173
5.4.1 Fluid Coupling 174
5.4.2 Torque Converter 176
5.5 Pump-Turbine, Power Storage System . . . . 178
E
6.4 Entropy . . . . . . . . . . . . . . . . . 191
6.5 Work . . . . . . . . . . . . . . . . . . 191
6.5.1 For a constant volume process 192
6.5.2 For a constant pressure process 192
6.5.3 For a constant temperature 193
process 193
6.5.4 For an adiabatic process 193
6.5.5 For polytropic process
6.6 First Law of Thermodynamics . . . . . . . . 194
6.7 Second Law of Thermodynamics . . . . . . . 194
6.8 Compression of Gases . . . . . . . . . . . 195
6.8.1 Adiabatic Compression 195
6.8.2 Isothermal Compression 196
6.8.3 Polytropic Compression 197
6.9 Plane Compressible Flow . . . . . . . . . . 198
6.10 Gothert's Rule . . . . . . . . . . . . . . 199
6.11 Prandtl-Glauert Rule . . . . . . . . . . . 201
F
CHAPTER (VIII) VOLUMETRIC MACHINES
1 General . . . . . . . . . . . . . . . . . . 231
1.1 Calculation of Loss Coefficient 233
2 Diffuser Types . . . . . . . . . . . . . . . . 236
2.1 Vaneless Diffuser 237
2.2 Vaned Diffuser 239
2.3 Volute Type Diffuser 241
2.3.1 Parallel Walls 241
2.3.2 Tapering Side Walls 242
2.3.3 Rectangular Cross Section 244
References “Appendix I” . . . . . . . . . . . 245
1 Introduction . . . . . . . . . . . . . . . . 246
2 Inception of Cavitation . . . . . . . . . . . . 248
3 Signs of Cavitation . . . . . . . . . . . . . 250
G
3.1 Noise and Vibration 250
3.2 Drop in Head-Capacity and Efficiency 250
Curves 253
3.3 Impeller Vane Pitting and Erosion
4 Mechanisms of Damage . . . . . . . . . . . 253
5 Thermodynamic Effects on Pump Cavitation . . 257
6 Net Positive Suction Head . . . . . . . . . . 260
7 Net Positive Suction Head Test . . . . . . . . 262
8 Thoma’s Cavitation Constant . . . . . . . . 263
9 Suction Specific Speed . . . . . . . . . . . . 264
10 Some Discussions Concerning the NPSH . . . 266
11 Cavitation Noise in Centrifugal Pumps . . . . 268
12 Cavitation Detection by Digital Acoustic
Emission Analysis . . . . . . . . . . . . . 286
13 How to Prevent Cavitation . . . . . . . . . 288
References “Appendix II” . . . . . . . . . 290
Chapter I . . . . . . . . . . . . . . . . . 292
Solved Examples 294
Problems 301
Chapter II . . . . . . . . . . . . . . . . . 303
Solved Examples 305
Problems 311
Chapter III . . . . . . . . . . . . . . . . 312
Solved Examples 314
Problems 319
Chapter IV . . . . . . . . . . . . . . . . 320
Solved Examples 324
Problems 333
Chapter V . . . . . . . . . . . . . . . . 336
Solved Examples 338
Problems 347
Chapter VII . . . . . . . . . . . . . . . . 349
Solved Examples 350
Problems 353
Chapter VIII . . . . . . . . . . . . . . . 354
H
Solved Examples 356
Problems 359
Mech.MuslimEngineer.Net
Types and shapes of turbomachines (adopted from Sayers)
BASIC THEORY
HISTORICAL REVIEW
Mech.MuslimEngineer.Net
Fig. A.1 Hero’s rotating sphere Fig. A.2 Giovanni de Branca's
of 120 B.C. turbine of 1629
Throughout this text, the emphasize has put on the practical aspects
of the machines without going deep inside the mathematical formulation.
Some important applications are treated separately as; cavitation
Mech.MuslimEngineer.Net
phenomena, pressure recovery devices and maintenance of
turbomachines.
The rotors also can be classified by the direction of flow in the wheel:
- Radial Wheel,
- Axial Wheel,
- Mixed Wheel.
Energy Conversion
Impulse Reaction
Principle
Energy transfer direction + ve + ve - ve + ve - ve
Flow Direction Radial Radial Axial Axial
Turbomachine Pelton Francis Centrifugal Kaplan Propeller
Wheel Turbine Pump Turbine Pump
→ → →
C =U+W
U =ω D/2
D : wheel diameter
→ →
α = angle between U, C
→ →
β = angle between U , W
Cu = C cos α
C r = C sin α
Fig. 1.2 Velocity triangle
→
Cu : tangential component of C
→
C r : radial component of C
Mech.MuslimEngineer.Net
At time t + d t the momentum in a unit volume could generally be written
as follows:
∂ (ρ C )
ρC + dt .................................………........................ (1.2)
∂t
∂ (ρ C )
J + d J = ∫∫∫ ρ C∂ V + ∫∫∫ dt. ∂ V
V V ∂t
∂ (ρ C )
+ ∫∫∫ ′′ ρ C∂ V ′′ + ∫∫∫ ′′ dt.∂ V ′′
δV δV ∂t
∂ (ρ C )
− ∫∫∫ ′ ρ C∂ V ′ − ∫∫∫ ′ dt.∂ V ′ ..........…….................. (1.3)
δV δV ∂t
δ V = d S.Cn . d t
∂ (ρ C )
i.e. d J = ∫∫∫ .d t.dV + ∫∫ C.ρ .d S .C n .d t
V ∂t Σ
∂ (ρ C )
= dt ∫∫∫ .dV + ∫∫ C.ρ .C n .d S …………............. (1.5)
dt Σ
V
The above equation contains two terms, the first refers to the local change
in the control volume and the second refers to the connective change
through the boundary of the control volume. The mass flow rate element
is given by:
∂J ∂ (ρ C )
= ∫∫∫ d V .......................................................… (1.7)
∂t V ∂t
dJ ∂ J
d t ∂ t ∫∫Σ
= + C dQm …...............................................…..... (1.8)
Mech.MuslimEngineer.Net
d J .r ∂ J .r
∂ t ∫∫Σ
T = F .r = = + C.r.dQm ................................. (1.9)
dt
Taking the time average of the above equation i.e. ∂ J.r/∂ t = 0, the
mean value of the local change of momentum vanishes and we have only
the average value
d Q = d E + dW …………................................................. (1.11)
where d E is the change in energy per unit mass, this term includes
internal, kinetic and potential energy.
(
d E = d E n + C 2 / 2 + gZ ) ...................................……........ (1.12)
∂ (E n + C 2 / 2 + gZ )
dE = dt ∫∫∫ ρ .dV +
V ∂t .....................…..... (1.13)
dt ∫∫ (E n + C 2 / 2 + gZ ).ρ .C n .dS
S
In the case of an open boundary, the term ∫∫S τ .Cτ .dS denotes a
friction energy E f which is sometimes can be neglected compared to the
other stresses involved, despite in the case of a closed but moving
boundary as a turbine rotor, this tangential stress delivers power E s .
One can combine the convective part of energy with the hydrostatic
pressure to give;
∫∫S [P.C n + (E n + C ) ]
/ 2 + gZ ρ .C n d S = ∫∫ io .d Qm
2
......... (1.17)
S
10
Mech.MuslimEngineer.Net
dE
+ E s + ∫∫ io .d Qm = 0 ................……........……...…...... (1.19)
dt S
In many cases the change in energy E with time is constant, thus we can
write the energy equation as follows:
E s = − ∫∫ io .d Qm .................................….....................….... (1.20)
S
The above equations are for steady flow, the dimensions of E s is length
per unit mass. From equations (1.14) and (1.20):
11
where T torque
Qm mass flow rate = W / g
W weight flow rate = γ Q with Q the volume flow rate
we have:
C u1 = C1 cos α 1 , Cu 2 = C 2 cos α 2 ,
U = ω .R , Power = T .ω
12
Mech.MuslimEngineer.Net
= (U 1 .Cu1 − U 2 .Cu 2 ) g per unit weight ………....... (1.24)
Remarks:
- From equation (1.22) we can see that the maximum power generated
by the turbine occurs when cos α 2 = 0 , i.e. when α 2 = 90° , radial
exit.
- Also, from equation (1.22) we see that Euler theory does not take into
consideration the geometric form of the blades, it considers the inlet
and outlet angles only.
- Applying the triangle relations to equation (1.24) we find another form
of equation (1.24):
W 2 = U 2 + C 2 − 2 U C. cos α
or:
Ho =
(C12 − C 22 ) + (W22 − W12 ) + (U 12 + U 22 )
..……......…....... (1.25)
2g
- It may be noted that equations (1.24) and (1.25) do not contain the
specific weight of the fluid handled, an impeller operating at a given
speed will develop the same head for any fluid handled.
The change of head due to the fluid flow through the wheel can be
written as following:
P * C2 P* C2
H o = H1 − H 2 = + − ρ g + 2 g ...................... (1.26)
ρ g 2 g 1 2
where
P * = P + ρ g h …....................................………................. (1.27)
13
The degree of reaction σ is the ratio of the change in static pressure to the
change in total pressure.
P1* − P2*
σ= Ho
ρg
............................... (1.28)
=
(U 1
2
) (
− U 22 + W22 − W12 )
(
C12 − C 22 ) + (U 2
1 ) (
− U 22 + W22 − W12 )
σ = 0.0 designate an impulse machine.
σ ≠ 0.0 designate a reaction machine.
To explain the actual flow inside the impeller we may consider the
case of a centrifugal pump impeller as follows:
14
Mech.MuslimEngineer.Net
Fig. 1.6 Pressure distribution Fig. 1.7 Exit flow shape in
between two blades radial blades
a. During the energy transfer process between the blades and the fluid, a
positive pressure acting on the driving face of blade and negative
pressure on the trilling face, (See Figures 1.6 and 1.7). If one consider
the total energy is constant, the increment of the static pressure energy
from the trailing edge to the driving edge is made on the dispense of
the dynamic pressure energy. The relative velocity magnitude, which is
a form of dynamic pressure energy, is inversely proportional to the
pressure difference. The relative velocity exit direction is tangential to
the high pressure side of the blade, and inclined opposite to the rotating
speed in the low pressure side of the blade. This fact has been observed
early in the 60’s by Dean et al., after a measurement of flow pattern at
the impeller exit, this flow pattern is known as jet-wake flow pattern.
15
16
Mech.MuslimEngineer.Net
Fig. 1.9 Inertial effect Fig. 1.10 The flow pattern inside
the impeller
From a, b we can see that for actual machine the relative velocity
mean exit angle does not equal to the ideal exit angle suggested by Euler
theory, usually this deviation is about 5 - 10°.
In practice the designer used empirical correction factor based on
test results and experiments, Figure 1.11.
The energy transfer process between fluid and rotor will occur only
if E s not equals zero.
17
∂C ∂C 1 ∂ P ∂ ( gZ )
+C =− − ......................................... (1.30)
∂t ∂l ρ ∂l ∂l
or
∂C ∂ C2 P
=− + gZ +
∂t ∂ l 2 ρ
d P
Isentropic process di = , and after equation (1.18) the following
ρ
relationship is reached:
∂C ∂ C2 ∂i
=− + gZ + i = − o ....................................... (1.31)
∂t ∂ l 2 ∂l
D io ∂ io ∂ C ∂i
= −C = .........…........................................ (1.32)
Dt ∂t ∂ t ∂t
The above equation indicates clearly that i or H must not equal zero
to extract or to supply energy from a turbomachine. In other words, the
flow must be unsteady to permit energy transfer. The physical meaning of
equation (1.32) is experimentally well presented by the jet-wake flow
pattern; a pressure measurement impeller's exit at a fixed point will show
the nature of the flow unsteadiness.
18
Mech.MuslimEngineer.Net
of the impeller ω, and finally Stodola assumes that the deviation in the
tangential component of the absolute velocity ∆C u is proportional to the
channel width and the rotating speed.
a
∆Cu ∝ ω ............................................…........................... (1.33)
2
The space between two successive blades at the outer radius is given by:
t = π D2 / Z n
19
20
Mech.MuslimEngineer.Net
1.6.3 Leakage:
η ∞ = Cu 2 / Cu 2
U 2 Cu 2
Ho = .η ∞
g
21
where Po = γ Q H o Constant .
delivered weight
η vol =
delivered weight + internal leakage
Q
=
Q + QL
22
Mech.MuslimEngineer.Net
The hydraulic efficiency is the ratio of work done on the rotor to the
available head of water or (energy actually supplied to the turbine).
η hyd . = H o / H av
23
DIMENSIONAL ANALYSIS
AND
SIMILITUDE OF TURBOMACHINES
2.1 Introduction:
24
Mech.MuslimEngineer.Net
a) Any equation could be reorganized by grouping the physical quantities
assumed in dimensionless groups.
Lmodel Amodel
= Lratio , = L2ratio
Lprototype Aprototype
25
E = f ( D, N , Q v , ρ ,ν ) or
................................................... (2.3)
f ( E , D, N , Q v , ρ ,ν ) = 0
π 1 = D x N y E z ρ = (L )x (T )− y (M z L− zT −2 z ) ML−3
x− z −3= 0
− y − 2z = 0
z +1 = 0
D2 N 2
i.e. π 1 = D 2 N 2 ( gH o )−1 =
gH o
and π 1′ = 1 / π
gH 0
i.e. π 1′ = .............................………....................... (2.5)
N 2D2
26
Mech.MuslimEngineer.Net
π 2 = D x N y E z Qv = Lx T − y M z T −2 z L3T −1 L− z
x− z+3= 0
− y − 2z − 1 = 0
z=0
∴ z = 0, x = −3, y = −1
Q
i.e. π2 = ................................................................. (2.6)
ND 3
π 3 = D x N y E zν = LxT − y M z L− z T −2 z L2T −1
x−z+2=0
z=0
− y − 2z − 1 = 0
∴ z = 0, y = −1, x = −2
i.e. π 3′ = 1 / π 3 = ND 2 /ν .........................….................... (2.7)
2.4.1 Discussion
a) π 1′ = gH o / D 2 N 2 Manometric Coefficient
The head H o is directly proportional to ( N D) 2 , thus for similar
flows and the same efficiency, the proportionality factor for one
machine equals that for other machine.
c) π 3′ = ND 2 /ν Reynolds Number.
π 3′ is the ratio of inertia force to viscous force.
27
28
Mech.MuslimEngineer.Net
So, in non-dimensional form, the relation between the power and the flow
will be a relation between P / ρ N 3 D 5 and Q / N D 3 , Figure 2.2, also for
a geometrically similar pump of different sizes, the same curve could be
obtained.
29
As it could be seen from Figure 2.3 that over Re > 107 their is no
required correction because the flow would be fully-developed turbulent
flow and the effect of viscosity will be minor. With decreasing Reynolds
number, the correction factor Cη increases to take into account the excess
of losses.
It may be suggested that departures from the scaling laws are due
to:
30
Mech.MuslimEngineer.Net
Geometrical dissimilarities due to tolerances
Clearance variations
Surface finish
Hydrodynamic effects (skin friction)
Testing errors
Installation effects.
These were examined for the prototype and four models, three of
which were in aluminum and the other in fiberglass, but being one-eighth
scale for comparison with the aluminum surface produced. Figure 2.4
illustrates the variations observed from the prototype machine when the
three aluminum models were tested.
31
1 − ηY = δ T = δM + δL + δD + δF + δI
Total mechanical leakage disc skin inertia
loss loss loss friction friction loss
The estimation of disc friction loss has been a subject for argument,
as the classical work was done on plain thin discs rotating in a close
fitting closed casing. Nixon used work by Necce and Daily (1960) and
Watabe (1958) for "smooth" and "rough" discs, and showed an error from
measured data of about 10%. Sutton (1968) studies this problem,
particularly the effect of leakage flow through wear rings and its relation
to disc friction. Osterwalder (1978) commented that there is little current
data of general applicability, but Kurokawa and Toyokura (1976) and
Wilson and Goulburn (1976) extended the database.
32
Mech.MuslimEngineer.Net
not propose a general correlation, but indicates that a computer-based in-
house study is relevant.
Table 2.1 A selection of model scale formulae (as quoted for example
by Nixon, 1965). Subscript m denotes ‘model’
Moody (1942):
n
1 − η Dm
= 0 ≤ n ≤ 0.26
1−ηm D
n
1 − η Re m
or = 0 ≤ n ≤ 0 .2
1 − η m Re
Moody:
0.01 0.25
1−η H m Dm
=
1−ηm H D
Anderson:
1−η 0.94 − Q −0.32
=
1 − η m 0.94 − Qm−0.32
Pfleiderer:
0.01 0.25
1 − η Re m D m
= valid between 1/12 < Re m / Re < 20
1 − η m Re D
Hutton:
0.2
1−η Re
= 0 .3 + 0 .7 m
1−ηm Re
Ackeret:
1−η Re 0.2
= 0.5 1 + m
1−ηm Re
33
when n varies from 0.1 to 0.25, the above relation is also used to take into
account the seal effect.
(1 − η A ) p Re
0.2
= 0 .3 + 0 .7 m (2.9)
(1 − η A )m Re
p
34
Mech.MuslimEngineer.Net
becomes 7.23:1. This is clearly important in a multistage compressor
where the stages are designed assuming a constant flow coefficient for
air, but the machine used for a heavier gas (K = 1.05) even though use of
the gas reduces power demand for dynamically similar conditions since
acoustic velocity goes down. The only technique proved satisfactory is to
blend gases to give the right K, as Csanady (1964) showed.
Illustrative examples
35
55 x 10 6
= 25 x 9.81 x 10 3 x Q
0.93
Therefore
Q = 241 m3/s
Q
= constant
ND 3
gH
= constant
N 2D2
The model efficiency must now be found using one of the equations in
Table 2.1. The well-known turbine equation due to Hutton will be used:
0.2
1−η 1894 x 0.3 2
= 0.3 + 0.7
2
1 − η model 94. 7 x 6
ηmodel = 89.77%
36
Mech.MuslimEngineer.Net
= 132.54 kW
m& (RT01 )
= constant
P01 D 2
and
ND
= constant
(RT01 )
m& = 17.24 kg / s
and
2900 D N x 0 .5 D
=
4.124 x 10 3 x 300 0.287 x 10 3 x 288
N = 1499.15 rpm
37
gH o gH o
from π1 : N 2 D 2 = N 2 D 2
1 2
H1 / H 2 = ( N1 / N 2 ) ........................….....................… (2.10)
2
i.e.
Q Q
from π2 : ND 3 = ND 3
1 2
i.e. Q1 / Q2 = N1 / N 2 …..........................................….…... (2.11)
and Q1 / Q2 = H1 / H 2
Remarks:
38
Mech.MuslimEngineer.Net
Fig 2.5 Head-capacity curve for a pump at different speeds
N Q
Ns = 3
(rad) N in r.p.s., .............…...…... (2.12)
( gH ) 4 Q in m3/s,
g in m/s2,
H in m.
39
2.7.1 Pumps:
N Q
Ns = 3
(rpm) * N in r.p.m., ...........…...….... (2.13)
H 4 Q in G.P.M.,
H in feet.
[ Hint: The conversion between the different units of the specific speed is as follows:
gpm ft 3 /s l/s m 3 /s
rpm = 21 . 19 rpm = 1 . 633 rpm = 51 . 64 rpm ]
(ft ) 3 / 4 (ft ) 3 / 4 (m) 3 / 4 (m) 3 / 4
*
The dimensions of Ns is rpm.((gpm)1/2ft-3/4) in the English system or
rpm.((m3/s)1/2m-3/4) in the French system but it is common practice to
omit the units in brackets
40
Mech.MuslimEngineer.Net
Fig. 2.6 Approximate relative impeller shapes and good average efficiencies
obtained for commercial pump as a function of specific speed,
N Q
Ns = 3/ 4
where N in (rpm), Q in (m3/h) and H in (m)
60 H
41
3/ 4 To obtain N s = N Q (gH ) 3 / 4
Ns = N Q H N Q H
(rad), divide by:
Pumps, Compressors, rpm G.P.M. ft 2730
Blowers, Ventilators rpm cfs (ft3/s) ft 129
rpm m3/s m 53
rpm lit/s m 1675
5/4
To obtain
N s = N B.H .P. H N B.H.P. H N s = N B .H .P . [ ρ 1 / 2 ( gH ) 5 / 4 ]
(rad), divide by:
Turbines rpm hp ft 42
rpm Metric hp m 187
rpm kW m 218
N s (rad) for centrifugal impellers varies from 0.0025 to 0.3 and reaches 1
for some axial impellers.
N Q
Ns = 3
(rpm) N in r.p.m., ...........…...….... (2.14)
H 4 Q in C. F. M. (ft3/min),
H in feet.
N B.H .P.
Ns = (rpm) .………….….............…...….... (2.15)
H 5/ 4
42
Mech.MuslimEngineer.Net
2.8 Pressure and Flow Coefficients:
*
Pressure coefficient φ = H /(U 2 2 / g ) …..….............…...….... (2.16)
Flow coefficient ψ = (Q / A) / U ...….….............…...….... (2.17)
Each kind of machine has design values of ψ and φ those fall in a
somewhat range and thus those coefficients are characteristics of
machine's type. For example, considering axial flow compressors φ
values are from 0.2 to 0.6 and ψ from 0.3 to 0.9, for radial flow
compressors φ values are from 0.03 to 0.1 and ψ from 1.1 to 1.5.
1
D( gH ) 4
Dsp = .………….………………..............…...…... (2.18)
Q
*
Sometimes φ is known as speed factor 1 φ = U 2 gH . It's important to notice
that there are many other definitions to those coefficients so, every coefficient will
be redefined for every type of turbine.
43
A comparison of test data is shown in Fig. 2.9 where the solid line
shows the location of efficient turbomachine designs (optimum D sp
values), this line is frequently referred to as the "Cordier" line and the
diagram is known as the Cordier diagram. The diagram indicates that
axial machines dominate the high specific speed region, whereas radial
machines are more efficient in the low specific speed region. From the
diagram also it is clear that N s and D sp values of efficient turbomachines
are close to the Cordier line, thus Cordier diagram is a good tool to permit
the designer to reach the optimum efficiency. Certainly, the performance
of turbomachines depends to a large extent on the compromises made by
the designer regarding the various parameters involved.
44
Mech.MuslimEngineer.Net
Fig. 2.9 "Cordier" line in Ns, Dsp diagram
45
CASCADE MECHANICS
"TWO-DIMENSIONAL APPROACH"
3.1 Introduction:
b. Radial cascade: is the stream surface in the cross sectional planes, this
may be considered as plane flow from source or to a sink.
47
Mech.MuslimEngineer.Net
Fig. 3.1 Straight cascade
The line connecting the two points of the intersection of the camber
line with the inlet and exit of the cascade is known as the chord, the
stagger angle ε is the angle which the chord forms with the perpendicular
to the cascade.
48
49
Mech.MuslimEngineer.Net
Applying the principles of continuity on the control volume of a
unit blade depth, yields:
Fx = t ( P2 − P1 ) ...............................................……................... (3.2)
The forces in the above equation are forces per unit depth exerted
by the blade on the fluid.
To consider the energy losses due to skin friction and other losses
of fluid mechanics, Bernoulli’s equation will be applied neglecting
potential energy.
∆ Po
ζ = ...................………………..………………........ (3.6)
1
ρ Ca2
2
Also similarly, one can define the pressure rise coefficient and the
tangential force coefficient as follows:
50
Fy
Cf = = 2(tan α 1 − tan α 2 ) ..….........……….............. (3.8)
1
ρ t C a2
2
After some manipulations, one can write the following relation between
C f , C p , and ζ :
The resultant force acting on the fluid, F, can be resolved into two
components: one perpendicular to the blade Lift, L, and the other parallel
to the blade Drag, D, Figure 3.3. The lift and drag forces may be then
written in terms of tangential and axial forces as follows:
51
Mech.MuslimEngineer.Net
From equation (3.5):
C m = C a / cosα m
L D
CL = , CD = .........……. (3.14)
ρ lC m2 / 2 ρ lC m2 / 2
t ∆Po cos α m t
CD = =ζ cos 3 α m .....…................................ (3.15)
ρ lC m2 /2 l
t
C L = 2 cos α m (tan α 1 − tan α 2 ) − C D tan α m .…..…….......... (3.16)
l
52
L = ρ ΓC .................................................…….................... (3.19)
After equation (3.13) neglecting the total pressure loss and after
some manipulations, one can obtain the following relation:
Γ = t (Cu1 − Cu 2 )
and L = ρ ΓC m ...................................…................………....... (3.21)
53
Mech.MuslimEngineer.Net
3.5 Cascades in Motion:
Fy .U C a .U
E= = .C f ............…………........................... (3.22)
ρ t Ca 2
54
Note that E has the dimension of energy per unit mass or = gH, from
equations (3.22) and (3.23):
φ = ψ C f / 8 .....................................................…….....…..... (3.24)
η = 1 − E L / E …...................................…….......................... (3.25)
∆ Po
where E L is the energy loss per unit mass and equals gH L = . From
ρ
equations (3.6) and (3.22) we may write:
ζ
η =1− ψ ........................................…............................ (3.26)
Cf
55
Mech.MuslimEngineer.Net
Actually, the performance of a cascade can not be completely
determined theoretically. Recourse to experiment is necessarily, there are
many data published in this field mainly extensive research undertaken by
the NACA.
Figure 3.5 shows the experimental results between lift and drag
coefficients, the deflection θ, and the lift-drag ratio as a function of the
stagger angle at constant fluid angle.
56
The difference between fluid and blade angles at exit is called the
deviation δ = α 2 − α 2′ where α 2 is the deviated flow angle at exit. The
deviation is usually positive.
57
Mech.MuslimEngineer.Net
Fig. 3.7 Variation of nominal deflection with nominal outlet angle
for several space-chord ratios (after Howell)
58
Also exists a kind of correlation between the design point and the
off-design operating points. The off-design performance of a compressor
cascade is published by Howell and is shown in Figure 3.8.
59
Mech.MuslimEngineer.Net
From Figure 3.9, we can remark that:
- The fluid outlet angle remains relatively constant over the whole
range of incidence.
- The reaction blades have a much wider range of low loss
performance than the impulse blades.
Fig. 3.9 Variation in profile loss with incidence for typical turbine
blades after Ainley (adapted from Dixon)
60
( M − M c ) /( M m − M c ) .......................................................... (3.32)
where M m denotes the maximum Mach number which can reach unity.
61
Mech.MuslimEngineer.Net
when calculating velocity field. In real case, fluid particles adhere to the
walls resulting in zero velocity on the walls.
Consider equation (3.27) and neglect the effect of Mach number and
viscosity, thus;
a. Zero-Lift Angle δ ; where the flow does not exert any force on the
blade. The inlet flow angle will be, Figure 3.11;
α1 = δ = α 2
62
tan δ − q. tan σ
tan α 1 = ..…...............................…….......... (3.34)
1+ q
tan δ + q. tan σ
tan α 2 = ..….……….................................... (3.35)
1+ q
with;
Ao Ao
tan δ = , tan σ =
1 − A1 1 + A1
63
Mech.MuslimEngineer.Net
3.9 The Head-Capacity Curve of a Straight Cascade:
ρ
Fy = Q.Γ ................................…...........................…….... (3.39)
t
ρb c
Fy = Q − Q.U ........................………..................... (3.40)
2
t a a
64
Fy
C f = 2 (tan λ1 − tan λ 2 ) =
ρ .t .C a2 / 2
and Q = t.Ca
1
A= (tan α1 − tan α 2 ) , B =1
t
65
Mech.MuslimEngineer.Net
l Zn r
= ln 2 .............................................………..... (3.42)
t 2π cos ε r1
r2
≥ e ( 2π / Z n ) cos ε
r1
The details of radial cascade are lying beyond the scope of these notes.
For further reading, the reader is invited to refer to Csanady.
66
Γ y − yo
U= ................................................……............. (3.43)
2π r 2
Γ x − xo
V =− .............................................…….........….. (3.44)
2π r 2
If the local strength of the vorticity (circulation per unit length) is γ (s) ,
the circulation in an airfoil of chord 1, as shown in Figure 3.15;
l
Γ = ∫ γ ( s ) ds ................................................……................. (3.45)
0
67
Mech.MuslimEngineer.Net
the mean camber line by a vortex at another point P(ξ ,η ) on this line is
approximately that which would be induced at the point on the X-axis
P(x,0) by the same vortex γ ( s ) ds = γ (ξ ) dξ at the point P(ξ ,0) .
1 l γ (ξ )
2π ∫0 ξ − x
V = dξ ............…......................…...................... (3.46)
dy V
= + α ......……............….…......................….......... (3.47)
dx C o
68
1 σ (ξ ).dξ
dU = ......................…….................................. (3.48)
2π x − ξ
and the total induced velocity at P (x,0) by all the source and sink
elements is,
1 l σ (ξ )
U (x) =
2π ∫0 x − ξ
dξ .....….........................……................. (3.49)
C o l d Z dξ
π ∫0 dξ x − ξ
U= ...……….……..........…......…............... (3.50)
l
ξ= (1 − cosθ )
2 ..................................................………….. (3.51)
l
x = (1 − cos τ )
2
dy 1 π sin θ
=α + ∫ γ (θ ) dθ .........……….....….. (3.52)
dx 2π .C o 0 cos τ − cosθ
69
Mech.MuslimEngineer.Net
Then the lift coefficient is defined:
2Γ 2 l
γ ( x ).dx ............................…........…...... (3.54)
l C o l C o ∫0
CL = =
P2 − P1 γ
=4 ………......................................…......... (3.55)
( ρ / 2).C o2 2.C o
and finally the inclination equation can be written on the following form;
∞
dy
= α − A + ∑ An . cos (nτ ) .........….................................. (3.56)
dx 1
γb ∞
= ∑ An . sin (nθ )
2.C o 1
.............…............…......................... (3.57)
∞
dy b
= ∑ An . cos(nτ )
dx 1
70
Γ = ∆φ ..........................….................................................. (3.60)
71
Mech.MuslimEngineer.Net
Fig. 3.17 Conformal transformation of a flat plate cascade
Since the flow far upstream and downstream are parallel, a sink of the
same strength -q and circulation of the same magnitude but of opposite
sign of rotation should be introduced at a point +R as shown in Figure
3.17.
72
q q q q
φq = ln (ξ + R ) + ln(ξ + 1 / R ) − ln (ξ − 1 / R ) − ln(R − ξ )
2π 2π 2π 2π
↓ ↓ ↓ ↓
source at source at sink at sink at
ξ=-R ξ = - 1/R ξ = + 1/R ξ=+R
Similarly, the total potential due to the system of vortex is given by;
Γ Γ
φ Γ = −i ln (ξ + R ) + i ln (ξ + 1 / R )
2π 2π
Γ Γ
− i ln(ξ − 1 / R ) + i ln(R − ξ ) ................................ (3.61)
2π 2π
W = φ q + φ Γ = φ + iψ = c∞ Z e −i ε
t.c∞ −iε R + ξ ξ + 1/ R
W = e ln + e 2iε . ln …………............ (3.62)
2π R −ξ ξ − 1 / R
dW
= c x − i c y = c ∞ (cos ε − i sin ε )
dZ
Putting c ∞ = 1 , we get;
73
Mech.MuslimEngineer.Net
t R +ξ 2 iε ξ + 1/ R
Z=
2π ln R − ξ + e ln ξ − 1 / R ....…...………............... (3.64)
The stagnation points E and A are located on the circle and the difference
in potential φ E − φ A gives a relation between t / l and R.
74
For flat plate cascade, using Eq. (3.16) with C D = 0 , it may be shown
that:
2 4 sin i
CL = cos α m (tan α 1 − tan α 2 ) =
σ σ cos ε
C L∞ = 2π . sin i
then:
CL 2 1
λ= = ..................................…......……....... (3.65)
C L∞ π σ cos ε
75
Mech.MuslimEngineer.Net
2 .m .σ . sin β v
Ct = ...…................................................ (3.66)
4 + m .σ . sin β v
This analysis has been extended to find the potential flow through
logarithmic spiral blades of radial cascade. The analysis of the rotor in
this case differs from that of the stator. Considering, first, the stationary
radial cascade as shown in Figure 3.10, the equation of a logarithmic
spiral is θ = A. log r , from which it can be seen that the angle β v between
the normal to the radius and the tangent to the curve is a constant all
along the curve. In this respect, the spiral blade in the radial cascade is
similar to the flat plate blade of a linear cascade and that the potential
flow through a cascade of flat plates can be transformed by conformal
76
ln (r2 / r1 )
σ = Zn
2π . sin β v
with this definition of σ, the coefficient C t is equal to the one for linear
cascade of flat plates and these values of C t are applicable directly to the
stationary radial cascade. As for rotating cascade, neglecting friction, the
absolute flow may be still considered irrotational, in spite of the fact that
the cascade is rotating. This unsteady potential flow problem has been
solved by Busemann and the total pressure rise through a rotor of this
type may be written as;
Fig. 3.20 Values of shut-off head coefficient ψo for various vane angles βv
and number of vanes for a radial cascade of logarithmic spiral blades
77
Mech.MuslimEngineer.Net
CHAPTER (IV)
INCOMPRESSIBLE FLOW
TURBOMACHINES
( PUMPS )
4.1 Introduction:
78
a. Volute type pump: (Figure 4.2.a) In this kind of pumps a volute type
diffuser is used. The volute diffuser is a spiral, which surround the
periphery of wheel having increasing cross-sectional area as it approaches
the discharge. The volute type is commonly used for single stage pumps,
and for the last stage of multistage pumps.
b. Diffuser type pump: (Figure 4.2.b) In this kind of pumps the diffuser
consists of a number of relatively short expanding passages surrounding
the periphery of the wheel. The diffuser is used usually in multistage
pumps and also in stages developing high heads.
79
Mech.MuslimEngineer.Net
a. Volute Type b. Diffuser Type
As shown from Figure 4.3 for forward blades impeller, the fluid
leaves the impeller with relatively high speed which means that the major
part of the energy gained is kinetic energy, this type of impeller requires a
very good diffuser to convert this kinetic energy to pressure energy. In
practice, it is difficult to construct this kind of diffuser, also it is usually
more efficient to convert pressure energy to kinetic energy rather than
converting kinetic energy to pressure energy.
80
81
Mech.MuslimEngineer.Net
U 22 U 2 .Q
i.e. Ho = + (4.1)
g g .π .D2 .b2 . tan β 2
i. Efficiencies:
a. Hydraulic losses:
1- Friction and turbulence.
2- Disk friction.
The overall efficiency equal the hydraulic efficiency when η mech = 1.0.
82
N s = N Q / H 3/ 4 Q in gpm, H in ft.
N s = N Q / ( gH ) Q in m3/s, H in m.
3/ 4
or
83
Mech.MuslimEngineer.Net
The third type of losses is turbulence losses. The losses due to
turbulence will be a minimum at design flow, for reduced or increased
capacity the turbulence losses increase.
H a = U 22 / 2 g
Remark: to obtain the value of the actual head we may consider that even
at no flow conditions exists a small amount of fluid in the impeller, and it
will rotate with the impeller creating a forced vortex.
P2* C 22 P1* C 12 U 22
γ + 2g − γ + 2g = g
P2* − P1* U 22 U 22
i.e. + = (4.2)
γ 2g g
From the above equation we can see that the gained energy consists of
two equal parts of energy;
i. Pressure energy,
ii. Kinetic energy.
H1 / H 2 = N12 / N 22 , Q1 / Q2 = N1 / N 2
85
Mech.MuslimEngineer.Net
η B.H.P. H
% (hp) (m w)
50 2 20
45 1.8 18 η
40 1.6 16 H
35 1.4 14
30 1.2 12
25 1 10
20 0.8 8 B.H.P.
15 0.6 6
10 0.4 4
5 0.2 2
0 0 0
0 0 0 20 40 60 80 100 1200.1 140 0.01
160 180
Q (lit/min)
Fig 4.7 Pump characteristic curve for a constant speed (N = 2920 rpm)
Q∝N
H ∝ N2
Power ∝ H .Q ∝ N 3
It could be said that if, instead of plotting values of H, Q and power P for
various speeds we may plot values of Q/N, H/N2 and P/N3, Figure 4.9,
then all the points would fall on a single set of curves representing the
performance of the pump at a speed of 1 rpm.
86
i. Viscosity Effect:
ii. Density:
Changes in the density of the liquid will not affect the head
discharge characteristics but they will affect the pressure generated
( P = ρ g H ), Figure 4.10. In steam power stations, centrifugal pump may
have to handle water at higher temperatures, which mean low density, as
for feed pump, which required to deliver a specified weight of water per
87
Mech.MuslimEngineer.Net
hour against a stipulated boiler pressure. The power input to the pump
increases as the temperature rises, for this reason it is more economical to
place the feed heaters on the delivery rather than on the suction side of
feed pump.
Q = V . A = Cv . A 2 gh1 (4.3)
where:
Cv = velocity coefficient
A = clearance area
U 22 − U 12
h1 = head across the orifice = 0 . 75
2g
Leakage will cause efficiency drop and increase the required B.H.P.
where K is a constant.
∆P
γ
=
1
2g
(
C 32 − C 42 ) (0.2 to 0.3) (4.5)
where C 3 is the velocity at the inlet of the diffuser and C 4 is the velocity
at the exit of the diffuser. C 4 is related to C 3 by continuity equation. The
last equation may be written as a function of C 3 only:
∆P C 32
= (0.15 to 0.25 ) (4.6)
γ 2g
2
∆ H loss 1 C 3u
= (0.15 to 0.25 ) (4.7)
H 2φ U
a. Single Seals
90
91
Mech.MuslimEngineer.Net
In cases where high-viscosity liquids (from about 120 mm2/s have
to be pumped, a large spring or metal bellows must always be employed
and the elbows must be relieved of torque transmission by suitable
design. Auxiliary measures are often necessary with single seals too
(Figure 4.13, Table 4.1).
92
Measure Purpose
Pumping liquid circulation to Plan 11 Leading off frictional heat and flushing out
impurities
Cyclone separator to Plan 31 Separating solid particles up to 10%
Requirement: solids heavier than fluid
fluid viscosity ≤ 20 mm2/s
Cooling jacket round shaft seal Lowering the product temperature in the
(neck bush between pump and seal chamber to improve the lubricant film
seal chamber) to Plan C on the contact faces, increasing ∆ t
Pumping liquid circulation via
flow-inducing ring and
water/air cooler to Plan 23
93
Mech.MuslimEngineer.Net
b. Tandem Seals
The tandem seal, functioning analogously to the single type, is employed where
the pumped liquid must be isolated from the atmosphere on account of safety or
environmental considerations. In the principal application area, the secondary seal
should be able to take over the sealing function in full if the primary seal fails. In
subordinate cases, however, the secondary seal serves merely to contain a sealing
fluid. As a general rule the secondary seal ought to have a flow-inducing ring to
initiate and boost the thermosiphon system usually employed. For pumping
liquefied petroleum gas (LPG), intermediate degassing is often provided as well
(Figure 4.14). More recent investigations have shown that this can be omitted if a
liquid seal at about 35% of the shaft sealing pressure is provided instead of the
pressureless seal. It is necessary to establish in situ the exact pressure at which no
leakage takes place on the primary seal. In this way, the flaring-off of liquefied
petroleum gases with its attendant environmental offence is eliminated.
Nevertheless the pressureless seal will remain the commonest method by far,
because providing overpressure always entails higher investments and ongoing
costs. At all events the seal containers must be tested and satisfy the regulations
governing pressure vessels. Non-pressurized sealing fluids do not have to meet any
special requirements except good lubricity (low viscosity, wide margin to boiling
point and explosion-proof). In principle the sealing fluids described under double
seals are applicable.
94
The auxiliary facilities for the primary seal are largely analogous to
those described for single seals. Prime consideration is given to self-
circulation in accordance with Plan 11, or to outside jacket cooling as in
Plan C (Figure 4.13).
A half-way solution between single and tandem seals is the single seal
with stuffing box, which normally consists of two packing rings and an
adjustable gland and is supplied with quenching liquid at a minimum rate
of 40 1/h. Gland leakage is adjusted to 0.25 1/h. With this provision the
leakage can be led off for certain, and there is an assurance that the
function of the main seal can be taken over for a short time in the event of
failure, until the installation has run down to a stop.
c. Double Seals
95
Mech.MuslimEngineer.Net
Fig. 4.15 Double mechanical seal comprising two single seals, type B
(axially-located seat ring on product side and seal with
flow-inducing ring on atmosphere side)
96
For the first class, the mechanical seals should be equipped with a
pumping system to assist the thermosiphon action. Flow-inducing rings
give satisfactory delivery only at a speed of N = 2900 rev/min and seal
sizes of 40 mm.
Due to leakage from the impeller tip, Figure 4.16, the fluid with a
pressure P2 will act on the disk, the same pressure will act on the cover.
As the disk area is larger than the cover, the resultant will be an axial
thrust opposite to the inlet.
97
Mech.MuslimEngineer.Net
Fig. 4.16 Hydrodynamic pressures
Figure 4.17 shows the different methods to eliminate the axial thrust:
a. A double suction impeller type.
b. Balance holes on the disk, but this method may increase the
leakage.
c. Balance disk connected to the suction pressure, the disk area
depends on the amount of the axial thrust and the pressure
difference.
98
99
Mech.MuslimEngineer.Net
a.3. Double Suction Impeller a.4. Two-Stage Back-to-Back Impeller
Arrangement Arrangement
The first step in designing the impeller is the selection of the speed;
the speed relation depends on the driver type mainly; specially for small
units. In large units, the speed should be calculated to meet the optimum
hydraulic efficiency.
100
16 T
Ds = 3 (4.9)
πSs
101
Mech.MuslimEngineer.Net
where T is the torque, S s is the allowable shear stress. The hub diameter
is usually larger than the shaft diameter by 0.8 to 1.5 cm. The inlet
velocity to the impeller eye V o should be kept as low as possible, and
slightly higher than that of the suction flange. V o varies between 3 - 5
m/s.
π
Q /Vo = (D 2
o − D H2 ) (4.10)
4
The inlet blade diameter D 1 is usually made about the same as the eye
diameter D o to ensure smooth flow. C r 1 is usually 10 - 15 % higher than
Vo .
Q
b1 = (4.11)
π D 1ε 1 C r 1
where ε 1 is the contraction ratio which takes into consideration the blade
thickness, for preliminary calculations this could be taken equal to 0.8 -
0.9.
φ = U 2 / 2 gH act (4.12)
Q
b2 = (4.13)
π D 2 .C r 2 .ε 2
ε 2 is the contraction ratio and taken for preliminary calculations from 0.9
to 0.95.
Blade Shape:
D 2 + D1
Z n = 6 .5 sin (180 − β n ) (4.14)
D 2 − D1
where β n = ( β 1 + β 2 ) / 2
π D − t .Z n
ε = (4.15)
πD
103
Mech.MuslimEngineer.Net
NUMERICAL EXAMPLE FOR IMPELLER DESIGN:
SOLUTION:
γ .Q.H 1000 × 10 × 50
Water horsepower = = = 111 hp
75 60 × 75
1700 10 / 60
Ns = = 6 .658 French system
(9 .81 * 50 )3 / 4
or = 2245 English system
It's difficult to predict the bending moment at this time but to care
of it and to keep the critical speed above the running speed, a safety factor
of 1.25 will be introduced.
i.e. Ds =1.25 * 4.8 = 6 cm
The hub diameter DH could be 1 cm larger, i.e. DH = 7 cm. Vo, the suction
velocity (inlet velocity to the impeller eye) will be taken equal to 4 m/s.
π
Q /Vo = (D 2
o − D H2 )
4
4.Q 4 × 10
+ (0.07 ) = 24 cm
2
i.e. Do = + DH2 =
π .Vo π × 4 × 60
Assume Do = D1
104
Q 10
b1 = =
π .D1 .ε 1 .C r 1 60 * π * 0 .9 * 0 .24 * 4 .3
= 0.057 m
where ε 1 is taken = 0.9
** Outlet diameter D2
U 2 = π .D 2 . N / 60 = φ . 2. g .H , φ = 1.0
then,
D2 = 2 × 9.81 × 50.60 / (π × 1700)
= 35.2 cm
U 2 = π .D 2 .N / 60 = 31 .33 m/s
Assume C r 2 = C r 1
Q 10
b2 = =
π . D 2 .ε 2 .C r 2 60 * π * 0 . 351 * 0 . 9 * 4 . 3
= 0 . 039 m
= 3 .9 cm
Assume β 2 = 160°
C u 2 = U 2 + C r 2 / tan β 2 = 19 . 5 m/s
H o = U 2 .C u 2 / g = ( 31 . 33 ) (19 . 5 ) / 9 . 81 = 62 . 28 m
105
Mech.MuslimEngineer.Net
η hyd or η overall = 50 / 62.28 = 0.8 (which is not equal to the assumed
efficiency and seems to be relatively high).
To reach the assumed efficiency, the theoretical head must be
increased. The theoretical head could be increased by means of the
following methods:
D 2 + D1 35 .2 + 24
Z n = 6 .5 sin β n = 6 .5 sin 163 ° = 10
D 2 − D1 35 .2 − 24
Usually, single stage type or two stages with double suction. They
develop a pressure of 100 psia, and have standard capacities of 500, 750,
1000, and 1500 G.P.M. The efficiency is not particularly important as
they are used only occasionally and the flows are not large.
106
W ( xL + t 1 − t 2 )
Flow (GPM) =
500 (T1 − T 2 )
where:
107
Mech.MuslimEngineer.Net
Fig. 4.20 Circulating pump arrangement
These pumps are placed before the boiler. Special care must be
paid to prevent cavitation hence water temperature could be increased.
These pumps are usually driven at a speed of 5500 rpm; the flow through
the pump is generally regulated by water level control that throttles the
flow to the boiler. A check valve is placed in the discharge line to prevent
back flow into the pump.
QG
qGs = QG = gas flowrate referred to suction conditions
QF
QF = liquid flowrate
In this type of pumps the flow is axial, Figure 4.21. The axial
pumps are usually used in irrigation purposes. It develops relatively low
head with large discharge. It has large specific speed.
109
Mech.MuslimEngineer.Net
Fig. 4.22 Velocity triangles
=
(C 2
2 ) (
− C12 + W12 − W 22 )
2g
(**)
From the velocity triangles, Figure 4.22; W12 − W22 = 2.∆C u .W∞u
i.e. σ = W ∞ u / U
(*) This relation could be found by applying Bernoulli's equation on the flow inside
the rotor, and then between the inlet and exit as follows;
( )
P2* − P1* = ρ / 2 W 1 2 − W 22
P − P = ρ / 2 (C − C )
3
*
2
* 2
2
2
3
(**) W 12 − W = (C + W ) − (C + W )
2
2 2
a
2
u1
2
a
2
u2
= (W u 1 + W u 2 )(W u 1 − W u 2 ) = 2W ∞ u .∆ C u
110
The definition is as section 2.8 and using the notations of Figure 4.23:
Pressure coefficient ( )
φ = P3* − P1* / ρ U 2
Flow coefficient ψ = Ca / U
Assume that: θ 1 = θ 3 , C1 = C3
i.e. P3* − P1* = ρ .U .∆ C u
∆C u C a
φ= = (tan θ 2 − tan θ 1 )
U U
C
= 1 − a (tan θ 1 + tan γ 2 )
U
φ = 1 −ψ .(tanθ1 + tanγ 2 ) (4.16)
111
Mech.MuslimEngineer.Net
Fig. 4.24 φ - ψ relation
For the study of the particle motion inside the axial machine, we
can consider the particle trajectory as a co-axial cylinder as shown in
Figure 4.25. So we can write:
C 2
= C a2 + C u2 (4.17)
∂P * C u2
= ρ (4.18)
∂r r
112
H o = P * / γ + C a2 / 2 g + C u2 / 2 g = Constant (4.19)
∂ P* ∂ Ca ∂ Cu
or + ρ Ca + ρ Cu =0 (4.20)
∂r ∂r ∂r
From equations (4.18) and (4.20) we can write the following equation;
∂ C a2 ∂ C u2 2 C u2
+ + =0 (4.21)
∂r ∂r r
∂ C a2 C ∂
+2 u. (r .C u ) = 0 (4.22)
∂r r ∂ r
∂
(C a2 ) + 12 . ∂ (r.C u )2 = 0 (4.23)
∂r r ∂r
The equations (4.18), (4.20), and (4.21) are not sufficient to determine the
flow, another relation could be written as following;
C u = f (r ) (4.24)
113
Mech.MuslimEngineer.Net
2 Cu2
C a2 = Constant − Cu2 −∫ dr (4.25)
r
In the case of free vortex flow, the last equation could be simplified;
C u .r = Constant
∂ C a2
or = 0 i.e. C a = Constant
∂ r
which gives that the axial velocity is constant along the axis.
From the last four equations the flow could be determined to keep
H o = Constant .
115
Mech.MuslimEngineer.Net
4.4 Pump Selection and Applications:
For any application, selecting the proper pump type from the
available different styles and sizes can be a difficult task.
Fig. 4.29 Specific speed chart for the preliminary selection of single speed
centrifugal pumps
Usually used when large capacities are required in such the total
capacity could be divided to two or more identical units. Figure 4.30
shows two pumps in parallel.
116
117
Mech.MuslimEngineer.Net
4.4.3 Economic Considerations:
The unit total cost is made up of the initial purchase price plus the
annual charges required to keep them in operation. The annual charges
include; insurance, power cost, taxes, interest on the investment,
depreciation, and maintenance. For comparison between different bids the
initial purchase price should be divided by the expected lifetime in years.
The running annual cost except the power cost could be taken as
percentage of the initial purchase price.
4.4.4.1 General
118
Fig. 4.32
119
Mech.MuslimEngineer.Net
Fig. 4.33 Minimum submergence Mreq as a function of the flowrate
In the flow range II, the minimum submergence must prevent the
formation of air vortices, which could continue into the pump and be
interrupted by the impeller vanes causing severe vibrations, which could
damage the pump. Here M req is a function of the flow velocity at the
pump inlet.
In the flow range III, the pump NPSH (NPSHreq) is the determining
parameter. The minimum submergence must ensure that cavitation does
not occur at any point inside the pump. If several pumps have to be
installed in one intake chamber, separate bays for the individual pumps
provide the best solution, Figure 4.34.
Fig. 4.34
120
If, for any reason, the extended supply channel length (1 ≥ 5 D),
which is required for trouble-free pump operation cannot be provided, an
alternative consists in fitting a sloping cover to the intake chamber. These
covers are very effective in reducing swirl. Approximate
recommendations for the principal dimensions can be taken from Figure
121
Mech.MuslimEngineer.Net
4.38. However, in each case the dimensions should be determined by the
pump manufacturer.
Fig. 4.39
122
The inlet cross-sectional area A1 of the elbow should be large so that the
velocity at entry to the elbow is insufficient to allow formation of air
entraining vortices or to allow air to be drawn into the pump. In each
individual case an economic assessment of the design has to be made to
determine whether the higher construction costs of an inlet elbow are
acceptable compared to the lower costs of a single intake chamber. The
design and construction of an inlet elbow is more complex than that of an
intake chamber, and in some cases deeper excavations may be necessary.
∆ p = ρ . a . ∆v
2l
t ≤ reflection time t r =
a
124
lv
K=
H
125
Mech.MuslimEngineer.Net
Fig. 4.42 Actual pump installation
The disposition and the drive of the pumps can have a significant
influence of the pipeline layout. Figure 4.43 shows a comparison between
several arrangements using multi-stage, horizontal pumps to obtain the
same delivery rate. Table 4.2 lists the advantages and disadvantages
individually and balances them against each other.
126
b. Two standard pumps in series (each two pumps have one motor in common)
127
Mech.MuslimEngineer.Net
Table 4.2 Advantages and disadvantages of various pipeline layouts (Fig. 4.43)
Efficiency:
Single- and multi-stage, double-flow water transport pumps with the same
delivery head and nq = 30 (specific speed) approximately are used as basis. The
efficiency rises with increasing impeller size. This is brought about by the fact that
greater pumping capacity is achieved by increased impeller diameter, which in turn
requires a greater flow cross section and reduce the effect of the boundary layer.
It was assumed that the average operating time per year of the pumps
amounted to 4000 h. The standby pump; which is not normally functioning was
included in the calculation.
128
G = Increased power consumption x operating hours per year x operating life x price
per kWh.
Value of penalty:
G (Sfr. / pump ) R n − 1
K=
n. Rn R − 1
Table 4.3 shows that the power penalty for larger pumps has to be set higher,
since such pumps have a longer operating life. Statistics concerning power penalties
assuming life. Statistics concerning power penalties assembled over a long period of
time from international invitations to tender confirm that nowadays a mean value of
Sfr. 3500/kW is assumed (Fig. 4.44). The pumping station designer can influence to
some extent the type of pump required by distributing the total delivery flow over
several pumps. The highest possible efficiency is to be obtained by suitable choice of
specific speed; i.e. by variation of the number of stages, of the number of flows and of
the speed of rotation. Figure 4.45 provides information regarding this. It is an
illustration of the peak efficiencies for several types of water transport pumps as a
function of the specific speed of rotation.
129
Mech.MuslimEngineer.Net
Table 4.3 Calculation of power penalties
Fig. 4.44 Hydraulic efficiency as a function of the specific speed and the pump
design. Valid for cold water
130
Variable speed drives are often used for water transport pumps, since the
consumption of drinking water varies according to the season, month, day and time of
day.
Notation:
132
133
Mech.MuslimEngineer.Net
Fig. 4.46 Example for auxiliary and main piping of a centrifugal pump
(Worthington Corp.)
135
136
Mech.MuslimEngineer.Net
Fig. 4.48 Irrigation system using axial pumps (Courtesy of Sulzer Brothers)
137
138
Mech.MuslimEngineer.Net
Fig. 4.52 Vertical motor-mount type of a centrifugal pump
139
Fig. 4.55 An offset-volute design of casing (left) and cover (right). The casing is
designed for top centerline discharge, self-venting, and back pull-out
140
Mech.MuslimEngineer.Net
Fig. 4.56 Centrifugal pump and motor placed on a large sub-base
and connected by a suitable coupling for direct-drive
141
INCOMPRESSIBLE FLOW
TURBINES
( Hydraulic Turbines )
General Introduction:
Many designs has existed but only one has been currently used
named by man who designed it first in California at 1810 Mr. A. Pelton.
Usually, Pelton Wheel, Figure 5.1, is used for high head about more than
300 m.
142
Mech.MuslimEngineer.Net
Fig. 5.1 Pelton wheel
γ Q
Power = (U 1 .C u1 − U 2 .C u 2 ) (5.1)
g
143
γQ
Power = {U [ C1 − (U + (C1 − U ). cos β ) ] } (5.5)
g
γQ
=
g
{U C (1 − cos β ) − U
1
2
(1 − cos β ) } (5.6)
The expression (5.6) shows that the power equal zero when U = 0
or when U = C1 , the maximum value of the power occurs at some
intermediate ratio of U to C1 which can be found by differentiating and
equating to zero;
d ( Power )
= 0 = C1 (1 − cos β ) − 2U (1 − cos β )
dU
then; U = C1 / 2
Thus, the maximum power generated from the Pelton wheel equals:
γQ
(C1 / 4).(1 − cos β )
2
Pmax =
g
η ideal =
(γ Q / g ).(C1 / 2)
2
1
= (1 − cos β )
2
144
Mech.MuslimEngineer.Net
Figure 5.3 shows the theoretical efficiency as a function of velocity ratio.
a. Nozzle Losses:
The flow rate must vary in proportion to the required power for
impulse wheel. This is done by the size of the jet, with a little change in
jet velocity as possible. This can be done by varying the position of the
needle. The velocity of the jet could be given by the following relation:
C1 = C v 2 gh .
At full load, when the nozzle is fully opened, C v is about 0.98 - 0.99, as
the needle is moved to decrease the nozzle opening, the C v coefficient
decreases but it is still above 0.9, Figure 5.4.
145
b. Bucket Losses:
c. Wheel Diameter:
The ratio D/d varies and there is no upper limit, but usually in
practice the ratio is varying between 6-12.
146
Mech.MuslimEngineer.Net
Speed Regulation:
147
148
Mech.MuslimEngineer.Net
Fig. 5.8 Relation between φ, Ns and D/d
5.2.1 General:
The reaction turbines are classified into three types according to the
flow direction, Figure 5.9:
- Radial: (Francis), low specific speed, usually used for medium and
high head installations.
- Mixed: (Francis), medium specific speed, usually used for medium
head installations.
- Axial: high specific speed, usually used for low head installations
(Aswan Dam).
149
5.2.2.1 General:
The Francis turbine original design was inward radial flow with
high flow to make a more compact runner, Figures 5.10 - 5.12. The
diameter was reduced and the water was discharged with a velocity
having an axial component as well as a radial one. This type of runners is
called a mixed flow runner.
150
Mech.MuslimEngineer.Net
Fig. 5.11 Typical radial inward flow turbine profile
151
Using the notations in Figures 5.14 and 5.15, the following expressions
could be written:
γ Q
Power = (U 1 .Cu1 − U 2 .Cu 2 )
g
φ : Pressure coefficient = U 1 2 g h
ψ : Flow coefficient = C r1 2gh
N s = N B.H .P. H 5 / 4
η hyd . = H o H available
η mech. = B.H .P. (γ Q H o / Const.) = H H o
η over . = η hyd . .η mech. = ( H o H av. )( B.H .P. (γ Q H o / Const.))
= B.H .P. (γ Q H av. / Const.) = B.H .P. W .H .P.
For ideal case where η mech. = 1, the hydraulic efficiency equals the overall
efficiency.
Discharge Q = B1 .π D1C r1
= B2 .π D2 C r 2
152
Mech.MuslimEngineer.Net
Fig. 5.14 Velocity triangles (inward turbine)
To enable the turbine to be set above the water level, Figure 5.16, a
kind of diffuser is provided and connected between the water exit and the
tail race level. This is called a draft tube.
153
2
V
The term (hL + e ) represents the tail losses; hL represents the
2g
2
V
hydraulic losses due to friction and enlargement and e represents the
2g
losses due to the kinetic energy rejected in the tail race.
154
Mech.MuslimEngineer.Net
The losses could be presented as follows:
2
Ve
H L = ( hL + )
2g
2
(V1 − Ve ) 2 Ve
= K' +
2g 2g
2 2 2
A V V
= K ' 1 − 1 . 1 + e
A2 2 g 2 g
A1 A1 V1
2 2 2
= K ' 1 − + .
A2 A2 2 g
2
V
= K. 1
2g
2
V
i.e. H L = K. 1
2g
The absolute pressure head at point 1 (turbine exit) could be presented as
follows:
2
P1 Pa V
= − Z − (1 − K ). 1
γ γ 2g
When K equals zero, there will be no energy losses and hence the
efficiency of the draft tube could be written as;
η D.T . = 1 − K
(
hg = V1 − Ve
2 2
) 2g
Introducing the draft tube efficiency η D.T . , the above expression will be;
hg =
(V 1
2
− Ve
2
)η
D.T .
2g
155
The water leaves the turbine with a whirl velocity forming a free vortex.
According to physical laws when the radial velocity of vortex equals to
zero, whirl approaches infinity, which is physically impossible. The
156
Mech.MuslimEngineer.Net
control core of a free vortex cannot follow the free vortex laws and this is
conductive to eddy losses, which are avoided by the solid core.
Unfortunately, for low head installations, the spaces between the turbine
and the water tail race are rarely large enough to permit a draft tube of
types II and III. In these cases, draft tube like type I is usually used. To
assure a complete guidance of the water which help in reducing hydraulic
losses a splitter is used.
2
H net = H av − h f − Ve 2 g
The resulting cavitation that occurs in the turbine not only impairs
the turbine performance but it may also damage the machine itself.
K = ( P / γ − Pv / γ ) (V 2 / 2 g )
where P and V are the pressure and velocity of flow at any point. The
velocity and pressure of the flow are not easily measured and cannot be
157
2
P / γ + V2 / 2 g − Pv / γ
K = 2
H
This coefficient is applied on the turbine exit on the top of the draft tube
or preferably to the highest point in the runner where cavitation might
occur and where H is the net head on the turbine as defined before.
Applying energy equation between the section (2) and the tail race,
Figure 5.16;
2
P2 V2 P
+ + Z = a + hL
γ 2g γ
hL is the head losses between section (2) and the water tail race level.
These losses include the kinetic energy rejected to the water and the
hydraulic losses.
P P
σ = a − v − Z H (5.10)
γ γ
158
Mech.MuslimEngineer.Net
Fig. 5.19 Critical cavitation factor
Pa Pv
Z= − − σ c .H
γ γ
N B.H .P.
Ns = English Units
H 5/4
The critical cavitation number depends upon the type of turbine and is a
function of the specific speed. For Francis turbine, it is given by:
2
N N B.H .P.
σ c = 0.625 s , ( Ns = metric units)
444 H 5/ 4
159
2.73
N N B.H .P.
σ c = 0.3 + 0.0024 s , ( Ns = metric units)
100 H 5/ 4
For Kaplan turbine, σ c is 10% higher than that for a similar propeller
turbine.
Figure 5.21 shows the relation between the head and the specific
speed for different types of turbines. The data has been taken from
installations successfully at work.
160
Mech.MuslimEngineer.Net
Fig. 5.21 Turbine classification according to specific speed and head
161
Also, the maximum value φ is greater than one for the performance
at constant speed. The efficiency curve is not as flat as that for impulse
wheel. In the impulse wheel (turbine), the velocity triangles are
independent of the quantity of the flow and so theoretically the hydraulic
efficiency should be constant at all loads at constant speed.
A. Propeller Turbine:
The blades are fixed on the rotor and the number of the blades is
from 3 to 8. In the simplest form of axial flow propeller turbine, the
blades are cast integrally with the hub. The propeller turbine is used in the
measurement of the flow, especially in petroleum industry to measure the
flow rate.
B. Kaplan Turbine:
162
Mech.MuslimEngineer.Net
Although the propeller turbine is almost adequate for high flow and
low head operation, it has one quite serious disadvantage is its part load
efficiency is unsatisfactory. In Kaplan turbine, this problem is overcome
by using movable blades, so that their angles of inclination may be
adjusted while the turbine is in motion in such a way that the turbine can
operate continuously at its maximum efficiency.
Efficiencies:
1/ 5
1 − η Dm
=
1 − ηm D
163
For Francis and propeller turbines, Figure 5.25 shows the average values
of the ratios φ , ψ and D / d suited for various specific speeds.
a. Impulse Turbine:
164
Mech.MuslimEngineer.Net
Fig. 5.26 Pelton wheel installation
b. Francis Turbine:
165
166
Mech.MuslimEngineer.Net
Cedillo, Spain: Kaplan turbines on the Tagus river
The Cedillo hydroelectric power station on the Tagus River is fitted with
four ESCHER WYSS Kaplan turbines and governors, together with
electricity generators supplied by General Electrica Espaňola, Figures
5.28 to 5.30. The main data of the turbine installation are:
Head 43 m
Output per turbine 110 MW
Speed 93.8 rev/min
Runner diameter 7600 mm
167
168
Mech.MuslimEngineer.Net
Fig. 5.30 Turbine unit supported on three bearings.
The runner blades can be dismantled via a port in the draft tube
169
Head 6.2 m
Output per turbine 24.6 MW
Speed 62.1 rev/min
Runner diameter 7700 mm
Weight per turbine 1130 t
170
Mech.MuslimEngineer.Net
Fig. 5.31 Workshop assembly of complete guide vane apparatus (two-part guide
vane casing, 18 adjustable guide vanes in closed position) and of the runner
Fig. 5.32 Sectional drawing of the Racine bulb turbine, which has a runner
diameter of 7700 mm. The bulb nose is fastened into the concrete pier of the inlet
171
172
Mech.MuslimEngineer.Net
Fig. 5.34 The principles of hydraulic transmission
P2 T2 .N 2
η= =
P1 T1 .N 1
The fluid coupling, Figure 5.35, transmits the power between the
driver and the driven unit. The power transmission could be at the same
speed but it permits also a variation in the speed. This is very practically
in the starting of a large machine. Thus, it could be act as a slipping
clutch.
173
P2 T2 .N 2
η= =
P1 T1 .N 1
174
Mech.MuslimEngineer.Net
Fig. 5.36 Fluid coupling efficiency curve
175
The stator or guide vanes, Figure 5.38, are so shaped as to change the
direction of flow and thus it could increase the angular momentum and
the fluid flowing through the turbine given up an angular momentum
equal to the sum of that from the stator and the pump. So, the torque of
the turbine, which is delivered to the output shaft, is the sum of the pump
torque and the stator torque. This could be easily seen from Figure 5.39.
176
Mech.MuslimEngineer.Net
Fig. 5.39 Torque converter performance
177
178
Mech.MuslimEngineer.Net
Fig. 5.41 Application ranges of different types of water turbines
180
181
Mech.MuslimEngineer.Net
Fig. 5.43 Horizontal section in the plane of the annular duct, through one of the
two vertical-shaft four-nozzle impulse turbines for the power station El-Salto II
(Colombia). Short nozzles, in conjunction with a close-fitting turbine casing and
internally situated servomotors, result in small machine dimensions and a
correspondingly reduced distance between centres. The nozzles are bolted
simultaneously to the pressure pipe and to the turbine casing by means of a
common flanged section and can be removed, in the same way as the runner,
through the discharge pit. P = 37,000 kW, H = 415.5 m, N = 400 rpm
182
183
Mech.MuslimEngineer.Net
Fig. 5.45 Francis turbine runner
184
Fig. 5.47 Four 260 MW Pelton turbine runners under a head of 1140 m
with the model runner inside one of them
185
Mech.MuslimEngineer.Net
Fig. 5.48 Kaplan runner for a bulb turbine (diameter 3.8 m,
output 11.3 MW, head 14 m)
186
187
Mech.MuslimEngineer.Net
Fig. 5.51 Pelton wheel
188
COMPRESSIBLE FLOW
TURBOMACHINES
( Thermodynamic Principles )
f ( P, v, T ) = 0
where
P = the pressure,
v = specific volume = V/M,
T = absolute temperature.
189
Mech.MuslimEngineer.Net
R = 1.9875 BTU/(lbmole .°R), kcal/ (kgmole .K)
The ideal gas obeys the above equation, for the real gases the
above equation could be applied with good approximation at low
pressures. There are other equations of state, which deal with real gases;
this equation lay out of scope of this text. This could be easily found in
any textbook on thermodynamics.
6.3 Enthalpy:
C2 C2
H 2 + gZ 2 + 2 − H 1 + gZ 1 + 1 = WT + QT (6.5)
2 2
H = E + P.dv (6.6)
190
H T = H + gZ + C 2 / 2 (6.7)
i.e.
H T 2 − H T 1 = WT + QT (6.8)
∆H = ∆E = C v .∆T (6.9)
6.4 Entropy:
2
dQ
ds = ∫ kcal / K (6.10)
1 T
6.5 Work:
v2
W .D = ∫ P.dv (6.12)
v1
191
Mech.MuslimEngineer.Net
Fig. 6.1 System displacement
W .D = P.dv = 0
192
v2
R.T
W .D = ∫ dv = R.T . ln(V2 / V1 ) (6.14)
v1 v
and
v2
W .D = ∫ P dv
v1
k −1
(6.16)
k P
P1 .V1 2 − 1
k
=
k −1 P
1
k −1
k P2 k
W .D = R.T1 − 1
k −1 P1
k
= (P2 .V2 − P1 .V1 )
k −1
P.V n = Constant
n −1
n P
P1 .V1 2 − 1
n
W .D = (6.17)
n −1 P
1
193
Mech.MuslimEngineer.Net
and;
ln (P2 / P1 )
n=
ln (V1 / V2 )
This law states that the energy change in a system is exactly equal
to the resultant energy transfer across its boundary. Resultant energy
transfer means the difference in energy transfer in and out of the system.
i.e. E 2 − E1 = Q − W (6.18)
dU = dQ − dW (6.19)
194
W T1 − T2
η= = (6.20)
Q T1
W .Disen H 2 − H 1
η ad = = (6.21)
W .Dreal H 2 '− H 1
and H = C p .∆ T , i.e.
∆Tisen
η ad = (6.22)
∆Treal '
195
Mech.MuslimEngineer.Net
Fig. 6.3 Adiabatic compression
k −1
k P2 k
W .Dad = P1 .V1 − 1
k −1 P1
196
W .Diso
η iso =
W .Dreal
W .Diso = R.T1 . ln (P2 / P1 )
W .Dreal = H 2 − H 1 − Q = C p (T2 − T1 ) − Q
R.T1 . ln( P2 / P1 )
ηiso = (6.23)
C p (T2 − T1 ) − Q
W .D poly
η poly = or
W .Dreal
η poly =
(k − 1) / k (6.24)
(n − 1) / n
197
Mech.MuslimEngineer.Net
To go deeper insight, we may consider the treatment of isentropic
flow in two dimensions. Since this treatment is highly complicated, the
so-called linearised compressible flow theory will be considered.
It is assumed that there are only small deviations in the flow field
from an undisturbed parallel flow of high speed. For a stream of mean
velocity C, with perturbation velocities u and v, the vorticity is
everywhere zero.
∂v ∂u
− =0 (6.25)
∂x ∂ y
∂ ( ρ .C + ρ .u ) ( ρ .v )
+ =0 (6.26)
∂x ∂ y
∂φ ∂φ
u= , v= (6.27)
∂x ∂ y
Using the small perturbation hypothesis, the density and the speed of
sound can be written as:
ρ = ρ + ρ'
a = a + a'
∂u ∂v ∂ ρ'
ρ +ρ +C =0 (6.28)
∂x ∂ y ∂x
1 dP
dC2 + =0 (6.29)
2 ρ
198
− ρ .C
dρ= du (6.31)
a2
After equations (6.28) and (6.31), one can obtain the following;
∂ 2φ ∂ 2φ
(1 − M 2 ) + =0 (6.32)
∂ x2 ∂ y2
ξ = x, η = y 1− M 2 , φ * = K .φ
u* = K u
K
v* = v (6.33)
B
v* K
tan θ * = = . tan θ
C B
199
Mech.MuslimEngineer.Net
with B = 1 − M 2
dy v v
= ≈ ≈ tan θ (6.34)
dx C + u C
dη v* v*
= ≈ ≈ tan θ * (6.35)
dξ C + u * C
dη dy
=B (6.36)
dξ dx
ξ = x, η = B y, φ * = B 2 .φ (6.38)
P' = ρ .C.u
(6.39)
P* = B 2 .P'
It arises from the above analysis that the lift is 1 / B 2 times greater
for the profile in compressible flow than for the corresponding profile in
incompressible flow.
200
201
Mech.MuslimEngineer.Net
CHAPTER (VII)
7.1 General:
7.1.1 Fans:
Fans are usually used to circulate air and usually axial flow type.
The pressure does not increase, one can assume it constant during the
process. So, the problem could be simplified to an incompressible flow
problem. The energy supplied by fan to fluid is mainly a kinetic energy.
202
7.1.2 Blowers:
Usually used when large capacity is required. It is also used to
circulate gases. The pressure is slightly increased but not more than
2 kp/cm2. The blowers could be axial or mixed flow. If the pressure
increment is not sensible, the problem could be treated as incompressible
flow problem. The blowers have other names according to the service in
which it operates. For example, in gas service blower used to remove
gasses from a coke oven is known as an exhauster.
7.1.3 Turbo-compressors:
Usually used when a compressed gas is required, the exhaust
pressure is not less than 2 kp/cm2. Centrifugal compressors, Figure 7.2,
are used when high pressure is used and relatively low volumetric flow is
required.
The axial flow compressors are used when large volumetric flow
and low exhaust pressures are required. In this chapter, we shall deal with
centrifugal compressors.
203
Mech.MuslimEngineer.Net
7.2 Head and Power:
From the work done equation, we can write the head equation if we
divide by the weight flow rate.
i.e.
n −1
n P2 n
H = R.T1 − 1 ...................................(7.1)
n − 1 P1
U 2 .Cu 2 − U1.Cu1
HO =
g
φ = U 2 / 2 g .H
φ ′ = g.H / U 2 .....................................(7.2)
204
γ QH
Power P = which usually called the gas h.p.
Const.η pol
η pol : the polytropic efficiency, sometimes called the hydraulic efficiency.
φ′ = ……................... (7.4)
2
U
205
Mech.MuslimEngineer.Net
φ s = Cu 2 /U 2 ...............................................(7.5)
206
n n −1
H poly = R.T1 (P2 / P1 ) n − 1 ...........................(7.6)
n −1
207
Mech.MuslimEngineer.Net
Fig. 7.6 Performance curve for centrifugal compressor. Solid lines are
characteristics at various speeds; dotted lines are corresponding efficiencies
208
The surge flow and maximum flow rates for multistage machine
reflect the surge flow and maximum flow rates of the particular stages in
the series which reach unstable flow (as flow decreases) and maximum
flow (as flow increases) first.
209
Mech.MuslimEngineer.Net
the pressure ratio-capacity curve steep while the low molecular weight
gases tend to flatten the curve.
The term stability means the stable flow range between surge and
normal flow. Stability with respect to normal flow is defined as the
percent of change in capacity between normal and surge point at normal
speed; this equals 100 % minus the percent ratio of the actual volume
flow rate of which surge begins (at normal speed) to the normal volume
flow rate. The term “turndown” means the mass flow reduction, which is
possible before encountering surge, recognizing the effects of the control
method. “Turndown” with respect to normal flow is defined as the
percent of change in capacity between normal and the surge point at rated
head when operating at design inlet temperature and gas composition; it
equal 100 % minus the percent ratio of the surge point mass flow at
normal head to normal mass flow.
210
Shock and flow separation are distinct possibilities even when the
Mach number is less than unity.
Consider now the problem at the impeller eye;
W1
i.e. M1 = ...............................................(7.7)
g .K .R.T1
It could be seen that the maximum Mach number occurs at the outer
radius of the eye at which W has its maximum value. To avoid shock and
separation, Mach number at inlet must be about 0.85. Thus even if it
reaches unity at impeller tip, diffuser will decelerate the flow.
7.6 Pre-Whirl:
211
Mech.MuslimEngineer.Net
7.7 Surging:
Surging occurs when the losses are equal to the energy gained, this
will cause the flow to stagnate causing backward pressure waves and will
then move backward. This will cause the damaging of the machine.
It is difficult to predict the surge point, only experience could
determine exactly the surge point. Usually, any operating point should lay
10 to 15 % away from the surge point.
16.T
Ds = 3 …............................................(7.8)
π .S s
π
Q1 / V = (D 2
o − DH2 ) ...............................................(7.9)
4
212
b1 = Q1 (π .D1 .C r1 .ε 1 ) ..............................................(7.10)
H = φ ′.U 22 / g
k −1
k
and H = R.T1 ( P / P )
1 .η p − 1
k ......…................(7.11)
(k − 1) η p 2
U 2 = π .D2 .N / 60
213
Mech.MuslimEngineer.Net
Now, the outlet velocity triangle could be drawn. The circulatory
flow effect reduces the tangential component C u 2 by an amount equal to
X, which could be calculated by the following formula;
π . sin (180 − β 2 )
X = ....................................(7.12)
Z
214
215
Mech.MuslimEngineer.Net
Fig. 7.13 Sectional view in a centrifugal compressor showing one stage for test
purposes (courtesy Nuovo Pignone)
216
Centrifugal
Foreign objects centrifuge toward heavy blade root
No blade clearance problem
High pressure ratio provided with 2 stage compressor
Blade fineness ratio-low (Higher relative strength and greater shock
resistance)
Efficiency not materially affected by film build-up or coating of foreign
material
Large foreign objects rejected at compressor face without significant
performance degradation
Axial
Foreign objects centrifuge and mass toward and against delicate blade tips
Minute clearances between blades and stators increases sensitivity to fod
Exposure to fod is high due to number of stages required for a given
pressure ratio
Blade fineness ratio-high (lower relative strength and shock resistance)
Significant performance degradation due to foreign material filming
Rotating clearances and design geometry will not accept ingestion of large
foreign objects without possibility of total compressor failure
184
Mech.MuslimEngineer.Net
CHAPTER (VIII)
VOLUMETRIC MACHINES
218
π 2r n
Q= D2 (8.1)
4 60
The actual volumetric flow rate is lower than the theoretical because of
leakage.
η v = Qa / Q (8.2)
219
Mech.MuslimEngineer.Net
Fig. 8.2 The instantaneous velocity in delivery pipe
2π .n.r D 2
VP = . sin θ (8.4)
60 d 2
To overcome this fluctuation in the flow, air vessel is used, Figure 8.3.
The air vessel is a container filled with compressed air at a certain
pressure. The container is placed on the delivery pipe, during the delivery
stroke all fluid in excess of the mean discharge will be stored in the air
vessel, this will be replaced during the suction stroke.
220
Some slight change in the air pressure must naturally take place as
the fluid surface in the vessel rises and falls. At high pressure, air is
slowly dissolved in liquid, the dissolved air must be replaced in large
units and air compressor is used.
221
Mech.MuslimEngineer.Net
Fig. 8.4 The damping of fluctuation
222
223
Mech.MuslimEngineer.Net
Fig. 8.6 Typical positive pump installation
224
225
Mech.MuslimEngineer.Net
Fig. 8.8 Parallel cylinder rotary pump
2.a.l.n.N
Q= (8.5)
60
where:
a = area enclosed between any two adjacent teeth and the casing.
l = axial length of teeth.
226
Q = K C (D − C) l (8.6)
where:
K is a constant = 2.983
D is the gear diameter.
C is the distance between the two gear centers.
The gear pump is usually used in the lubrication of petrol engines, as lube
oil pump.
The performances of positive pumps are generally the same for all
kinds of this pump, Figures 8.10 and 8.11. This kind of pumps discharge
227
Mech.MuslimEngineer.Net
a definite quantity of liquid for each revolution, provided their suction is
flooded and the discharge head is zero.
228
229
Mech.MuslimEngineer.Net
Consider a short time interval, thus substitute dV p / dt for V p / t , the
acceleration of the water in the horizontal pipe = dV p / dt ;
dV p / dt = (dV p / dθ )(dθ / dt ) = ω .r.( D 2 / d 2 ). cos θ .(dθ / dt )
and dθ / dt is the angular speed of crank shaft (ω), thus;
dV p / dt = ω 2 .r.( D 2 / d 2 ). cos θ
It is clear from Figure 8.12 that the mean pressure in the delivery line is
AD. If fd is smaller than the vapor pressure, cavitation will occur. It is
important to notice that the work done is constant, hence the area fADg
equals the area fBCg.
230
1. General:
The energy of the fluid leaving the blade consists of kinetic energy
C2 P
and pressure energy ( + ). Usually, any pump or blower is used to
2g γ
raise the fluid pressure, so the efficiency of the machine will be increased
as much as one could convert the kinetic energy to pressure energy. The
pressure recovery devices, as a general, consist of gradually expanded
passages, which provide a gradual deceleration to the fluid motion. As a
general, the ratio of the kinetic energy with reference to static pressure
energy is higher in forward curved impeller than backward. This process
of pressure recovery could be responsible of fairly high losses.
where A1 and A2 are the inlet and exit areas, respectively. The efficiency
of the diffuser is defined as the ratio of the static pressure gain to the
kinetic energy.
(∆P) actual
ηd = ............................. (2)
C12 (1 − A12 / A22 )(γ / 2 g )
231
Mech.MuslimEngineer.Net
Fig. 1 Diffuser. (a) Boundary-layer increase in diffuser. (b) and (c) Stall zones.
(d) Short vane dividers. (e) Boundary-layer absorption.
(f) Diffuser with standing vortices.
232
233
Mech.MuslimEngineer.Net
Fig. 3 Boundary layer effects at the exit of a cascade
+∞
W1 − W
δd = ∫ W
dn ................................................. (3)
−∞ 1
+∞
W W1 − W
δm = ∫ W1 .W
dn ................................................ (4)
−∞ 1
1 t/2
C a = ∫ W cos α 1 .dy ................................................ (5)
t −t / 2
234
t/2
ρ sin α 1 cos α 1 ∫ W 2 dy = ρ C a C r t …………………............ (7)
−t / 2
with
W 2 = W12 − (WW1 − W 2 ) − W1 (W1 − W ) ….......................... (8)
δd δm
∆= , θ = ................................................. (9)
t cos α 1 t cos α 1
then
Ca
W1 cos α 1 = ............................................... (10)
1− ∆
The physical meaning of equation (10) is that the through flow velocity is
increased by 1 /(1 − ∆) . After equations (7), (8), (9), and (10), one can
write;
or
1− ∆
W1 sin α 1 = C r .....................................................….. (12)
1−θ − ∆
Apply now the momentum balance in X-direction between (1) and (2);
235
Mech.MuslimEngineer.Net
P2 − P1 cos 2 α 1 t 2
∫
2
= W dy − C a ...........................…......…… (13)
ρ t 0
P2 − P1 1 −θ − ∆
= C a2 − 1 ..................................…...…....... (13*)
ρ (1 − ∆ ) 2
W12 P2 W22
P1
EL = + − − ......................................……......... (14)
ρ 2 ρ 2
with
C r2 (1 − ∆ )
2
C a2
W =2
+
1
(1 − ∆ )2 (1 − θ − ∆ )2
W22 = C a2 + C r2
EL
ζ = or
C a2 / 2
2θ 2δ m
ζ = = ..........................…........................... (15)
cos 2 α 2 t cos 3 α 2
2. Diffuser Types:
236
237
Mech.MuslimEngineer.Net
Fig. 4 Vaneless diffuser
238
1) Continuity equation:
∂ Ca ∂ C a Cu ∂ C a Cu2 ∂P
ρ
. + Ca + . − =− .......................... (17)
∂ t ∂ r r ∂ θ r ∂ r
∂ Cu ∂ Cu Cu ∂ Cu C a Cu 1 ∂P
ρ . + Ca + . + =− . ................. (18)
∂t ∂r r ∂θ r r ∂θ
The main dimensions of the vaneless diffuser are its outer diameter D4 ,
impeller diameter D2 , the distance between the two parallel walls b3 , and
the impeller width b2 . Following are the ratio between these dimensions:
D4 / D2 = 1.8 , b3 / b2 = 0.8
239
Mech.MuslimEngineer.Net
Fig. 6 Vaned diffuser
Cu 3 .r2 C r 2 .r2
Cu = , Cr = ……………....... (19)
r r
where Cu and C r are the tangential and the radial velocities of fluid,
respectively. Consequently:
C r Cr 2
tan α = = = Const. ..................................... (20)
Cu Cu 3
2π
where φ corresponds to the blade log spiral pitch. φ = where Z is
Z
equal to number of blades. Z should not be equal to the number of
impeller blades.
240
241
Mech.MuslimEngineer.Net
2.3.1 Parallel Walls:
ro
Cr = Cro .................................................................. (26)
r
C r C ro
i.e. tan α = = = Const. .........................………......... (27)
Cu Cuo
242
r C
ln = φ . ro = φ . tan α o ..…...………......................... (29)
ro C uo
The above equation is the outside boundary curve for the volute. For the
outline curvature, the following equation could be used;
r
R= ...............................………...……............ (30)
cos α
The last solution is for constant width diffuser. In some cases, bo does not
equal b. This would result in an additional factor bo / b , resulted from the
continuity equation as follows:
ro bo
C r = C ro ....................................................... (31)
rb
r b
ln = φ . tan α o . o .................……….…........... (32)
ro b
C ro ro bo
Cr = ..........................……... (33)
[bo + 2 (r − ro ) tan δ ] r
and;
Cr dr C ro bo
tan α = = = ................ (34)
C u r dφ C uo [bo + 2 (r − ro ) tan δ ]
243
Mech.MuslimEngineer.Net
1 ro r ro r
φ=
1 − 2 tan δ ln
r + 2 − 1 tan δ ...... (35)
tan α o bo o bo ro
and;
φ = f (r ) ..................................... (36)
ro
1− 2 tan δ = 0 ............................................... (37)
bo
which means:
bo
tan δ =
2 ro
δ is usually from 8° to 15°.
This type consists of gradually increasing rectangles, Figure 10. Thus, the
area expands in axial and radial directions. We can put h/b = Const. At
C r
position x, C u = uo o . Thus, we can consider a deviation from the line
ro + x
244
φo h
Cuo ro h
V =∫ dx .............................................. (38)
360 o ro + x c
360 Cuo ro h h
φ= ln 1 + ..................................... (39)
Vc ro
245
Mech.MuslimEngineer.Net
REFERENCES “Appendix I”
246
THEORY OF CAVITATION IN
CENTRIFUGAL PUMPS
1. INTRODUCTION:
246
Mech.MuslimEngineer.Net
Fig. 1 A condition whereby an air bubble can remain stable
indefinitely, without solution on a hydrophobic particle
possessing a concave surface (after Liebermann)
247
e- From airfoil theory (Fig. 2), when a flow around an airfoil occurs a
region of low pressure is generated on the upper part and high
pressure on the lower part. The inlet edge of impeller vanes acts in
a manner similar to an airfoil, giving a local pressure rise on the
leading face and a pressure reduction on the trailing face.
248
Mech.MuslimEngineer.Net
Fig. 2 (a) Actual pressure distribution on the airfoil designed
249
3. SIGNS OF CAVITATION:
250
Mech.MuslimEngineer.Net
the suction pressure, increasing for higher specific speed and lower
suction pressure. The difference in the behavior of pumps of different
specific speeds results from the difference in the impeller design. Low
specific speed impeller vanes form a definite channel, the length of which
depends on the vane angles, the number of vanes, and the ratio of the
impeller eye diameter D1 to the impeller outside diameter D2 (Fig. 5 (a)).
When the pressure at the impeller eye reaches the vapor pressure, usually
on the backside of the vane entrance tips, it extends very rapidly across
the whole width of the channel, A-B, Fig. 5 (a)-(a), with a small increase
in capacity and decrease in head.
251
A further drop in the discharge pressure does not produce any more
flow because the pressure differential moving water to the impeller eye
cannot be increased anymore. This differential is fixed by the suction
pressure outside the pump, and the vapor pressure across the whole
channel between any two vanes at the impeller entrance.
With high specific speed impellers, the channel between two vanes
is wider and shorter, Fig. 5 (a)-(b). More drop in head and a greater
increase in capacity are required to extend the vapor pressure zone across
the whole channel. Therefore, the drop in the head-capacity curve extends
through a wider range before the sudden break-off occurs. With propeller
pumps the vanes do not overlap, Fig. 5 (a)-(c).
252
Mech.MuslimEngineer.Net
constant speed and capacity, Fig. 5 (b). The drop in efficiency and head
curves indicates the beginning of cavitation.
4. MECHANISMS OF DAMAGE:
253
Fig. 7 (a) The jet collapse model for damage due to cavitation
254
Mech.MuslimEngineer.Net
Fig. 7 (b) Schematic representation of successive stages of nonsymmetrical
cavity collapse with microjet impingement against a metallic surface
Tillner et al. (1990) and Pearsall (1978) who show ways that damage can
be produced in the laboratory for test and development purposes describe
typical damage. The extent of the damage suffered depends on the fluid,
the materials and the hydrodynamic system, and it has been found that
even with advanced material loss the machine has developed the duty
255
Material Relative
volumetric loss
Stellite 1
Cast stainless steel: 12.88% Cr, 0.17% Ni, 0.43% Mn, 0.38%Si 7
Stainless steel 18:8 Cr:Ni 5
Monel 16
0.33% carbon steel 37
14% Cr stainless steel (forged or drawn) 98
Manganese bronze 118
Gun metal 230
Cast iron (as cast without skin) 374
Typical cast aluminum alloy 1176
256
Mech.MuslimEngineer.Net
Fig. 9 (a) Schematic representation of successive ranges of growth, collapse,
and rebound of a single traveling cavity.
(b) Graph of cavity diameter as a function of time for the cavity in (a)
The cavitation performance of a pump varies with the fluid condition. For
example, when pumping water the required net positive suction energy
(NPSE) is highest when passing cold water and decreases as the
temperature rises. A similar effect is noted when other liquids are being
pumped. The change is too large to be explained by the so-called
Reynolds number effects. An empirical approach to this problem is
known as the B-factor or β-factor method outlined by Knapp and Daily
(1970) and Stahl and Stepannof (1956); the technique correlates suction
hydraulic behavior with vapor volume in the cavitating region.
When examining bubbles in cold water it has always been assumed with
some justification that all energy terms involving the vapor in the cavities
are negligible compared with those of the surrounding liquid. When the
temperature increases such an assumption may not be applied since the
latent heat required to supply vapor to the cavity can no longer be
ignored, and neither can the energy exchange during the expansion and
contraction of the bubble. If they are neglected and the NPSE calculated
257
The heat required by vaporization must come from the liquid surrounding
the cavity, thus casing a drop in temperature and vapor pressure in the
immediate vicinity of the bubble. This has the effect of reducing the
bubble size from that which would apply in cold liquid, thus reducing the
effect on the flow. This reasoning has led to the approach of Stahl and
Stepannof (1956). Figure 10 shows the conventional NPSE variation
curves for a cold and a hot liquid. The two critical points where the 3%
head drop applies are B and C. Point B is known for the cold test, and the
NPSE reduction can be found to estimate point C.
volume of vapor
B= (1)
volume of liquid
258
Mech.MuslimEngineer.Net
2
v
B = ∆( NPSE )C p T v (2)
v h
1 fg
or
2
B v
B`= = C pT v (3)
∆ ( NPSE ) v h
1 fg
Figure 11, taken from Stahl and Stepannof (1956), plots B` for a number
of fluids based on refinery pumps of the double suction design and a 3%
fall in gH. Also shown are lines of NPSE adjustment in the relation:
− NPSE adjustment
σ corrected = σ (4)
gH pump
The method is based on the assumptions that the cavities are uniformly
distributed across the flow cross-section and that there is the 3% drop
criterion. There are many other approaches but these indicate that there is
a considerable difference in opinion (see for example Hutton and Furness,
1974).
259
T.E.L.
H.G.
Hmd
hsd
(2)
Hst
Hm
hss
Hms h1
Pat (1)
T.E.L. h3 Datum
H.G. hf
Vs2/2g
H m = Manometric head
H ms = Manometric suction head
H md = Manometric delivery head
hss = Static suction head
h f = Friction head loss (included all losses such as entry losses
2 2
Vs / 2 g , so h f = h fs + Vs / 2 g ).
Apply Bernoulli’s equation between (1) and (2), (i.e. between the
suction pipe flange and pump center assuming that point (2) is on
impeller entrance).
2
P2 C2 Pat C1 2 C
2
+ + h1 = + + h3 − h f ( 1 =0)
γ 2g γ 2g 2g
260
Mech.MuslimEngineer.Net
2
P2 C2 P
+ = at − hss − h f (5)
γ 2g γ
P2 C 2 2
+ − hv ≤ 0 (6)
γ 2 g
The above term is called “Net Positive Suction Head” 'NPSH' or H sv , i.e.
Pat
H sv = ± hss − h f − hv (7)
γ
The available suction head must be larger than the required suction head.
Example: Find the available NPSH for pump pumping water, if the pump
is above the water surface by 3 m and h f = 0.7 m and the atmospheric
pressure is 1.03 kp/cm2, vapor pressure is 0.0355 kp/cm2 at 20°C,
γ = ρ g = 1000 kp/m 3 .
1.03 x 10 4 0.0355 x 10 4
H sv = − 0 .7 − 3 − = 6.24 m
1000 1000
261
From the above points, all these losses occur between the pump
suction flange and impeller entrance and are proportional to W 2 and
could be equal to K .W 2 / 2 g where K is a coefficient. If the required
suction head is larger than the available suction head cavitation will
occur.
30( H − h) sg
Ps = Pb − Pv + (8)
34
where:
Ps vacuum to be maintained at pump suction flange, inches of
mercury.
262
Mech.MuslimEngineer.Net
Pb barometric pressure, inches of mercury.
Pv absolute vapor pressure of water at the temperature of test
water, inches of mercury.
H NPSH, in feet of water as quoted by vendor.
h suction velocity head, feet of water at test conditions.
sg specific gravity of water at test temperature.
Thoma has suggested that the dynamic depression, including the velocity
head at the impeller eye, can be expressed as a fraction of the total head
H m or:
NPSH = σ H m (9)
Pat
± hss − h f − hv
γ
σ= (10)
Hm
σ 1 − σ c H m2
=
σ 2 − σ c H m1
263
Q
3
= Const
De N
4
NPSH .De
= Const
Q2
NPSH
= Const
N 2Q 2
264
Mech.MuslimEngineer.Net
N Q
S= (11)
NPSH 3 / 4
Since the two parameters σ and S are likely to have parallel use in
the same field of application, it becomes necessary to state an analytic
expression for the relation between them. This relation which may be
derived from the explicit expressions for S and N s , is:
4/3
N
σ = s (12)
S
σ = Const.N s 4 / 3 (13)
265
266
Mech.MuslimEngineer.Net
sewing too far. This was resulted in large pumps operates at low speed.
Also one of the recommendations usually used to avoid cavitation is the
use of a larger pump than the size that would normally be chosen, one
merely offers a “4 inches” pump instead of “3 inches” or one offers a
double-suction impeller type pump instead of a single-suction, Figure 14.
This practice thereby reduces the required NPSH to a significantly lower
value and in this relatively painless manner pump application engineers
have been falling into a trap lately.
267
268
Mech.MuslimEngineer.Net
in head of 3 percent is usually accepted as evidence that cavitation is
present”. It is important to note that 3 percent head drop point is not
always the start-of-cavitation point, neither is a 1 percent head drop point.
At these head drop points, cavitation has already progressed sufficiently
to affect the flow in the impeller passages.
Nomenclature
* Dimensionless
c Velocity of sound in the fluid, (m/s)
CL Velocity of sound in casing material E / ρ p , (m/s)
CNL Acoustic cavitation pressure, (N/m2)
CV Solid-borne noise acceleration (effective value), (m/s2)
CV* Dimensionless solid-borne noise acceleration, CV* = CV ( u12 / D1 )
269
270
Mech.MuslimEngineer.Net
and to cavitation on the pressure side. All parameters, which influence the
hydrodynamic cavitation intensity, have also an influence on the
cavitation noise.
If the inlet pressure is lowered during a suction test with constant flow
quantity and constant speed of rotation, signals can be measured using a
pressure transducer at the pump inlet. These signals are graphed in
Figure 15.
271
The lowering of the noise level in the range B can be explained by the
following mechanisms:
- When the bubble volume at the impeller inlet is large, a great part
of the sound in these zones is absorbed with two phase flow: the
bubbles can, for example, implode within the impeller ducts, whereby
areas with two-phase flow between implosion zone and pressure
transducer form an effective barrier.
- At low inlet pressure, air can be precipitated in the pump inlet (or
in the upstream throttle valve), so that a two-phase flow is again
present, which makes the cavitation noise measurement completely
unusable.
272
Mech.MuslimEngineer.Net
Basing on the fact that the noise in the range B sinks under the
background noise level, the conclusion can be drawn that the main source
of the background noise level is not in the inlet but in the impeller, i.e. in
the unstable blade forces. Were the background noise level to be induced
by flow processes in the inlet, it could not be blocked off from the sound
pressure detectors by the two-phase flow in the impeller.
CNL R =
Z2, R
Z2
(
NL2 − NL20 )
2 CNL
CNL* =
ρ u12
273
Fig. 16 Influence of the air content and of the frequency range on the cavitation
noise (high-pass filter set in each case at the given frequency)
274
Mech.MuslimEngineer.Net
Moreover, tests show that the gas content has a significant influence
(particularly at frequencies of less than 10 kHz). A bubble implosion
gives rise to a sharp pressure peak and lower pressure peaks in the
rebound phase.
1 p − p∞
fm =
2 Ro ρ
The bubble radius Ro and the local pressure p are unknown, but
proportionality can be derived using the following equation [11]:
275
The frequency with the highest energy density is thus strongly dependent
on the cavitation coefficient. As well as being indicated by the
observations of Figure 16, this is also confirmed in the literature [10, 15].
Summarizing it can be established:
1.463
Fcor I ac
PER = 8,8 .10 −8
FMat IR
CNL2R
I ac =
ρc
276
Mech.MuslimEngineer.Net
the erosion rate can be estimated from the fluid-borne noise
measurements. There is, however, a basic difficulty in that the inlet
pressure can generally not be varied during measurements carried out in
pumping installation. Hence it is also not possible to determine the
background noise level. For approximation purposes, the total noise level
can be utilized. This, however, implies the presence of developed
cavitation.
In the case of large impeller inlet angles and unsuitable shaping of the
runner blade inlet edges, a coefficient of 2 is no rarity. At an installation
cavitation coefficient value of about 0.5, cavitation is generally to be
expected. If the coefficient is between 0.2 and 0.3, it is possible that the
range B (Figure 15) has already been attained and thus the measured
acoustic pressure is no longer unquestionably a measure of the danger of
cavitation.
277
If the noise is measured at a defined spatial point, the sound power of the
source can be calculated basing on the acoustic characteristics of the
space concerned. By estimations of the acoustic efficiency (ratio of the
radiated acoustic power to the mechanical power of the implosion
process), the implosion pressure can be calculated approximately. Then
the material stressing can be compared with the fatigue strength of the
material.
Using the data given in [27], a surprisingly good agreement with Figure
18 is found.
278
Mech.MuslimEngineer.Net
Fig. 18 Specific erosion power as a function of implosion pressure
279
280
Mech.MuslimEngineer.Net
into the casing wall and the transmission of solid-borne noise in the
casing are dependent on the geometry and on the material (Figure 21).
However, for a given system, there is a correlation between erosion and
solid-borne noise.
281
ρ pC L R hD η
NL = CVρ h 1 + . 1 +
3π ρ C L h η 21
ρ Cσ rad
η 21 =
ρ phω
When using this method, the specific erosion power as function of the
solid-borne sound (Figure 20) is transferred by calculation into the
specific erosion power as a function of the fluid-borne sound (Figure 22),
it is clear that the groups of data from the solid-borne measurements now
group themselves round the erosion correlation as a function of the fluid-
borne sound. This is particularly to be noted for the tests carried out on
airfoils [4] (for these tests, it was not considered appropriate to normalize
the CNL to the reference blade number).
282
Mech.MuslimEngineer.Net
Fig. 22 Specific erosion power as a function of the fluid-borne noise
Since the acoustic models utilized are approximate and sound absorption
as well as sound transmission are dependent on many parameters
including pipeline connections and foundations, the degree of uncertainty
and the amount of scatter associated with this calculation method are
considerable. In addition, the measured data have been obtained from
283
The use of an acoustic model for the calculation of the fluid-borne noise
from the solid-borne noise also has the advantage that the implosion
pressure can be estimated and thus a value of absolute material loading is
derived.
PER ~ CNLx ~ CV y
284
Mech.MuslimEngineer.Net
measured previously) without the necessity of having to calculate to
absolute values.
285
286
Mech.MuslimEngineer.Net
Fig. 25 Pump impeller showing cavitation erosion
287
288
Mech.MuslimEngineer.Net
Fig. 27 AEA facility for cavitation investigation
289
7- The impeller velocities, impeller vane load, and head per stage
should be low for minimum suction head. All these factors lead to
a bigger pump operated at a low speed, and possibly to location of
the operating point to the left of the b.e.p.
290
Mech.MuslimEngineer.Net
Fig. 28 Effect of suction nozzle on NPSH requirements. 3-in. pump, 3500 rpm.
Full lines: straight end suction;
Dotted lines: flat elbow suction nozzle
291
292
Mech.MuslimEngineer.Net
32. TURTON, R.K., “Principles of Turbomachinery”, 2nd edition, Chapman & Hall,
London, 1995.
33. VARGA, J.J., SEBESTYEN, G. and FAY, A., “Detection of Cavitation by
Acoustic and Vibration-Measurement Methods”, La Houille Blanche, 1969,
No. 2, pp. 137-149.
34. WISLICENUS, G.F., “Fluid Mechanics of Turbomachinery”, Dover Publication
Inc., New York, 1965.
293
CHAPTER I
BASIC THEORY
The velocity of the fluid particle moving through an impeller channel must
satisfy the vectorial relationship.
→ → →
C = W +U
where:
πDN
2 ω.r = C
60 W Cr
β
where:
α
292
Mech.MuslimEngineer.Net
The power generated from or by turbomachine rotor is defined by:
→
Power = ∫∫ r .ω C . d Q m
U 2 C u2 − U 1 C u1
The theoretical head generated Ho =
g
(U 2
1 −U 2
2 ) + (W
2
2
− W12 )
The degree of reaction σ =
(C
2
1 − C22 ) + (U 12 − U 22 ) + (W22 − W12 )
η∞ = C u / Cu
2 2
γ .QH o
η mech. = / BHP (Pump)
Const.
or
γ .QH o
η mech. = BHP / (Turbine)
Const .
η hyd . = H o / H av
293
Derive an expression for the head created by the centrifugal forces only in a
rotating element containing fluid such as pump impeller.
Solution
F = M.a = M r ω2
dF = dM r ω2
i.e.
∂P = ρω 2 r∂ r (2)
∂P
P.r.dr.dφ − P + dz r.d r.dφ − ρ g.r.d r.dφ .dz = 0
∂ z
or
∂P
= −ρ g (3)
∂z
∂P ∂P
dP = dr + dz (4)
∂r ∂z
294
Mech.MuslimEngineer.Net
From equations (2), (3) and (4):
dP = ρω 2 r.dr − ρg.dz
Integrate from r1 to r2
ρω 2
P2 − P1 =
2
(r
2
2
)
− r12 + ρ g (Z 1 − Z 2 )
P2 − P1 represents the increment in total pressure due to centrifugal forces only. The
last formula is important to the study of turbomachine. If Bernoulli's equation is
applied, ( P2 − P1 ) equals the difference in potential kinetic energy.
Example (2)
Starting from Bernoulli's equation for a steady flow, determine the head
developed by a turbomachine.
Solution
Consider also that the axis of reference rotate with the rotating element in such way,
that we can consider that there is no power gained from the machine and vice-versa
no power induced to the machine.
P1 W 12 P2 W 22
+ + V = + +V (1)
ρg 2g 1 ρg 2g 2
i.e.
P1 P2 W22 − W12 ω 12 r12 − ω 22 r22
− = + + Z 2 − Z1
ρg ρg 2g 2g (2)
U = ωr
i.e.
P C2 P1 C 12
Ho = 2 + 2 + Z2 − + + Z1 (3)
ρg 2 g ρg 2 g
295
which is the general equation for the power delivered by the turbomachine.
W 2 = U 2 + C 2 − 2UC cos α 2
Example (3)
A pump impeller rotating at 1400 rpm has an outside radius of 21 cm, the vane
outlet angle β2 is 158° and the radial velocity at the outlet Cr2 is 4 m/s. Assuming
radial flow at inlet, draw the theoretical outlet velocity diagram and calculate the
various velocities and angles. What is the theoretical head Ho assuming that the
circulatory flow coefficient η∞ = 1.
Solution
β 2 = 158°
2π N
ω= = 146.6 rad/s
60
U 2 = ω r2 = 146.6 × 0.21 = 30.8 m/s
C r2 4
W2 = = = 10 .66 m/s
sin (180 − β 2 ) sin 22 °
C r2
Cu2 = U 2 −
tan (180 − β 2 )
4
= 30.8 − = 20.9 m/s
tan 22°
2 2
C 2 = C r 2 + Cu 2
i.e. C2 = (4)2 + (20.9)2 = 21.3 m/s
C r2
α 2 = tan −1 = 10 ° 5 0 ′
C u2
U 2 C u 2 − U 1 C u1 U 2 C u2
Ho = =
g g
296
Mech.MuslimEngineer.Net
30 .8 × 20 .9
Ho = = 65 .6 m
9 .81
Example (4)
Assuming that the mechanical efficiency ηmech = 0.95, the hydraulic (or
manometric) efficiency ηman = 0.8, find the required power to drive the pump. Also,
calculate the water horsepower. The flow is 30 lit/s.
Solution
C u2
η∞ =
C u2
U 2 C u2 30.8 × 16.7
Ho = = = 52.5 m
g 9.81
γ Q Ho
η mech =
75 B.H .P.
i.e.
1000 × 0.030 × 52.5
B.H .P. = = 22.1 hp
75 × 0.95
Example (5)
A turbine rotates at 150 rpm and discharges 0.8 m3/s. The radial velocity at
inlet Cr1 = 2 m/s and equals 6 m/s at exit. The physical data are: r1 = 0.5 m, r2 = 0.2 m,
297
a) The head, power delivered by the turbine, (no draft tube are used) neglect
circulatory flow and hydraulic losses.
b) The pressure head at entrance.
c) The degree of reaction.
Solution
a)
2π N
ω= = 15.7 rad/s
60
U 1 = ω r1 = 7.85 m/s
U 2 = ω r2 = 3.14 m/s
C r2
Cu 2 = −U 2
tan (180 − β 2 )
= 6 − 3.14 = 2.86 m/s
C U − C u 2U 2
H o = u1 1 = 5.05 m
g
γ Q Ho
Power = = 53.9 hp
75
P1 C2 P C2
+ Z 1 + 1 − H loss − H o = 2 + Z 2 + 2 (Z 1 = Z 2 )
ρg 2g ρg 2g
2 2
C12 = C u1 + C r1 = 57.65 m 2 /s 2
2 2
C 22 = C u 2 + C r 2 = 44.18 m 2 /s 2
P1 57.65 44.18
+ − 2 − 5.05 = 6 +
ρ g 2 (9.81) 2 (9.81)
P1 P
+ Z 1 − 2 + Z 2
ρg ρg
c) The degree of reaction =
Ho
In this case, we must subtract the entry losses from P1 / ρg to obtain the real
pressure head at entrance.
(12.36 − 2 ) − 6
σ = = 0.864
5.05
298
Mech.MuslimEngineer.Net
The turbine is reaction type.
N.B.: We can find the same value for σ using the following form of σ:
σ=
(U 1− U 22 ) + (W22 − W12 )
2
W12 = U 12 + C12 − 2U 1C u1
= 61.62 + 57.65 − 2 × 7.46 × 78.5 = 2.146
W22 = U 22 + C 22 − 2U 2 C u 2 = 36.08
W2 = 6 m/s
σ=
(61.62 − 9.86) + (36 − 2.1)
(61.62 − 9.86 ) + (36 − 2.1) + (57.6 − 44.18)
σ = 0.864
Example (6)
It is desired to pump 100 lit/s of water to the top of a cooling tower, the
required manometric head was 19 meters. A radial centrifugal pump type was selected
to give the required manometric head. The pump technical data are as following:
impeller outside diameter D2 = 21 cm, vane outlet angle β2 = 158°, impeller tip width
= 5 cm and rotating speed N = 1750 rpm. Find the following:
299
U 1 C u1 − U 2 C u 2 U 2 C u2
a) Ho = =
g g
2π N
ω= = 183.25 rad/s
60
U 2 = ω r2 = 183.25 * 0.21 / 2 = 19.24 m/s
C r2
Cu2 = U2 −
tan (180 − β 2 )
Q 0.1
Cr2 = = = 3.03 m/s
π D2 b2 π × 0.21 × 0.05
19.24 × 11.74
Ho = = 23 m
9.81
C u2
b) η∞ =
C u2
β 2′ = 162° 54′
δ = β 2′ − β 2 = 4° 54′
300
Mech.MuslimEngineer.Net
Problems
1- Determine the expression of the degree of reaction and show that the propeller
pump is a reaction machine.
2- A turbine rotates at 100 rpm and discharges 0.84 m3/s, the hydraulic efficiency
under these conditions is 75.5 %. The physical data are: r1 = 0.46 m, r2 = 0.22 m,
α1 = 15°, β2 = 135°, A1 = 0.12 m2, A2 = 0.078 m2. Neglect the circulatory flow
coefficient and take the mechanical efficiency ηmech = 0.95. Determine the power
delivered by the turbine.
( 39 HP )
3- A centrifugal pump rotates at 600 rpm. The following data are taken: r1 = 5.08 cm,
r2 = 20.3 cm, radial area A1 = 769 cm2, radial area A2 = 295 cm2, β1 = 135°,
β2 = 120° and assume radial flow at entrance to blades. Neglecting friction,
calculate the relative velocities at entrance and exit and the power transmitted to
the water.
4- A reaction turbine is working under a head of 25 meters running at 300 rpm. The
velocity of the periphery of the wheel is 30 m/s and the radial velocity Cr1 is 4 m/s.
If the hydraulic losses are 20 % of the available head and the discharge is radial,
find:
i) the inlet angle β1 and guide blade angle at inlet α1.
ii) the rotor diameter.
(170° 18′, 31° 27′, 1.9 m)
301
8- An inward flow reaction turbine discharges radially and the velocity of flow is
constant and equal to the velocity of discharge from the suction tube. Show that the
hydraulic efficiency can be expressed by:
1
ηh =
0.5 tan 2 α 1
1−
tan α 1
− 1
tan β 1
******
302
Mech.MuslimEngineer.Net
CHAPTER II
DIMENSIONAL ANALYSIS
L model
a- Geometric similitude: = L ratio
L prototype
303
N Q
- Specific Speed: Universal Ns =
( gH )3 / 4
- Pressure Coefficient φ = H / (U22 / 2g)
- Flow Coefficient ψ = (Q / A) / U
D ( gH )
1/ 4
- Specific Diameter Dsp =
Q
304
Mech.MuslimEngineer.Net
Solved Examples
Example (1)
N Q
Show that the specific speed expression N s = depends upon the
( gH )3/ 4
proportions of the wheel rather than the operating speed, consider pump impeller and
radial flow.
Solution
N ( Q)
1/2
Ns = (1)
( gH ) 3/ 4
Q = A1 C r1 = π . D1b1 C r1
= π . D1b1U 1 tan(180 − β1 )
C u2 U 2
H= radial flow at inlet.
g
U2 Cr2
= U 2 −
g tan(180 − β2 )
Ns =
[
N π 2 D12 tan (180 − β 1 ) b1 N ]1/ 2
3/ 4
D tan (180 − β 1 )
2
b1
(π 2 2
D N )
2 3/ 4
1 − 1
tan (180 − β 2 )
2
D 2 b2
305
Example (2)
Solution
Q Q
π1 = =
ND 3 1 ND 3 2
0.240 0.01
= N = 1350 rpm
1200 (3 )
3
N
gH gH
π2 = 2 2
= 2 2
N D 1 N D 2
30 H2
= H 2 = 3.75 m
(1200) ( 3)
2 2
(1350) 2
γ Q1 H 1 1000 × 0.240 × 30
P1 = = = 120 hp
75 η 75 × 0.8
0 . 25
1 − η1 D 2
=
1 − η 2 D1
0 . 25
1 − 0.80 1
= η2 = 0.736
1 − η2 3
306
Mech.MuslimEngineer.Net
Example (3)
Solution
2πN
i) ω = = 11.5 rad/s
60
U 1 = ω R1 = 25.9 m/s
C r1 2
C u1 = = = 7.46 m/s
tan α 1 tan 15
Cu U 1 − C u2 U 2
Ho = 1
g
7.46 × 25.9
Ho = = 19.7 m
9.81
1000 × 11 × 19.7
Power = = 2889 hp
75
N B. H . P.
Ns = = 142
H 5/4
P1 C12 P2 C 22
ii) + + Z1 = + + Z 2 + hloss + ho
ρg 2 g ρg 2 g
2 2
C12 = C u1 + C r1 = 7.7 m/s
C 2 = C r 2 = 9.5 m/s
P1
+
(7 .7 )
2
= 6+
(9 .5 )
2
+ 2 + 19 .7
ρ g 2 * 9.81 2 * 9.81
307
H H
=
N D22
1 N D2 2
2
19.7 5 .5 × 0 .8
= ∴ N = 780 rpm
(110) (4.5)2 2
N 2 (0.3)
2
Q Q
=
N. D3 1 N. D3 2
11 Q
= ∴ Q = 0.025 m 3 /s
110(4.5) 780(0.3)
3 3
N B. H . P.
Ns = = 142
H 5/ 4
Example (4)
A turbine develops 144 HP running at 100 rpm under a head of 7.7 meters.
What power would be developed under a head of 11 meters? At what speed should
the turbine run?
Solution
H H
π2 = =
N 2D2 1 N 2D2 2
H H
For the same machine: =
N2 1 N2 2
7 .7 11
2 = N 2 = 120 rpm
(100 ) N 22
P α ρgQH Q α ND3
i.e. P α ND3H
P1 N H
i.e. = 1 1
P2 N 2 H2
308
Mech.MuslimEngineer.Net
144 100 × 7 .7
=
P2 120 × 11
∴ P2 = 247 HP
Example (5)
Solution
γ QH P1 Q1 H 1
P= =
Const . P2 Q2 H 2
Q1 P1 H 2 48,000 11
i.e. = = = 1.05
Q2 P2 H 1 36,000 14
Q Q
=
N. D3 1 N. D3 2
Q2
N 2 D 23 = N 1 D13
Q1
= 0 . 95 × 82 × (7 ) = 26720
3
H H
=
N D2 2
1 N D2
2
2
H 2 2 2 11
N1 D1 = (82 ) (7 ) = 258874
2 2
N 22 D22 =
H1 14
N 2 D2 = 509
26720
D2 = = 7.25 m
509
309
310
Mech.MuslimEngineer.Net
Problems
1- An impulse wheel at best speed produced 125 HP under a head of 85 m. By what
percent should the speed be increased for a 90 m head? Assuming equal
efficiencies, what power would result?
(2.9 %, 136.1 hp)
2- Assuming the power delivered to a pump is a function of the specific weight of the
fluid, the flow in lit per minute, and the head delivered, establish an equation by
dimensional analysis [power equation].
3- A centrifugal pump discharged 65 lit/min against a head of 17 meters when the
speed was 1500 rpm. The diameter of the impeller was 30 cm and the brake
horsepower was 6 HP. A geometrically similar pump 38 cm in diameter is run at
1750 rpm. Assuming equal efficiencies, what head will be developed? How much
water will be pumped and what brake horsepower will be required?
(37.1 m, 154 1it/min, 31.7 hp)
4- A turbine model, built to a scale of 1:5 was found to develop 4.5 BHP at a speed of
400 rpm under a head of 2 meters. Assuming the overall efficiency of the full-size
turbine = 0.8, find the speed and the power of the full-size turbine under a head of
9 meters.
(169.7 rpm, 1073.8 hp)
5- Seawater of specific gravity 1.03 is to be circulated through condensers by a
propeller pump 120 cm in diameter, it is found that a scale model of the pump
25 cm in diameter gives its best efficiency when pumping 97 lit/s of fresh water
against a head of 4 meters at a speed of 2060 rpm. What should be the speed of the
full size pump to deliver 90 tons per minute and what pressure difference would it
generate?
(279.6 rpm, 17177 Pa)
6- A centrifugal pump of diameter 30 cm (throat diameter) running at 1450 rpm
delivers 74 lit/s of water against a head of H = 46 m at optimum efficiency of 81
%. Investigate the operation of the pump handling:
(a) oil of sp. gr. 0.86, viscosity 0.857*10-5 m2/s,
(b) fuel oil of sp. gr. 0.94, viscosity 0.232*10-5 m2/s.
7- The same model pump of Problem (6) is now built to a scale 2.5 times larger and
used to pump oil (b) of the previous problem. Find at a speed of 1450 rpm head
capacity, power, and efficiency at the optimum point.
8- A centrifugal pump of impeller 25 cm diameter revolving at 1450 rpm delivers at
its best efficiency point 29 lit/s of water (ν = 0.9*10-6 m2/s) against a head of 20
meters at an efficiency of 80 percent. A geometrically similar pump is to be built to
produce a head of 25 m when handling oil of specific gravity 0.89 and kinematic
viscosity 4.5*10-5 m2/s at the best efficiency point, at the same speed, what will the
diameter of the pump have to be and what will be the quantity delivered, power
consumption, and efficiency at the same optimum operating point?
(Mans. Univ., B.Sc., 1978)
******
311
CASCADE MECHANICS
312
Mech.MuslimEngineer.Net
The head-capacity curve:
H = B U 2 − AU Q
l
A = (tan α 1 − tan α 2 )
t
B=1
∆ Po
Total pressure loss coefficient ζ =
1
ρCa 2
2
∆P
Pressure rise coefficient Cp =
1
ρ Ca 2
2
Tangential force coefficient Cf = 2 ( tan α1 − tan α2 )
Mean angle tan αm = ½ ( tan α1 + tan α2 )
Lift coefficient CL = 2 ( t/l ) cos αm ( tan α1 − tan α2 ) − CD tan αm
Drag coefficient CD = ζ ( t/l ) cos3 αm
Discharge coefficient ψ = Ca / U
Cascade efficiency η = 1 − ζ ψ/Cf
Deviation from ideal direction δ = α2 − α2'
Nominal deviation δ* = m θ' ( t/l )n
313
A compressor cascade has a space-chord ratio of unity and blade inlet angle of
45°, stagger angle ε = 28°, using the NACA results, Fig (3.5), find the loss coefficient
ζ, Cf, the drag coefficient, the lift coefficient and pressure rise coefficient. Also, find
the pressure rise across the cascade for ρ = 9.5 x 10-6 kg/m3, Ca = 47 m/s. If the above
cascade is for an axial air compressor of tangential velocity U = 150 m/s, find its
efficiency and construct the velocity diagrams (Adapted from Csanady).
Solution
It is easy to notice that the term CD tan αm can be easily neglected since CD is small
and equals zero for ideal lift.
314
Mech.MuslimEngineer.Net
α1
α2 W2 Ca = 47 m / s
C1 α1 W1
C2
∆Cu
U
Example (2)
A compressor cascade has a space-chord ratio of unity and blade inlet and
outlet angles of 45° and 15°, respectively. If the blade camber line is a circular arc
(i.e. a/l = 0.5) where a is the distance of maximum camber from leading edge. The
cascade is designed to operate at Howell’s nominal condition. Determine the fluid
deflection, incidence and ideal lift coefficient at the design point (Adapted from
Dixon).
Solution
315
Experimental compressor cascade results suggest that the stalling lift coefficient of a
3
C1
cascade blade may be expressed as C L = 18 . where C1 and C2 are the entry
C2
and exit velocities. Find the stalling inlet angle for a compressor cascade of
space/chord ratio of unity if the outlet air angle is 30°. If the total pressure loss
coefficient at zero incidence i = 0 (corresponding to minimum pressure loss) ζ = 0.02,
find the corresponding pressure coefficient.
Solution
NOTE: Stall point is arbitrarily specified as the incidence at which the total pressure
loss is twice the minimum loss in total pressure.
C1 cos α1 = C2 cos α2 = Ca
C1 cos α 2
i.e. =
C 2 cos α 1
3
C
C L = 1.8 2 (2)
C1
3
cos α 1
C L = 1.8 = 1.51
cos α 2
Assume α1 = 48°
tan αm = 0.84 and αm = 40° 06'
From equation (1) CL = 0.82 and from equation (2) CL = 0.83 which is acceptable.
316
Mech.MuslimEngineer.Net
Example (4)
The cascade shown in Figure 3.5, is to be used for an axial flow propeller pump
impeller. The cascade has a space-chord ratio of unity and blade inlet angle of 45° and
stagger angle ε = 28°, find the pressure rise coefficient if the net positive suction head
(total pressure at inlet minus vapor pressure) is 7 meters of water. Specify the
maximum allowable axial velocity with regard to cavitation danger. The maximum
surface velocity may be taken as 25 percent above inlet velocity C1. Also, find Euler
head and the efficiency.
Solution
317
318
Mech.MuslimEngineer.Net
Problems
1- For the cascade data shown in Example 1, choose the optimum ( ψ / Cf a minimum)
operating point, and find the corresponding pressure rise and tangential force per
centimeter height of blade, if axial velocity of 80 m/s, density of 9.5*10-6 kg/m3,
blade spacing 10 cm.
(0.022 Pa, 3.04 x 10-5 N/cm)
2- Explore one or two other possible choices of stagger angle in Example 1 with a
view to ascertaining the most favorable design possibility with regard to efficiency.
3- A compressor cascade has the following data of t/l = 1.0, α1' = 45°, α2' = 15°. The
nominal conditions were θ* = 24° and i* = 2.3, the off-design performance of this
cascade is required at an incidence i = 3.8° ( referring to Fig. 3.8 ).
4- A compressor cascade has a space-chord ratio of unity, blade angle of 45° and
stagger angle ε = 28°, find the loss coefficient, the drag coefficient, the lift
coefficient, and the pressure rise coefficient. Also, find the pressure rise across the
cascade for ρ = 0.00001 kg/m3, Ca = 47 m/s (using NACA data: θ = α1 − α2 = 30°
and CD cos2 α1 / cos2 αm = 0.0065).
( B.Sc., Mans. Univ., 1982 )
******
319
INCOMPRESSIBLE FLOW
TURBOMACHINES
Theoretical head:
(U C cosα
1 1 1 − U 2 C2 cosα 2 )
Ho =
g
Assuming radial inlet and substitute C 2 cos α 2 by U 2 + (C r 2 / tan β 2 )
U 22 U 2 Cr2
Ho = + and Q = π . D2 b2 C r2
g g tan β2
Then:
320
Mech.MuslimEngineer.Net
U 22 U2 Q
Ho = +
g g π . D2 b2 tan β2
H a (measured head )
ηm = Manometric efficiency =
Ho
γ QH o
ηmech. = Mechanical efficiency =
Const. × BHP
ηoverall = ηmech. . ηm
η mech ηm
ηoverall
ψ = Discharge coefficient =
Q/A
U
or Q D 2 . 2 gH( )
N Q
Universal specific speed N s = in SI units
( gH ) 3/ 4
The net positive suction head NPSH
Pat
NPSH = ± hss − h f − hv
ρg
where:
hss : Static suction head
h f : Friction loss
hv : Vapour pressure
NPSH
Thoma cavitation factor σ =
Hm
321
- Leakage Calculation
Q L = C v A p 2 gh
where:
C v = Velocity coefficient
A p = Clearance area
U 22 − U 12
h = Head across the orifice = 0.75
2g
- Disk Friction
16 T
- Shaft diameter Ds = 3
π Ss
where:
T is the torque,
Ss is the allowable shear stress.
- Eye Diameter Do
4 Q
Do2 = + D H2
π Uo
Q
b1 =
π .D1 .ε .C r1
U2
φ= to find D2 and assume β2 is slightly lower than β1
2 gH
322
Mech.MuslimEngineer.Net
D2 + D1
The number of blades Z n = 6.5 sin (180 − βm )
D2 − D1
where βm = ( β1 + β2 ) 2
π D − tZn
ε=
πD
323
Test results on a single stage single suction centrifugal mixed flow type pump
operating at 375 rpm designed to deliver 2.4 m3/min of water are given in curve form
as follows:
Flow (m3/min) 0 4 8 12 16 20 24 28 34
Head (m) 12 11.3 10.6 9.9 9.1 8.2 7.2 5.9 4.9
η (%) 0 23 45 62 75.5 84 88 85 81
Solution
γ QH 1000 Q H
(a) The horsepower = =
η . Const 60 * 75 * η
Thus the power will be:
(b) If the fluid density changes, this should not affect the head-discharge curve, so the
head-discharge will remain constant, the delivery pressure will change only, (P =
ρgh).
324
Mech.MuslimEngineer.Net
The power should be reduced as it could be seen from the power relation.
γ QH
Power =
η.Const .
The efficiency is a relation between the head delivered and the power consumed
which practically will alter (The hydraulic losses will change, as it is a function
of the density and the viscosity of water), but it could be assumed constant.
(c) To draw the performance of the pump at different speeds, one can make use of the
affinity laws as follows:
Q1 N 1 N2 500
= i.e. Q2 = Q1 ∴ Q2 = Q1
Q2 N 2 N1 375
2
H1 N 12 N 500
2
= i.e. H 2 = 2 H 1 ∴H2 = H1
H 2 N 22 N1 375
Example (2)
is interposed in a pipe in which the loss in suction pipe = 1/3 loss in delivery line
when the static lift was 6 m, the maximum possible discharge was found to be
300 lit/s. Find the highest possible position of pump above sump level, given that
hv = 0.3 m, hat = 10.3 m, cavitation factor = 0.5. Neglect the effect of kinetic energy
loss in pipe.
Solution
H sv
σ= where Hsv is the NPSH
Hm
i.e. Hsv = 0.5 Hm
The manometric head at Q = 300 lit/s could be taken from the performance curve,
i.e. Hm = 14 m
i.e. Hsv = 14 x 0.5 = 7 m
325
1
The losses in suction pipe = ×6 = 2 m
3
From equation (1) hss = 1 m
i.e. the highest possible position of the pump above the sump level would be 1 meter.
Example (3)
326
Mech.MuslimEngineer.Net
Solution
H. G.
Hst
Hm d
hss
c
Hms
H. G.
hf
327
Example (4)
Solution
a) Q = AVo =
π
4
(D o
2
− DH
2
)V
o
328
Mech.MuslimEngineer.Net
Q
C r1 =
π D1 b1ε 1
40
= = 4.6 m/s
1000 × π × 0.133 × 0.023 × 0.9
4 .6
tan (180 − β 1 ) = = 0.3770
π × 0.133 × 1750
60
i.e. β 1 = 159° 18'
U2
b) φ=
2 gH
π D2 N
U2 = = 1.08 2 × 9.81 × 19 = 20.85 m/s
60
D2 = 0.228 m or 22.8 cm
C r2
ψ= = 0.21
2 gH
C r2 = 4.05 m/s
C r2 4.05
Cu2 = U 2 − = 20.85 − = 10.9 m/s
tan (180 − β 2 ) tan 22°
U 2 Cu2
20.85 × 10.9
Ho = = = 23.2 m
g 9.81
γ Q H o 1000 × 40 × 23.2
Power = = = 12.37 hp
Const. 1000 × 75
19
η= = 0.815
23.2
C u2 gH 2 U 2
c) η∞ = = . = 0.815 (Neglecting all mechanical losses and friction
C u2 U 2 gH o
losses).
Thus, we can say that η ∞ = η h only if we neglect the mechanical losses and
hydraulic losses in the pump. This is a theoretical case which does not exist
actually.
C u 2 = 8.924 m/s
C r2
tan(180 − β 2′ ) =
U 2 − C u2
i.e. β 2′ = 161°15'
δ = β 2′ − β 2 = 161°15'−158° = 3°15'
329
NPSH − 0 .7
σ = = = − 0 .0368
Hm 19
Example (5)
Competitive bids from three companies for a water circulating pump to handle
2000 lit/min of water with a head of 50 meters and operates 8 hours/day, 300
day/year, are as follows:
Company A B C
Price in £ 230 410 260
Efficiency % 71 77 73
Estimated life in years 10 15 10
The annual fixed charges are 20 % of the initial price and the power costs 2 piasters
per kW.hr, which pump could be selected?
Solution
γ QH
Water horsepower =
Const.
1000 × 2000 × 50
= = 22.2 hp
1000 × 60 × 75
Company A B C
B.H.P. 31.2 28.8 30.4
B.H.P.-hr/year 74,800 69,100 73,000
Power cost/year 1115 1030 1080
Fixed charges/year 46 82 52
Depreciation 23 27.3 26
Total annual cost 1184 1139 1158
It may be seen from the above example that the total annual cost are about the
same for the bids A, B, and C, despite the differences between their initial prices, but
if the operation time is increased, the bid B will be in a better position.
330
Mech.MuslimEngineer.Net
Example (6)
Solution
Example (7)
A centrifugal pump, which runs at 16.6 rev/s, is mounted so that its centre is
2.4 m above the water level in the suction sump. It delivers water to a point 19 m
above its centre; the friction loss in the suction pipe is 68 Q2 meter and that in the
delivery pipe is 650 Q2 meter where Q in m3/s is the rate of flow. The impeller of the
pump is 350 mm diameter, and the width of the blade passages at outlet is 18 mm.
The blades themselves occupy 5 % of the circumference and are backward facing at
35° to the tangent. At inlet, the flow is radial and the radial component of velocity
remains unchanged through the impeller. Assuming that 50 % of the velocity head of
the water leaving the impeller is converted to pressure head in the volute, and that
friction and shock losses in the pump, the velocity heads in the suction and delivery
pipes are negligible, calculate the rate of flow and the manometric efficiency of the
pump.
Solution
N =16.6 rev/s
Hs = 2.4 m, Hd = 19 m
Hfs = 68 Q2, Hfd = 650 Q2 19 m
D2 = 35 cm
b2 = 1.8 cm
2.4 m
ε2 = 0.95
β2 = 145°
Cu1 = 0
Cr1 = Cr2
331
P1 C12 P C2
+ + Z 1 + H o = 2 + 2 + Z 2 + hL ( hL = 0 , Z 1 = Z 2 )
γ 2g γ 2g
P2 − P1 C u 2 U 2 C 2 − C1
2 2
= − (6)
γ g 2g
Total pressure rise = Pressure rise in impeller + Gain of pressure in the volute
P − P1 C2
Hm = 2 + 0 .5 2
γ 2g
2
C U C 0.5 C 22
H m = u2 2 + 1 −
g 2g 2g
332
Mech.MuslimEngineer.Net
Problems
1- A centrifugal pump lifts water against a static head of 37 m of which 4 m is the
static suction lift. The suction and delivery pipes are both 15 cm diameter, the
head loss in the suction pipe is 2.2 m and in the delivery pipe is 7.25 m. The
impeller is 38 cm diameter, and 2.5 cm wide at the exit. It revolves at 1200 rpm
and its effective exit blade angle is 145°. If the manometric efficiency of the pump
is 82 % and the overall efficiency is 70 %, what discharge would you expect, and
what horsepower would be needed to drive the pump? What would be the pressure
head indicated at the suction and delivery branches of the pump?
(12.8 lit/s, 11.3 hp)
2- If the static lift for a centrifugal pump is hs in meter, speed of rotation N, rpm, and
the exit diameter of the impeller is D meters, show that N = 84.6 hs / D for the
speed at which the pumping begins.
(Mans. Univ., B.Sc., 1980)
333
7- A centrifugal pump, having 4 stages in parallel, delivers 180 lit/s of liquid against a
head of 25 m, the diameter of impellers D2 = 23 cm and the speed is 1700 rpm. A
pump is to be made up with a number of identical stages in series to run at
1250 rpm and to deliver 240 lit/s against a head of 250 m. Find the diameter of
the impellers and the number of stages required if these impellers are of similar
construction to those of the first pump.
(5 stages, 44.5 cm)
9- For a specified duty, 55 tons of water per minute are required to be pumped against
a head of 3 meters, it is desired to compare the probable speeds, overall
dimensions of a centrifugal pump and a propeller pump. Give an estimate of these
values assuming that the maximum diameter of the centrifugal pump casing is 2.5
times the impeller diameter and that the maximum diameter of the propeller pump
casing is 1.25 times the propeller diameter.
(lit/s) 0 4 8 12 16 20 24
Pump I m (m) 50.0 51.8 50.8 48.0 47.5 32.5 18.3
Pump II m (m) 46.5 45.9 44.2 40.3 40.3 26 17
Both pumps are installed together and are required to pump water through a pipe
15 cm diameter having f = 0.02. Calculate the heads under which pumps are
working and discharges in lit/s pumped by them if:
(a) The pumps are connected in series, static lift is 65 m and suction and delivery
pipes are 800 m long.
(b) The pumps are connected in parallel, static lift is 15 m and suction and
delivery pipes are 360 m long.
(39 m, 31 m, 17.8 l/s; 28 m, 21 l/s, 19 l/s)
11- A pump operating continuously (8760 hours a year) has an efficiency of 70
percent and circulates 8000 lit/min against a head of 100 m. A new pump costing
1000 LE and having an efficiency of 82 % can be purchased. If power costs 2
334
Mech.MuslimEngineer.Net
piasters per kW-hr and the annual fixed charges are 20 percent of the initial cost,
would you advise purchasing the pump?
When this pump is used with a long pipeline 50 cm diameter, it gave a maximum
discharge as 150 lit/s with static lift 10 meters. If another pipeline 40 cm has the
same length as the first one is used in parallel with the first line, calculate the
approximate pump discharge and SHP. Neglect the effect of suction pipe K.E.,
(f = 0.03 for all pipes).
(B.Sc., Alex. Univ., 1982)
13- A centrifugal pump is used to deliver water from a water main at point A to the
atmosphere at a height of 23 m above pump level. The suction pipe, which is
horizontal is of length 5 m and diameter 5 cm, the delivery pipe has 3 bends
(C = 0.8) and is of length 45 m and diameter 5 cm. If the water is collected in a
tank of dimensions 3 m x 2 m x 2 m, find the minimum time to fill the tank given
that the water main pressure at point A is 2 kg/cm2, coefficient of friction of both
pipes = 0.03, and the performance of the pump is given by:
Q (lit/s) 0 3 6 9 12
H (m) 22 22.5 20 16 10
ηm (%) 0 40 70 85 65
14- A centrifugal pump having 3 stages and delivers 100 lit/s of liquid against a head
of 30 m, the diameter of impellers are De = 10 cm, D2 = 22 cm, the speed is
1750 rpm. Find the diameter of the balance drum if it has to be connected to the
suction pipe, what would be the drum diameter if the impellers are to be arranged
in a way back to back? (The third impeller is in opposite direction to the 1st and
2nd).
******
335
γ Q cos β
Power = 1 − (UC1 − U 2 )
g 1+ k
where:
k is the friction coefficient
C1 is the jet velocity = C v 2 gh
γQ
The maximum power =
g
(C1
2
/ 4)(1 − cos β )
γQ
Power =
g
(U C
1 u1 − U 2 Cu2 )
Pressure coefficient φ = U 1 / 2 gh
Flow coefficient ψ = C r1 / 2 gh
η hyd = H o / H available
η mech = BHP / (γ Q H o / Const.)
η overall = η hyd .η mech = BHP / WHP
336
Mech.MuslimEngineer.Net
hg
Draft tube efficiency =
V − Ve2
1
2
2g
where:
hg is the actual head gained,
Subscripts 1 and e denotes inlet and exit of draft tube.
Z = Pa / γ − Pv / γ − σ g H
T2 N 2
η=
T1 N 1
337
Solution
Wheel velocity:
U = φ 2 gh = 0.45 2 × 9.81 × 264 = 32.4 m/s
U = π D N / 60
Example (2)
338
Mech.MuslimEngineer.Net
(b) The required number of jets, comment if the ratio of jet diameter to wheel
diameter = 1/10, also calculate the head coefficient φ and the specific speed
Ns .
(Mansoura University, 1977)
Solution
a) BHP = η h × WHP
BHP = 0.8 x 1000 x 3.5 x 290/75 = 10827 HP
γ Q cos β 2
BHP = γ Q H / 75 = 1 −
(C1U 1 − U 12 )
75 g 1+ k
− b ± b 2 − 4c
U=
2a
(U − 48.8) or (U − 25.1) = 0
U = 48.8 m/s or 25.1 m/s
U × 60
D=
π ×N
25.1× 60
i.e. D= = 1.6 m
π × 300
0.0474 × 4
Diameter of jet dj = = 0.14 m = 14 cm
3×π
339
N Q 300 3.5
Ns = = = 0.0239
( gH )
3/ 4
60(9.8 × 290 )
3/ 4
Example (3)
Solution
(a) P = γ Q H / 75
η overall = η mη h = 0.85 × 0.95 = 0.808
P = 1000 × 10 × 57 / 75 = 7600 HP
η overall = BHP / WHP
N BHP
Ns = (metric units)
H 5/4
= 200 6140 / 57 5 / 4 = 100
340
Mech.MuslimEngineer.Net
∴ Z = 7.52 meters
gH gH
(b) =
D2N 2 1 D2N 2 2
2
D2H 5 6
N = N × 12 2 = 40000
2
2 1
2
D2 H 1 0.5 57
N 2 = 649 rpm
Q Q
=
ND 3 1 ND 3 2
3
Q 10 × 649 0.5
Q2 = 3
N 2 D 23 =
ND 1 200 5
649 2.1267
Ns = = 100
65 / 4
Example (4)
The following data refer to a Pelton wheel: the bucket angle β = 165°, the
coefficient of discharge for the nozzle Cv = 0.98, the friction coefficient k = 0.2.
Assume that the optimum speed ratio differs from 0.5 as a result of losses due to
windage and bearing friction, which are proportional to the square of the rotational
speed (take the proportionality constant equal to 0.2), obtain a formula for the
optimum speed ratio, hence calculate it for the above given data and also calculate the
maximum overall efficiency.
( B.Sc., Mansoura University, 1981 )
Solution
1
Constant = J = i.e. W 2 = J W1
1+ k
γQ
Power =
g
(C u1 )
U 1 − C u2 U 2 − Losses in bearing, .. etc.
341
γQ
i.e. Power =
g
(C1 cosα1U 1 − C 2 cosα 2U 2 − PU 2 ) (1)
U 1 = U 2 = U , W2 = J W1 (2)
C1 cos α 1 = C
Cu2 = C 2 cos α 2 = U − W2 cos (180 − β )
W2 = J W1 = J (C − U )
C u 2 = U + J (C − U ) cos β (3)
γQ
Power =
g
{
U [C − (U + J (C − U ) cos β )] − PU 2 }
γQ
Power =
g
[
UC (1 − J cos β ) − U 2 (1 − J cos β ) − PU 2 ]
d Power
= 0 = C (1 − J cos β ) − 2 U (1 − J cos β ) − 2 UP
dU
U 1 − J cos β 1
= 0.5 , J= = 0.913
C 1 − J cos β + P 1 + 0.2
U 1 − 0.913 × cos163°
= 0.5 = 0.4519
C 1.2 − 0.913 × cos163°
B.H .P.
η overall =
W .H .P.
γQ
B.H .P. =
g
[
U C (1 − J cos β ) − U 2 (1 − J cos β ) − PU 2 ]
When U / C = 0.45 = φ i.e. U = 0.45 C
B.H .P.
η overall =
W .H .P.
η max =
1
gH
[
φ C 2 (1 − J cos β ) − φ 2 C 2 (1 − J cos β ) − Pφ 2 C 2 ]
=
C2
gH
[φ (1 − J cos β ) − φ 2 (1 − J cos β ) − Pφ 2 ]
2
=
C v .2 gH
gH
[ (
(1 − J cos β ) φ − φ 2 − Pφ 2 ) ]
[ (
= 2 C v (1 − J cos β ) φ − φ 2 − Pφ 2
2
) ]
[
= 2 × (0.98) 1.88 × 0.2475 − 0.2 × (0.45) = 0.832
2 2
]
342
Mech.MuslimEngineer.Net
Example (5)
Solution
(BHP )1 (BHP )2
a) η= =
γ Q1 H 75 γ Q 2 H 75
950 500
i.e. = i.e. Q2 / Q1 = 0.526
Q1 Q2
(BHP )1 (BHP )2
b) η= =
γ Q1 H 75 γ Q (H − H e ) 75
950 500
= (1)
Q1 H Q 2 (H − H e )
Q 2 (200 − H e ) = 50.6
Q2 / Q1 = 0.807
343
In a vertical shaft inward flow reaction turbine, the sum of the pressure and
kinetic heads at entrance to the spiral casing is 120 m and the vertical distance
between this section and the tail race level is 3 m. The peripheral velocity of the
runner at entry is 30 m/s, the radial velocity of the water is constant at 9 m/s and the
discharge from the runner is without whirl. The estimated hydraulic losses are:
(1) between turbine entrance and exit from the guide vanes 4.8 m.
(2) in the runner 8.8 m.
(3) in the draft tube 790 mm.
(4) kinetic head rejected to the tail race 460 mm.
Calculate the guide vane angle, the runner blade angle at inlet and the pressure head at
entry to the runner.
Solution
U1 = 30 m/s
Cr3 = Cre = 9 m/s
e
P1 C12
+ = 120 m
γ 2g
Losses1-2 = 4.8 m
Losses3-e = 8.8 m
Lossese-4 = 0.79 m
Kinetic energy rejected = 0.46 m
P1 C12
Ho = + + Z − Losses
γ 2g
H o = 120 + 3 − Losses
= 123 − (4.8 + 8.8 + 0.79 + 0.46)
= 108.15 m
U C − U e C ue
H o = 3 u3 (The second term = zero)
g
U .C
108.15 = 3 u 3
g
i.e. C u 3 = 35.36 m/s
C
tan α 3 = r 3 = 0.2545 ∴ α 3 = 14° 17'
Cu 3
Cr3
tan β 3 = = 1.677 ∴ β 3 = 59° 12'
Cu 3 − U 3
2 2
C3 = C u 3 + C r 3 = 36.49 m/s
344
Mech.MuslimEngineer.Net
P1 C12 P3 C 32
+ = + + Losses1-2
γ 2g γ 2g
P (36.49) 2
120 = 3 + + 4.8
γ 2(9.81)
P3 (36.49) 2
Pressure head = 120 − − 4.8
γ 2(9.81)
= 47.3 m
Example (7)
Solution
A1 = 0.2 m2 c
H = 30 m N = 6.25 rev/s
α 1 = 15° β 1 = 75°
D = 0.76 m Ddraft = 45 cm
V 22
hloss = h f + = 0.03 Q 2L
2g d
345
P1 V12 P V2
+ + Z = 2 + 2 + hf
ρ g 2g ρ g 2g
Q2
3.6 + 2
+ L = 10.3 + 0.03 Q 2 L
2 g A1
3.6 + 1.467 + L = 10.3 + 0.03 (0.853) 2 L i.e. L = 5.35 m
346
Mech.MuslimEngineer.Net
Problems
1- For the following two-nozzles Pelton wheel, calculate its maximum efficiency and
output. Pitch circle diameter of runner is 2.5 m, nozzle diameter is 0.15 m, head
on the nozzle is 1600 m, Cv for nozzle is 0.98, bucket angle at exit is 173°,
horsepower lost in disk friction is 6.6 x 10-6 N3, horsepower lost in bearing
friction is 2.1 x 10-3 N2, where N is the speed of wheel in rpm, speed ratio is
0.45, the relative velocity is reduced by 10 % of its magnitude before leaving the
buckets due to friction in buckets.
(88.5 %, 115981 hp) (Alex. Univ., B.Sc., 1970)
2- An inward Francis turbine rotates at 150 rpm and discharges 0.8 m3/s. The radial
velocity at inlet Cr1 = 2 m/s and equal 6 m/s at exit, the physical data are:
r1 = 0.5 m, r2 = 0.2 m, α1 = 15°, β2 = 135°, the pressure head at exit is 6 m.
Assuming a loss of head of 2 m, find:
(a) The head, power delivered by the turbine (no draft tube), neglect hydraulic losses.
(b) The pressure head at entrance.
3- A Kaplan turbine develops 2000 H.P. under a head of 6 m. The turbine is set 2.5 m
above the tail water level. A vacuum gauge inserted at the turbine outlet records
a suction head of 3.1 m. If the turbine efficiency is 85 %, what will be the
efficiency of the draft tube, having inlet diameter of 3 m ? Neglect hydraulic
losses (Neglect hydraulic losses in the draft tube).
4- The hydro-power station comprises three propeller turbine sets, each develops
10,000 HP at full load and has the following characteristics at constant speed:
Head available is 20 m
Load in % of maximum load 20 40 60 80 100
Overall efficiency (%) 35 58 73 85 92
If the average monthly load on the station during a certain year is:
Months 1,2,3,4 5,6,7 8,9,10,11,12
Average load H.P. 16,000 8,000 22,000
Determine the total amount of water passing in station during this year in the
above case, and in the case when one turbine only of 30,000 HP capacity is used
having the same above characteristics, which case would you prefer and why?
347
10- A Francis turbine has a runner diameter 2 m, outlet diameter 1 m, outer width
0.5 m, running at 310 rpm under a net head of 90 m. If at full load, flow velocity
is 6 m/s, mechanical and hydraulic efficiencies are 80 and 86 %, respectively,
calculate the output and speed ratio.
If the velocity of flow is reduced by the use of guide vanes to 3 m/s, draw
carefully the inlet and outlet velocity triangles for the above two cases. Estimate
approximately the values of mechanical and hydraulic efficiencies.
( Alex. Univ., B.Sc., 1982 )
******
348
Mech.MuslimEngineer.Net
CHAPTER VII
FANS, BLOWERS
AND
TURBO-COMPRESSORS
R.T1 P
m
The polytropic head H poly = 2 − 1
m P1
n −1 k −1
with m= =
n kηp
U 2 C u 2 − U 1C u1
Ho =
g
Nominal meter cube (N m3) is the volume of a gas at 0°C and 760 mm Hg.
A standard cubic foot per minute SCFM is the volume of a gas at 66°F (520°R) and
1497 psia.
349
Solution
R.T1 P
m
(a) The polytropic head: H poly = 2 − 1
m P1
k −1
where m=
kηp
P2 remains constant, thus,
H poly 1 Pv T
= 1 1 = 1, Pv = RT
H poly 2 P2 v 2 T2
and R = R / MW = 848/29 = 29.25 m/K
1.4 −1
29.25 × (20 + 273) 1.56 1.4×0.7
H1 = − 1
1.4 − 1 1.03
1.4 × 0.7
H 1 = 21005 (1.1845) = 3876 m
i.e. H 2 = H 1 (T2 T1 ) = 3876(316 293) = 4180 m
The nominal meter cube is the volume of a gas at atmospheric pressure and
0°C
Q = v o Q2 / v 2
350
Mech.MuslimEngineer.Net
vo is the specific volume of air at 0°C = R.To / Pa
Po = P1 inlet pressure, then:
4180 × 35311
i.e. Power = = 869 H .P.
3600 × 0.899 × 0.7 × 75
869
The brake horsepower = ∴ B.H.P. = 965 H.P.
0.9
Example (2)
The impeller of a centrifugal fan has an inner radius of 250 mm and width of
187.5 mm; the values at exit are 375 mm and 125 mm, respectively. There is no whirl
at inlet, and at outlet the blades are backward facing at 70° to the tangent. In the
impeller there is a loss by friction of 0.4 times the kinetic head corresponding to the
relative outlet velocity, and in the volute there is a gain equivalent to 0.5 times the
kinetic head corresponding to the absolute velocity at exit from the runner. The
discharge of air is 5.7 m3/s when the speed is 13.5 rev/s. Neglecting the thickness of
the blades and whirl slip, determine the head across the fan and the power required to
drive it if the density of the air is sensibly constant at 1.25 kg/m3 throughout and
mechanical losses account for 220 W.
Solution
r1 = 250 mm
b1 = 187.5 mm
r2 = 375 mm
b2 = 125 mm
N = 13.5 rev/s
180° − β2 = 70°
W2 C2
Losses: 0.4 2 , 0.5 2
2g 2g
3
Q = 5.7 m /s
ρ air = 1.25 kg/m 3
Mechanical Losses = 220 W
351
352
Mech.MuslimEngineer.Net
Problems
1- In a centrifugal compressor the inlet absolute flow to the runner is radial, and the
exit relative flow from the runner is radial and the runner operates at 22000 rpm
and has an outside diameter of 28 cm, the used gas is air and the inlet temperature
is 18°C. If the pressure coefficient φ = 0.5, find the pressure ratio.
2- A centrifugal fan runner consists of some blades, two parallel circular disks and a
short length of circular pipe in an axial direction at inlet. Assume an
incompressible fluid with a certain density in slugs per cubic meter. The runner
speed is N (rpm), the outside diameter of the runner is D (meter), and the inlet
diameter is 0.4 D, assume that the area of the inlet pipe equals the circumferential
area (between the parallel disks) at the inlet diameter, the relative velocity of exit
is radial, with a magnitude 0.3 of the velocity. Develop an expression for the
torque in meter-kg of the runner on the fluid in terms of density, N, D, and some
numerical constants.
******
353
VOLUMETRIC MACHINES
π N
Q= D 2 2r Stroke = 2 r
4 60
π
The swept volume = D 2 2r
4
η v = Qa / Q
2π n.r D 2
The velocity in the delivery pipe V p = sin θ , where d is the pipe diameter.
60 d 2
P1 l 2 D 2
= ω r 2 cos θ
ρ g g d
H ma = H ss + H sd +
2
(h fs + h fd )
3
354
Mech.MuslimEngineer.Net
8.4 Gear Wheel Pump
2a..l.n.N
Q=
60
where a area enclosed between any two teeth and the casing,
l axial length of teeth,
n number of teeth in each gear,
N speed in rpm.
355
Solution
Discharge per revolution = 350 x 1000 /750 = 466.6 cm3 / revolution
Q = 466.6 / 0.9 = 518.4 cm3 / revolution
Q = KC ( D − C ) l
Q = 2.98 C ( D − C ) l
D D
Q = 2.98 D− (15
. D)
.
118 .
118
Q is the discharge per revolution, from which the rotor diameter D = 9.65 say 10 cm.
L = Length = 15 cm
Example (2)
A single-acting reciprocating water pump, with a bore and stroke of 150 mm
and 300 mm respectively, runs at 0.4 rev/s. Suction and delivery pipes are each
75 mm diameter. The former is 7.5 m long and the suction lift is 3 m. There is no air
vessel on the suction side. The delivery pipe is 300 m long, the outlet (at atmospheric
pressure) being 13.5 m above the level of the pump, and a large air vessel is
connected to the delivery pipe at a point 15 m from the pump. Calculate the absolute
pressure head in the cylinder at beginning, middle and end of each stroke. Assume
that the motion of the piston is simple harmonic, that losses at inlet and outlet of each
pipe are negligible, that the slip is 2 % and that f for both pipes is constant at 0.01.
(Atmospheric pressure is 10.34 m water head).
356
Mech.MuslimEngineer.Net
Solution
I. Suction Stroke
357
Absolute Pressure = Hat + Hsd + hfd1 − hfd2 − Had2 + Vd2/2g (The last term = 0)
ALN
= (π / 4) (0.15) × 0.30 × 0.4 = 2.12 × 10 −3 m 3 /s
2
Q=
60
Q − Qa
Q slip = 2 % = × 100
Q
∴ Qa = 0.98 Q = 0.98 (2.12 × 10 −3 ) = 2.0776 × 10 −3 m 3 /s
Qa 2.0776 × 10 −3
Vd = = = 0.47 m/s
A p (π / 4) (0.075)2
hf d = f
l d1 Vd 2
= 0.01 ×
(300 − 15) (0.47 )
×
2
1 d d 2g 0.075 2 × 9.81
= 0.428 m
2
ld2 A
hf d2 = f (ω r. sin θ )2
2 g d d a d
4
15 0.15
× (2π × 0.4 × 0.15) sin θ
2
= 0.01 ×
2
- At the end θ = 2 π
Absolute Pressure = 24.268 − 0 − 5.795 = 18.473 m water
358
Mech.MuslimEngineer.Net
Problems
1- A plunger pump works against a total static head of 96 m and when running at
42 rpm. It is required to force 45 lit/s of water along a delivery pipe 25 cm
diameter and 130 m long, there are no air vessels, the stroke of plungers is twice
the diameter. If the number of cylinders chosen were 1, 2, 3 and 4, calculate in
each case:
i) the plunger diameter,
ii) the maximum pressure in the pipe.
2- A gear wheel pump is required to deliver 4 lit/s of oil of specific gravity of 0.94
when running at 700 rpm. The suction pressure is 0.2 kg/cm2 and the delivery
pressure is 6 kg/cm2, the overall efficiency is 90 %. The length of the gear wheels
or rotors is 2 x maximum diameter, and in shape they are geometrically similar to
those shown in the figure (Chapter VIII), what would be their outside diameter
and what would be the power input to the pumps?
3- A reciprocating pump has two double acting cylinders each 200 mm bore, 450 mm
stroke, the cranks being at 90° to each other and rotating at 20 rev/min. The
delivery pipe is 100 mm diameter, 60 m long and there are no air vessels.
Assuming simple harmonic motion for the piston, determine the maximum and
mean water velocities in the delivery pipe and the inertia pressure in the delivery
pipe near the cylinders at the instant of minimum water velocity in the pipe.
(2.666 m/s, 2.4 m/s, ± 236.9 kPa)
************
******
359