Turbine Broadband Noise Predictions Using Linearised Frequency Domain Navier-Stokes Solvers †

: A linear frequency domain Navier-Stokes solver is used to retain the influence of turning, thickness, and main geometric parameters on turbine broadband noise. The methodology has been applied to predict the broadband interaction noise produced by a representative low-speed low-pressure turbine section. The differences in the spectra with respect to those yielded by state-of-the-art flat plate based methodologies are up to 6 dB. The differences are caused by multiple effects that semi-analytical methodologies do not account for. The most important are blade thickness and turning, which have been studied separately to quantify their impact on the broadband noise footprint. The influence of changing the turbine operating conditions has been discussed as well. The outlet sound pressure level scales with the third and second power of the inlet and outlet Mach number, respectively, for constant turbulence intensity, within most of the frequency range considered.


Introduction
Since the last decade, there has been an increasing concern regarding the impact of turbine-generated broadband noise in the context of civil aviation. Nesbitt [1] described what he considered a 'turbine noise storm' caused by the confluence of several design choices. Namely, their reduced blade-count and solidity, the removal of acoustic liner treatments to reduce engine weight, and the reduction of noise emissions by other engine components. Nesbitt claimed that it is necessary to develop tools that can predict turbine noise with sufficient accuracy to understand its sources and reduce them in the future. Regarding turbine noise sources break-down, Serrano et al. [2] quantified experimentally the relative importance of turbine broadband noise, which can result in being dominant in many operating points and designs and hence, there is a need of industrial methods to predict it accurately. Turbulent wakes generated by a turbine row (see Figure 1) impinge on the downstream one, giving rise to acoustic waves with a broadband energy content. While industrial predictions of fan broadband noise have been based traditionally on the flat plate response to incoming vortical disturbances, this is much less justified for turbines. This approach can be valid for plate-like shapes such as fans, but it is unclear how to define an equivalent flat plate for a turbine, and more importantly, what the accuracy of such predictions would be. Several physical effects are not modelled by flat-plate methodologies. Among them, the influence of blade thickness and blade turning stands out, which entails a relevant flow acceleration.
The effect of blade thickness on broadband noise has been studied analytically by Gershfeld [3] and Evers and Peake [4], numerically by Gill et al. [5] and Gea-Aguilera et al. [6], and experimentally by Devenport et al. [7] within the past years, to cite some examples. The most extended conclusions are that thick airfoils lead to weaker broadband noise footprints, especially at high frequencies and in the downstream arc, as numerically shown by Gea-Aguilera et al. [6].
On the other hand, the impact of blade turning has been analysed analytically by Myers and Kerschen [8] and Evers and Peake [4], numerically by Gea-Aguilera et al. [6] and Paruchuri et al. [9], and experimentally by Devenport et al. [7], among other authors. There is a general agreement that camber has a very limited impact on the broadband noise footprint, impacting only the azimuthal modal decompositions, i.e., directivity, as shown by Myers and Kerschen [8] and Paruchuri et al. [9].
All these works, and some others not mentioned here, are either asymptotic studies or are applied to geometries with moderate thickness and low camber as those found in Fan/OGV interaction. However, for turbine geometries, thickness and camber can be very important, and the conclusions extracted from the past may not be applicable.
To shed light on the influence of the turning, thickness, and main geometric parameters on turbine broadband noise, the use of a computationally efficient linear frequency domain Navier-Stokes solver [10] is proposed. The solver runs on commodity GPUs [11], enabling the computation of the broadband noise spectra within an industrial design loop. The method has been validated previously for Fan/OGV interaction against experimental data and in a numerical benchmark in the context of the TurboNoiseBB EU project [12,13].
The objective of the present work is to assess quantitatively and qualitatively the effect of the airfoil geometry on turbine broadband noise, compare the results to the flat plate simplifications, and finally, investigate the impact of the operating point. The comparison of the present methodology to experimental data is postponed for the future since it requires other building blocks such as accurate turbulence modelling, and transmission effects through the turbine stages.

Methodology
The methodology has been thoroughly described for multi-stage applications [13] however, for completeness, it will be briefly described herein. Synthetic turbulence methods aim at reproducing a given turbulent spectrum by explicitly introducing vortical content into the simulation domain. They consist of three well-differentiated steps, namely incoming turbulence modelling, computation of the blade's acoustic response to the synthetic turbulence, and post-processing of the radiated acoustic power. The original methodology can retain certain 3D effects by using several strips at different radial positions. However, the analyses will be restricted here to a single strip for simplicity. For more information about three-dimensional effects, please refer to Blázquez-Navarro and Corral [13].

Turbulence Modelling
When turbulent wakes impact a turbine row, they give rise to broadband sound generation. These wakes can be characterised by their velocity power spectral density (PSD). Synthetic turbulence methods aim at reproducing the turbulence spectral characteristics through the summation of individual vortical gusts [14]. Their interaction with the turbine cascade is modelled under the Rapid Distortion Theory (RDT) hypothesis [15], which allows linearising their propagation through the airfoil if the fluctuations are small compared to the mean flow and the eddies remain coherent through the blade passage.
Since usually experimental data or results from scale-resolved simulations are not available, certain assumptions shall be made to build the turbulent spectrum. In the present work, the synthetic turbulent spectrum is built using the Liepmann isotropic model Liepmann [16]. Its upwash component (x 2 direction in Figure 1), which is the most relevant one for noise generated by thin bodies [17], results: is the dimensionless velocity wave-number spectrum, which is a function of the wave-vector, k, non-dimensionalised by the turbulence longitudinal integral length-scale, Λ f :k = kΛ f = (k 1 ,k 2 ,k 3 ) andk 2 =k 2 1 +k 2 2 +k 2 3 . The turbulence integral length-scale, Λ f , and the velocity root mean square value, u rms , are usually characterised by the RANS simulation of the previous row [18]. Since the purpose of this paper is not to compare numerical predictions against empirical data, reference data from the literature will be used.
The non-dimensional wave-vectork = (k 1 ,k 2 ,k 3 ) of the perturbation is not a degree of freedom that can be chosen when performing a numerical simulation. Instead, the excitation frequency, ω, and the azimuthal, β, and spanwise, γ, wave-numbers have to be explicitly imposed. The former is related to the streamwise wave-number using the hypothesis of Frozen Turbulence, i.e., the turbulent waves are assumed to be convected at the mean flow speed. The rest are computed using trigonometric relationships: where V 0 is the mean flow velocity and θ is the velocity mean flow angle, depicted in Figure 1. The spectrum is discretised into upwash vortical gusts of prescribed frequency and wavelength, w =ŵ k e i(k.x+ωt) . Its amplitude,ŵ k , is a function of the wavenumber and the bandwidth used in the discretisation, (∆k) 3 = ∆k 1 .∆k 2 .∆k 3 . We can define a dimensionless perturbation,ŵ, which is directly related to the turbulent spectrum: To keep a simpler formulation than in a previous publication [13], and as previously done by Ju et al. [19], the vortical gusts, which are sub-critical, i.e., do not produce any cut-on response, are a priori filtered out. This is performed to avoid an over-prediction of the number of cut-on modes that would occur in a purely 2D case leading to large over-predictions of the resulting spectra. Thus: After this pre-filtering process, the problem becomes fully bi-dimensional.

Response Computation
The modelling of the interaction of the turbulent spectrum with the turbine gives rise to the resolution of a large number of gusts with different ω, β which build-up the unsteady response of the OGV cascade, u(x, t): where ω(n) = n∆ω, β(µ) = µ∆β are a discrete set of frequencies and inlet wavenumbers. The summation limit µ ± n is chosen at each frequency to sample a certain fraction of the TKE, which in this work is 90% of the total. The explicit dependence of the wave-numbers on the summation indices will be omitted whenever possible to avoid unnecessary overloading of the nomenclature. The interaction of the airfoil cascade with each individual gust is solved in the frequency domain. Thus, thousands of linearised Navier-Stokes (LNS) simulations are required to compute the noise spectrum. A time instant solution can be reconstructed through the summation of all the simulations as shown in Equation (6) and displayed in Figure 2 for illustrative purposes. As mentioned above, the acoustic response is computed using a frequency-domain linearised Navier-Stokes Solver [10]. The methodology is explained with the Linearised Euler Equations (LEE) for simplicity even though the problem actually solved retains viscous effects. The extension of the ideas presented here to the solution of high Reynolds number viscous flows is straight forward. The flow is decomposed into a mean base flow, U 0 (x) and a small perturbation u(x, t) U 0 , i.e., U(x, t) = U 0 (x) + u(x, t). The base flow is obtained solving the nonlinear problem A(∂U 0 /∂x) + B(∂U 0 /∂y) = F, which in practice is done using a RANS code known as Mu 2 s 2 T [20]. The perturbations are solved with the linear counter-part of the same code [10]. If the problem is assumed harmonic in time, i.e., f =f ω e iωt and u =û ω e iωt , the linearised Euler equations can be Fourier-transformed, leading to: where A 0 (U 0 ) and B 0 (U 0 ) are the Jacobian matrices particularised at U = U 0 = (ρ 0 u 0 v 0 p 0 ) T , u ω = (ρ ωûωvωpω ) T is the vector of perturbed variables, andf ω contains just the small unsteady perturbations injected through the inlet. I is the identity matrix and L a linear operator which depends on the base flow and the frequency. The latter is actually discretised using a second-order accurate scheme to obtain the discrete operator L D , but the superindex is omitted for clarity. The harmonic solution can be formally solved asû ω =L −1f ω , though in practice the inverse of L is never computed. The forcing term associated with the synthetic turbulent spectrum,f ω , in practice also depends on β and is expressed as: where e R w is the right eigenvector associated with a vortical wave andv ω,β is the Fourier component of the perturbed velocity along the inlet, which is proportional to the upwash velocity component,ŵ ω,β . The readers are referred to Corral et al. [10], Burgos et al. [20] and Blázquez-Navarro and Corral [13], respectively, for detailed information about the solver numerics and the mesh criteria established to perform broadband noise predictions.

Post-Processing
The incoming vortical perturbation with azimuthal wavenumber β is scattered by the airfoils, giving rise to acoustic modes that conform the solutionû ω (β). In order to post-processû ω (β), the solution at the inlet and the outlet (x = x − S and x = x + S in Figure 1, respectively) is Fourier transformed in the azimuthal direction: where ϕ( ) = ∆ϕ is an azimuthal wave-number used to compute the Fourier Transform. Nonetheless, the only non-null azimuthal modes that result from the Fourier Transform are the Tyler and Sofrin Modes,β(r is the r th Tyler-Sofrin mode of the solutionû ω (β) at x ± S , obtained applying the spatial Fourier operator in the circumferential direction, i.e.,û ± ω,β = F θûω (β, x ± S ). Next, the pressure associated with the outgoing waves, p (a) ω,β (x ± S ), is obtained assuming a uniform flow-field in these regions.
where e L p ± is the left eigenvector of the 2D LEEs of either the downstream (+) or upstream (−) running pressure waves. Within the same wave-splitting process, the cut-off acoustic waves are filtered out and do not contribute to the pressure spectrum. Thus, the acoustic perturbations upstream and downstream of the turbine row,p ω,β (x S± ), are the result of applying a concatenation of different linear operators over each harmonic of the turbulent spectrum. These can be more conveniently expressed in non-dimensional form, thus obtaining:p where L −1 is the dimensionless linearised Navier-Stokes operator. It is readily seen that the acoustic response is proportional to the incoming gust strength,v ω,β . The pressure PSD can be then computed as: where S pp ω,β contains the squared summation of all the vortical waves with impinging azimuthal wave-number β(µ) = µ∆β that produce an acoustic response with wave-numbeȓ β(µ, r) = β(µ) + rV/R. This actually means that sources with different β are assumed to be uncorrelated. This result can be expressed more generally by using the non-dimensional variablep, as defined in Equation (11), leading to a non-dimensional pressure PSD: where the subindices have been omitted for clarity. Hence: where the identities ∆k 2 cos θ = ∆β and ∆k = ∆kΛ f have been used. Finally, upon integration along the azimuthal wave number and expressing S pp ω in non-dimensional form it is readily obtained: where β( ) = ∆β and ∆β = ∆k 2 cos θ is the bandwidth employed to discretise the azimuthal wave-numbers.
The cross contribution to the acoustic power of pairs of waves with either different, ω, or β is as well null since they are considered uncorrelated sources with random phases between them. As usually in aeroacoustics applications, the sound loudness is represented in logarithmic scale, by means of a Normalised Sound Pressure Level, NSPL = 10 log(S pp ω /10 −11 ), where no reference pressure is needed asS pp ω is already nondimensional, but an arbitrary factor is introduced to obtain a positive spectrum.

Results
The presented methodology has been applied to a representative section of a modern aeronautical low-pressure turbine. The geometry selected corresponds to a section that has been analysed numerically [21] and tested [22] for aeroelastic purposes. The total turning is slightly higher than 100 • , and its maximum relative thickness is about 12%. Several Operating Points (OP) have been considered however, as they all share the following parameters: where V represents the number of blades of the cascade, s/c is the pitch to chord ratio and Re is the Reynolds number. Note that the sub-indices 1 and 2 indicate respectively properties at the inlet and outlet axial stations, i.e., θ(x − S ) = θ 1 . A baseline operating point, OP1, which can be considered as the On-Design conditions, has been run. Several modifications have been performed to the boundary conditions, leading to the rest of the OPs described in Table 1.

Comparison with Flat Plate Based Methods
The first part of the analysis assesses the suitability of industry state-of-the-art flat plate methodologies to predict broadband noise of low-pressure turbine airfoils. The impact of the airfoil geometry on broadband noise predictions has been studied thoroughly in the literature with a focus on the fan outlet-guide-vanes. For these configurations, there is an extended agreement that the airfoil geometry impact on broadband noise is small, up to the highest frequencies of interest [4,6,23], since the airfoils are thin and feature a low camber. Recently, the present methodology has been applied to a modern Fan, concluding that the effect of the OGV detailed geometry on sound generation is, in general, small [12,13]. This section compares the results obtained accounting for the airfoil geometry with those obtained by replacing it with a flat plate cascade. Since which one is the best-suited equivalent flat plate is not apparent, several approaches are compared herein.
A broadband noise prediction has been performed at OP1, for reduced frequencies between f red = 2π f c/V 0 = 0.75 and 25. The maximum reduced frequency corresponds approximately to 10 4 Hz. Figure 3 compares the NSPL obtained accounting for the actual airfoil geometry and substituting it by two different flat plates. The flat plate A has been constructed using the turbine inlet mean flow properties, whereas the flat plate B is defined by the outlet conditions (see Table 2 for their definition and flow conditions).  The comparison of the flat plate A with the LPT airfoil yields the following conclusions. The pressure spectra obtained with the actual geometry in the inlet and outlet are up to 6 dB higher, at reduced frequencies lower than 12. At higher frequencies, the predictions using a flat plate cascade result in higher noise levels (about 4 dB), especially in the outlet. If the spectra are integrated along the frequency range, both effects are somehow compensated, as shown in Figure 4. The inlet prediction accounting for the airfoil geometry is overall 0.75 dB higher than the corresponding flat plate approximation while the outlet is 0.25 dB lower. Despite this similarity, the spectral distributions are very different, which can also lead to different perceived noisiness because the low-mid frequency octave bands are penalised by noise regulations.  Alternatively, the spectra retrieved by the flat plate B overestimate the Airfoil results by up to 10 dB. At low frequencies, i.e., f red < 10, there is a reasonable agreement between the airfoil's results and flat plate B. Nonetheless, at higher frequencies, the differences between both approaches become large. Consequently, flat plate A yields a better agreement with the airfoil's result than flat plate B, as shown in Figure 4. This highlights the role of the leading edge as the major contributor to sound generation, especially at f red > 10. The following analyses will therefore focus on further comparing the flat plate A results with the airfoil.
To provide a more detailed broadband noise footprint, the NSPL modal decomposition is depicted in Figure 5. The y axis corresponds to the reduced frequency, while the x-axis represents the azimuthal mode index: m = βVs/(2π). Figure 5a,c show the results obtained at the inlet with the actual and flat plate geometry. The latter exhibits a straight line of extremely low pressure, caused by the dipole behaviour of the flat plate, corresponding to a directivity −θ, because the plate's dipole cannot radiate in directions parallel to itself. The highest noise levels arise close to the cut-off condition for positive m because the regions of maximum radiation are normal to the plate. On the other hand, the simulations accounting for the actual geometry do not show such preferred directions. At low frequency, f red < 10, noise is radiated fairly homogeneously along the whole cut-on region. This is likely to be produced by the presence of additional noise sources that do not arise in flat plates and also by the reflection of acoustic waves between neighbouring blades.
Regarding the downstream noise (Figure 5b,d), the differences are even greater. The results obtained with the equivalent flat plate, (d), have the same cut-on region as in the inlet since the flow properties are uniform. The observed behaviour is analogue to that in the inlet, but with opposite directivities. Contrarily, the cut-on region exhibited by the actual geometry is broader and tilted towards positive directivities due to the effect of the swirl. At low frequencies, f red < 10, the directivity pattern is quite uniform while at higher frequencies there are very low-pressure levels close to the cut-on line, for positive m. It turns out that the cut-off transition of upstream-travelling waves has an effect on the downstream radiated noise because the pattern of the inlet cut-on line at positive m is seen also in Figure 5b. At f red 5, the pressure level clearly increases through the cut-on-cut-off transition of the inlet (m 20). This is, upstream travelling modes are reflected towards the outlet, analogously as trapped modes [24], leading to an increase of downstream noise. The information depicted in Figure 5 certainly shows that the resemblance of the integrated noise spectra is caused by statistical averaging and not by the accurate modelling of the physics. The reasons underlying the aforementioned qualitative differences will be explored next. Figure 6 shows the pressure root mean square value (RMS) at three different frequencies for both the airfoil and flat plate geometries, where a different colour scale is used in each column for enhanced visualisation. This figure highlights the different acoustic sources present in both geometries. However, a note of caution is sounded as a high pressure RMS region is not directly equivalent to a sound source because the hydrodynamic pressure does not lead to sound generation. Anyhow, the pressure RMS can provide valuable information to understand the effect of the airfoil on noise generation, compared to flat plate geometries. At low reduced frequency, f red = 3, the flat plate clearly presents a dipole pattern in the leading-edge, justifying the observations performed in Figure 5c,d. The regions with high-pressure RMS are located in the leading-edge, as predicted by quasi-steady aerodynamic theories. On the other hand, the turbine airfoil presents several regions of high levels of unsteady pressure, namely the leading-edge, the suction side peak Mach number region, and the pressure side where the flow starts to accelerate. These may correspond to volumetric noise sources, associated with the stretching of vorticity, which are known to produce sound [4]. Moreover, as the acoustic wave-length normalised with the chordλ ac = a 0 /( f · c) is about 0.2, the case can be considered acoustically compact and we may relate the pressure directly with the radiated noise, as deduced from Curle's analogy [25]. This explains the larger noise levels exhibited by the airfoil at low frequencies, compared with the flat plate case (see Figure 3). This noise mechanism is analogue to that described by Marble and Candel [26] for entropy waves, which lead to sound generation when they propagate through a non-uniform mean flow. At higher frequencies, f red = 15 and f red = 21, the behaviour observed on the airfoil is rather similar to the previous one, although the large pressure regions on the suction side shift slightly forward. On the other hand, the flat plate cascade exhibits new pressure lobes which are not present in the low frequency case. The plates do not behave any more as compact acoustic sources and the interference of multiple leading and trailing edge scatterings lead to the more complex pressure pattern observed in Figure 6e,f. It is also seen that the high pressure regions at the plate's leading-edge are more prominent than those of the airfoil. The reason for that is the finite blade thickness and leading edge radius of the airfoil, which according to Gill et al. [5] at high frequencies ( f red .t/c) > 1.5 reduces the radiated noise compared to the flat plate case. This limit corresponds in this case to approximately f red = 15. The effect of blade thickness and leading-edge radius on broadband noise have been studied by several authors previously and there is a general agreement that it diminishes the radiated noise, especially the downstream component, which is in agreement with Figure 3 where the NSPL roll-off increases at f red > 10. This point will be further explored in the next section.

Separation of Geometrical Features
The objective of this section is to highlight the origin of the discrepancies between the results produced by the flat plate approximation and the actual geometry. For that aim, the effect of blade turning and thickness will be separated as follows. The blade camber line and thickness profile are computed, so that a new uncambered airfoil with the same thickness profile as the original turbine is built. This process is sketched in Figure 7. The base flow around the cambered plate is computed using the same boundary conditions as those used with the original geometry. Alternatively, the base flow around the uncambered airfoil is computed so that the inlet Mach number is the same as in the airfoil. The selected approach to separate the thickness and turning effects relies on the assumption that the impact of the geometry on the flow is linear. Therefore the results obtained with the uncambered airfoil with a flat plate cascade will be compared first to derive a thickness correction. Next, we will compare the results obtained with the cambered plate with those of the original airfoil. We will apply the thickness correction to the cambered plate results to assess if the combined effect explains the results seen in the previous section. Figure 8 compares the spectra retrieved by the uncambered airfoil cascade to the equivalent flat plate. In agreement with other authors, it is observed that the thickness decreases the radiated noise at high frequencies. As proposed by Gershfeld [3], the thickness-based reduced frequency ( f t max /V ∞ ) is displayed, so the results may be compared with others in the literature. Figure 8a,b compare the absolute spectra obtained with the airfoil and flat plate geometries at the inlet and outlet, respectively. Clearly, the airfoil thickness has a stronger effect on the downstream radiated noise. This point can be better observed in Figure 8c, which displays the difference between the flat plate and the airfoil. In spite of some large oscillations, most likely attributed to the effect of resonances and the 2D simplification, a clear trend can be observed. A least-squares fit has been computed, which highlights that at very low frequencies, thickness appears to increase the noise, especially in the upstream direction. On the other hand, the slope of the outlet least-squares fitting is larger than that of the inlet, denoting a stronger effect of thickness on the downstream arc. These conclusions are well in agreement with the literature; Gea-Aguilera et al. [6] show remarkably similar trend lines for NACA0012 profiles for Mach Numbers between 0.3 and 0.5. The present analysis has also been performed with a cascade with a stagger angle of 44 • , similarly as the inlet flow angle at OP1, yielding remarkably similar results.
The impact of flow turning is assessed by comparing the results obtained with the original airfoil with those obtained by the cambered plate. The spectra retrieved are displayed in Figure 9. Figure 9a compares the spectra at the inlet. Three curves are displayed, namely, the spectrum produced by the original airfoil, that of the cambered plate, and the latter modified by the previously derived thickness correction law. The cambered plate features a significantly similar shape compared to the spectra obtained with the original airfoil, but the levels are up to 3 dB higher at some frequencies. The thickness correction improves the agreement with the airfoil results but is insufficient to obtain a good matching. Moreover, the cambered plate results exhibit a resonant behaviour at f red ∼ 9, which is not seen when the actual airfoil geometry is included. On the other hand, the outlet results, displayed in Figure 9b, show a better agreement between the two approaches. At the highest frequencies, the cambered plate results overpredict those of the airfoil by ∼3 dB. When the thickness correction law previously derived is applied, the matching between the curves improves. Overall, the outlet spectra are remarkably well reproduced by the cambered plate with a thickness correction.
To further explore the similarities and differences between the different spectra, their modal decomposition is analysed next. Figure 10 displays the NSPL modal decomposition for the original airfoil and the cambered plate. At first glance, their modal content exhibits a remarkable qualitative agreement. The inlet and outlet cut-on regions are very similar, as the OP1 boundary conditions were attempted to be matched. There are slight differences because a perfect agreement in the inlet cannot be achieved because the geometries are different. The same high-pressure structures appear in both geometries. At low reduced frequencies ( f red < 10) the modal content is fairly uniform, while at higher frequencies, the pressure decreases substantially and shifts closer to the cut-on lines, especially at the inlet.
At the inlet, higher pressure levels are observed close to the positive cut-on line, denoting an important contribution of positive modes to the overall footprint. This suggests strong leading-edge radiation at those frequencies, similarly as it occurred in the flat plate case (see Figure 5c). On the other hand, the inlet NSPL obtained with the cambered plate exhibits some high-pressure structures at 5 < f red < 10 and at 15 < f red < 20, which are responsible for the mismatching of the spectra observed in Figure 9a. The origin of such a high-pressure structure is yet unknown. At the outlet, the high-pressure structures occurring at f red < 10 seem peakier in the case of the cambered plate but are qualitatively similar. Moreover, for both geometries, there are very low-pressure levels at f red > 15 and m > 150 and at m ∼ −20. Unlike at the inlet, there is not a clear preferred direction of radiation at high frequency. The reason for that resides most likely in the multiple reflections suffered by the acoustic waves when they propagate from the leading edge through the inter-blade channel.
As shown by this analysis, the similarity among the spectra shown in Figure 9 is due to a similar noise generation mechanism and not to an averaging process. Turbine broadband noise generation, in comparison with flat plates, can be better understood through the separation of the two main geometrical features: thickness and camber. At the highest frequencies considered, tthe hickness reduces the radiated noise in low-speed low-pressure turbines by up to 1 and 4 dB, in the upstream and downstream direction, respectively. On the other hand, highly cambered lines, similar to those featured by turbine airfoils, increase the noise levels at low and middle frequencies ( f red 12) strongly.

Effect of the Operating Conditions
Six different Operating Points (OPs) have been run to unravel the effect of varying the Mach number on the spectra (see Table 1). Their isentropic Mach number distributions over the airfoil are shown in Figure 11. The fundamental changes between the different OPs have been the incidence and Mach number. Some changes can be adressed using the flat plate theory, but others require a detailed simulation of the flow. There are two different groups, the high Mach number cases, i.e., OPs 1-3 (see Figure 11a), and the low Mach number ones, i.e., OPs 4-6 (see Figure 11b). Among them, OP2 and OP5 increase the incidence, i, by 10 • , which leads to an increase of the inlet Mach number and the airfoil loading. On the other hand, OP3 and OP6 keep the higher incidence but reduces the pressure ratio to retrieve the same inlet Mach number as OP1 and OP4, respectively with a significant decrease in the exit Mach number.  Figures 12 and 13 show, respectively, the inlet and outlet NSPL at all the OPs considered obtained. Figure 12a,b show respectively the results of OPs 1-3 and 4-6, which have been segregated into two sub-figures for clarity. Figure 12c shows the spectra of OP1-3 and OP6 scaled as explained in the lines below; the spectra of OP4-5 are not included for readability purposes. The information depicted in these Figures is complemented by Figure 14, where the spectra modal decompositions are shown, and the cut-on regions at all the OP are seen. The inlet NSPL is weakly affected by changes in the outlet Mach number, as long as the inlet Mach number remains constant. Figure 12a shows differences between OP1 and OP3 smaller than 2.5 dB, while Figure 12b depicts an even better agreement between OP4 and OP6. Their modal decomposition exhibits significant similarities as well (See Figure 14a,e,g,k), denoting very similar physics and not just an artifact of statistical averaging.
Interestingly, the inlet NSPL scales as NSPL ∼ M 2 1 due to the increase of the cut-on region, which has the same scaling. The cut-on regions displayed in Figure 14 scale as M, but there is an additional M factor whose origin is the spanwise wavenumber integration limit of Equation (4). In other words, the total number of modes is a combination of circumferential and span-wise modes whose cut-on regions scale with M. The consistency of this scaling is demonstrated by the rather good collapse of the curves shown in Figure 12c. On the other hand, the outlet noise levels do not vary significantly with the inlet Mach number, except at low frequency ( f red < 5). OP1 and OP2 (see Figure 13a) and OP4 and OP5 (see Figure 13b) exhibit a close agreement. Figure 14b,d,h,j respectively show as well similar NSPL modal decompositions. The outlet Mach number widens or shortens the downstream cut-on region, as seen in Figure 14, and therefore controls the noise radiated in that direction through the number of modes that are cut-on or cut-off. Changes in the inlet Mach number seem to have a small impact on the outlet NSPL. However, it must be kept in mind that the physical noise is changed as it is scaled with the inlet speed, as clearly reflected by Equation (15). At low reduced frequency, it appears that loading increases by about 3 dB the outlet spectra, which may be related to a stronger stretching of the vortical waves by the flow field.
Following the same argument as for the inlet NSPL, the Outlet noise scales as NSPL ∼ M 2 2 due to the increase of cut-on modes. The OP1 and OP3 curves in Figure 13a and the OP4 and OP6 curves in Figure 13b can be fairly well collapsed when these are scaled with M 2 2 , which supposes, respectively a +3.3 and +6.9 dB shift on OP3 and OP6 outlet spectra. The collapsed curves are shown in Figure 13c. The agreement between the OP3 and OP1 is remarkable. The scaled OP6 slightly over-estimates the NSPL up to f red < 17. However, the matching is still quite good considering the high difference in mean flow properties.
The inlet Mach number is the key parameter for noise generation as it modifies the spectral shapes. The outlet Mach number merely acts as a scaling factor on the noise due to the cut-on region variation. This is justified by the fact that noise is generated mainly at the leading-edge, within most of the frequency range. Therefore, the outlet physical noise, i.e., not normalised with the inlet conditions, scales as SPL ∝ M 3 1 M 2 2 . The first factor arises from assuming constant turbulence intensity and the response scaling with the inlet velocity (see Equation (15)). The second arises from the effect of the Mach number on the cut-on condition, as previously discussed. Only at low reduced frequencies does the loading seem to produce relevant changes in the radiated noise. The conclusions obtained herein are not a by-product of statistical averaging as they are also sustained by the NSPL modal decompositions. The trends found in the non-dimensional noise spectra and modal decompositions are consistent in the OP noise analysed, which explore changes in loading at low and medium Mach number levels.

Conclusions
A methodology to include the effect of real geometry on the prediction of turbine broadband noise has been presented. The methodology has been applied to a representative low-speed low-pressure turbine airfoil and compared to state-of-the-art industrial techniques, based on the flat plate cascade theory response. Unlike in Fan/OGV noise, the airfoil geometry becomes relevant. The differences in the spectra between the actual airfoil and the equivalent flat plate analysis are up to 6 dB. At low reduced frequency ( f red < 12), the turbine airfoil produces substantially higher noise levels than the flat plate, while at higher frequencies, the flat plate becomes louder. Furthermore, the spectrum modal decomposition using both approaches has been shown as well, exhibiting big discrepancies resulting from different physics. On the other hand, when the spectra are integrated in frequency, the differences become much smaller, of the order of 1dB, due to the compensation of the aforementioned effects. The reason for such a small discrepancy is the statistical averaging of the response, which could become much larger under anisotropic turbulence, as pointed out by Devenport et al. [7].
The geometrical effects, namely blade camber, and thickness have been studied separately to understand their impact on the broadband noise footprint. For a typical lowpressure turbine airfoil, the thickness decreases the downstream radiated noise by up to 4 dB at the highest frequencies considered. The blade camber produces a large change of the mean flow properties through the turbine. As described by Evers and Peake [4], amongst other authors, the stretching of vortical waves by a non-uniform flow gives rise to sound generation. This noise generation mechanism seems to appear in highly cambered airfoils up to f red < 12, giving rise to the aforementioned 6 dB extra noise compared to flat plate methodologies.
The effect of the operating conditions of the turbine airfoil has been assessed by running six different operating points with different Mach numbers and incidences. The inlet Mach number does not significantly affect the outlet NSPL within the OP considered. Conversely, the outlet Mach number acts as a scaling factor of M 2 2 on the outlet noise through the cut-on condition. The loading increase due to the incidence increases the downstream dimensionless noise at low reduced frequencies ( f red < 5) but seems to have little effect at higher frequencies. Since the results are non-dimensionalised with the inlet properties, the outlet physical noise scales with M 3 1 M 2 2 , assuming constant turbulence intensity, within most of the frequency range considered. On the other hand, the inlet noise remains relatively unaffected by changes in the incidence that lead only to a change in the outlet Mach number if the pressure ratio is kept. These conclusions are also supported by the NSPL modal decomposition, which shows that the rather good collapse of the spectra is not just a product of the statistical averaging. This is the first time that the impact of the turbine airfoil geometry has been included in a broadband noise prediction and compared against state-of-the-art methodologies.