Computation of Airfoil Trailing Edge Noise Based ...

3 downloads 0 Views 3MB Size Report
Fan Noise Simulation. Source modelling (LES) and experiments (Christophe, Moreau,. Roger, Wang) for simplified fan geometry: Valeo Controlled Diffusion.
Computation of Airfoil Trailing Edge Noise Based on Steady Flow Data Sebastian Remmler von Karman Institute for Fluid Dynamics, Rhode-St-Gen` ese (Belgium) Dresden University of Technology (Germany)

EPFDC Presentation, 15th July 2009

Outline

Introduction

Theory

Validation

Application

Outline Computation of Airfoil Trailing Edge Noise Based on Steady Flow Data 1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Conclusion

Outline

Introduction

Theory

Validation

Outline

1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

Fan noise generation mechanisms Tonal noise:

Broadband noise:

thickness noise (for M  0)

turbulence-interaction noise (leading edge)

load noise

boundary layer noise (trailing edge)

rotor-stator interaction noise

Conclusion

Outline

Introduction

Theory

Validation

Application

Fan noise generation mechanisms Tonal noise:

Broadband noise:

thickness noise (for M  0)

turbulence-interaction noise (leading edge)

load noise

boundary layer noise (trailing edge)

rotor-stator interaction noise

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Fan Noise Simulation Source modelling (LES) and experiments (Christophe, Moreau, Roger, Wang) for simplified fan geometry: Valeo Controlled Diffusion airfoil in uniform flow Analogies: extended Amiet, FW-Hall, FW-Hawkings

Objective: complex geometries → simpler methods necessary

Outline

Introduction

Theory

Validation

Application

Acoustic Source Modelling

DNS

LES

U-RANS

RANS stochastic or statistical model

pressure fluctuations = acoustic sources

broadband and tonal noise

tonal noise

broadband modelled

Conclusion

Outline

Introduction

Theory

Validation

Outline

1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

Statistical Model: General Idea Assumption: streamwise and spanwise homogeneity ∂ hu1 i (x2 )/∂x2 , u02 (x2 ), Λ(x2 ) (from RANS computation) R22 (˜ r1 , r˜3 ), Uc , α (modelled) ⇓ Statistical model explicit expression computed by Monte Carlo integration ⇓ Φpp (ω): surface pressure spectrum equivalent acoustic source

Conclusion

Outline

Introduction

Theory

Validation

Application

Statistical Model of Panton & Linebarger [5] Wavenumber spectrum of pressure: ∞ 2 RRR k1 −k(x2 +ˆ x2 ) x2 ) S ∂hu1 i(x2 ) ∂hu1 i(ˆ φpp (k1 ) = 8ρ2 e dx2 dˆ x2 dk3 22 ∂x2 ∂x ˆ2 k2 0

Fourier transform of vertical velocity covariance: RR∞ u0 (x )u0 (ˆ x ) R22 cos(αk˜1 r˜1 ) cos(k˜3 r˜3 )d˜ r1 d˜ r3 S22 = 2 2π2 2 2 Λ2 0

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Statistical Model of Panton & Linebarger [5] Wavenumber spectrum of pressure: ∞ 2 RRR k1 −k(x2 +ˆ x2 ) x2 ) S ∂hu1 i(x2 ) ∂hu1 i(ˆ φpp (k1 ) = 8ρ2 e dx2 dˆ x2 dk3 22 ∂x2 ∂x ˆ2 k2 0

Fourier transform of vertical velocity covariance: RR∞ u0 (x )u0 (ˆ x ) R22 cos(αk˜1 r˜1 ) cos(k˜3 r˜3 )d˜ r1 d˜ r3 S22 = 2 2π2 2 2 Λ2 0

Modelled terms Velocity fluctuation correlation coefficient:  √ 2 2 ˜ 2 R22 = 1 − √ 2 r˜ ˜ 2 e− r˜ +(˜x2 −xˆ2 ) 2

r˜ +(˜ x2 −x ˆ2 )

Ratio of streamwise to spanwise turbulence scales: 1 . α . 3 Convecion velocity (Φpp (ω) = φpp (ω/Uc )/Uc ): 0.5 . Uc /U∞ . 1

Outline

Introduction

Theory

Validation

Application

Amiet’s theory Analytical formula for the sound generated by the interaction of boundary layer turbulence with the trailing edge of an airfoil, developed by Amiet [1] and extended by Roger & Moreau [6] Sound radiated sound in the midspan plane (x = [R, θ, 0]):   L sin θ 2 (k c)2 |I|2 Φpp (ω) ly (ω) Spp (x, ω) = 2π R 2 Φpp (ω) power spectral density of surface pressure fluctuations near trailing edge ly (ω) spanwise correlation length near trailing edge I transfer function for airfoil geometry c chord length L span

Conclusion

Outline

Introduction

Theory

Validation

Outline

1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Test Case I: BL with Adverse Pressure Gradient Boundary layer profiles available from Bradshaw [2], no RANS computation required. Experimental data directly used as input for statistical model: 1

0.12

〈u1〉 / U∞

0.08 0.6 0.06 0.4 0.04 〈u1〉 / U∞ u'2 / U∞ lm / δ

0.2 0 0

0.2

0.4

0.02 0 0.6

x2 / δ

0.8

1

u'2 / U∞ , lm / δ

0.1

0.8

Outline

Introduction

Theory

Validation

Application

Test Case I: BL with Adverse Pressure Gradient Scale anisotropy factor α: -40 -45

10 log[Φ(ω) / U3∞ ρ2 δ*]

-50 -55 -60 -65 -70 experiment (Bradshaw) spectral model (α = 1) spectral model (α = 2) spectral model (α = 3) spectral model (α = α(k1δ))

-75 -80 -85 0.1

1 ω δ* / U∞

10

α = 3 for low wavenumbers, α = 1 for high wavenumbers, linear blending inbetween

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Test Case I: BL with Adverse Pressure Gradient Convection velocity Uc :

-40

-40

-45

-45

-50

-50 10 log[Φ(ω) / U3∞ ρ2 δ*]

10 log[Φ(ω) / U3∞ ρ2 δ*]

Scale anisotropy factor α:

-55 -60 -65 -70 experiment (Bradshaw) spectral model (α = 1) spectral model (α = 2) spectral model (α = 3) spectral model (α = α(k1δ))

-75 -80 -85 0.1

1 ω δ* / U∞

-55 -60 -65 -70 -75 experiment (Bradshaw) spectral model (a) Uc / U∞ = 0.9 spectral model (b) Uc / U∞ = f(k1δ) spectral model (c) Uc / Uτ = f(k1δ)

-80 -85 10

α = 3 for low wavenumbers, α = 1 for high wavenumbers, linear blending inbetween

0.1

1 ω δ* / U∞

Uc /U∞ = 0.9 for low wn, Uc /U∞ = 0.55 for high wn, linear blending inbetween

10

Outline

Introduction

Theory

Validation

Application

Conclusion

Test Case II: BL with Zero Pressure Gradient mean velocity

velocity fluctuations

25

mixing length

5 experiment (Schewe) DNS (Spalart) RANS (K-ω) RANS (v2f) RANS (RSTM)

20

0.14

4.5 0.12 4 0.1

3.5

10

3

0.08 lm / δ

u+

100 u'2 / U∞

15

2.5

0.06

2 1.5

5

0.04

1

experiment (Schewe) RANS (K-ω) RANS (v2f) RANS (RSTM)

0.5 0

0 1

10

100 +

y

1000

lm (Prandtl) lm = 1.9 L (K-ω) lm = 0.65 L (v2f) lm = 1.9 L (RSTM)

0.02

0 0

0.2

0.4 0.6 x2 / δ

0.8

1

0

0.2

0.4 0.6 x2 / δ

0.8

⇒ rather similar RANS results for all used turbulence models (K-ω SST, v 2 f , RSTM)

1

Outline

Introduction

Theory

Validation

Application

Conclusion

Test Case II: BL with Zero Pressure Gradient -45

-55

2 3

10 lg[Φ(ω) / ρ U∞δ*]

-50

-60

-65 experiment (Schewe) statistical model; RANS (K-ω) statistical model; RANS (v2f) statistical model; RANS (RSTM)

-70

-75 0.1

1 ω δ* / U∞

10

⇒ spectral results better with v 2 f and RSTM, compared to K-ω

Outline

Introduction

Theory

Validation

Outline

1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

CD Airfoil in Jet Wind Tunnel

Experimental data: hot-wire and wall pressure spectra (ECL and MSU) [4] LES computations: Christophe [3] and Wang [7] BL data extracted on suction side near trailing edge Trailing edge noise: measured and computed by Amiet’s theory

Outline

Introduction

Theory

Validation

Application

Conclusion

CD Airfoil: RANS Results mean velocity

velocity fluctuations

18

mixing length

3 experiment RANS (v2f) RANS (K-ω) RANS (RSTM)

16 2.5

RANS (v2f) RANS (K-ω) Prandtl RANS (RSTM)

1

14 0.8

10 8 6

lm [mm]

2 u'2 [m/s]

〈u1〉 [m/s]

12

1.5

0.6

0.4

1

4 2

experiment RANS (v2f) RANS (K-ω) RANS (RSTM)

0

0.2

0.5

0 0 1 2 3 4 5 6 7 8 9 10 x2 [mm]

0 0

1

2

3 4 5 x2 [mm]

6

7

8

0

1

2

3 4 x2 [mm]

5

⇒ similar results with K-ω and v 2 f , strong scattering of mixing length results, RSTM fails predicting the mean flow field

6

7

Outline

Introduction

Theory

Validation

Application

CD Airfoil: Wall Pressure Spectra 75 70 65

Φ [dB re 2⋅10-5 Pa]

60 55 50 45 40 35 30 25 102

experiment statistical model; experimental flowfield statistical model; RANS (K-ω) statistical model; RANS (v2f) 103 f [Hz]

⇒ good spectral result with v 2 f , acceptable with K-ω

104

Conclusion

Outline

Introduction

Theory

Validation

Application

CD Airfoil: Wall Pressure Spectra 75 70 65

Φ [dB re 2⋅10-5 Pa]

60 55 50 45 40 35 30 25 102

experiment LES (Christophe) LES (Wang) statistical model; experimental flowfield statistical model; RANS (K-ω) statistical model; RANS (v2f) 103 f [Hz]

104

⇒ good spectral result with v 2 f , acceptable with K-ω ⇒ differences due to turbulence model comparable to LES

Conclusion

Outline

Introduction

Theory

Validation

Application

CD Airfoil: Far-field Noise 40 30

SPL [dB re 2⋅10-5 Pa]

20 10 0 -10 -20 -30 -40 102

experiment statistical model; experimental flowfield statistical model; RANS (K-ω) statistical model; RANS (v2f) 103 f [Hz]

104

⇒ differences in wall pressure spectra are translated directly into similar differences in calculated far-field noise spectra

Conclusion

Outline

Introduction

Theory

Validation

Application

CD Airfoil: Far-field Noise 40 30

SPL [dB re 2⋅10-5 Pa]

20 10 0 -10 -20 -30 -40 102

experiment LES (Christophe) LES (Wang) statistical model; experimental flowfield statistical model; RANS (K-ω) statistical model; RANS (v2f) 103 f [Hz]

104

⇒ differences in wall pressure spectra are translated directly into similar differences in calculated far-field noise spectra

Conclusion

Outline

Introduction

Theory

Validation

Outline

1

Introduction: Aeroacoustics of fans and airfoils

2

Theoretical Background

3

Validation

4

Application to CD Airfoil

5

Conclusion and Perspective

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Conclusion Spectral model was implemented and validated with literature test cases (flat plate boundary layer with and without adverse pressure gradient) Calculated spectrum strongly depends on RANS turbulence model, optimal choice of turbulence model is case-dependent

Outline

Introduction

Theory

Validation

Application

Conclusion

Conclusion Spectral model was implemented and validated with literature test cases (flat plate boundary layer with and without adverse pressure gradient) Calculated spectrum strongly depends on RANS turbulence model, optimal choice of turbulence model is case-dependent Application to CD airfoil: Good agreement of wall pressure spectra with v 2 f turbulence model, K-ω also acceptable, RSTM did not work Noise computation with Amiet’s theory: agreement as good as wall pressure spectra

Outline

Introduction

Theory

Validation

Application

Conclusion

Conclusion Spectral model was implemented and validated with literature test cases (flat plate boundary layer with and without adverse pressure gradient) Calculated spectrum strongly depends on RANS turbulence model, optimal choice of turbulence model is case-dependent Application to CD airfoil: Good agreement of wall pressure spectra with v 2 f turbulence model, K-ω also acceptable, RSTM did not work Noise computation with Amiet’s theory: agreement as good as wall pressure spectra

+ Quality of calculated spectra is not inferior to LES results, computational costs are lower by orders of magnitude - preliminary knowledge about flow field or validation of RANS required, restricted applicability due to assumptions

Outline

Introduction

Theory

Validation

Application

Conclusion

Perspectives Extension from one-point method to surface distributed pressure fluctuations ⇒ possible use of other analogies (e. g. Curle, Ffowcs Williams & Hawkings...) ⇒ more complex geometries like 3D fan blades ⇒ acoustical post-processing with standard tools, such as SYSNOISE

Outline

Introduction

Theory

Validation

Application

Conclusion

Perspectives Extension from one-point method to surface distributed pressure fluctuations ⇒ possible use of other analogies (e. g. Curle, Ffowcs Williams & Hawkings...) ⇒ more complex geometries like 3D fan blades ⇒ acoustical post-processing with standard tools, such as SYSNOISE

Trend verification ⇒ even if absolute noise levels are not always precisely predicted, at least trends must be reproduced by the statistical model ⇒ parametrical studies for reference velocity, angle of attack, incoming turbulence ⇒ useful for design of aeroacoustic applications

Outline

Introduction

Theory

Validation

Application

Conclusion

Perspectives Extension from one-point method to surface distributed pressure fluctuations ⇒ possible use of other analogies (e. g. Curle, Ffowcs Williams & Hawkings...) ⇒ more complex geometries like 3D fan blades ⇒ acoustical post-processing with standard tools, such as SYSNOISE

Trend verification ⇒ even if absolute noise levels are not always precisely predicted, at least trends must be reproduced by the statistical model ⇒ parametrical studies for reference velocity, angle of attack, incoming turbulence ⇒ useful for design of aeroacoustic applications

Include wall functions to make the statistical model work also with coarse RANS simulations

Outline

Introduction

Theory

Validation

Thank you for your attention. Questions?

Application

Conclusion

Outline

Introduction

Theory

Validation

Application

References R. K. Amiet. Acoustic radiation from an airfoil in a turbulent stream. J. Sound Vib., 4(4):407–420, 1975. P. Bradshaw. The turbulence structure of equilibrium boundary layers. J. Fluid Mech., 29(4):625–645, 1967. J. Christophe, J. Anthoine, and S. Moreau. Trailing edge noise of a controlled-diffusion airfoil at moderate and high angle of attack. In 15th AIAA/CEAS Aeroacoustics Conference, 2009. S. Moreau, D. Neal, Y. Khalighi, M. Wang, and G. Iaccarino. Validation of unstructured-mesh LES of the trailing-edge flow and noise of a controlled-diffusion airfoil. In Proceedings of the Summer Program 2006, pages 1–14. Center for Turbulence Research, 2006. R. L. Panton and J. H. Linebarger. Wall pressure spectra calculations for equilibrium boundary layers. J. Fluid Mech., 65(2):261–287, 1974. M. Roger and S. Moreau. Back-scattering correction and further extensions of Amiet’s trailing-edge noise model. Part 1: Theory. J. Sound Vib., 286:477–506, 2005. M. Wang, S. Moreau, G. Iaccarino, and M. Roger. LES prediction of pressure fluctuations on a low speed airfoil. In Annual Research Briefs. Center for Turbulence Research, 2004.

Conclusion

Outline

Introduction

Theory

Validation

Application

Appendix: Turbulence data from RANS solution Turbulence length scale Λ general:

Λ = 1.5 lm = 1.5 Cm L

K-ε, K-ω, RSTM:

L = Cµ K ε

v2f :

L = CL max

3/2

= 

Cµ K 1/2 β ω

K 3/2 ε ;





ν3 ε

1/4 

Wall-normal velocity fluctuations u2 √ K-ω: u2 = 2 α2 K α2 - empirical anisotropy factor √ v 2 - wall-normal Reynolds-stress v2f : u2 = v 2 RSTM:

u2 ← hui uj i

hui uj i - Reynolds-stress tensor

Conclusion

Outline

Introduction

Theory

Validation

Application

Conclusion

Appendix: Monte Carlo Method Monte Carlo integral for one dimension (N random input points xi ): b

Z

f (x)dx = (b − a) hf i =

I=

i=1

1

1

0.8

0.8

0.6

0.6

f(x, y)

f(x)

a

N →∞ (b − a) X f (xi ) N

0.4 f(x) theoretical f(x) average f(xi)

0.2 0 0

0.2

0.4

0.6 x

0.4 f(x, ymax) theoretical f(x, ymin) theoretical f(x, y) average f(xi, yi)

0.2 0 0.8

1

0

0.2

0.4

0.6 x

convergence: ∆I ∝ √σN ⇒ well suited for multidimensional integration (coordinate transform necessary to reduce σ).

0.8

1