Acessibilidade / Reportar erro

Effect of Boundary Layer and Rotor Speed on Broadband Noise from Wind Turbines

ABSTRACT:

Trailing edge surface of aerofoil is an important source of broadband aerodynamic noise production. In this paper, three aerofoil self-noise mechanisms from turbulent boundary layer near trailing edge surface are studied. Numerical computations were performed for a three bladed 2 MW horizontal axis upwind turbine of blade length 37 m and source height of 80 m, for wind speeds of 8-15 m/s. A weighted 1/3rd octave band sound power levels (SPL) are evaluated for receiver located at distance of total turbine height and at 2 m above ground. The results obtained for sound power level using baseline models showed maximum values occurring between 300 Hz and 1 kHz region of spectrum. The trends for BPM model showed a reduction of ~2 dBA near 1 kHz region of spectrum at 10 m/s, but Grosveld’s and Lowson model were identical and agreed over the entire spectrum. The effect of rotational speed on sound power levels using three baseline models are illustrated at a wind speed of 8 m/s for 2 MW turbine. Results showed that for a change of ±10% rotor speed from the rated value, there is an increase of 2 to 6 dBA over the entire sound spectrum due to differences in blade tip speed.

KEYWORDS:
Boundary layer; Wind turbine; Sound power; Rotor speed; Blade

INTRODUCTION

Sound is composed of mechanical waves that are produced when hydrodynamic pressure fluctuations in fluid vary above or below the reference atmospheric pressure values (20 µPa for air). It results in noise when this pressure amplitude exceeds the human hearing response and cause effects like sleep disturbances or annoyance (O’Neal et al. 2009O’Neal R, Hellweg Jr. RD, Lampeter RM (2009) A study of low frequency noise and infrasound from wind turbines. Maynard: Epsilon Associates Inc. ). For free field and unbounded flows, these pressure disturbances depend on velocity fluctuations that convect along free stream direction of flow. However, in diffuse sound field the aerodynamic sound production becomes significant when the turbulent boundary layer flow interacts with the fixed or moving surfaces such as an aerofoil, flat plate or corrugated geometry undergoing edge scattering (Doolan et al. 2012C Doolan, D J Moreau, L Brooks (2012) Wind turbine noise mechanisms and some concepts for its control. Acoustics Australia 40(1):7-13.; Fink 1979Fink MR (1979) Noise component method for airframe noise. Journal of aircraft 16(10):659-665. https://doi.org/10.2514/3.58586
https://doi.org/10.2514/3.58586...
).

For a wind turbine, noise production from mechanical equipment is not critical (Son et al. 2009Son E, Kim Hyunjung, Kim Hogeon, Choi W, Lee S (2009) Integrated numerical method for the prediction of wind turbine noise and long range sound propagation. Current Applied Physics 10(2):S316-S319. https://doi.org/10.1016/j.cap.2009.11.034
https://doi.org/10.1016/j.cap.2009.11.03...
). Instead, the dominant noise mechanism results from rotating blades due to boundary layer near the trailing edge surface of an aerofoil that varies with flow conditions. Another important factor responsible for noise production is rotational speed of turbine. For subsonic flows, aerodynamic noise produced from moving blades exhibits tonal and broadband character. The broadband nature is due to scattering of incident turbulence field or random velocity fluctuations in flow field, and tonal nature from blade passing frequency or its multiples. Further, it has been found from several experiment studies that most of noise occurs during the downward motion of blade and radiates in symmetrical manner with respect to the chord line (Moriarty and Migliore 2003Moriarty P, Migliore P (2003) Semi-empirical aeroacoustic noise prediction code for wind turbines. (TP-500-34478). NREL Technical Report. https://doi.org/10.2172/15006098
https://doi.org/10.2172/15006098...
; Oerlemans 2011Oerlemans S (2011) Wind turbine noise: primary noise sources. (NLR-TP-2011-066). NLR Technical Report.; Schepers et al. 2007Schepers JG, Curvers APWM, Oerlemans S, Braun K, Lutz T, Herrig A, Koegler K (2007) SIROCCO-Silent rotors by acoustic optimization. Presented at: 2nd International Meeting on Wind Turbine Noise; Lyon, France.).

In this work, empirical trailing edge noise models proposed by Brooks et al. (1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.) (BPM), Grosveld (1985) for broadband noise prediction are numerically evaluated for different flow and operating conditions for a 2 MW horizontal axis wind turbine of blade length 37 m at 0° downwind position. Inputs to simulation consist of geometric and aerodynamic data for given blade and the pitch angle of blade is set to 3.2°. In the analysis of sound pressure level, for each model, 2D incompressible quasi-uniform flow was assumed along blade span. Blade elements are considered as point source and boundary element momentum method was used to calculate aerodynamic data (Hansen et al. 2008Hansen MOL (2008) Aerodynamics of wind turbines. 2nd ed. UK: Earthscan Publishers.). For prediction of single source sound power level, octave band attenuation due to geometric divergence, also known as spherical spreading, is taken into account ignoring other factors due to atmospheric absorption, ground surface impedance and barrier or screen effects. A-weighted filter corrections were used since it is most suitable to filter the low-mid frequency range of sound spectrum (Bastasch 2011Bastasch M (2011) Summary of international wind turbine noise regulations. Technical Report. Renewable Northwest Project.).

METHODS

BPM (BROOKS, POPE, MAROCOLINI) MODEL -TBL-TEN

Brooks et al. (1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.) derived this model using data from wind tunnel experiments with NACA0012 aerofoil with varying chord lengths in anechoic chamber to measure sound pressure levels. The turbulent boundary layer noise occurs due to suction and pressure sides of aerofoil and is considered as common source of noise from wind turbine blade. Aerofoil is assumed as half-infinite flat plate. The thickness of turbulent boundary layer (δ), local Mach number (M), length of blade segment (L), as well as the distance between the observer and source (re) are important parameters to predict the acoustic field from 2D lifting surfaces (Brooks et al. 1989; Amiet and Schlinker 1981Amiet R.K, Schlinker R (1981) Helicopter rotor trailing edge noise. Presented at: 7th Aeronautics Conference; Palo Alto, USA. https://doi.org/10.2514/6.1981-2001
https://doi.org/10.2514/6.1981-2001...
). Equations 1 and 2 are used to calculate the sound pressure levels involving spectral functions, while Eq. 3 is used for angle dependent noise source, which is caused at moderate to high angle of attack (AOA) (Brooks et al. 1989; Zhu et al. 2005Zhu WJ, Heilskov N, Shen WZ, Sørensen JN (2005) Modeling of aerodynamic generated noise from wind turbines. Journal of Solar Energy Engineering 127(4):517-528. https://doi.org/10.1115/1.2035700
https://doi.org/10.1115/1.2035700...
; Moriarty and Migliore 2003Moriarty P, Migliore P (2003) Semi-empirical aeroacoustic noise prediction code for wind turbines. (TP-500-34478). NREL Technical Report. https://doi.org/10.2172/15006098
https://doi.org/10.2172/15006098...
; Dijkstra 2015Djikstra P (2015) Rotor noise and aero acoustic optimization of wind turbine aerofoils (Master Thesis). Delft: Technology University Delft.). Noise is produced due to interaction of turbulent boundary layer with trailing edge surface of an aerofoil for an incident hydrodynamic pressure field. This causes local pressure in flow higher than atmospheric pressure, resulting in sound. The sound pressure levels are therefore obtained by adding all the components logarithmically for each blade section and given by Eq. 4:

(1) SPL p = 10 . log 10 δ p * M 5 L D h r e 2 + A St p St 1 + K 1 3 + K 1

(2) SPL S = 10 . log 10 δ s * M 5 L D h r e 2 + A St s St 1 + K 1 3

(3) SPL α = 10 . log 10 δ s * M 5 L D h r e 2 + B St s St 2 + K 2

(4) SPL Total = 10 . log 10 10 SPL α 10 + 10 SPL p 10 + 10 SPL s 10

The Strouhal number (St) is used for describing oscillating flows, which involve center frequency as well as characteristic dimension of source. For flow over aerofoil it is calculated using pressure and suction side displacement thickness, given by Eq. 5:

(5) SPL Total = 10 . log 10 10 SPL α 10 + 10 SPL p 10 + 10 SPL s 10

It can be seen that St1 is related to turbulent boundary layer trailing edge noise and varies with free stream velocity (U) and Mach number (M), while St2 is related to separation noise in spectra (Grosveld 1985; Zhu et al. 2005Zhu WJ, Heilskov N, Shen WZ, Sørensen JN (2005) Modeling of aerodynamic generated noise from wind turbines. Journal of Solar Energy Engineering 127(4):517-528. https://doi.org/10.1115/1.2035700
https://doi.org/10.1115/1.2035700...
). These parameters are set between 0.01 and 10 for low Mach number flows and vary along span direction. The Reynolds number expresses the relation between inertial and viscous forces in flow, measured along the chord direction of blade, given by Eq. 6. This parameter is function of pressure side displacement thickness and also chord length of aerofoil for calculating level adjustment function, ΔK1, and amplitude functions, K1 and K2 respectively. For wind turbines the blade experiences moderate to high Reynolds number (Re) of order, Re from 3.5 × 106 to 1.2 × 107 flows and vary along the blade span. This type of source uses high and low frequency directivity functions, given by Eq. 7:

(6) Re p = δ p * U ϑ ; Re c = Uc ϑ ;

(7) D h θ , = 2 sin 2 1 2 θ sin 2 1 + M cos θ . 1 + M M c cos θ 2 ; D L θ , = sin 2 θ sin 2 1 + M cos θ 4

The boundary layer thickness (δ) and displacement thickness (δ*) are calculated as function of the chord length of aerofoil (c) and angle of attack (AOA) for both pressure and suction sides of aerofoil, as given in Brooks et al. (1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.). For all components of this source, the acoustic pressure produced near trailing edge varies as fifth power of Mach number dependence or M 5 and exhibits broadband characteristics. Further, for low Mach number flows, M ~0.2, noise radiated from pressure and suction sides of aerofoil depend upon the spectral functions A and B, which are correlated with the aerodynamic and boundary layer properties (Brooks et al. 1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.; Grosveld 1985FW Grosveld (1985) Prediction of broadband noise from horizontal axis wind turbines. Journal of Propulsion and Power 1(4):292-299. https://doi.org/10.2514/3.22796
https://doi.org/10.2514/3.22796...
; Lee and Lee 2013Lee Seunghoon, Lee Seungmin (2013) Numerical modeling of wind turbine aerodynamic noise in time domain. Journal of Acoustical Society of America 133(2):EL94-EL100. https://doi.org/10.1121/1.4774072
https://doi.org/10.1121/1.4774072...
; Moriarty and Migliore 2003Moriarty P, Migliore P (2003) Semi-empirical aeroacoustic noise prediction code for wind turbines. (TP-500-34478). NREL Technical Report. https://doi.org/10.2172/15006098
https://doi.org/10.2172/15006098...
). The spectral function A is related to turbulent boundary layer-trailing edge noise and function B with flow separation noise. For non-compact sources and attached flows, high frequency directivity function (Dh) is used and given by sin2(θ/2). The combined Doppler shift and convective amplification terms use (1 + Mcosθ) [1 + (M - Mccosθ)]2 show a cardioid pattern of sound field near trailing edge. For stalled or fully separated flow on blade, the trailing edge noise source reduces to complex dipole pattern and becomes a function of M6. Therefore, for such case, angles of attack usually exceed 12.50 for which low frequency directivity, DL is also used. The directivity angles, θ and Ø, are aligned in the azimuth and polar directions of rotor plane for turbine and shown by coordinate reference system in ref (Brooks et al. 1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.; Grosveld 1985FW Grosveld (1985) Prediction of broadband noise from horizontal axis wind turbines. Journal of Propulsion and Power 1(4):292-299. https://doi.org/10.2514/3.22796
https://doi.org/10.2514/3.22796...
; Moriarty and Migliore 2003Moriarty P, Migliore P (2003) Semi-empirical aeroacoustic noise prediction code for wind turbines. (TP-500-34478). NREL Technical Report. https://doi.org/10.2172/15006098
https://doi.org/10.2172/15006098...
).

GROSVELD MODEL - TBL-TEN

In this model, noise is produced due to turbulent boundary layer attached to trailing edge surface of an aerofoil. It is similar to BPM model and based on scaling law prediction for trailing edge noise spectrum (Grosveld 1985). Since a wind turbine blade is tapered and twisted, the individual blade segments of length (l), encounter different local free stream velocities and experience unsteady blade lift and drag forces. However, this model assumes linearly tapered rotor and ignores twist, hence does account for corrections in the angle of attack over individual blade segments as in the case of BPM. For uniform 2D lifting flows, the trailing edge noise for an isolated aerofoil can be approached with Eq. 8 (Grosveld 1985; Zidan et al. 2014Zidan E, Elnady T, Elsabbagh A (2014) Comparison of Sound Power Prediction Models of Wind Turbines. Presented at: International conference on Advances in Agricultural, Biological & Environmental Sciences; Dubai, UAE.):

(8) SPL total = 10 . log 10 B δ U 5 L D h r e 2 . K K 2 + c

where KK2 is the frequency dependent scaling function, given by Eq. 9:

(9) KK 2 f = 10 . log 10 S t St max 4 S t St max 1 . 5 + 0 . 5 4

where Stmax is the peak Strouhal number set to ~0.1; and St’ is the Strouhal number as function of boundary layer thickness ; B is the number of blades; and U is the free stream velocity scaled to U 5, Empirical constant C -5.44 dB.

Grosveld (1985) model uses the directivity function (Dh) proposed by Fink (1979Fink MR (1979) Noise component method for airframe noise. Journal of aircraft 16(10):659-665. https://doi.org/10.2514/3.58586
https://doi.org/10.2514/3.58586...
), which is similar to Eq. 7. For turbulent flows, the thickness and displacement thickness of turbulent boundary layer for aerofoils are approached using the empirical curve fitting function dependent on chord length (c) of aerofoil and angle of attack (AOA) and given by Eqs. 23 and 33 in Grosveld (1985). It must be noted that no blades for a wind turbine determine the rotor solidity. Hence this model also considers the rotor solidity factor in order to estimate the sound pressure level from a turbine. The Doppler shift and convective amplification term as discussed in BPM model is set equal to 0.8 times the Mach number to include the effects of dipole radiation of sound along the blade chord (Grosveld 1985; Schepers et al. 2007Schepers JG, Curvers APWM, Oerlemans S, Braun K, Lutz T, Herrig A, Koegler K (2007) SIROCCO-Silent rotors by acoustic optimization. Presented at: 2nd International Meeting on Wind Turbine Noise; Lyon, France.). The overall 1/3rd octave band far field sound pressure level is obtained by logarithmic integration of individual contributions of blade elements for each blade.

LOWSON MODEL - TBL-TEN

In this model, the far field sound pressure level is predicted based on empirical equations similar to BPM method; however, the directional nature of sound is ignored (Brooks et al. 1989Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.; Zidan et al. 2014Zidan E, Elnady T, Elsabbagh A (2014) Comparison of Sound Power Prediction Models of Wind Turbines. Presented at: International conference on Advances in Agricultural, Biological & Environmental Sciences; Dubai, UAE.). Sound pressure calculation involves the spectral function G6, which depends on the Strouhal number, Mach number, boundary layer thickness, distance between source and receiver, as well as empirical constant. The empirical constant changes according to flow turbulence and evaluated using mathematical gamma function. The frequency dependent scaling factors are also used to determine the sound pressure level, given by Eqs. 10 and 11:

(10) SPL total = 10 . log 10 δ M 5 L r e 2 . G 6 f + 128 . 5

(11) G 6 f = 4 f f max 2 . 5 f f max 2 . 5 + 1 2

where fmax is dependent on free stream Mach number M and boundary layer thickness δ. The boundary layer thickness is also approximated according to the Eqs. 23 and 33 given in Grosveld (1985)FW Grosveld (1985) Prediction of broadband noise from horizontal axis wind turbines. Journal of Propulsion and Power 1(4):292-299. https://doi.org/10.2514/3.22796
https://doi.org/10.2514/3.22796...
. L-span segment length (re) is the distance between the source and receiver, in m.

Figure 1a shows the illustration of the major noise generation mechanism from a wind turbine blade. It has been known that flow along blade chord is predominantly turbulent during turbine operation. The turbulent flow remains non-uniform along the blade span and varies with turbulence integral length scale, turbulence intensity level; Therefore trailing edge noise levels are found to vary with the thickness of turbulence boundary layer as well as the slope of trailing edge surfaces. Figure 1b shows streamlined flow over aerofoil surface. Near the trailing edge of aerofoil, the acoustic waves are emitted normal to chord line and caused due to interaction of turbulent boundary layer with aerofoils surface (Dijkstra 2015; Oerlemans et al. 2007Oerlemans S, Sijtsma P, López BM (2007) Location and quantification of noise sources on a wind turbine. Journal of Sound and Vibration 229(4-5):869-883. https://doi.org/10.1016/j.jsv.2006.07.032
https://doi.org/10.1016/j.jsv.2006.07.03...
; Oerlemans 2011Oerlemans S (2011) Wind turbine noise: primary noise sources. (NLR-TP-2011-066). NLR Technical Report.; Schepers et al. 2007Schepers JG, Curvers APWM, Oerlemans S, Braun K, Lutz T, Herrig A, Koegler K (2007) SIROCCO-Silent rotors by acoustic optimization. Presented at: 2nd International Meeting on Wind Turbine Noise; Lyon, France.). A feedback mechanism can also be observed between the emitted acoustic waves and the incoming flow at leading edge to cause instabilities which propagate downstream along the chord. It must be noted that for all three trailing edge noise models, the boundary layer thickness or displacement thickness is essential for predicting the overall sound pressure levels. Figure 1c shows the 37 m blade developed using NuMAD software using the blade geometric properties twist and chord length. The maximum twist and chord length of each blade for the present work is 13° and 3.22 m respectively.

Figure 1
(a) Schematic of major noise mechanisms from wind turbine blade, adapted from (Grosveld 1985FW Grosveld (1985) Prediction of broadband noise from horizontal axis wind turbines. Journal of Propulsion and Power 1(4):292-299. https://doi.org/10.2514/3.22796
https://doi.org/10.2514/3.22796...
; Doolan et al. 2012C Doolan, D J Moreau, L Brooks (2012) Wind turbine noise mechanisms and some concepts for its control. Acoustics Australia 40(1):7-13.); (b) noise generation mechanism near the trailing edge of aerofoil; (c) illustration of a 37 m turbine blade developed using NuMAD software.

RESULTS AND DISCUSSION

In the analysis of turbulent boundary layer trailing edge noise, 2D incompressible flow over aerofoil surface was assumed. The programs for three models were implemented in MATLAB 2017b software and numerical simulation performed on 12GB RAM Intel core i5 8th gen computing environment. For each wind speed, the clock time for computation took ~3/4 h. Minimum three aerofoils, NACA 0012, NACA 63215 and NACA 6320 were chosen to approximate the entire blade shape. The values for intermediate aerofoils have been linearly interpolated along the blade span. Boundary layer data for each aerofoil was obtained using XFOIL computations. It can be seen from Fig. 2a that at constant wind speed of 10 m/s the BPM model predicts peak sound power level of ~90 dBA near 800 Hz in frequency spectra; while for Grosveld shows a peak near low frequency region of spectra 200 Hz with a ~80 dBA when displacement thickness (δ*) is considered. However, pressure amplitudes for Grosveld’s model are found to agree more closely with those predicted from BPM when thickness (δ) is considered. Similar trends can be observed for Lowson’s model, in which the directivity function is ignored in prediction of sound pressure level.

Figure 2
(a) Comparison of frequency spectra of TBL-TE, A-weighted sound power level using BPM, Grosveld and Lowson models at free stream velocity, U = 10 m/s; (b) Δ dBA, between thickness δ, displacement thickness δ* at free stream velocity, U = 10 m/s, receiver height of 2 m and at distance of ~120 m, i.e. total turbine height (HH + D/2).

Figure 2b shows the difference in amplitudes across the spectra due to boundary layer thickness and displacement thickness parameters. A minimum difference of 7.5 dBA can be seen for whole spectrum and step change is observed for all three models near 1 kHz equal to 8.5 dBA. From Figs. 3a and 3b, results for computed 1/3rd octave A-weighted frequency spectra for BPM, Lowson and Grosveld baseline models at wind speeds of 8m/s and 10m/s were validated using experiment data obtained for GE 1.5sle and Siemens SWT-2.3MW turbines. It must be noted that a penalty up to 10 dB consisting of wind shear and tone was added to simulated results due to the fact that experiment results take account up to 5 dB tonal penalty and miscellaneous effects (O’Neal et al. 2009O’Neal R, Hellweg Jr. RD, Lampeter RM (2009) A study of low frequency noise and infrasound from wind turbines. Maynard: Epsilon Associates Inc. ). Good correlation between experiment data and Lowson model has been observed for 100 Hz and 1 kHz region, while data agrees closely within 5 dBA with respect to BPM and Grosveld between 1 kHz and 2 kHz region of spectrum. Further, the difference in low frequency region can be caused due to other noise mechanisms such as inflow turbulence, which are taken into account during field measurements. Figures 4a, 4b and 4c show the contour plot of frequency spectra for BPM, Grosveld and Lowson baseline TBL-TEN models discussed before. The SPL levels are plotted with varying wind speeds, between 8 m/s to 15 m/s and between 20 Hz to 3.15 kHz region of spectrum. The amplitude of power spectra is found to reach peak values of 100 dBA and 107 dBA near 500 Hz for all three models.

Figure 3
Validation of BPM, Grosveld and Lowson baseline models with OASPL experiment results of Dijkstra (2015)Djikstra P (2015) Rotor noise and aero acoustic optimization of wind turbine aerofoils (Master Thesis). Delft: Technology University Delft. and O’Neal et al. (2009)O’Neal R, Hellweg Jr. RD, Lampeter RM (2009) A study of low frequency noise and infrasound from wind turbines. Maynard: Epsilon Associates Inc. for 2.3 MW turbine and GE 1.5sle turbine (O’Neal et al. 2009O’Neal R, Hellweg Jr. RD, Lampeter RM (2009) A study of low frequency noise and infrasound from wind turbines. Maynard: Epsilon Associates Inc. ) (a) at wind speed of 8 m/s, receiver height of ~2 m located at 120 m from turbine; (b) at wind speed of 10 m/s, receiver height of ~2 m located at 120 m from turbine.

Figure 4
Illustration of TBL-TE contours, Sound Power Level (dBA) prediction from (a) BPM; (b) Lowson; and (c) Grosveld baseline models at fixed rotor speed of 15 RPM for wind speeds 8-15 m/s, receiver height of 2 m, at distance of ~120 m, i.e. total turbine height (HH+D/2).

It can be noted that the A-weighted apparent sound power level can be evaluated based on the equivalent continuous sound pressure level for a receiver located in far field and given by Eq. 12:

(12) L wA = L Aeq + 10 log 10 4 π R 2 S 0 6

The second term in Eq. 12 represents the geometric divergence of sound waves. R is the slant distance from hub height of turbine to the receiver location; S0 is the reference area, 1 m2; LAeq is the far field equivalent continuous A-weighted sound pressure level evaluated using the noise prediction methods discussed earlier.

From Fig. 4a, it is evident that the sound power levels (SPL) increase with free stream velocities (U) between 8 and 15 m/s at constant rotational speed of 15 RPM. The maximum values are obtained for Lowson model, which range between 100 and 107 dBA and agree closely with those obtained from Grosveld’s model. However, sound power level from BPM agrees closer with Grosveld’s model predictions than with Lowson’s model between 2 kHz and 3 kHz, as depicted in the Figs. 4a and 4c. For all the cases, the SPL values are considered when the observer is in downwind or 0° position. In Fig. 5a all the baseline TBL-TEN models are compared with increasing rotational speed of 14, 15 and 17 RPM, respectively. The influence of rotational speed can be seen higher near 500 Hz region. For frequencies f < 100 Hz, a difference of 10 dBA can be observed for each model. The level change at each frequency in spectrum is shown in Fig. 5b. A change of 3 dBA was found for BPM’s and Lowson’s models, while for Grosveld’s model it is ~2 dBA. A step change can also be seen near 1 kHz region for BPM’s and Lowson’s models while for Grosveld’s model it is observed near 100 Hz. For each of models, the peak is found to occur near 1 kHz, caused by turbulent boundary layer interaction from suction side of aerofoils; while for BPM’s model a secondary peak can also be seen between 1 kHz and 2 kHz, caused by turbulent boundary layer from the pressure side (Oerlemans et al. 2007Oerlemans S, Sijtsma P, López BM (2007) Location and quantification of noise sources on a wind turbine. Journal of Sound and Vibration 229(4-5):869-883. https://doi.org/10.1016/j.jsv.2006.07.032
https://doi.org/10.1016/j.jsv.2006.07.03...
; Schepers et al. 2007Schepers JG, Curvers APWM, Oerlemans S, Braun K, Lutz T, Herrig A, Koegler K (2007) SIROCCO-Silent rotors by acoustic optimization. Presented at: 2nd International Meeting on Wind Turbine Noise; Lyon, France.). Hence, model sensitivity on account of varying boundary layer parameters and rotor speed are demonstrated well. Fig 6a shows the orientation of the observer at 60° and 240° positions with respect to the rotor azimuth of machine. The maximum differences in sound power levels are computed numerically for those positions and for downward moving blade azimuth locations using one single blade element located at 75% r/R. The results in Fig. 6b show that peak difference of ~10 dBA is observed between 80° and 170° blade azimuth angle during downward motion and from outboard sections of blade. This can be attributed to increasing relative velocity field over blade span during downward motion.

Figure 5
Frequency spectra of BPM, Grosveld and Lowson baseline models at wind speed of 8 m/s, receiver height of 2 m, located at distance of 120 m (a) for 14, 15 and 17 RPM; (b) Δ dBA, between 14 and 17 RPM for baseline models.

Figure 6
(a) Illustration of plan and side views of observer azimuth position and distance of receiver location from turbine according to IEC61400-11 standards. Horizontal bar represents the rotor plane and dashed line across rotor plane shows location of 60° and 240° observer positions; (b) difference in sound power level (Δ dBA) between 60° and 240° observer positions for an isolated blade element 75% R of single blade using BPM model at fixed speed of 15 RPM, at a wind speed 8 m/s for 2 MW turbine.

Further, the influence of receiver heights on perceived noise levels (LwA), is also evaluated for 2 MW turbine. From Table 1 it can be seen that a constant difference of ~0.3 dB is obtained between 63 Hz and 8 kHz for both BPM and Grosveld methods when the receiver is shifted from 10 m to 2 m height above ground level. However, this trend is not seen for Lowson’s method and increase of ~0.1 dB in sound power level is found when the receiver location is at 2 m above ground level. It must be noted that for all simulation runs, distance between source and receiver and source heights are kept constant and in accordance with IEC-61400-11 regulations for measurement of acoustic emissions from wind turbines. Hence the receiver height difference exhibits negligible change in computed sound power levels for constant operating conditions of the turbine. However, it is not known clearly whether receiver height difference affects the measured data of sound power levels given the same operating conditions of the machine.

Table 1
Comparison of 1/3rd octave A-weighted sound power level for trailing edge noise source of a 2 MW turbine due to change in receiver height from 10 m to 2 m between 63 Hz and 8 kHz using BPM, Grosveld and Lowson method. Receivers are located at ~120 m and downwind position, for both cases at wind speed of 8 m/s and at 15 RPM.

CONCLUSIONS

Numerical computations of Grosveld, BPM and Lowson empirical trailing edge noise models were done for 2 MW wind turbine at different operating and flow conditions. All models are based on scaling laws and utilize boundary layer, geometric data of blade to predict the 1/3rd octave band sound pressure levels. The A-weighted sound power level was computed for 2 MW wind turbine using all three models and found to increase with the free stream velocity (U) and rotor RPM by ~2-10 dBA between 500 Hz and 1 kHz in sound spectrum. The results of BPM, Grosveld and Lowson baseline models are in better agreement at mid-high frequency region with existing experiment values obtained for GE1.5sle, Siemens SWT-2.3 MW turbines. A difference of ~7.5 dBA was found for all trailing edge noise models due to boundary layer thickness and displacement thickness effect. Model sensitivity results demonstrated a maximum change of 6 dBA for entire sound spectrum due to increase of rotor speed from 14 to 17 RPM. Most of the noise produced from blade occurred during the downward motion. A maximum difference of 10 dBA due to change in observer position was found between blade azimuth angles of 80° and 170°, respectively.

  • FUNDING
    There are no funders to report.

REFERENCES

  • Amiet R.K, Schlinker R (1981) Helicopter rotor trailing edge noise. Presented at: 7th Aeronautics Conference; Palo Alto, USA. https://doi.org/10.2514/6.1981-2001
    » https://doi.org/10.2514/6.1981-2001
  • Bastasch M (2011) Summary of international wind turbine noise regulations. Technical Report. Renewable Northwest Project.
  • Brooks TF, Pope DS, Marcolini MA (1989) Aerofoil self-noise prediction. (RP-1218). NASA Technical Report.
  • Djikstra P (2015) Rotor noise and aero acoustic optimization of wind turbine aerofoils (Master Thesis). Delft: Technology University Delft.
  • C Doolan, D J Moreau, L Brooks (2012) Wind turbine noise mechanisms and some concepts for its control. Acoustics Australia 40(1):7-13.
  • Fink MR (1979) Noise component method for airframe noise. Journal of aircraft 16(10):659-665. https://doi.org/10.2514/3.58586
    » https://doi.org/10.2514/3.58586
  • FW Grosveld (1985) Prediction of broadband noise from horizontal axis wind turbines. Journal of Propulsion and Power 1(4):292-299. https://doi.org/10.2514/3.22796
    » https://doi.org/10.2514/3.22796
  • Hansen MOL (2008) Aerodynamics of wind turbines. 2nd ed. UK: Earthscan Publishers.
  • Lee Seunghoon, Lee Seungmin (2013) Numerical modeling of wind turbine aerodynamic noise in time domain. Journal of Acoustical Society of America 133(2):EL94-EL100. https://doi.org/10.1121/1.4774072
    » https://doi.org/10.1121/1.4774072
  • Moriarty P, Migliore P (2003) Semi-empirical aeroacoustic noise prediction code for wind turbines. (TP-500-34478). NREL Technical Report. https://doi.org/10.2172/15006098
    » https://doi.org/10.2172/15006098
  • Oerlemans S, Sijtsma P, López BM (2007) Location and quantification of noise sources on a wind turbine. Journal of Sound and Vibration 229(4-5):869-883. https://doi.org/10.1016/j.jsv.2006.07.032
    » https://doi.org/10.1016/j.jsv.2006.07.032
  • Oerlemans S (2011) Wind turbine noise: primary noise sources. (NLR-TP-2011-066). NLR Technical Report.
  • O’Neal R, Hellweg Jr. RD, Lampeter RM (2009) A study of low frequency noise and infrasound from wind turbines. Maynard: Epsilon Associates Inc.
  • Schepers JG, Curvers APWM, Oerlemans S, Braun K, Lutz T, Herrig A, Koegler K (2007) SIROCCO-Silent rotors by acoustic optimization. Presented at: 2nd International Meeting on Wind Turbine Noise; Lyon, France.
  • Son E, Kim Hyunjung, Kim Hogeon, Choi W, Lee S (2009) Integrated numerical method for the prediction of wind turbine noise and long range sound propagation. Current Applied Physics 10(2):S316-S319. https://doi.org/10.1016/j.cap.2009.11.034
    » https://doi.org/10.1016/j.cap.2009.11.034
  • Zhu WJ, Heilskov N, Shen WZ, Sørensen JN (2005) Modeling of aerodynamic generated noise from wind turbines. Journal of Solar Energy Engineering 127(4):517-528. https://doi.org/10.1115/1.2035700
    » https://doi.org/10.1115/1.2035700
  • Zidan E, Elnady T, Elsabbagh A (2014) Comparison of Sound Power Prediction Models of Wind Turbines. Presented at: International conference on Advances in Agricultural, Biological & Environmental Sciences; Dubai, UAE.

Edited by

Section Editor: Paulo Greco

Publication Dates

  • Publication in this collection
    03 June 2019
  • Date of issue
    2019

History

  • Received
    13 May 2018
  • Accepted
    21 Sept 2018
Departamento de Ciência e Tecnologia Aeroespacial Instituto de Aeronáutica e Espaço. Praça Marechal do Ar Eduardo Gomes, 50. Vila das Acácias, CEP: 12 228-901, tel (55) 12 99162 5609 - São José dos Campos - SP - Brazil
E-mail: submission.jatm@gmail.com