Theory of turbo machinery Chapter 3
D cascades Let us first understand the facts and then we may seek the causes. (Aristotle)
D cascades High hub-tip ratio (of radii) negligible radial velocities D cascades directly applicable r hub r tip 1 Low hub-tip ratio Blade speed varying Blades twisted from hub to tip
Twist
Generic airfoils NACA airfoil Double circular arcs (DCA) C4 prescribed pressure distribution
D cascades FIG. 3.1. Compressor cascade wind tunnels. (a) Conventional low-speed, continuous running cascade tunnel (adapted from Carter et al. 1950). (b) Transonic/supersonic cascade tunnel (adapted from Sieverding 1985).
D cascades FIG. 3.6. Cascade of high-speed turbine airfoils
D cascades FIG. 3.. Streamline flow through cascades (adapted from Carter et al. 1950).
D cascades How long must the infinite direction be to make derivatives negligible?
D cascades Camber line y( x) Profile thickness tx ( ) Max camber y ( a ) Position of Max camber a / l % b t y x a FIG. 3.4. Compressor cascade and blade notation.
Profile families Notice how the maximum thickness point differs C4 at 30%, NACA65 at 40% and DCA at 50% Maximum thickness close to leading edge => wider operating range but poorer high speed performance
http://airfoiltools.com/airfoil/naca4digit
Maximum t max /l for a modern high speed compressor is below 0.05 Wind turbine blade section (close to hub) requiring a very wide operating range. Low speed with extremely large t max /l
D cascades Spacing Stagger angle Camber angle Change in angle of the camber line Blade entry angle Blade exit angle Inlet flow angle Incidence s ' 1 ' 1 i FIG. 3.4. Compressor cascade and blade notation.
D cascades Fluid deviation Incidence is chosen by designer With limited number of blades: ' So that the deviation may be defined as ' FIG. 3.4. Compressor cascade and blade notation.
Aspect ratio Two aspect ratios exist: H/l, where H = blade height = r tip -r hub and l is the corda H/b, where b is the axial corda
Terminology
D cascades Cascades are used to measure losses generated through blade boundary layers flow separation shock waves And to estimate deviation As stated, the flow does not exactly follow the blade, but it is underturned ' FIG. 3.7. The flow through a blade cascade and the formation of wakes (Johnsson and Bullock, 1965).
D cascades performance parameters Y p, compressor p, turbine c, is p p c 01 01 c, is p p 0 1, Y p, turbine p p 01 01 p p 0 Pressure loss coefficient Energy loss coefficient
D cascades performance parameters The purpose is to establish functional relationships where Y p f f f 3 1 ( M ( M ( M 1 1 1,,Re),,Re),,Re) 1 1 1 (Deviation) Pressure loss coeff. Energy loss coeff.
D cascades Experimental Techniques in separate lecture Experiments should help determining Blade shape (thickness, max camber, position ) Space chord ratio Deviation.. Generalized experiments
Cascade flow characteristics Cascade measurements are D => properties such as deviation and loss coefficients have to be determined by integration: Incoming flow Slots for traversing instruments m Y s 0 tan c x p, compressor Hdy s 0 s 0 c c c c Hdy c s dy x 0 dm c c y x x y 0 s x xhdy cx 0 s m d p p 01 0 cxdy p p 0 01 1 s dy dy,is Turntable As probe is traversed downstream of cascade a periodic flow is observed Stagnation pressure change as probe is traversed through one wake
D cascades (incompressible) Continuity: c1cos1 ccos cx Forces (x and y): 1 X p p s Y sc c c or Y sc x x y1 y tan1tan FIG. 3.10. Forces and velocities in a blade cascade. Forces per unit depth!
D cascades Energy losses Loss in total pressure from skin friction Δp0 Δp0 p1 p 1 c c 1 1 X p p s y cx y cx y1 y y y c c c c c c c c 1 1 1 Y sc c c x y1 y Δp0 X Y tan1 tan X Y tan m Def of m s s s s
D cascades Energy losses Dimensionless forms are obtained normalizing with axial or absolute velocity : Stagnation pressure loss coefficients Δp c 0 x Δp c 1 0 Pressure rise coefficient and tangential force coefficient are C C p f p p X 1 cx scx Y sc x C C tan p f m
D cascades Lift and drag Corrected lift and drag coefficients versus angle of attack, A, from wind tunnel measurements for the profiles NACA 441 to 444. Adapted from Paulsen, U.S.: Aerodynamics of a full-scale, non rotating wind turbine blade under natural wind conditions, Risø National Laboratory, Roskilde Denmark 1989 Note: different definitions of angle of attack exist! C C L D L c l m D c l m
D cascades Lift and drag FIG. 3.11. Lift and drag forces exerted by a cascade blade (of unit span) upon the fluid. c c cos m x m Lift and drag forces are same as Y and X, but in the coordinates of the blades FIG. 3.1. Axial and tangential forces exerted by unit span of a blade upon the fluid. L X sin Ycos DYsin X cos m m m m
D cascades Lift and drag Rearranging previous equations: L sc tan tan sec sδp sin D sδp x 1 m 0 m 0 cos m sec 1 cos Dimension less forms are C C L D L D cl sc tan tan sec sδp sin L x 1 m 0 m m cl m D cl sδp cosm cl 0 m m CL sec C m f tan1 tan sec m C (3.6 b) D
D cascades Circulation and lift Bases in potential theory Kutta-Joukowski Theorem: L c
D cascades Circulation and lift
D cascades Velocity on suction and pressure surfaces Low pressure (suction side) is obtained through high velocity At the trailing edge the velocities must match This requires a large velocity decrease (diffusion) with potential for flow separation A diffusion factor may be defined (Lieblein): DF loc c max, s c c max, s
D cascades Incidence changes velocity, and thereby, pressure distribution FIG. 3.14. Effect of incidence on surface Mach number distribution around a compressor blade cascade
D cascades Velocity on suction and pressure surfaces What does the velocity distribution look like? (Inviscid computation) Please note: Fluid elements passing pressure and suction side do NOT pass the foil in the same amount of time
D cascades Efficiency of a compressor cascade Compressor blade cascade efficiency defined as diffuser efficiency: D p c p 1 1 c so that Δp0 0 when D 1 D 1 cx tanmtan1 tan p p -Δp c c Δp c c c c 1 0 1 0 1 1 In terms of drag and lift coefficients, this can be shown to become: D 1 C L CD sin m
D cascades Efficiency of a compressor cascade Assuming constant ratio between lift and drag C D C const. L An optimum of D 1 C L CD sin m may be found by differentiation: 4C cos 0 45deg D D m m CLsin m mopt, And the corresponding efficiency becomes D,max C C 1 D L
D cascades Efficiency of a compressor cascade Efficiency variation with average flow angle (adapted from Howell 1945).
D cascades FIG. 3.. Streamline flow through cascades (adapted from Carter et al. 1950).
D cascades FIG. 3.8. Contraction of streamlines due to boundary layer thickening (adapted from Carter et al. 1950).
Profile loss (Ainley & Mathieson correlation) Axial flow turbine Determine Y p at zero incidence (Y p,i=0 ). Y p at any other incidence is predicted as a function of the ratio i/i s where i s is the stalling incidence. Y p p p 01 01 p p 0 (3.6) The stalling incidence i s is defined as the incidence where the profile loss has doubled from the zero incidence loss. So, what about Y p,i=0?
Profile loss (Ainley & Mathieson correlation) Y p,i=0 is determined by an interpolation procedure between two extreme turbine design concepts: a) Nozzle blades: large amount of acceleration in blade row b) Impulse blades: no acceleration in relative frame From what do you know about diffusion why is the profile loss coefficient Y p much higher for the impulse blades than for the nozzle blades? The expression used for Y p,i=0 is: Y p( i0) Y p( 0) 1 1 1 tmax / l ( ) [ Yp( ] ( ) 1 ) Yp( 10) 0. (3.49) Use β for α 1 and β 3 for α if rotor blades are being considered
Reynolds number correction (Ainley & Mathieson correlation) Correct around the nominal Re-number = 10 5 according to: 1 1 tt,corrected tt,nominal Re 10 5 0.0 Correction is valid down to Re-number 10 4. Renumber is computed on chord at blade mean and exit flow conditions from the turbine.
Simplified correlation (Soderberg) Turbines If turbines are designed to operate at optimal space chord ratio (defined through the Zweifel criterion), and at zero incidence turbine blade losses can be correlated according to: tmax f (,, AR,Re) l where ε = fluid deflection, t max /l= maximum thickness to chord ratio, AR is blade aspect ratio and Re is Reynolds number.
Simplified correlation (Soderberg) For AR = 3.0 and RE=10 5 a nominal loss coefficient is predicted from: * 0.040.06( ) 100 where the deflection ε is in degrees. If the aspect ratio deviates from the nominal aspect ratio AR=3.0, the following correction is used (b=axial chord, H, blade height): * 0.01b 1 1 (1 )(0.993 ) ( nozzle H * 0.075b 1 1 (1 )(0.975 ) ( rotor H row) row)
Simplified correlation (Soderberg) If the Re-number deviates from Re=10 5, a correction must be predicted from: 5 10 Re 1/4 1 Where the Re-number is based on the hydraulic diameter according to: D h shcos scos H
D cascades Incidence: i ' 1 1 Deflection: 1 Compressor cascade characteristics (Howell 194). (By courtesy of the Controller of H.M.S.O., Crown copyright reserved).
D cascades Generalizing experimental results Deviation by Howell: Nominal deviation a function of camber and space chord ratio: ' m s l * n with the following constants for compressor cascades n 0.5 * m0.3 a l a 500
D cascades Generalizing experimental results FIG. 3.16. Variation of nominal deflection with nominal outlet angle for several space/chord ratios (adapted from Howell 1945). Example 3. recommended
D cascades Optimum space chord ratio of turbine blades (Zweifel) Note angle def! FIG. 3.8. Pressure distribution around a turbine cascade blade (after Zweifel 1945).
D cascades Optimum space chord ratio of turbine blades (Zweifel) Tangential load (force) Y m c y c 1 y Maximum tangential load (force) Y id p 01 p bh ( b is passage width ) Incompressible loss free flow p p 01 c Ratio of real to ideal load for minimum losses is around 0.8, after some manipulations: Y T sbcos tan1tan0.8 Y id For specified inlet and outlet angles sb or sl may be determined
Turbine limit load FIG. 3.3. Schlieren photograph of flow in a highly loaded transonic turbine cascade with an exit Mach number of 1.15 (from Xu, 1985).