Turbomachinery Slides PDF
Turbomachinery Slides PDF
1
Lect-1
Turbomachinery Aerodynamics
A Video course by
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Course Outline
Introduction to Turbomachineries:
Axial flow compressors and Fans: Introduction;
Aero-Thermodynamics of flow thru axial flow compressor
stage; Losses in axial flow compressor stage; Losses and
Blade performance estimation; Secondary flows (3-D);
Tip leakage flow and scrubbing; 3-D flow analysis; Radial
Equilibrium Equation; Axial compressor characteristics;
Design of compressor blades-2-D blade designs; Airfoil
Data; Axial Flow Track Design; Multi-staging of
compressor characteristics; Transonic Compressors;
Shock Structure Models in Transonic Blades; Transonic
Compressor Characteristics; 3-D Blade shapes of Rotors
and Stators; Instability in Axial Compressors; Loss of
Pressure Rise; Loss of Stability Margin; Noise problems
in Axial Compressors and Fans
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Course Outline
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Course Outline
Centrifugal Compressors :Introduction;
Elements of centrifugal compressor/ fan; Inlet
Duct ; Impeller flow; Effect of Slip factor; Concept
of Rothalpy; Ideal and real work done; Incidence
and lag angles; Diffuser ; Centrifugal Compressor
Characteristics ; Surging and Rotating stall;
Design variants of modern centrifugal compressors
Radial Turbine: Introduction; Thermodynamics
and Aerodynamics of radial turbines; Radial
Turbine Characteristics; Losses and efficiency;
Design of radial turbine
Use of CFD for Turbomachinery analysis and
design
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Course Pre-requisites
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Text/References
Nicholas Cumpsty, Compressor Aerodynamics,
2004, Kreiger Publications, USA
Johnson I.A., Bullock R.O. NASA-SP-36, Axial Flow
Compressors, 2002 (re-release), NTIS
NASA-SP-290, Axial Flow turbines, 2002 (re-
release), NTIS, USA.
J H Horlock, Axial flow compressors, Butterworths,
1958, UK
J H Horlock, Axial Flow Turbines, Butterworths,
1965, UK
B Lakshminarayana; Fluid Mechanics and Heat
Transfer in turbomachineries, 1995, USA
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Lecture schedule
No. Topic Speak
er
1 Introduction to Turbomachineries : Syllabus, BR /
References and schedules AMP
2 Axial flow compressors and Fans : Introduction to AMP
compressor aerothermodynamics
3 A two dimensional analytical model (Cascade) AMP
4 2-D Losses in axial flow compressor stage AMP
primary losses
5 Tutorial-1 solved examples and tutorial problems AMP
6 3-D flows in Blade passages, Secondary flows, BR
Tip leakage flow, Scrubbing
7 Three dimensional flow analysis Radial BR
Equilibrium concept
8 Classical blade design laws Free vortex and other BR
Laws
9 Tutorial-2 - solved examples and tutorial BR
problems and Quiz 1
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
A Brief Introduction
to
Compressors and Turbines
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
Fan
Turbine
Compressor
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-1
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 2
1
Lect-2
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Introduction
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Thermodynamics of compression
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Thermodynamics of compressors
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Thermodynamics of compressors
The compression
process is usually
expressed in H-s or T-
s diagrams.
The ideal compression
process is assumed to
be isentropic.
Deviation from this is
expressed as
isentropic efficiency.
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Thermodynamics of compressors
X1 , X2 are the
losses in the rotor
and the stator
respectively
Compression in terms
of static parameters
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Thermodynamics of compressors
Compression in terms
of total parameters
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Velocity triangles
Elementary analysis of axial compressors begins
with velocity triangles.
The analysis will be carried out at the mean height
of the blade, where the peripheral velocity or the
blade speed is, U.
The absolute component of velocity will be
denoted by, C and the relative component by, V.
The axial velocity (absolute) will be denoted by Ca
and the tangential components will be denoted by
subscript w (for eg, Cw or Vw)
denotes the angle between the absolute velocity
with the axial direction and the corresponding
angle for the relative velocity.
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Velocity triangles
1 2 3
V2
2
V1
1 2
C2
V1 3
V2 C3
U
U
C1 C2
Rotor Stator
C = U +V
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Velocity triangles
V2
2 Vw2
2
V1 Vw1
Cw
U
1 C2 Cw2
1
C1 Cw1
Ca
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
Total enthalpy
h01 h02 h03
Absolute velocity
C1 C2 C3
Static pressure
P1 P2 P3
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
In this lecture...
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-2
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 3
1
Lect-3
In this lecture...
Design parameters
Two dimensional analysis: Cascade
aerodynamics
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Design parameters
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Degree of reaction
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Degree of reaction
Static enthalpy rise in the rotor
Rx =
Stagnation enthalpy rise in the stage
h h1 h h1
= 2 2
h03 h01 h02 h01
For a nearly incompressible flow,
(P03 P01 )
1
and for the stage, h03 h01
h2 h1 P2 P1
Rx =
h02 h01 P02 P01
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Degree of reaction
From the steady flow energy equation,
V12 V22
h1 + = h2 +
2 2
h2 h1 V12 V22
Rx = =
h03 h01 2U(Cw 2 Cw1 )
For constant axial velocity, V12 V22 = Vw21 Vw2 2
And, Vw1 Vw 2 = Cw1 Cw 2
1 Ca
On simplification, R x = (tan 1 tan 2 )
2 2U
Ca
or, R x = (tan 1 + tan 2 )
2U
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Degree of reaction
Special cases of Rx
Rx=0, 2 = 1 , There is no pressure rise in the
rotor, the entire pressure rise is due to the
stator, the rotor merely deflects the incoming
flow: impulse blading
Rx=0.5, gives 1 = 2 and 2 = 1 , the velocity
triangles are symmetric, equal pressure rise in
the rotor and the stator
Rx=1.0, 2 = 1 , entire pressure rise takes
place in the rotor while the stator has no
contribution.
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Degree of reaction
2
2
2 = 1
V2 V2 V2
2 2
2 2
V1 V1
U U U
1 C2 1 C2 C2
1 1 V1
C1 C1
C1
2 = 1 1 = 2 and 2 = 1 1
1
Diffusion factor
Fluid deflection (2-1)is an important
parameter that affects the stage pressure rise.
Excessive deflection, which means high rate of
diffusion, will lead to blade stall.
Diffusion factor is a parameter that associates
blade stall with deceleration on the suction
surface of the airfoil section.
Diffusion factor, D*, is defined as
Vmax V2
D = Where, Vmax is the ideal surface velocity at
V1
the minimum pressure point and V2 is the ideal velocity
at the trailing edge and V1 is the velocity at the leading edge.
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Diffusion factor
Suction
surface
V1
Vmax
V2
Velocity
Pressure
surface
0 50 100
Percent chord
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Diffusion factor
Lieblein (1953) proposed an empirical
parameter for diffusion factor.
It is expressed entirely in terms of known or
measured quantities.
It depends strongly upon solidity (C/s).
It has been proven to be a dependable indicator of
approach to separation for a variety of blade
shapes.
D* is usually kept around 0.5.
V2 Vw1 Vw 2
D = 1 +
V1 C
2 V1
s
Where, C is the chord of the blade and s is the spacing
between the blades.
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Cascade aerodynamics
A cascade is a stationary array of blades.
Cascade is constructed for measurement of
performance similar to that used in axial
compressors.
Cascade usually has porous end-walls to
remove boundary layer for a two-dimensional
flow.
Radial variations in the velocity field can
therefore be excluded.
Cascade analysis relates the fluid turning
angles to blading geometry and measure
losses in the stagnation pressure.
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Cascade aerodynamics
The cascade is mounted on a turntable so
that its angular direction relative to the
inlet can be set at different incidence
angles.
Measurement usually consist of pressures,
velocities and flow angles downstream of
the cascade.
Probe traverse at the trailing edge of the
blades for measurement.
Blade surface static pressure using static
pressure taps: cp distribution.
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Cascade nomenclature
C
C = Chord
s = spacing / pitch
t = thickness
= camber C
= stagger
i = incidence angle
= deflection angle
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Cascade aerodynamics
The cascade is mounted on a turntable so
that its angular direction relative to the
inlet can be set at different incidence
angles.
Measurement usually consist of pressures,
velocities and flow angles downstream of
the cascade.
Special nulling type probes (cylindrical,
claw or cobra type) are used in the
measurements.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Performance parameters
Measurements from cascade: velocities,
pressures, flow angles ...
Loss in total pressure expressed as total
pressure loss coefficient
P01 P02
WPLC = 1
2 V1
2
Performance parameters
Blade performance/loading can be
assessed using static pressure coefficient:
Plocal Pref
CP = 1
2 V1
2
Performance parameters
Deflection, degrees
Total pressure loss coefficient
Location of the
blade trailing edge
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Performance parameters
C2
C2
2
2
1 Stalled or
C1
1 separated flow
C1
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
Performance parameters
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
In this lecture...
Design parameters
Two dimensional analysis: Cascade
aerodynamics
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-3
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 4
1
Lect-4
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
Performance parameters
Measurements from cascade: velocities,
pressures, flow angles ...
Loss in total pressure expressed as total
pressure loss coefficient
P01 P02
PLC = 1
2 V1
2
Performance parameters
Blade performance/loading can be
assessed using static pressure coefficient:
Plocal Pref
CP = 1
2 V1
2
Performance parameters
Deflection, degrees
Total pressure loss coefficient
Location of the
blade trailing edge
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
Performance parameters
C2
C2
2
2
1 Stalled or
C1
1 separated flow
C1
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
Performance parameters
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
= P + sh + s + L + E
Where, P : profile losses
sh : shock losses
s : secondary flow loss
L : tip leakage loss
E : Endwall losses
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
[ ]
+ 2.6 + 0.18(1' 650 ) 10 2 (1 1' )
where, 1' is the blade inlet angle.
This is valid for an incidence angle upto 5o.
These empirical correlations are however,
derived using the 2-D assumption.
Actual flows are seldom 2-D in nature.
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
In this lecture...
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-4
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 5
1
Lect-5
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Problem 1
Air at 1.0 bar and 288 K enters an axial flow
compressor with an axial velocity of 150 m/s. There are
no inlet guide vanes. The rotor stage has a tip diameter
of 60 cm and a hub diameter of 50 cm and rotates at
100rps. The air enters the rotor and leaves the stator in
the axial direction with no change in velocity or radius.
The air is turned through 30.2 degree as it passes
through the rotor. Assume a stage pressure ratio of 1.2
and overall pressure ratio of 6. Find a) the mass flow
rate of air, b) the power required to drive the
compressor, c) the degree of reaction at the mean
diameter, d) the number of compressor stages required
if the isentropic efficiency is 0.85.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution : Problem # 1
V2
2
U
C2
V1
U
1
Ca
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 1
dt + dh 0.6 + 0.5
U= N = 100 = 172.76m / s
2 2
U
1 = tan = 49.2
1
Ca
2 = 49.2 30.2 = 19
U Ca tan 2
tan 2 = = 80.75
Ca
2 = 38.92
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 1
4
2
(
m = dt dh Ca 2
2
) & T1 = T01
Ca2
2Cp
= 276.8 K
1
P 02
T 02 = T 01 T 02 = 303.41K
P 01
C2 Ca
2
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 1
192.79 2
T 2 = 303.41 = 284.91 K
2010
P 2 = 1.216 bar
1.216 101325
2 = = 1.507 Kg / m3
287 284.9
m = 19.53 Kg / s
P = U Ca m (tan 1 tan 2 )
= 172.76 150 19.53 (tan 49.2 tan 19 ) = 412 KW
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 1
Ca
RX = 1 (tan 1 + tan 2 )
2U
(tan 49.2 + tan 19 ) = 1 0.65
150
=1
2 172.76
= 0.35
U Ca
T0s = (tan 1 tan 2 )
Cp
172.76 150
= (tan 49.2 tan 19 ) = 20.99 K
1005
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 1
T1
1
T0overall = C 1
C
=
288
0.85
(
6 0.286 1 = 226.5K )
226.5
n= = 10.79 11
20.99
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Problem # 2
An axial flow compressor is to be designed to
generate a total pressure ratio of 4.0 with an
overall isentropic efficiency of 0.85. The inlet and
outlet blade angles of the rotor blades are 45
degree & 10 degree respectively and the
compressor stage has a degree of reaction of 50
percent. If the blade speed is 220 m/s and the
work done factor is 0.86, find the number of
stages required. Is it likely that the compressor
will suffer from shock losses? The ambient air
static temperature is 290 K and the air enters the
compressor through guide vanes.
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 2
V2
2
2
V1
U
1 C2
1
C1
Ca
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 2
U
Axial velocity, Ca = = 187 m / s
tan 1 + tan 2
Ca
Absolute velocity at C1 = = 190m / s
cos 1
inlet,
Solution: Problem # 2
T 03 s = T 02 c
= 493.9 K
Solution: Problem # 2
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 2
To determine whether the compressor
will suffer from shock losses, we need to
find the relative Mach number
V1
M rel =
RT2
Ca
V1 = = 264.5m / s
cos 1
M rel = 0.77
Since relative Mach number is less than
unity, the compressor is not likely to suffer
from shock losses.
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Problem # 3
The conditions of air at the entry of an axial
compressor stage are P1=1 bar and T1=314
K. The air angles are 1=51o, 2=9o, 1= 3=7o.
The mean diameter and peripheral speed are
50 cm and 100 m/s respectively. Given that
the work done factor is 0.95, stage
efficiency is 0.88, mechanical efficiency is
0.92 and the mass flow rate is 25 kg/s,
Determine a) air angle at stator entry, b)
blade height at entry and hubtip diameter
ratio, c)Stage loading coefficient, d) Power
required to drive the stage.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 3
a) U
= tan 1 + tan 1
Ca
100
= tan 7 + tan 51 Ca = 73.65m / s
Ca
U
tan 2 + tan 2 =
Ca
100
tan 2 + tan 9 = 2 = 50.18
73.65
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 3
b)
m = Ca ( d h),
Substituti ng known values in the above, h = 0.19 m
dt = 50 + 19 = 69cm,
dh = 50 19 = 31cm
dh
The hub - tip ratio is = 0.449
dt
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
Solution: Problem # 3
d) P = m
w
= 204.75KW is the power required.
m
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-5
In this lecture...
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
1
Lect - 6
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
Fabricated Blades
Top View
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
In certain
blade shapes
the flow, in
passing
through the
blades, develop
two passage
vortices
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
Boundary layer
development at
casing and hub
(due to adverse
pressure
gradient of
main flow)
further
contributes to
3-D flow
development
End-wall Boundary layer development
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
Scrubbing
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
Passage
vortex
development
across blade
passage
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 6
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
1
Lect - 7
Three Dimensional
Flow Analysis in
Axial Compressor
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
It may be recalled
that this element is
also executing a
path through the
curved diffusing
passage between
the rotor blades.
w2
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
Subscript w refers
to tangential / whirl
component of the
flow
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
(p+dp)(r+dr).d.1 p.r.1.d
2(p+dp/2).dr.(d /2).1
= . dr. r. Cw2 /r
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
1 dp = 1 .C2
dr r w
This is called the
Simple Radial Equilibrium Equation
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
dH dCa dCw 1 dp p dp
= Ca +Cw + . - 2
dr drdr drdr -1
9 10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
d dp
=
dr .p dr
Substituting this in the new energy equation
we get
dH dCa dCw 1 dp
= Ca +Cw +
dr dr
dr dr
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
1 dp = 1 .C2
dr r w
We get
2
dH dC a +C dC C
w+ w
= Ca w dr
dr dr r
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
dCw dr
=-
dr r
This yields, on integration
Cw . r = constant.
This condition is commonly known
as the Free Vortex Law
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 7
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
1
Lect - 8
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
Cw . r = constant.
12 3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
w1 w1
w2
w2
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 8
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
1
Lect - 9
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Assume that a
fluid particle, p
is moving in an
arbitrary path
O
within the two
coordinate
systems
o
In axial flow compressors rotors need a
rotating co-ordinate system, whereas the
stator may use a static co-ordinate system
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
d'
Velocity of P w.r.to small xyz, v xyz =
dt xyz
Vectorially, motion of P is summation of the
motion of the moving system w.r.t the fixed
r = R + '
system and particle P w.r.t the moving system
dr ' dR d
Or, = +
dt XYZ dt XYZ dt XyZ
Velocity of P w.r.t .
V XYZ R V xyz + .'
the fixed system = +
XYZ
9 9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
dV XYZ
a XYZ =
dt XYZ
And, acceleration of P w.r.t. Rotating coordinates
xyz,
dV xyz
a xyz =
dt xyz
10 10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
d V XYZ =
a =
dt XYZ
XYZ
dV xyz ** ( '
d )
dt + R + dt
xyz XYZ
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
dV xyz dV xyz
dt = + .V xyz
XYZ dt xyz
12 12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
aXYZ = a vxyz +
xyz + R + 2 . . ' + '
For motion with constant angular velocity
a =a
a xyz+ 2 v xyz + '
XYZ
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Final accn.
a XYZ= a xyz +2 v xyz + ( ' )
14 14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
1p
Hence, =a
s XYZ
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
1 Dv
.p = r +
- 2 Vx
2
x
Dt
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
17 17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Dt Dt r
w r r
using, V = V r i r +Vwt i wt +V a i a
By D( ) = V D ( ) + ds s is length
definition in any
Dt Ds
dt direction
0
as
ds For Steady
=0 State flow
dt
19 19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
DV
DV r d DVw
DV a
d = i r + i t +i a
t
-V w +V r
Dt Dt dt w Dt dt Dt
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
( )
2
1 D V r V tw+r
- =V - ---------- (a)
p r
1 D V a V tw.V r
- . =V + + 2.
r. p
r -------- (b)
.V r
1 D Va
- =V -------------- (c)
p
21 21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Using the
coordinate
systems the
vw flow velocity
in relative
frame is V
and its
components
may be
shown as Va ,
Vr , Vw
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
C wt = V wt +
Then equation (a) and (b) from slide 21 can be
rewritten as
( )
2
1p DVr Cwt
- =V - ---------- (d)
r Ds r
1 p V D (r.Ctw)
- = ----------- (e)
Dr. s r
23 23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
( )
Now, equation 2
(d) from the last 1 DV r C wt
slide can be re- - =V m -
written as : p r
24 24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Now, by our D( )= D( )
earlier definition V Ds V m
Dm
25 25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
D sin. D D 1
Now, = cos. and =-
Dm Dm Dm rm
Where rm is the radius of curvature of the meridional flow
The negative sign is arbitrary. But, for axial flow compressor the
flow track inside generally moves towards lesser or higher
rm , i.e. the flow later on flattens out. Hence,
2 2
1 p C V t D m V
r = + cos. w
V
m
r rDm r
- m
This is the full radial equilibrium Equation for
circumferentially averaged (blade to blade) flow properties
inside of a turbo machine blade row
21 27
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
30
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
31
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
32
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 9
Next Class
33
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
1
Lect-10
Lecture-10
Tutorial -2
Solved Problems and Tutorial Problems
On
Three Dimensional flow in Axial Flow
Compressor
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Solution 1 :
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Example -2 .
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Solution 2
Rx-hub = 0.29
Rx-tip = 0.744
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Using the prescribed law --- in front and behind the rotor;
At the hub :
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Using the prescribed law --- in front and behind the rotor;
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-10
Next Class :
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 11
1
Lect-11
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
V1
V2
U
U
C1 C2
3
C3
Rotor Stator
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
V2
2 Vw2
2
V1 Vw1
Cw
U
1 C2 Cw2
1
C1 Cw1
Ca
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
Ca Ca
U design U
1.0
st
Measured
Stage efficiency
Ca Ca
U design U
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
V1 V1 V1
1 1 1
V1 V1
V1
U 1 U 1 U
1 1
1 1 C1
C1 C1
Ca Ca
Ca
Design condition : Off - design condition : Off - design condition :
Normal operation Positive incidence flow separation Negative incidence flow separation
Ca Ca Ca Ca Ca Ca
= < >
U U design U U design U U design
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
P02
P01
m
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
P02
P01
m
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
Ca Ca
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
P02
P01
m
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
P0
U 2
Stage characteristics
Throttle characteristics
Ca
U
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-11
In this lecture...
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 12
1
Lect-12
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
P02
P01
m
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Definition of Instability
Stability is related to the response of a
compressor to a disturbance which
perturbs compressor operation from a
steady operating point.
If the disturbance is transient, the
performance is considered stable if the
system returns to the original point of
equilibrium.
If the response is to drive the operation
away from the original point, the
performance is unstable.
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Definition of Instability
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Operational Stability
In this form of stability analysis, the complete
compressor system including inlet ducting,
guide vanes, rotors, stators and any pressure
recovery exhaust system are considered.
The operational stability depends upon the
rate of change of pressure rise or pressure
drop as mass flow rate varies.
System stability is established when the rate
of change of compressor pressure rise with
mass flow is algebraically less than the rate of
change of throttle pressure drop with mass
flow.
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Operational Stability
Compressor
Flow device
Matching characteristics
between the compressor and a
downstream flow device
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Aerodynamic Stability
Aerodynamic stability is the ability of the
entire compressor system to maintain or
increase the delivery pressure when the
compressor operation has been perturbed
to a lower flow.
The part of the compressor characteristic
with a positive slope is a region containing
subsystem stall or complete instability
resulting in surge.
In this region, operational stability may be
theoretically possible, but not aerodynamic
stability.
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Aerodynamic Stability
Pressure ratio
Flow
Stall
There are two types of stall that have been
experimentally identified.
Individual blade stall
Rotating stall
Individual blade stall occurs when the
entire blade row stall at once. This type of
stall would be expected in airfoil or
cascade data.
If all the airfoils are identical, all the blades
would experience the stall incidence angle
and hence the complete row of blades
stall.
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
Rotating stall is the most common type of
stall.
Rotating stall: progression around the
blade annulus of a stall pattern, in which
one or more adjacent blade passages are
instantaneously stalled, then are cleared
for unstalled flow as the stall cell
progresses.
Rotating stall causes alternate loading and
unloading of the blades: fatigue failure.
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
Rotating Stall
Stall cells
Rotating Stall
Rotating stall often precedes surge.
The stall patterns move in a direction
opposite to that of the rotor revolution.
The stall frequency can be as high as
50% of the rotor frequency.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
Rotating stall may be initiated due to a
variety of reasons: off-design operation,
inflow distortion, blade stagger/profile
mismatch etc.
If allowed to propagate, rotating stall may
lead to surge of the compressor.
The number of stall cells can be as high as
9 or more or as low as one.
The number of stall cells is associated with
the type of stall.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
There are two types of stall: progressive and
abrupt.
Progressive stall has a gradual reduction in
total pressure ratio after initiation of stall.
Abrupt stall has a sharp discontinuity in the
pressure ratio characteristic.
Progressive stall usually has multiple stall cell,
whereas abrupt stall seems to always occur
with a single cell.
Some of the recent researchers have termed
progressive and abrupt stall as modal and
spike type of stall.
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
Pressure rise
Pressure rise
Stall
Stall
Flow Flow
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Rotating Stall
Tip stall
Part-span
stall
Full-span
Pressure rise
stall
Deep stall
Flow
Variation of stall pattern during progressive stall
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Surge
Surge is a state of operation during which the
complete compression system is unstable.
Surge is characterised by fluctuations in the
average flow throughout the compressor.
The net flow through the compressor can be
positive or negative in this highly transient
state.
This differs from rotating stall wherein the
average flow through the compressor system
is constant with time.
During surge, the system is unable to attain
a stable match with a throttle.
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Surge
Surge line denotes the locus of unstable
operation of the compressor.
Surge is characterised by violent, periodic
oscillations in the flow.
Surge might lead to flame blow-out in the
combustion chamber.
Surge can lead to substantial damage to
compressors and must be avoided.
The operating line of the compressor is
therefore kept slightly away from the surge
line: surge margin.
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Surge
P02
P01
m
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
Surge
Hysteresis
Surge
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
In this lecture...
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-12
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 13
1
Lect-13
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Inlet distortion
Engine performance significantly depends
upon the quality.
Air inlets are required to provide the
necessary quantity of good quality air to
the engine.
The exit flow may become non-uniform
under a variety of circumstances:
manoeuvre, geometry of the intake,
boundary layer ingestion, wakes/jet plume
from freestream, cross-wind etc.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Inlet distortion
Intakes of civil and military combat aircraft
have very different geometries.
Combat aircraft intakes can have very complex
geometries leading to inherent problem of flow
non-uniformity.
Inflow non-uniformity or distortion is
detrimental to engine operation.
Several aircraft in the past, that were
operating with engines not designed for
distortion have had serious operational issues
including several engine failures.
Some of these are F100 (1954), F101 (1954),
Hunter (1955), Britannia (1956), F111 (1966)
etc.
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Inlet distortion
Inlet distortion
Inflow distortion can lead to earlier
initiation of instabilities: rotating stall and
surge.
Flow distortion causes local change in
incidence angles.
If these angles exceed the critical angles,
one or more adjacent blades falling in this
zone of distortion, stalls.
The stall cell(s) propagate if these are able
to withstand the system dynamics. Else,
they dissipate and the system does not
undergo any instabilities.
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Inlet distortion
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Distortion coefficient
Quantification of distortion: total pressure
non-uniformity
Most commonly used measure: distortion
coefficient.
Inlet distortion
Distortion coefficient
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Swirl
Many of the military aircraft have engines
that are offset from the intake centerline.
Such intakes referred to as S-type or Y-
type intake, inherently suffer from strong
secondary flows.
In the absence of guide vanes, the flow
entering the compressor is likely to have
some amount of swirl.
This swirl may get amplified under certain
operating conditions, leading to severe
inflow distortion.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Swirl
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Swirl
Control of instabilities
Compressor instabilities limit the operating
range of an engine.
Operating the compressor too much away
from the surge line compromises the
efficiency.
Ability to operate the compressor close to high
efficiency points (and possibly closer to the
surge line) is of immense interest.
This would require methods of preventing or
controlling the occurrence of instabilities.
The other way of preventing this altogether is
to control the inflow from the inlet by flow
control methodologies.
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
There are several methods that have been
proposed by researchers over the past 50
years or so.
These can be broadly classified as Passive
and Active control techniques.
Passive control
Does not involve any external energy addition.
Control scheme incorporated by design changes
on the compressor blade and/or the compressor
casing.
Simpler to design and implement.
Disadvantage: cannot be controlled, may lead to
performance penalties when the control is not
required.
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
Active control
Involves addition of energy external to
the system.
Separate control scheme and associated
components need to be designed and
integrated with the compressor system.
More complex, difficult to design and
implement.
Can be controlled, switched-off when
not required, minimal performance
penalties.
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
Passive control methods
Casing treatments
Proposed in late 40s
Involves making grooves/slots on the casing
above the rotor.
Affects the tip flow behaviour.
Delays stall and therefore offers better stall
margin.
However reduces the efficiency.
Area of active research to develop casing
treatments that improve stall margin without
efficiency penalty.
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
Control of instabilities
Passive control methods
Blade shape modifications
Sweep and dihedral
Non-radial blade stacking methods.
Depending upon the orientation, can
significantly alter the rotor tip flow
characteristics.
Envisaged to improve the stability
characteristics as well as the efficiency.
Currently under research and
development.
Other methods: tandem blading, vortex
generators, fins etc.
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
Control of instabilities
Active control methods
Compressed air injection from casing
Energises the blade tip region, making the tip
flow more resistant to adverse pressure
gradients.
Use air from later stages of the compressor
for injection.
Expected to improve the stability margin and
possibly efficiency.
Variants of tip injection scheme: steady
injection, pulsed injection, injection at
varying angles (skew and pitch).
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
28
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
Control of instabilities
Other active control methods
Variable IGVs
Bleed valves
Typically used during starting to
prevent stall due to front and rear
stage mismatch
Plasma actuators and synthetic jets
Are in premature state of research
Seem to show promise under certain
operating conditions.
29
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-13
In this lecture...
30
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
1
Lect 14
Transonic Compressors
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
V1 or C2
may be
supersonic
at tip or at
Normally a transonic many
compressor has either sections of
V1 or C2 supersonic the blade
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
1. Supersonic flow
enters the rotor.
2. Shock is
contained within
the rotor blade
and leaves the
rotor subsonically
3. Flow in the stator
is subsonic
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
1. Rotor performs a
large flow turning
subsonically
2. Very large energy
transfer in rotor
3. Rotor exit flow has
large K.E.
4. Large diffusion
needs to be done
in the stator
5. Thus, stator needs
to be supersonic
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
1. High stage
performance
2. Flow enter
both rotor
and stator
supersonically
3. And exit
subsonically
4. Both rotor
and stator
blades are
highly loaded
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
M = 1.3 Developed in
60s 2 arcs
CDA Blades :
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
DCA blades
Flow through a
transonic
blading would
diffuse through
the shocks
before further
diffusing and
exiting as
subsonic flow
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
Multiple
Circular Arc
M = 1.5
M = 1.8
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
DCA Blades :
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
MCA Blades :
MCA blades are used for compressors /fans with
low solidity and higher Mach number (>1.4).
This shape was created for greater control of the
blade profile by using multiple arcs.
These blade shapes create a bow shock.
These MCA blades, used near the tips, are set at
high stagger, due to which the inflow experiences
a mildly converging (virtual) passage. The suction
surface of the blade is convexly curved resulting in
a series of mild shock fans.
The entry flow through the shock fans is, thus,
supersonically diffused till the passage shock ,
through which it finally becomes subsonic.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
S-type Blades:
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 14
Next Class :
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
1
Lect-15
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
PT
m T N
02
, 0C , 0
=f 01
,R e,
P01 T01 P01 T
01
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
Stall
2
margin
N/T01
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
4 14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
0.7 - 0.8
W
W C,m = C,1
For transonic stages
0.8 - 0.9
WC
number of stages Z=
WC,m
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
Next Lecture
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-15
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
1
Lect - 16
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
rm
Cw,1-r = Cw,1m .
r
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
C1-r = (C 2
a,1r +C 2
w,1r ) Absolute Vel
r2
U 2,r = U 2,m . If, dm = constant, U1,m=U2,m
r2,m
r2m Cw,rm
Cw,2r = Cw,2m . Check DR = 1 -
r 2.U.rm
Degree of Reaction, Rx should never be
zero anywhere on the rotor blade
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
11 6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
V 2
> V 1 transonic fan design possibility
2,r = tan
-1
Cw,2r
2,r = tan -1
(
U -C
2,r w,2r )
C Ca,2r
a,2r
= 2,r - 1,r
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
c/s = , solidity
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
Mu = U/a
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
= 1,r + i r
1,r
2,r = 2,r - r
At any radius
s
Deviation,
r
- 2,r
= 2,r = mr ..
c
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
- i r
- 1,r
r = 2,r =
s
1+ mr
c
2
ai 90 - 2,r
Where, m = 0.23 2 * + 0.1
ci 50
a
ai
= 0.4 to 0.5
ci
c
17 12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
19 14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
Transonic Compressor
Basic
Characteristics
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 16
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
1
Lect 17
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Airfoils at various
radial stations are
arranged in a manner
such that the airfoils
are at constant radius
(from LE to TE).
The centroids of the
airfoils are in a radial
locus
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Flat Tip
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Casing
Smaller tip gap
Large tip gap at mid
chord section
Tip rounding
LE TE
Hub
A Flat tip creates a divergent-convergent tip gap (from LE to
TE) which is meridionally arranged as per tip-airfoil stagger
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
3D Airfoil
Isometric View Top View
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Cp distribution of a
standard airfoil shall
change when it is set
on a meridional plane
and in cylindrical
coordinate system
Restoration of the
original Cp
distribution shall
require the original
airfoil shape to be
altered
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Trajectory of
sweep
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Swept
blade
portion
LE TE
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Straight
Original leaned
radial airfoils
stacking line
Lean
stacking line
Hub
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Swept
Blades
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 17
Next Lecture
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
1
Lect - 18
2
Prof Bhaskar Roy, Prof A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
3
Prof Bhaskar Roy, Prof A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Fundamentals of acoustics:
4
Prof Bhaskar Roy, Prof A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Fundamentals of acoustics:
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Measure of Sound
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Measure of Sound
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Measure of Sound
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Noise meters
are designed
to measure
either in 1/3
octave spectra
Or
In Full Spectra
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Noise due to
interaction between
the rotor and the
stator may propagate
or decay
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Rotor stator
spacing plays an
important role in
noise generation.
Higher spacing
produces less
noise but also
reduces
compressor stage
efficiency and
increases engine
size and weight
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Noise created
during
different
operations
are different
in nature
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 18
Next Class
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
1
Lect 19
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
Introduction
Gas Turbine engine derives its name from the turbine,
which is at the heart of the work producing
mechanism of the engine.
Principle
A fluid with large kinetic energy content is allowed to
hit a freely rotating set of blades, certain amount of
energy can be extracted from the passing fluid as shaft
power
3 3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
03 is the
delivery
from C.C.
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
Impulse turbines
High energy flow is accelerated in a stator and made
to impinge on the rotor with high momentum and
then made to take huge turn through the passage
between the blades. The work transfer is through
large angular momentum change through the blades.
Reaction turbines
The flow is accelerated through the rotor blade
passage which is an converging curved nozzle
passage; Jet effect creates a reaction force as per
Newtons 3rd law of motion. This work is in addition
to the work done by large turning.
8 8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
Impulse turbines
V2=V3
reaction turbines
V3>V2
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
; Isentropic Efficiencies
T0T
Total-to-total 0T = /
efficiency, T0T
Static-to-static TT TT
efficiency,
= =
TT TStator + TRotor
T / / /
12 12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
Next Class :
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 19
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 20
1
Lect-20
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Disc
1 2 3
Velocity triangles
Elementary analysis of axial turbines too begins
with velocity triangles.
The analysis will be carried out at the mean height
of the blade, where the peripheral velocity or the
blade speed is, U.
The absolute component of velocity will be
denoted by, C and the relative component by, V.
The axial velocity (absolute) will be denoted by Ca
and the tangential components will be denoted by
subscript w (for eg, Cw or Vw)
denotes the angle between the absolute velocity
with the axial direction and the corresponding
angle for the relative velocity.
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Velocity triangles
C1 1
1
Stator/Nozzle
2
2 V2
C2 2
U
Rotor
3
V3 3 3
C3
U
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
C3
Cw3 C2
3 U
C
Cw2
3 2 2
Vw3
V3 V2 3
V2 2 V3
Vw2
2 2
U
Ca
Stator/Nozzle Rotor
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
C3
C2
U 2
V3
V2
C2
2
V2 3
U
V3
Stator/Nozzle Rotor
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
A cascade is a stationary array of blades.
Cascade is constructed for measurement of
performance similar to that used in axial
turbines.
Cascade usually has porous end-walls to
remove boundary layer for a two-dimensional
flow.
Radial variations in the velocity field can
therefore be excluded.
Cascade analysis relates the fluid turning
angles to blading geometry and measure
losses in the stagnation pressure.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
Turbine cascades are tested in wind tunnels
similar to what was discussed for compressors.
However, turbines operate in an accelerating
flow and therefore, the wind tunnel flow driver
needs to develop sufficient pressure to cause
this acceleration.
Turbine blades have much higher camber and
are set at a negative stagger unlike
compressor blades.
Cascade analysis provides the blade loading
from the surface static pressure distribution
and the total pressure loss across the cascade.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
From elementary analysis of the flow through
a cascade, we can determine the lift and drag
forces acting on the blades.
This analysis could be done using inviscid or
potential flow assumption or considering
viscous effects (in a simple manner).
Let us consider Vm as the mean velocity that
makes and angle m with the axial direction.
We shall determine the circulation developed
on the blade and subsequently the lift force.
In the inviscid analysis, lift is the only force.
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
Turbine Cascade
Circulatio n, = S(Vw 2 Vw1 )
and lift, L = Vm = VmS(Vw 2 Vw1 )
Expressing lift in a non - dimensiona l form,
L VmS(Vw 2 Vw1 )
Lift coefficien t, CL = 1 =
2
V 2
m C 1
2
VmC
2
S
=2 (tan 2 tan 1 ) cos m
C
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
Viscous effects manifest themselves in the
form to total pressure losses.
Wakes from the blade trailing edge lead to
non-uniform velocity leaving the blades.
In addition to lift, drag is another force that
will be considered in the analysis.
The component of drag actually contributes to
the effective lift.
We define total pressure loss coefficient as:
P01 P02
= 1
2
V2
2
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
Turbine Cascade
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
Turbine Cascade
Based on the calculation of the lift and drag
coefficients, it is possible to determine the
blade efficiency.
Blade efficiency is defined as the ratio of ideal
static pressure drop to obtain a certain change
in KE to the actual static pressure drop to
produce the same change in KE.
CD
1 CL tan m
b = CD
1+ CL cot m
If we neglect the CDterm in the lift definition,
1
b =
2CD
1+
CL sin 2m
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
In this lecture...
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-20
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 21
1
Lect-21
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Degree of reaction
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Velocity triangles
C1 1
1
Stator/Nozzle
2
2 V2
C2 2
U
Rotor
3
V3 3 3
C3
U
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Degree of reaction
Static enthalpy drop in the rotor
Rx =
Stagnation enthalpy drop in the stage
h h3
= 2
h01 h03
Since, in a coordinate system fixed to the rotor, the
apparent stagnation enthalpy is constant,
V32 V22
h2 h3 =
2 2
If the axial velocity is the same upstream and downstream
of the rotor, this becomes,
h2 h3 = 1
2
(Vw2 3 Vw2 2 ) = 1
2
(Vw3 Vw 2 )(Vw3 + Vw 2 )
Also, sin ce h01 h03 = U(Cw 2 Cw3 )
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Degree of reaction
(Vw3 Vw 2 )(Vw3 + Vw 2 )
RX =
2U(Cw 2 Cw3 )
Since, (Vw3 Vw 2 ) = (Cw3 Cw 2 )
(Vw3 + Vw 2 )
Therefore, R X =
2U
We know that, Vw3 = Ca tan 3
and Vw 2 = Ca tan 2 U
Ca
so that R X = 1
2 1 U (tan 2 + tan 3 )
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Degree of reaction
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
C3
Cw3 C2
3 U
C
Cw2
3 2 2
Vw3
V3 V2 3
V2 2 V3
Vw2
2 2
U
Ca
Stator/Nozzle Rotor
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
C3
C2
U 2
V3
V2
C2
2
V2 3
U
V3
Stator/Nozzle Rotor
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
We noted that the aerodynamic losses in
the turbine differ with the stage
configuration, or the degree of reaction.
Improved efficiency is associated with
higher reaction, which implies less work
per stage and therefore a higher number of
stages for a given overall pressure ratio.
The understanding of losses is important to
design, not only in the choice of the
configuration, but also on methods to
control these losses.
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
There are two commonly used turbine
efficiency definitions.
Total-to-static efficiency
Total-to-total efficiency
The usage of the efficiency definition depends
upon the application.
In land-based power plants, the useful turbine
output is in the form of shaft power and
exhaust KE is a loss.
In this case the ideal turbine process would be
isentropic such that there is no exhaust KE.
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
P01
P1
T
01
P2
1
2 P03
03 P3
C 32
2s 2c p
03s 3
3s
Efficiency
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
In many applications (turbojets), the exhaust KE is
not considered a loss as this is converted to thrust
in such machines.
The ideal turbine work in such cases would be
WT, ideal = cP (T01 T03s )
The total - to - total efficiency is defined as
T01 T03
ts =
T01 T03s
T01 T03 1 (T03 / T01 )
= =
[
T01 1 (P03 / P01 )( 1) /
] [
1 (P03 / P01 )( 1) / ]
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
We can compare the two definition s of efficiency by making
an approximation :
T03s T3s T03s T3 = C 32 / 2c p
ts
Therefore, tt =
1 C 32 [2c p (T01 T3s )]
We can see that, tt > ts
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Efficiency
Losses in a turbine
Nature of losses in an axial turbine
Viscous losses
3-D effects like tip leakage flows, secondary
flows etc.
Shock losses
Mixing losses
Estimating the losses crucial designing loss
control mechanisms.
However isolating these losses not easy and
often done through empirical correlations.
Total losses in a turbine is the sum of the
above losses.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Losses in a turbine
Viscous losses
Profile losses: on account of the profile or nature of
the airfoil cross-sections
Annulus losses: growth of boundary layer along the
axis
Endwall losses: boundary layer effects in the corner
(junction between the blade surface and the
casing/hub)
3-D effects:
Secondary flows: flow through curved blade passages
Tip leakage flows: flow from pressure surface to
suction surface at the blade tip
3-D effects are likely to be stronger in a turbine blade
as compared to compressor blade due to high camber
and flow turning
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Losses in a turbine
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
= P + sh + s + L + E
Where, P : profile losses
sh : shock losses
s : secondary flow loss
L : tip leakage loss
E : Endwall losses
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Deviation
Flow at the exit of the rotor does not leave
at exactly the blade exit angle.
It has been found from experience that the
actual exit angle at the design pressure
ratio is well approximated by
2 = cos 1 (d / s)
This is true as long as the nozzle is not
choked.
Under choked condition, a supersonic
expansion may alter the flow direction at
the exit.
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
2
s
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
Reattachment shock
2
Flow
In this lecture...
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-21
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 22
1
Lect-22
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
P02
P01
N / T01
m
T01
P01
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
1.0
0.8
1.0
0.8 0.6
0.4
0.6
t 1.0 N / T01
0.8
Multistage
3-stage
P02
2-stage
P01 1-stage
m
T01
P01
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
1 P2 P1
C1 1
C12
2
SS
PS
SS PS
2 1
C 22
P2 2 P0
C2
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
Turbine
m
To m
To3 m
To Nozzle
Po Po3 Po
m
To 2
m
To 2 Po 2
Po 2
a P02/P03 c b a P03/P0
c b
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
In this lecture...
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-22
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 23
1
Lect-23
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 1
A single stage gas turbine operates at its
design condition with an axial absolute flow at
entry and exit from the stage. The absolute
flow angle at the nozzle exit is 70 deg. At
stage entry, the total pressure and
temperature are 311 kPa and 850oC
respectively. The exhaust static pressure is
100 kPa, the total to static efficiency is 0.87
and mean blade speed is 500 m/s. Assuming
constant axial velocity through the stage,
determine (a) the specific work done (b) the
Mach number leaving the nozzle (c) the axial
velocity (d) total to total efficiency (e) stage
reaction.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 1
C1 1
1
Stator/Nozzle
2
2 V2
C2 2
U
Rotor
3
V3 3 3
C3
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 1
We know that total to static efficiency,
wt
ts =
[ ]
c p T01 1 (P3 / P01 )( 1) /
[
Specific work is, w t = ts c p T01 1 (P3 / P01 )( 1) / ]
[
= 0.87 1148 1123 1 (1 / 3.11)0.248 ]
= 276 kJ / kg
Solution: Problem # 1
C 2 = C w 2 / sin 2 = 588 m / s
We know that T 2 = T01 12 C 22 / c p = 973 K
Hence, M2 = 588 / 1.33 287 973 = 0.97
tt = 0.93
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 1
P01
P1
T
01
P2
1
2 P03
03 P3
C 32
2s 2c p
03s 3
3s
Solution: Problem # 1
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 2
Combustion gases enter the first stage of a gas
turbine at a stagnation temperature and
pressure of 1200 K and 4.0 bar. The rotor blade
tip diameter is 0.75m, the blade height is 0.12
m and the shaft speed is 10,500 rpm. At the
mean radius the stage operates with a reaction
of 50%, a flow coefficient of 0.7 and a stage
loading coefficient of 2.5. Determine (a) the
relative and absolute flow angles for the stage; (b)
the velocity at nozzle exit; (c) the static
temperature and pressure at nozzle exit assuming a
nozzle efficiency of 0.96 and the mass flow.
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 2
1 2 3
C3
C2
U 2
V3
V2
C2
2
V2 3
U
V3
Stator/Nozzle Rotor
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 2
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 2
(b) At the mean radius, rm = (0.75 0.12) / 2 = 0.315m
the blade speed, Um = (10500 / 30) 0.315 = 346.36 m / s
The axial velocity, C a = Um = 242.45 m / s and
Therefore, velocity at the nozzle exit,
C 2 = C a / cos 2 = 242.45 / cos 68.2 = 652.86 m / s
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 2
h01 h 2 1 T2 / T01
The nozzle efficiency, n = =
h01 h 2s 1 (P2 / P01 )( 1) /
1 T2 / T01
(P2 / P01 )( 1) / = 1 = 0.84052
n
P2 = 4 0.84052 4.0303 = 1.986 bar
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 3
A single stage axial flow turbine operates with
an inlet temperature of 1100 K and total
pressure of 3.4 bar. The total temperature drop
across the stage is 144 K and the isentropic
efficiency of the turbine is 0.9. The mean blade
speed is 298 m/s and the mass flow rate is
18.75 kg/s. The turbine operates with a
rotational speed of 12000 rpm. If the
convergent nozzle is operating under choked
condition determine (a) blade-loading coefficient
(b) pressure ratio of the stage and (c) flow
angles.
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 3
C1 1
1
Stator/Nozzle
2
2 V2
C2 2
U
Rotor
3
V3 3 3
C3
U
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 3
(a) The blade loading is defined as
c p T0 1148 144
= = = 1.8615
U2 298 2
Problem # 3
(c) Since the nozzle is choked, the exit Mach number is unity.
Therefore, C 2 = RT2
T02 +1
and = = 1.165
T2 2
The static temperature at the nozzle exit is T2 = 944.2 K.
The absolute velocity of the gases leaving the choked nozzle is
therefore, C 2 = 600.3 m / s.
The axial velocity C a = U = 298 0.95 = 283 m / s.
From the velocity triangles,
cos 2 = C a / C 2 = 283 / 600 = 0.4716 and 2 = 62 o
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 3
U 1
= tan 2 tan 2 =
Ca
1
tan 2 = tan 2 = 0.828 or 2 = 39.6 o
The turbine specific work, w t = c p To = UC a (tan 2 + tan 3 )
c p To 1148 144
or tan 3 = tan 2 = 1.8807 = 0.0793
UC a 298 283
or 3 = 4.54 o
U 1
= tan 3 tan 3 or tan 3 = + tan 3 = 1.132
Ca
and therefore, 3 = 48.54 o
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 4
A multi-stage axial turbine is to be designed
with impulse stages and is to operate with
an inlet pressure and temperature of 6 bar
and 900 K and outlet pressure of 1 bar. The
isentropic efficiency of the turbine is 85 %.
All the stages are to have a nozzle outlet
angle of 75o and equal inlet and outlet rotor
blade angles. Mean blade speed is 250 m/s
and the axial velocity is 150 m/s and is a
constant across the turbine. Estimate the
number for stages required for this turbine.
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Problem # 4
1 2 3
C3
Cw3 C2
3 U
C
Cw2
3 2 2
Vw3
V3 V2 3
V2 2 V3
Vw2
2 2
U
Ca
Stator/Nozzle Rotor
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Solution: Problem # 4
Since the overall pressure ratio is known,
( 1) /
T01 P
= 01 = 6 0.33 / 1.33
T0es P0e
T0es = 576.9 K
Hence, T0overall = t (T01 T0es ) = 0.85(900 576.9) = 274.6 K
From the velocity triangles, C 2 = C a / cos 2 = 150 / cos 75 = 579.5 m / s
T02 = T2 + C 22 / 2c p T2 = T02 C 22 / 2c p
Since, there is no change in stagnation temperature in the nozzle,
T2 = T01 C 22 / 2c p = 900 579.5 2 / 2 1148 = 753.7 K
Solution: Problem # 4
From the velocity triangles at rotor entry,
tan 2 = (C 2 sin 2 U) / C a = (579.5 sin 75 250) / 150 = 2.065
2 = 64.16 o
V2 = C a / cos 2 = 344.15 m / s
We can see that V2 = V3 = C 3 for cons tan t axial velocity.
Therefore, T03 = T2 + C 32 / 2c p = 753.7 + 344.14 2 / 2 1148
= 805.28 K
The temperature drop per stage is
T01 T03 = 900 805.28 = 94.7 K
The number of stages required for the turbine is
T0overall /(T01 T03 ) = 274.6 / 94.7 = 2.89 3 stages.
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Exercise Problem # 1
An axial flow turbine operating with an
overall stagnation pressure of 8 to 1 has a
polytropic efficiency of 0.85. Determine the
total-to-total efficiency of the turbine. If the
exhaust Mach number of the turbine is 0.3,
determine the total-to-static efficiency. If, in
addition, the exhaust velocity of the turbine
is 160 m/s, determine the inlet total
temperature.
Exercise Problem # 2
The mean blade radii of the rotor of a mixed
flow turbine are 0.3 m at inlet and 0.1 m at
outlet. The rotor rotates at 20,000 rev/min and
the turbine is required to produce 430kW. The
flow velocity at nozzle exit is 700 m/s and the
flow direction is at 70 to the meridional plane.
Determine the absolute and relative flow angles
and the absolute exit velocity if the gas flow is 1
kg/s and the velocity of the through-flow is
constant through the rotor.
Exercise Problem # 3
An axial flow gas turbine stage develops 3.36MW at
a mass flow rate of 27.2 kg/s. At the stage entry
the stagnation pressure and temperature are 772
kPa and 727C, respectively. The static pressure at
exit from the nozzle is 482 kPa and the
corresponding absolute flow direction is 72 to the
axial direction. Assuming the axial velocity is
constant across the stage and the gas enters and
leaves the stage without any absolute swirl velocity,
determine (a) the nozzle exit velocity; (b) the blade
speed; (c) the total-to-static efficiency; (d) the
stage reaction.
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-23
Exercise Problem # 4
A single stage axial turbine has a mean radius
of 30 cm and a blade height at the stator inlet
of 6 cm. The gases enter the turbine stage at
1900 kPa and 1200 K and the absolute velocity
leaving the stator is 600 m/s and inclined at an
angle of 65 deg to the axial direction. The
relative angles at the inlet and outlet of the
rotor are 25 deg and 60 deg respectively. If the
stage efficiency is 0.88, calculate (a) the rotor
rotational speed, (b) stage pressure ratio (c)
flow coefficient (d) degree of reaction and (e)
the power delivered by the turbine.
Ans: 13550 rpm, 2.346, 0.6, 0.41, 34.6 MW
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
1
Lect 24
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
Multi-staging
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
Multi-spooling of Turbines
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
C.C
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
IP Turbines
HP Turbines LP Turbines
3-spool
Civil
aircraft
engine
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
h2 - h1
tan =
2S
S =S + S i= S.Z p+ S .Z
i
tan =
h2 - h1
=
(
Si h2 - h1
=
)
h2 1 -
h2
2Z p ( SS
i + ) 2Z S 1+ S 2Z 1+ S
p i p
S i S i
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
Multi-staging of Turbine
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
TWO-SPOOL TURBOJET
ENGINE
16 13
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
Compressor-Turbine spool-by-spool
matching needs to be augmented
with various engine controls that
facilitates better matching and safer
operation at all operating conditions
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 24
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
1
Lect 25
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
rm
simplified radial equilibrium equation is valid
1 dp = 1 .C2
dr r w
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
2 = constant
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
C a2
cot 2 = = const
Cw2
; dC a2 dCw2
C a2 = Cw2 .cot 2 which yields = .cot 2
dr dr
Now invoking the radial equilibrium equation in
energy eqn
2
dH dC a +C dC C
w+ w and, dH = 0
= Ca w dr
dr dr r dr
dCa dCw Cw 2
We get, Ca +Cw + =0
dr dr r
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
2
dCw2 dCw2 Cw2
Cw2 .cot 2 .
2
+ Cw2 . + =0
dr dr r
2
dCw2 Cw2
Cw2( 1 + cot 2 ) .
2
+ =0
dr r
dCw2 dr
= sin 2 .
2
dr r
which, on integratio n yields
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
sin 2
sin
2 rm 2
Cw2 .r = const ; and then C =C
w2m r
2
w2
sin
2
alternate l y, C a2 .r 2
= const
sin 2
rm 2
and then C =C
a2 a2m r
and finally in terms of absolute velocity ,
sin 2
rm 2
C =C
2 2m r
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
if 2 = constant
then,
C C C
w2 = a2 = 2 = r
C C C rm
w2m a2m 2m
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
Now,
There are three possibilities:
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 25
Problem Solving on
Turbine 3-D flow theories
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
1
Lect 26
Tutorial - 4
Solved Problems
And
Exercise Problems
On
3-D flows in Axial Flow
Turbine
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
Example 1.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
Solution 1 :
At the rotor inlet station we know,
sin 2
C C C r 2
w2 = a2 = 2 =
C C C r
w2m a2m 2m m
And, at the rotor exit
cos2
r 2
Ca 3 = Ca 3m + 2U mCw2m 1
2 2
r
m
and
rm /rt = 0.875, and rm /rh = 1.166
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
sin 2
rm 2
Ca2t C
=
a2m r = 257 m/s
sin 2
rm 2
C
w2t
=C
w2m r = 447 m/s
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
cos2
2
Ca 3 = C 2 + 2U C r
m w2 m r
1
a 3m
m
From which we can calculate the axial velocities,
Ca3t = 262 m/s
Ca3h = 306 m/s; Cw3 =is constant radially
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
Example 2
It is proposed that for design of an axial flow
turbine two design methods are to be explored :
Solution 2 :
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
Case (B)
sin 2
r 2
C = C h
Prescribed condition is a2t a2h r
t
C C C
a2t = a2t = a2m
Which essentially means : C C C
a2h a2m 2h
sin 2
For constant nozzle angle: rm 2
C =C
a2 a2m r
sin 2
2 sin 2
rm rm 2
C =C C =C
w2 w2m r ; 2
2m r
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
3 = 0 ; Cw3 = 0
cos2
r 2
Ca 3 = Ca 3m + 2U mCw2m 1
2 2
r
m
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 26
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 27
1
Lect-27
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
Combustion products
Stator Rotor
Average radial
temperature profile
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
Possibility of
shock-boundary
layer interaction
Unsteady wake
Stagnation point flow
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-27
In this lecture...
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
1
Lect 28
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
2 Where,
T01 +T02 U mean
T0 -bl = - (1 - 2.DR) DR= Rx
2 2.c p-gas
Heat transfer coefficient = Where,
t time ;
Quantity of heat transferred T Temp
surface area x t x T between hot gas & surface difference
Heat Transfer
The heat transfer is mainly by surface convection,
conduction and then internally mainly by forced convection.
Radiation heat transfer is negligible.
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
Internal
Cooling
air
Hot spot
Hot spot at LE passages
at TE
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
CFD
simulation
Internal
Cooling
air
passages
Blade
Fabrication
Distributed
Internal
Blade cooling
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Turbine Cooling
Lect 28
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
Need for
high turbine
temperature
was reduced
due to high
compressor
pressure
ratio
Advanced
cooling has
extended
both TET and
Compr ratio
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 28
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
1
Lect - 29
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Selection of
1) Design requirements *
T , T*g-max, *exit, M*2
2) Design constraints : for both HPT & LPT
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Important
Aerodynamic
Parameters for
blade airfoil
selection
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Airfoil Geometry
Selection
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Airfoil pressure
distribution
applied as per
requirement to
deviate from a
starting profile
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Supersonic
Turbine airfoils
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 29
Next Class
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
1
Lect - 30
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
Selection of parameters
Design at mean diameter, Dm
Radial variation of parameters
Profiling of stator and rotor blades
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
1. Design at Dm
a) Parameters to be selected from the preliminary cycle
(thermodynamic) calculations are: (@ flt velocity & altitude),
Tg, g , Pa , Ta at design point Pg , ,TP T
* * *
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
2. Design at Dm
Matching : Peripheral velocity, u is selected from
turbine-compressor matching criteria. The outlet
velocity from the nozzle may be supersonic, but inlet
relative velocity to rotor is generally brought down
below sonic speeds. The respective sonic speeds are:
At rotor exit
a3 =RT 3
5 5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
-1
Twt
T0t + 1
(
= 1 - . ..sin
2. u -2 2
2
u )
where
2-abs cos2 Twt 2-abs =C2 /acr
and = 2-rel =V2/acr
2-rel cos2 T0t u =U/acr
6 6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
M
2-rel
=J
1 ( M 2 - M 2u ) , where
1 - sin 1
J1 = 1 +
M1 M1u
+
M1u M1
-1
2
= 1.2 to 1.4, is a flow coefficient
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
For moderate pressure drops ( <2.0) the flow in the rotor may
be entirely subsonic. However for high pressure drops ( >2.5),
the flow becomes transonic at the stator trailing edge.
Stator Exit flow conditions:
From continuity, for an unit length of the blade, at throat
t .Vt .st .1 = s.1.3 .V3 .sin3eff
s = blade spacing or pitch
st = O = Throat width, Subscript t for throat
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
st
sin3-eff = Throat area ratio
s
The exact relationship between and st / s can be found
2
experimentally by accurate cascade analysis
st
3 = sin k 2 -1 Initially assume K2 = 1
s
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
a) Actual Rx ,
hrotor hrotor V 2 -V 2
Rx = = = 3 2
hrotor + hstator h 2h
T T
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
12 12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
v) m g
A2 =
.C2 a
From which A2
Blade height h2bl =
Dm2
h2
passage
r= h2
bl
( 1+ )
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
vi)
( D m2 h2
bl
)-1
(D )+ 1
d=
m2 h2
bl
C 2w - U
vii) 2 = tan -1
C 2a
Find V2 and 2-rel
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect - 30
A modern turbine
blade obtained through
design, optimization
and stress calculations
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 31
1
Lect-31
In this lecture...
Centrifugal compressors
Thermodynamics of centrifugal
compressors
Components of a centrifugal
compressor
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Centrifugal compressors
Centrifugal compressors were used in the
first jet engines developed independently
by Frank Whittle and Hans Ohain.
Centrifugal compressors still find use in
smaller gas turbine engines.
For larger engines, axial compressors
need lesser frontal area and are more
efficient.
Centrifugal compressors can develop
higher per stage pressure ratios.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Centrifugal compressors
Besides small aero engines, centrifugal
compressors are used in the auxiliary
power units (APUs) in many aircraft.
Some of the aircraft air conditioning
systems employ centrifugal compressors.
In a few engines, centrifugal compressors
are used as the final stage of the
compression process downstream of a
multi-stage axial compressor. Eg. GE T
700, P&W PT6, Honeywell T53.
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
2
1
T01 C12 01
2c p
T1
1
s
T-s diagram for a centrifugal compressor
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
b U2
U1
r2
r1
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Conservation of Rothalpy
If we were to assume steady, viscous flow
without heat transfer
C12 C 22
h1 + U1C w1 = h 2 + U2 C w 2 = I
2 2
Here, I, is the rotational enthalpy or rothalpy.
It is now known that rothalpy is conserved for
the flow through the impeller.
Any change in rothalpy is due to the fluid
friction acting on the stationary shroud (if
considered in the analysis).
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Impeller
Impeller draws in the working fluid. It is
the rotating component of the centrifugal
compressor.
The diverging passages of the impeller
diffuses the flow to a lower relative
velocity and higher static pressure.
Impellers may be single-sided or double-
sided, shrouded or un-shrouded.
In the impeller, the working fluid also
experiences centripetal forces due to the
rotation.
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Impeller
Impeller
2 2
U2 U2 U2
V2 V2 V2
C2 2 C2 C2
2
Forward leaning blades Straight radial Backward leaning blades
(2 is negative) (2 is positive)
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Inducer
Inducer
m m
rt
rm
rh
Inducer Ut Um Uh
t Ct
Vt
h
m
Leading edge velocity triangles
U1
1 C1
V1
Section m-m
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Inducer
The diffuser
High impeller speed results in a high
absolute Mach number leaving the impeller.
This high velocity is reduced (with an
increase in pressure) in a diffuser.
Diffuser represents the fixed or stationary
part of the compressor.
The diffuser decelerates the flow exiting the
impeller and thus reduces the absolute
velocity of the working fluid.
The amount of deceleration depends upon
the efficiency of the diffusion process.
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
The diffuser
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
The diffuser
Diffuser vanes
Vaneless space
r3
r3>r2>r1
r2
r1
Impeller
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
The diffuser
Logarithmic spiral
Radial direction
Cr
C
CW
The diffuser
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
Centrifugal compressors
Coriolis acceleration
Slip factor
Performance characteristics
Stall and surge
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-31
In this lecture...
Centrifugal compressors
Thermodynamics of centrifugal
compressors
Components of a centrifugal
compressor
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 32
1
Lect-32
In this lecture...
Centrifugal compressors
Coriolis acceleration
Slip factor
Performance characteristics
Stall and surge
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Coriolis acceleration
Coriolis acceleration
dCw
dr Vd
dC V
C C
d
r
(r+dr)
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Coriolis acceleration
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Coriolis acceleration
The existence of the tangential pressure
gradient means that there will be a positive
gradient of V in the tangential direction.
1 dP
=
d V2 / 2
=
(
V dV )
rd rd r d
1 dV
Therefore, = 2
r d
This means that there will be a tangential
variation in relative velocity.
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Coriolis acceleration
Cw2
C2
V2
U2
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Slip factor
Towards the outlet of the impeller, as the
Coriolis pressure gradient disappears, there
will be a difference between Cw2 and U2.
This difference in the velocities is expressed
as slip factor, s = C w 2 / U2
The slip factor is approximately related to
the number of blades of the impeller.
For a straight radial blade, the slip factor is
empirically expressed as s 1 2 / N, where
N is the number of blades.
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Slip factor
As the number of blades increases, the slip
factor also increases and thus the slip lag at
the tip of the impeller reduces.
The effect of slip is to reduce the magnitude
of swirl velocity and therefore the pressure
ratio.
The presence of slip means that to deliver
the same pressure ratio, either the impeller
diameter or the rotational must be
increased.
This in turn may lead to either increase in
frictional losses or stresses on the impeller.
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
The centrifugal compressor performance
characteristics can be derived in the
same way as an axial compressor.
Performance is evaluated based on the
dependence of pressure ratio and
efficiency on the mass flow at different
operating speeds.
Centrifugal compressors also suffer from
instability problems like surge and
rotating stall.
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
The compressor outlet pressure, P02, and
the isentropic efficiency, C, depend upon
several physical variables
P02 , C = f(m
, P01 , T01 , , , R, , design, D)
In terms of non - dimensionl ess parameters,
P02 m RT01 D D 2
, C = f , , , , design
P01 P01D 2 RT
01
P02 m T01 N
The above reduces to , C = f ,
P01 P01 T
01
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
Surging limit
B
D
Constant speed line
Pressure ratio
A E
Choking limit
Mass flow C
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
Surge line
P02 Locus of points of
P01 maximum efficiency
m
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
There are two limits to the operation of
the compressor.
Operation between A and B are limited
due to occurrence of surge.
Surging: sudden drop in delivery pressure
and violent aerodynamic pulsations.
Operation on the positive slope of the
performance characteristics: unstable
Surging usually starts to occur in the
diffuser passages.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
The pressure ratio or the temperature rise in
a centrifugal compressor also depends upon
the blade shaping.
There are three possible types of blade
shapes: forward leaning, straight radial and
backward leaning.
Theoretically, the forward leaning blading
produces higher pressure ratio for a given
flow coefficient.
However such a blading has inherent dynamic
instability.
Therefore, straight radial or backward leaning
blades are popularly used.
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Performance characteristics
Backward leaning
or Cr 2
m U
2
Performance characteristics for different
blade geometries
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Rotating stall
Rotating stall might also affect the
compressor performance.
In this case a stall cell (that might cover
one or more adjacent blades) rotates
within the annulus.
Full annulus rotating stall may eventually
lead to surge.
Rotating stall may also lead to
aerodynamically induced vibrations and
fatigue failure of the compressor
components.
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
Rotating stall
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
In this lecture...
Centrifugal compressors
Coriolis acceleration
Slip factor
Performance characteristics
Stall and surge
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-32
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 33
1
Lect-33
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 1
At the inlet of a centrifugal compressor
eye, the relative Mach number is to be
limited to 0.97. The hub-tip radius ratio of
the inducer is 0.4. The eye tip diameter is
20 cm. If the inlet velocity is axial,
determine, (a) the maximum mass flow
rate for a rotational speed of 29160 rpm,
(b) the blade angle at the inducer tip for
this mass flow. The inlet conditions can be
taken as 101.3 kPa and 288 K.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 1
U1
1 C1
rt V1
rm
rh
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 1
The rotational speed at the inducer tip is
U1 = dN / 60 = 0.2 29160 / 60 = 305.36m / s
From the velocity traingle, we can see that
V1 C12 + U12
M1rel = =
RT1 RT1
T1 = T01 C12 / 2c P = 288 C12 / 2010
C12 + U12
M1rel =
R(288 C12 / 2010)
C 2
+ 305 .63 2
0.97 2 = 1
115718.4 0.2C12
Simplifyin g, C1 = 114.62m / s
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 1
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 1
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 2
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 2
C1 25 V1
Cw1 U1
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 2
Exit stagnation temperature is
T02 = T01 ( c ) = 293(4)
( 1) / (1.4 1) / 1.4
= 435.56K
Therefore the isentropic temperature rise,
T0s = 435.56 293 = 142.56K
The actual temperature rise, T0 = T0s / c
T0 = 178.2K
Work done per unit mass is, w = c P T0
w = 1.005 178.2 = 179 kJ/kg
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 2
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 2
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 3
Air at a stagnation temperature of 22C enters
the impeller of a centrifugal compressor in the
axial direction. The rotor, which has 17 radial
vanes, rotates at 15,000 rpm. The stagnation
pressure ratio between diffuser outlet and
impeller inlet is 4.2 and total-to-total efficiency
is 83%. Determine the impeller tip radius.
Assume the air density at impeller outlet is
2kg/m3 and the axial width at entrance to the
diffuser is 11mm, determine the absolute Mach
number at that point. Assume that the slip factor
= 1 - 2/N, where N is the number of vanes.
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 3
The specific work required is
w c = U2 C w 2 U1C w1
Since C w1 = 0, w = U2 C w 2 = U22
Expres sin g U2 in terms of efficiency and pressure ratio,
c p T01 ( (c 1) / 1)
U22 =
c tt
= 1 2 / N = 1 2 / 17 = 0.8824
Substituti ng all other values, U2 = 452 m / s
Since, = 15000 2 / 60 = 1570 rad / s
Therefore, the impeller radius is
rt = U2 / = 0.288 m
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 3
Mach number, M2 = C 2 / a2 = C 2 / RT2
where, C 2 = C w2 2 + Cr22
Cr 2 = m
/( 2 2rtb 2 ) = 2 /(2 2 0.288 0.011) = 50.3 m / s
C w 2 = U2 = 400 m / s
C2 = 50.3 2 + 400 2 = 402.5 m / s
We know that h02 = h01 + w c = h01 + U22
or, h2 = h01 + U22 12 C 22
or, T2 = T01 + (U22 12 C 22 ) / c p
= 394.5 K
Therefore, M2 = 402.5 / 1.4 287 394.5 = 1.01
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 4
A centrifugal compressor with backward leaning
blades develops a pressure ratio of 5:1 with an
isentropic efficiency of 83 percent. The
compressor runs at 15000 rpm. Inducers are
provided at the inlet of the compressor so that
air enters at an absolute velocity of 120 m/s.
The inlet stagnation temperature is 250 K and
the inlet air is given a pre-whirl 22o to the axial
direction at all radii. The mean diameter of the
eye of the impeller is 250 mm and the impeller
tip diameter is 600 mm. Determine the slip
factor and the relative Mach number at the
impeller tip.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 4
U2
2 V2 Ca
C2
Inducer
U1
1
V1 C1
Fixed inlet
guide vane
Solution: Problem # 4
T01 = 300 K
T02s = T01 ( c )
( 1) /
= 250 (5)0.4 / 1.4 = 395.95 K
T0s = 395.95 300 = 95.95 K
Actual temperature rise, T0actual = T0s / c = 95.95 / 0.83
= 115.6 K
The specific work required, w c = c p T0actual = 1005 115.6
= 116.186 kJ / kg
Given that C1 = 150 m / s, C w1 = C1 sin 1 = 150 sin 22
= 56.2 m / s
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Solution: Problem # 4
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-33
Problem # 4
U2
2 V2 Ca
C2
Inducer
U1
1
V1 C1
Fixed inlet
guide vane
Solution: Problem # 4
From the impeller exit velocity triangle,
V2 = C a2 + (U2 C w 2 )2 = (C1 cos 1 )2 + (U2 C w 2 )2
= 222.9 m / s
Mrel = V2 / RT2
T2 = T02 C 22 / 2c p
T02s T01
T02 = T01 + = 365.61K
c
and C 2 = C w2 2 + C a2 = 269.9 2 + 139.08 2 = 303.68 m / s
T2 = 365.61 303.68 2 / 2 1005 = 319.73 K
The relative Mach number at the impeller tip is
Mrel = 222.9 / (1.4 287 319.73) = 0.62
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
1
Lect 34
Centrifugal Compressors
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
a) Impeller
b) Diffuser
Vanes
c) Vaneless
diffuser
e) Inlet Guide
vanes
f) Volute
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
Other
important
issues to be
designed
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
The velocity
Diagrams
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
Slip factor
In a real compressor relative velocity vector V2 is at
angle 2 because of non-radial exit from the impeller tip
as the real viscous flow detaches near the tip from the
impeller vane (trailing) surface
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
In highly
forward and
highly backward
curved
(2 >-600 )
impellers
slip factor
looses its
meaning
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
Diffusion Limit :
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
Vaneless diffuser
The theoretically
obtained points to
the right of b are
considered
choked, i.e. the
compressor
cannot process
greater mass
flows. The
compressor is
said to go in to
stall at m a , this
happens when
high pressure rise
is attempted at
low mass flow With Mass Flow Control only
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
In aircraft engines,
rotating speed is
variable during
actual running.
Thus the zone of
operation is bounded
between the points
a,b,c and d .
The min lines and
the speed lines, nmax
and nmin, define the
boundaries (shaded
area) of operation.
With Speed Control and Flow Control
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 34
Next Lecture
Radial Turbines
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 35
1
Lect-35
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Radial turbines
Radial turbines
There are two types of inward flow radial
turbines
Cantilever turbine
90o IFR turbine
Cantilever turbine
Similar to the impulse type turbine
Little change in relative velocity across the rotor
Aerodynamically very similar to the axial impulse
turbine
Can be designed in a similar manner as axial
turbines.
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Radial turbines
1
Nozzle blades
2
Rotor blades
3
Flow
C2 V3
V2 Ca3
U2 U3
Radial turbines
90o IFR turbine
This turbine has a striking similarity with a
centrifugal compressor.
The flow direction and blade motion are
reversed.
The flow enters the turbine radially and exits the
turbine axially.
Straight radial blades are generally preferred as
curved blades would incur additional stresses.
The rotor or impeller ends with an exducer.
Usually the flow exiting the rotor passes through
a diffuser to recover KE, which would otherwise
be wasted.
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Radial turbines
Volute/Scroll
1 Nozzle blades
2 Rotor blades
3 4 Flow
Diffuser
C2 2 3 V 3
V2 Ca3
U2 U3
s
T-s diagram for an IFR turbine
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Radial turbines
Volute/Scroll
1 Nozzle blades
2 Rotor blades
3 4 Flow
Diffuser
C2 2 3 V 3
V2 Ca3
U2 U3
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Nominal design
Defined by a relative flow of zero
incidence at the rotor inlet (V2=Cr2).
An absolute flow at the rotor exit that is
axial (C3=Ca3).
Therefore, with Cw3=0 and Cw2=U2, the
specific work for nominal design is
W=U22 .
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
s
T-s diagram for an IFR turbine
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
Incidence losses
At off-design conditions, the fluid is likely to
enter the rotor at a relative flow angle different
from the optimum angle.
This leads to an additional loss component due
to incidence angles.
Often defined as equal to the kinetic energy
corresponding to the component of velocity
normal to the rotor vane at inlet.
There is an increase in entropy and hence a
corresponding loss in enthalpy due to incidence.
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
2: after 4 P3 C 3
V2 2, opt 2c p
entry to
rotor 3
3s
3ss
s
T-s diagram for an IFR turbine
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-35
In this lecture...
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 36
1
Lect-36
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Problem # 1
The rotor of an IFR turbine, which is
designed to operate at the nominal
condition, is 23.76 cm in diameter and
rotates at 38,140 rpm. At the design
point the absolute flow angle at rotor
entry is 72o. The rotor mean exit
diameter is one half of the rotor
diameter and the relative velocity at
rotor exit is twice the relative velocity at
rotor inlet. Determine the specific work
done.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 1
Volute/Scroll
1 Nozzle blades
2 Rotor blades
3 4 Flow
Diffuser
C2 2 3 V 3
V2 Ca3
U2 U3
Solution: Problem # 1
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 1
We can sum up the three terms and divide
by 2 to get the specific work as
Problem # 2
A radial inflow turbine develops 60 kW
power when running at 60,000 rpm. The
pressure ratio (P01/P3)of the turbine is
2.0. The inlet total temperature is 1200
K. The rotor has an inlet tip diameter of
12 cm and an exit tip diameter of 7.5 cm.
The hub-tip ratio at exit is 0.3. The mass
flow rate is 0.35 kg/s. The nozzle angle is
70o and the rotor exit blade angle is 40o.
If the nozzle loss coefficient is 0.07,
determine the total-to-static efficiency of
the turbine and the rotor loss coefficient.
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 2
Volute/Scroll
1 Nozzle blades
2 Rotor blades
3 4 Flow
Diffuser
C2 2 3 V 3
V2 Ca3
U2 U3
Solution: Problem # 2
The rotor tip rotational speed is
U2=D2N/60=377 m/s
From the velocity triangle at the rotor
inlet, 2=0, therefore,
sin 2 =U2/C2
C2=U2 cosec 2 = 401.185 m/s
T2=T02-(C22/2cp)=1130 K
To calculate the stagnation temperature
drop (isentropic) across the turbine, we
shall use the pressure ratio.
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 2
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 2
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 2
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Problem # 3
An IFR turbine with 12 vanes is required to
develop 230 kW at an inlet stagnation
temperature of 1050 K and a flow rate of 1
kg/s. Using the optimum efficiency design
method and assuming a total-to-static
efficiency of 0.81, determine (i) the
absolute and relative flow angles at rotor
inlet; (ii) the overall pressure ratio, P01/P3;
(iii) the rotor tip speed and the inlet
absolute Mach number.
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 3
For optimum design, we use the Whitfields
equation, cos22=1/N, where, N is the
number of vanes.
2=2(90- 2)=33.56o
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 3
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 3
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 3
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Problem # 4
Compare the specific power output of axial
and radial turbines in the following caseL
Axial turbine : 2=3=60o and 3=2=0o
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 4
U2=U3=U
2 3
C3 V2
C2
V3
Axial turbine
C2 2 3 V3
V2 Ca3
U2 U3
Radial turbine
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-36
Solution: Problem # 4
Axial flow turbine:
Since 2=3=60o and 3=2=0o
The specific work is
Waxial= U(Cw2+Cw3)=U2
Radial flow turbine:
2=60o and 3=3=2=0o
The specific work is
Wradial= U2Cw2-U3Cw3)= U2U2-U3x0
=U22
Therefore, the specific work done in both the
turbine configurations are the same, given the
conditions of operation.
20
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
1
Lect 37
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
Thermodynamic
parameters of
Radial Flow Turbine
(a)
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
Aerodynamics of
Radial flow Turbine
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
Losses in
Radial
turbine
components
Radial Turbine
Characteristics
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
Compressor
Turbine
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 37
Next class
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 38
1
Lect-38
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Grid/mesh
Structured, unstructured and hybrid grids
Structured grid
More suited for well-defined geometries
More difficult to generate
Easier to control near-wall clustering of cells
Unstructured grid
Primarily intended for complex geometries
Easier to generate
Not much control over the near-wall clustering of
cells
Easily automated
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Blocks
Grid Generation
Unstructured grid
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Multi-block structured grid
In order to generate structured grid
over curved surfaces, multiple blocks
need to be defined.
Interface of the blocks need to be
carefully managed.
Grid topology needs to be
appropriately defined.
The Grid topology also needs to
account for the change in geometry of
the blade from hub to tip.
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Topology
Is a structure off blocks that acts as a
framework for placing mesh elements.
Blocks are laid out without gaps with
shared edges and corners.
Blocks contain same number of
elements along each side.
Is usually invariant from hub to tip.
Can be edited on 2-D layers from hub to
tip sections.
Number of blocks will dictate the
skewness of the grid elements.
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Grid topology schemes
O-grid:
Usually used around the blade by
forming a continuous loop around it
Yields excellent boundary layer
resolution
gives good control over the y+
values that needs to be tightly
monitored
Provides near orthogonal elements
on the blades
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
O-grid topology
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
J-grid:
Usually used near leading and trailing
edges
Wraps up in opposite directions at the
leading and trailing edges
H-grid:
Tends to complete the meshing by
adding some blocks in an unstructured
manner
The structured blocks extend from
upstream of the LE, downstream of
the TE and between the blades and
the periodic surfaces
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Grid Generation
Other topology options include C-grid
and L-grid.
These are also often used at the
leading and trailing edges.
All the above grid topologies are used
along with an O-grid for proper
resolution of the boundary layer.
Proper resolution of the leading and
trailing edge radii are important.
Establishing grid-independence or
grid-insensitivity of the results is now
a standard practice.
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Grid Generation
Boundary conditions
To capture the flow physics correctly, the
boundary conditions must be set
appropriately.
Quality of the solutions is a strong
function of the boundary conditions.
Turbomachinery flows
Inlet boundary
Exit boundary
Periodic boundary
Walls or surfaces
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Boundary conditions
Boundary conditions
Inlet boundary conditions
Depends upon the application
Flow conditions (incompressible or
compressible)
Total pressure, total temperature,
velocity components/profile (most
commonly used)
There are other forms of specifying
the inlet boundary conditions: velocity
inlet, mass flow inlet etc.: not
commonly used due to several
limitations.
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Boundary conditions
Exit boundary conditions
Exit static pressure to achieve the
required mass flow
It is also possible to specify a static
pressure distribution at the exit domain.
Alternatively, mass flow can be directly
specified at the exit.
For incompressible flows, using either of
the two does not affect the results.
However, for compressible flows, static
pressure outlet condition yields better
results.
26
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Boundary conditions
For single passage simulations, periodic
boundary conditions are used for
simulating the effect of a blade row.
The domain must be appropriately
chosen to ensure that periodic boundary
conditions are indeed valid.
On surfaces (blade, hub and shroud),
no-slip and adiabatic conditions are
usually used.
In turbines with hot gases present, the
adiabatic condition may be replaced by
constant heat flux condition.
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-38
Boundary conditions
In this lecture...
29
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect- 39
1
Lect-39
In this lecture...
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence models
Reynolds averaging of the Navier-Stokes
equation can be expressed as time, space
or ensemble.
The main idea of Reynolds averaging is to
decompose the mean and fluctuating
components: Reynolds decomposition.
However this introduces additional
variables for which there are no available
relations: closure problem.
Modelling this is a major challenge in CFD
even today.
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence Models
Zero equation or algebraic eddy viscosity
model
Use an algebraic form for the turbulent
stresses
Valid for simple 2D shear flows
Mild pressure gradient
3D boundary layers with small cross
flows
Not accurate for flows with pressure or
turbulence driven secondary flows
Cannot predict shock-induced separated
flows
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence Models
One equation model: Spalart Allmaras
Employ an additional PDE for a
turbulence velocity scale
Usually used in design-iteration type
simulations
Popular in recent times due to
inherent problems with more refined
models
Very robust models, rarely produce
completely unphysical results
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence Models
Two equation models:-,- models
One PDE for turbulence length scale and
one PDE for velocity scale
Good for 2D flows with moderate pressure
gradients
Not satisfactory for flows with rotation,
strong swirl, and separated flows
Modified two equation models: improved
results
Modified via ARSM
Coupled -/ARSM
Realizable -
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence Models
Near wall treatment
On-design flows without large separated
regions, wall function model close to the
wall
Off-design, low Re model, over
production of turbulent KE must be
checked
Menters SST - and Durbins v2f
Works well for adverse pressure
gradients and separating flows
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Turbulence Models
Transition prediction
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Case studies
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Measured
Passage vortex
Passage vortex
Passage vortex
Tip flow
Significant effect on flow stability,
pressure rise, and efficiency
Self induced unsteadiness related to
spike initiated stall?
Role of shock wave in the flow physics
Experimentally capturing the tip flow
dynamics very challenging
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Secondary flows
Secondary flow losses: significant portion of the
total losses
Accurate prediction of secondary flow: reduction
in total losses
Literature on secondary flows: mostly cascades
Case 2: Yu (2004)
Turbine flow simulations using structured and
unstructured grids
In-house code (Penn State)
Secondary flow structures well captured
Validation with experimental data
Tran. ASME, Journal of turbomachinery, 2004
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
Secondary flows
Secondary flows
Computed Experimental
Aeroelasticity
Computing requirements
Sources
Type of simulation Systematic
procedure for
Geometry errors
estimating these
Modeling errors errors like say the
Boundary ASME / AIAA
conditions standards for
Numerical errors experimental
Discretization uncertainty
Round-off analysis
Convergence
27
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect-39
In this lecture...
29
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
1
Lect 40
Fundamentals of CFD
for use in
Turbomachinery Analysis
2
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
3
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
4
4
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
2 2 2 2
A 2 +B +C 2 + D 2
x x .y x y
A,B,C,D,E,F,G are functions of x,y & +E + F + G =
0
y
Assume that f = f (x,y) is a solution of the above
differential equation.
This solution, typically is a surface in space, and
the solutions produce space curves called
characteristics.
2nd order derivatives along the characteristics are
often indeterminate and may be discontinuous
across the characteristics.
The 1st order derivatives are continuous.
5
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
2
dy B B - 4AC
dx =
2A
6
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
7
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
8
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
9
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Computational space
11 10
10
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Domain
Transformation
11
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Orthogonal
Grid
Grid without
Orthogonality
12
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Unstructured Grid
generation
13 13
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
14
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
15
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Three-Dimensional Flow
Analysis Full Blade Structural and
Aero-elastic analysis
Blade Construction
16
16
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
17
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
18
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
19
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Blade-to-Blade program
Input : Blade geometry
Inlet and Exit Velocity
distribution
Streamline Definition
Output : Surface velocity distribution
Profile and loss distribution
Section Stacking Program
Blade-to-Blade program
2D MISES code for
Cascade Analysis
19 21
21
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
22
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Three-Dimensional
Flow Analysis
Solution of :
1) Continuity Equation
2) Momentum equation
3) Energy Equation
4) Radial Equilibrium
Equation
23
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
1800 Speed
Lines
P1600
0 (Pa) 90%
100%
80%
1400
1200
1000
Mass Flow rate kg/s
1.6 1.8 2 2.2 2.4 2.6 2.8 3
24
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay
Lect 40
Thank you
for
participating in this
NPTEL course
25
Prof. Bhaskar Roy, Prof. A M Pradeep, Department of Aerospace, IIT Bombay