Axial Flow Machines
Axial Flow Machines
Axial Flow Machines
further concepts
Euler Turbine Equations Axial Flow Gas Turbines Axial Flow Compressors Losses / Zweifel Free-Vortex Flow Blade Stresses Re: Peng Ch. #8 & #7
Stator
Rotor
Stator
TURBINE STAGE
STATOR
STAGE
V1
W1
ROTOR
V2
W2
STATOR
AXIAL TURBINE
V2
W2
V1
W1
100
percent reaction
U
stator rotor
TURBINE STAGE FLOW ANGLE TERMINOLOGY ALPHA () ~ STATOR ANGLES BETA () ~ ROTOR ANGLES
AXIAL DIRECTION
STATOR
Negative Angle V1
U
W1
Positive Angle 0o
ROTOR
V2
W2
STATOR
W = m ( VU1U1 VU2U2) / gc
Energy per unit mass
Right end: free wheeling , runaway condition. No torque on rotor.no power no use. Left end: Lots of torque but rotor velocity goes to zero. Since power = torque x angular velocity = 0 7
V2
W2
V1
W1
U
V2
W2
V1
W1
MAXIMUM UTILIZATION Work per unit mass = W / m = U (VU1 - VU2 ) For maximum work per unit mass: i.e. V2 is in axial direction only VU2 = 0
MAXIMUM UTILIZATION
Work per unit mass = W / m = U (VU1 - VU2 )
W2 V2 = Vx
V1
W1
U S R
10
V2
W2
V1
W1
U
MAXIMUM UTILIZATION AXIAL TURBINE
11
EXAMPLE
An axial-flow turbine has a mean diameter of 40 cm and a flow coefficient ( = Va/U) of one-half. Flow enters the stator in the axial direction with a velocity of 50 m/s and exits the stator with a nozzle angle (, measured between V1 and U) of 17 degrees. The fluid flows through the stage with constant axial velocity and produces an output torque of 1600 Nm. a) Sketch the velocity vector diagram for this stage. b) Calculate the power output of the turbine stage. c) What is the degree of reaction?
12
V2 = 50
W2
= 17o
V1
W1
13
SOLVING FOR VELOCITY VECTORS: = Va/U = VU1 = 50 tan 73o = 164 m/s WU1 = VU1 U U = 100 m/s VU1 = 164 m/s V1 = 171 m/s VU2 = 0 WU2 = U = 100 m/s W2 = 112 m/s W1 = 80.85 m/s EULER EQNS: Torque = m (VU1 R1 VU2R2) m = 1600 / (164)(0.2) = 49 kg/s
14
COMPRESSOR STAGE
STATOR
STAGE
V1 W1 U ROTOR
W2 V2 STATOR V1
15
AXIAL COMPRESSOR
V1
100
V2 W1
percent reaction
W2
0
U
stator rotor
16
V1
100
V2
W1
percent reaction
W2
0
U
stator rotor
17
V2 V1
100
W2 W2
0
V2
W1
percent reaction
U
stator rotor
18
Allied Signal ASE120 Gas Turbine Engine Phoenix, AZ Designed for cogeneration, emergency power and mechanical drive. ~35% thermal efficiency ~10 MW m-dot~34kg/s Texhaust ~500oC NOx ~ 10-25 ppm Weight ~ 5,000 lb
19
16,000
Power
8,000
-20
50
(NOT RELATED TO UNRECOVERABLE PRESSURE LOSSES LECTURE) ASE120 engine performance as a function of inlet air temperature Notice that maximum output power decreases about 50% when the air let increases from ~ -20C to +50C (~16,000hp to 8,000 hp)
20
21
IDEAL TURBINE
V2
V1
W2
W1
U S R
22
IDEAL
ACTUAL
CU
Actual performance is always less than the ideal theoretical performance due to un-recoverable losses
23
P = C2 / gc
m = A 2Pgc
discharge coefficient
CD
CD =
CD
Turbine performance follows the general behavior of ideal flow through a nozzleTot max delP ~ den x sq vel (Eqn W 2.12) Incompressible flow, conservation of energy-ss, no change in PE P = Po Pst = total static = dynamic pressure Dont forget gc when using USCS units. Notice the total maximum pressure drop ~ density x velocity2 Mass flow ~ flow area x rho x delP CD ~ discharge coefficient: correction coefficient for real fluids with viscosity, friction losses. Usually empirically based, rooted in experiment.
24
CD =
2FD V2A
In addition to classic profile / skin friction drag.. VIBRATIONS: fluttering / vortex shedding esp. problem w/ thin bladesfans, some compressors
25
FLOW
High Reaction
(no separation)
Low Reaction
(separated flow)
Separation leads to COMPRESSOR LOSSES STALL losses in efficiency. Separation leads to losses in pressure, enthalpy It is very complicated to accurately predict the conditions in which B.L separation occurs. Therefore, empirically-based rules-of-thumb have been developed. As well see, flow separation can be minimized with a cascade of blades..i.e. form flow passages
26
CD ~ 0.02 (s/h)
Due to boundary layer formation within the annular flow passage. Empirical correlation.
27
End wall boundary layers are convected inward along the suction side of the blades. Secondary currents are set up in a plane transverse to the flow, dissipating energy. Results in un-favorable flow redistribution, competing vortices, and loss in stagnation pressure. Can be a significant fraction of total P-loss. These currents also take place in the wake regions down stream from the trailing edges of the blades.
28
CD = 0.018CL2
The constant 0.018 stems from a complex relationship between flow acceleration, aspect ratio (h/c), spacing & pitch. RE: Logan Ch 6 Because the trailing vortices are similar to wing vortices, it is expected that the corresponding drag is proportional to the LIFT COEFFICIENT CL.
29
30
TIP LEAKAGE
TIP LEAKAGE
CD = 0.29 (k/h)CL3/2
The difference in pressure on the two sides of the moving blades results in leakage around the tip. Shroud reduces leakage ~50% k=clearance gap = f(T, rpm, loading) h= blade length CL=lift coefficient
31
32
GENERATOR
LPT1
LPT2
HPT
Pressure changes across casing boundaries are another source of losses. Each place where the rotor penetrates a pressure boundary needs a seal. Typically LABYRINTH SEALS.
33
LABYRINTH SEALS
34
35
LABYRINTH SEALS
CD = f(P, h, , D, .)
36
Clearances change with thermal expansion, radial stresses, axial stresses, fluid loadings, rotor flex, rotor unbalance, startup/shutdown transients,. Rotor alignment is critical for rotor seals.
37
stage
= U2
= (R - ) / (+ R) + (1 R - ) / ( + (1 R))
rotor
stator
= flow coefficient = Cx / U = blade loading coefficient = Cu / U R =degree of reaction = drag to lift ratio U = rotor velocity P = pressure gain
LOGAN38 P115
39
EXAMPLE: Axial compressor with Rn=50% and the velocity vectors shown. The lift and drag coefficients are the same for the rotor and stator and are 1.4 and 0.1, respectively. CD includes all losses. Calculate the stage pressure rise and required energy input if the entering air has a density of 1.2 kg/m3. Each square corresponds to 50 m/s x 50 m/s
50 m/s
V1
V2
W1
W2
U stator rotor
40
EXAMPLE (cont)
stage
= U2 = Pactual
Vx = 150 m/s VU2 = 250 m/s
= (R - ) / (+ R) + (1 R - ) / ( + (1 R)) =
(0.5-0.5x0.071)/(0.5 + 0.071x0.5) + (1-0.5-0.5x0.071)/ (0.5+0.071(1-0.5))
41
EXAMPLE (cont)
stage
= U2 = Pactual
2 2
42
FLOW SEPARATION
de Haller Number: dH
dH = Wexit / Win = W2 / W1 dH > 0.72 No Separation
Swiss engineer Established criterion for flow separation in compressor cascades (for compressor cascades, RE: Wilson, Ch#4 p.182) Butalso sometimes applied to turbine blades as well using V1 & V2
43
V1
V2
W1
W2
U S R
44
V1
V2
W1
W2
U S R
45
Consider each stage of a compressor as a compressor The pressure rise at each stage can only increase so much, Limits on pressure rise Cpr ~< 0.5
47
DEGREE OF REACTION Rn = hrotor / hstage = (h2 h3) / (h1 h3) = Protor / Pstage = (P2 P3) / (P1 P3) = 1 (Vu1 + Vu2) / 2U
(Text eqn 5.8)
Vu2
Vu1
V2
V1
W2
W1
48
Stator
Rotor
Stator
49
FL = CL
(AV2) / 2gc
FD = CD
(AV2) / 2gc
Comments on LIFT and DRAG forces: Lift and drag forces are difficult to model exactly due to complex flow phenomena including: viscous effects, boundary layer separation, etc. Hence, empirically determined correction coefficient is typically used. CL and CD are lift and drag coefficients, usually determined empirically. CD typically accounts for both pressure and skin friction drag. A is some characteristic area. V is some representative velocity: free stream, average or mean velocity, etc.
50
Cascade of Blades
CLcas=K(CL)
S C K=cascade coefficient = f(beta, s/c) S= pitch or spacing C = chord
51
52
c
Low solidity Losses by separation High solidity ~No separation Losses by skin friction
53
54
55
1 2
56
CL,op ~ 0.8
57
EXAMPLE The spacing (pitch) between turbine or compressor blades influences the machines performance. Shown below is a cascade of blades for an axial compressor.
- 30o
Flow in
c = 12 cm b = 10 cm
Optimum Spacing
= + 60o
Zweifel Correlation
(b/s)
opt
58
AXIAL TURBINES NUMBER OF BLADES Wind: Water: Gas / Steam: 1-24 3-30 11-110
Rotors usually have an even number of blades (for balance) Stators usually have an odd number of blades (for vibration)
59
VORTEX FLOW
Radial Equilibrium
Tornado
A tornado is natures efficient vortex flow. Characteristic of air / gas at higher velocities (200+ mph ~ 300 ft/s ~ 100 m/s) It pumps air from the higher pressure ground to overhead cloud. We can capture that flow characteristic in the design of axial compressor (and turbine) blade design that encourage free vortex flow.
60
Velocity Distribution
in a Tornado
R = eye of tornado
Vu = r Vurn = constant
61
62
Patm
R -P
Minimum pressure is in center of eye. The entire field pressure is sub-atmospheric pressure (AKA Under-pressure)
63
Fc
Fc = mac = mv2/r
Fc = centripetal force (in +) ac = centripetal acceleration v = tangential speed = vu
F = ma
= angular position (rad) = angular velocity (rad/s) ac = centripetal acceleration (m/s2) 2 radians = d/dt = time for one revolution 2 = 2r/v = v/r
ac = d2/dt2 = v2/r
64
P + dP P r d
P + dP dr
Fc
For Radial Equilibrium: The radial components of the pressure forces (F=PA) must be balanced by the centrifugal force on the fluid mass.
65
Free-Vortex Flow:
Vur = constant
66
TIP
HUB
67
TIP
Vu1
HUB
Constant axial velocity Axial V2 = maximum utilization Cur = const rVu = const
68
69
TIP (RXN)
HUB (IMP)
70
Station A is entry to nozzle / stator Station B is entry to blade / rotor Station C is exit from blade Nozzle blade is twisted to give overall free vortex flow. The Blade is straight. Nozzle sets up flow in radial decreasing velocity profile.upon passage over uniform rotor blade, reestablishes uniform vel profile. Illustrates other losses: Wall friction losses (skin friction) and tipclearance leakage (PsB > PsC)
71
W = m ( Vu1U1 - Vu2U2) / gc
Energy per unit mass
W / m = ( Vu1U1 - Vu2U2 ) / gc
72
TURBINE
V2
V1 W2 W1
U S R
W ~ Vu1U1 - Vu2U2
V2
V1
W2
W1
73
frozen isochromatic fringes on three-dimensional photoelastic model of gas turbine blade hub.
74
Fluid forces: Lift & Drag Centrifugal forces Thermal forces / stresses
Fluid forces: Lift & Drag Of special interest are the tangential components of the lift and drag forces since these forces directly affect the power and efficiency. Produce bending forces / stresses due to changes in fluid pressure and momentum. Total Drag force = pressure drag (form drag, shape drag) + Viscous Drag (viscous drag) Due to complex flows, losses occur. Excitation and vibrations may result from fluctuating pressure gradients. Centrifugal forces Produces radial & bending stresses (bending when centroids of all cross sections do not lie along a radial line) Limits the design length and rpm Easiest force to model Thermal forces / stresses When its not at a uniform temperature, thermal stresses arise.
75
LIFT
AXIAL
DRAG V1
Vm
V2
76
FL = CL
(AV2) / 2gc
FD = CD
(AV2) / 2gc
Comments on LIFT and DRAG forces: Lift and drag forces are difficult to model exactly due to complex flow phenomena including: viscous effects, boundary layer separation, etc. Hence, empirically determined correction coefficient is typically used. CL and CD are lift and drag coefficients, usually determined empirically. CD typically accounts for both pressure and skin friction drag. A is some characteristic area. V is some representative velocity: free stream, average or mean velocity, etc.
77
Cascade of Blades
CLcas=K(CL)
S C K=cascade coefficient = f(beta, s/c) S= pitch or spacing C = chord
78
79
80
Upstream variations in the flow can lead to blade vibrations. Can lead to high frequency, high cycle fatigue.
81
Campbell diagram illustrating example of natural frequencies of a rotor blade. Centrifugal forces on the blade causes the blade to stiffen, centrifugal stiffening and increase the natural vibration frequency with rotor rpm.
82
h R
HUB HUB
83
84
h= length of blade = 20 inches = 1.67 ft N= rate of rotation = 1800 rpm r= hub radius = 1.5 ft
85
CENTRIFUGAL FORCE: Liberated Blade If a blade broke loose - liberated, how high could it go? N = 1800 rpm
K.E
P.E.
m V2 = mgH
V = centroid velocity
86
Thermal Stresses
=E L/L= T
88
BLADE COATINGS
thermal & corrosion protection
Materials: NiCrAlY alloy ( nickel, chromium, aluminum, yttium) Zirconia (ceramic, ZrO2) Alumina (Al2O3)
Silicon carbide
89
90
BLADE COATINGS
SUBSTRATE
BOND COAT
TOP COAT HOT COMBUSTION GASES OPERATING TEMPERATURE SUBSTRATE MELTING TEMPERATURE
TEMPERATURE
91
~MELTING TEMPERATURE o o F C
2650 2100 2300 2400 2550 3035 3000 1450 1150 1260 1315 1400 1668 1650
For most metals used in blades, creep becomes significant at about one-half the melting point.
92
93