The Effect of Blade Thickness and Angle of Attack on Broadband Fan Noise Stewart Glegg Florida Atlantic University and William Devenport Virginia Tech Work Supported by ONR, Prog. Manager: Dr. Ron Joslin Authors would like to thank Ed Envia for providing additional slides 1
Sources of Broadband Fan Noise Rotor Alone Rotor/Stator Interaction Trailing Edge Noise Blade Tip Gap Noise Duct Wall BL Inflow Noise Inflow Noise Rotor Turbulence Rotor Tip Vortices Duct Wall BL Self Noise AIAA Papers No: 93-4402, 98-2316 2
Comparison of Predicted R/S Interaction Noise and Measurements for the SDT Fan Rig (Envia et al, 2008) Envia et al AIAA paper no 2008-2991 3
Prediction of R/S Interaction Noise with Rotor Alone Noise Added for the SDT Fan Rig (Envia et al, 2008) Envia et al AIAA paper no 2008-2991 4
Outline of Talk Rotor Alone Noise Predicting Trailing Edge Noise Sources Blade Surface Pressure Spectra and AoA Recent Developments Rotor Stator Interaction Noise Review of Existing Studies Physics of Leading Edge Noise AoA and Thickness Effects Wind Tunnel Measurements Of Leading Edge Noise Conclusions 5
Rotor Alone Noise 6
Sources of Rotor Alone Noise Tip Flows may dominate for low speed fans with large tip gaps and low AR Wall BL Inflow Noise -typically BL is very thin and source is not important Trailing Edge Noise -assumed to dominate for moderate AR blades with small tip gaps AIAA Paper No: 98-2316 7
Predicting TE Noise 8
Trailing Edge Noise Single Blade BL Turbulence Sound Waves Sharp Trailing Edge Scattering of BL pressure fluctuations causes propagating acoustic waves Kutta condition is assumed at the trailing edge (Amiet (1977), Brooks(1984), Howe (1998) 9
Trailing Edge Noise Blade Row Sound Waves BL Turbulence In a blade row there is acoustic scattering from multiple blades AIAA Jnl Vol 36(9) 1998 10
TE Noise Prediction Blade is broken down into sections Surface pressure wavenumber spectrum defined at each blade section Coupled Duct/ 3D TE Scattering Response defined for each blade section Spanwise wavenumber varies with radius Assumes small spanwise correlation lengthscale Blade Surface Pressure Blade Row Scattering Duct Mode Amplitudes AIAA Jnl Vol 36(9) 1998 11
Properties of the Surface Pressure Assume that the boundary layer on a single blade has the same surface pressure spectrum as a blade in a cascade if U and δ* are the same AIAA Jnl Vol 36(9) 1998 12
Deducing the Surface Pressure Spectrum for a given U and δ* Single Blade Blade Surface Pressure * TE Scattering = Far Field Sound Blade Row Blade Surface Pressure * Blade Row = Scattering Duct Mode Amplitudes AIAA Jnl Vol 36(9) 1998 13
Source Spectra at Different Radii Source Power Spectrum db re 20e-6Pa/sqrt(Hz) 70 60 50 40 30 Constant AoA of 7 deg blade radius 20 0 5 10 15 20 25 30 35 40 frequency khz AIAA Jnl Vol 36(9) 1998 14
Source Spectra at Different Radii Source Power Spectrum db re 20e-6Pa/sqrt(Hz) 70 60 50 40 30 blade radius 20 0 5 10 15 20 25 30 35 40 frequency khz AoA is 8 deg at hub and 4 deg at tip AIAA Jnl Vol 36(9) 1998 15
Effect of Increasing AoA on Sound Power Output Sound Power Spectrum db re 1e-12W/Hz 90 85 80 75 70 65 60 55 9 8 7 6 5 Increasing AoA 50 0 5 10 15 20 frequency khz AoA constant from hub to tip AIAA Jnl Vol 36(9) 1998 16
Effect of AoA on Total Sound Power Sound Power Output 130 125 120 115 110 105 100-2 -1 0 1 2 angle of attack re case A and B Sound power increases as ~2.4 db per degree AIAA Jnl Vol 36(9) 1998 17
Recent Developments In collaboration with Bruce Morin, Oliver Attasi and Ramon Reba AIAA paper no 2008-2993 18
Using RANS to get surface pressure spectrum Surface pressure spectrum Both models give tke Local Linearized Unsteady Flow RANS Mean Flow Dissipation Scale AIAA paper no 2008-2993 19
SDT Test Comparison Rotor Alone Noise 115 110 SDT Test Case Nominal VAN Dnstrm Data Dnstrm Pred RMG Dnstrm Pred Brooks Brooks Input 105 PWL (db) 100 95 Data RANS Input 90 10 2 10 3 10 4 F (Hz) The downstream sound power radiated from the fan blade versus frequency. The blue curve corresponds to the measured data, the green curve is the prediction based on the current model and the red curve is the prediction using the Brooks et al (1989) self noise database. AIAA paper no 2008-2993 20
Rotor Stator Interaction Noise 21
Rotor/Stator Interaction Noise Hanson and Horan (AIAA 98-2319) Glegg and Walker (AIAA 99-1888) Evers and Peake (JFM v.43, 2002) Nallasamy and Envia (JSV v. 283, 2005) Cheong et al (JASA 119, 2006) Attasi and Vinogradov (AIAA 2007-3691) u.n U s b 22
Rotor/Stator Interaction Noise Effect of AoA Evers and Peake (JFM v.43, 2002) b Gave a high frequency analysis of loaded cascade, Showed effects of AoA for tonal noise from blade wakes was large ~10dB u.n U s However only a small effect ~2dB was found on broadband noise for a homogeneous turbulent inflow 23
Rotor/Stator Interaction Noise Fully Loaded Rotor Attasi and Vinogradov (AIAA 2007-3691) Gave a complete analysis of loaded set of stator vanes in a circular duct Showed up to 3-4dB of real blade effects on broadband noise Predictions were similar to flat plate results at high frequencies Verified predictions using experimental data Concluded that real blade effects were small, but what is the physics?? 24
Leading Edge Noise 25
Phase of a Vortical Gust AIAA paper no 2006-2424 26
Streamwise Component of Goldstein s Gust AIAA paper no 2006-2424 27
Time Domain Methods Calculate the unsteady loading from a blade vortex interaction and hence the noise Arbitrary gust is modeled by multiple vortices Point vortex Works well for problems which can be reduced to 2D Glegg and Devenport JSV 2008 28
Unsteady Loading from a BVI unsteady lift 2 1 0-1 Unsteady Lift 10 deg 0 deg 5 deg -2-5 -4-3 -2-1 0 1 2 3 4 5 Ut/a 0.5 Unsteady Drag unsteady drag 0 0 deg 5 deg 10 deg -0.5-5 -4-3 -2-1 0 1 2 3 4 5 Ut/a Time Glegg and Devenport JSV 2008 29
Response to Step Upwash Gust F(t) U Step Upwash Gust Glegg and Devenport JSV 2008 30
Unsteady Loading from a Step Gust for 0,5,10 deg AoA 5 Net Unsteady Loading Unsteady Loading unsteady lift 0-5 -10-8 -6-4 -2 0 2 4 6 8 10 Ut/a Time drag 5 Glegg and Devenport JSV 2008 31
Unsteady Loading from a Step Gust for thickness to chord ratios 0.001,0.06,0.15 Net Unsteady Loading 5 Unsteady Loading unsteady lift 0-5 -10-8 -6-4 -2 0 2 4 6 8 10 Ut/a Time 5 32 ag Glegg and Devenport JSV 2008
Conformal Mapping (a) Streamlines at an AoA α deg (b) Streamlines at zero AoA (a) Flow at an angle of attack For 2D incompressible flow the airfoil can be mapped onto a circle 3 2 1 0-1 -2 2 on the distance of the 3 0-2 (b) Flow at zero a Unsteady loading depends vortex from the singular 1 point Location of the singular point depends Overlay of on (b) airfoil -1 onto thickness (a) rotated clockwise by 2α deg -3-3 -2-1 0 1 2 3 3-3 -3-2 -1 0 33 Glegg and Devenport JSV 2008
Conformal Mapping (a) Streamlines (a) Flow at at an an angle AoAof α attack(b) deg (b) Flow Streamlines at zero angle at zero of attack AoA 3 3 2 2 1 1 0 0-1 -1-2 -2-3 -3-2 -1 0 1 2 3-3 -3-2 -1 0 1 2 3 After rotation streamlines are overlaid Increased levels from vortices in upper plane cancels reduced level from vortices in lower plane 3 2 1 0-1 -2-3 -3-2 -1 0 1 2 3 Overlay of (b) onto (a) rotated clockwise by 2α deg Glegg and Devenport JSV 2008 34
Experimental Verification Experiments carried out at Virginia Tech in collaboration with William Devenport and Joshua Staubs 35
VT Stability Wind Tunnel Air exchange tower 0.5MW Drive and Fan AIAA-2008-2911 Removeable Test section 6'x6'x24' test section 15' dia., 600 h.p. fan produces flow speeds to 80m/s (Re > 5000000 per meter) Control room Anechoic system installed in control room area 36
AIAA paper no 2008-3018 From Inside Loss (db) 8 7 6 5 4 3 2 1 0 0 5000 10000 15000 20000 25000 Frequency (Hz) 37
Turbulence Grid AIAA paper no 2008-3018 5cm bars on 30cm centers 69.4% open area ratio In contraction (area 32% larger than test section) 19.5 mesh sizes ahead of airfoils 10 0 Longitudinal spectrum functions, large grid, 30m/s lf=1.28 1 0.8 0.6 0.4 0.2 Transverse corrrelation functions, large grid, 30m/s R uu (y) R ww (y) R uu (z) R vv (z) von Karman 10-1 10-2 10-3 10-4 10-5 u'/u = 3.8% v'/u = 3.9% w'/u = 4.0% L = 0.082m E uu E uu (single wire) von Karman modified von Karman 0-0.2 0 2 4 6 8 10 12 Separation (inches) 10-6 10-1 10 0 10 1 10 2 10 3 k/k e 38
Measurement locations and conditions c α o NACA 0012 0.2m -12,-8,-4,0,4,8,12 1.8m Kevlar acoustic window NACA 0015 0.6m 0,2,4,6,8,10,12 NACA 0012 0.9m 0,2,4,6,8,10,12 α 1.83 DU96 0.9m 0,2,4,6,8,10 S831 0.9m 0,2,4,6,8,10 Foam transition All measurements at 30m/s (additional at 40m/s not shown) AIAA paper no 2008-3018 39
WJD1 0.2-m NACA 0012 Effects of Angle of Attack 69 68 67 90 80 SPL 1/3 (db) 66 65 64 63 70 62 L SPL 1/3 (db) 60 50 40 30 0 o 4 o 8 o 12 o 61 60 59 300 400 500 600 20 10 2 Frequency (Hz) 10 3 2.1 Reduced Frequency ω r 21 AIAA paper no 2008-3018 40
Diapositiva 40 WJD1 Sync animation change Gershfeld to Howe William Devenport; 06/05/2008
0.2-m NACA 0012 AIAA paper no 2008-3018 Spp (db) 50 40 30 α eff 0 o 1.8 o 3.6 o 5.4 o 20 2.3 7.8 23 Reduced Frequency Moreau et al. (2005) NACA 0012 h/c=1.3 (jet) L=6.6%c u'/u=5% U=40m/s SPL (db) 55 45 35 25 α eff 0 o 12 o Reduced Frequency NACA 0012, 0.2-m h/c=9 L=40%c u'/u=3.9% U=30m/s 0 o 4 o 8 o 12 o 4 6 8 10 20 30 41
Noise Calculations Glegg, Devenport and Staubs (2008) Theory. Von Karman S pp 2 2 πω UR3 Lˆ( ω, k2) = Ω33( ω / U, k 2 2 2,0) 2 r c0 2R 3 ( ω) dk Sound from a flat plate 80 75 Panel Method k 2 70 65 Drift coordinates ω Wavenumber space SPL1/3(dB) 60 55 50 45 40 35 Amiet Analytical Solution 30 1 10 100 42 Reduced Frequency ω r
WJD2 0.2-m NACA 0012 Effects of Angle of Attack with predictions 90 80 70 Flat plate L SPL 1/3 (db) 60 50 40 30 0 o 4 o 8 o 12 o 20 10 2 Frequency (Hz) 10 3 2.1 Reduced Frequency ω r 21 AIAA paper no 2008-3018 43
Diapositiva 43 WJD2 Sync animation change Gershfeld to Howe William Devenport; 06/05/2008
0.6-m NACA 0015 Effects of Angle of Attack 90 80 0.203-m 0012, 12 o 0 o 2 o 4 o 70 6 o o 8 L SPL 1/3 (db) 60 50 0.203-m 0012, 0 o 10 o o 12 40 30 20 10 2 Frequency (Hz) 10 3 6.4 Reduced Frequency ω r 64 AIAA paper no 2008-3018 NACA 0015, 0.6-m h/c=3 L=13%c u'/u=3.9% U=30m/s 44
0.6-m NACA 0015 Effects of Angle of Attack with predictions 90 80 Flat Plate 0 o 2 o 4 o L SPL 1/3 (db) 70 60 50 10 db 6 o o 8 10 o o 12 40 30 20 10 2 Frequency (Hz) 10 3 6.4 Reduced Frequency ω r 64 AIAA paper no 2008-3018 45
Conclusions Rotor alone noise and rotor stator interaction noise are equally important Trailing edge noise can be predicted from calculations of the blade surface pressure spectrum Blade stall dramatically increases trailing edge noise Angle of attack effects on blades in isotropic turbulence are small (but not necessarily for anisotropic turbulence and asymmetric blades) Blade thickness can significantly reduce high frequency noise (but the details of the blade shape are important) 46