The increasing importance of aerodynamic noise is a result of recent emergence of, amongst others, wind turbines and multirotor unmanned (and potentially manned) vehicles. Prospects for incorporation of multirotor transport in densely populated urban areas [1] require the effective solutions to mitigate the aerodynamic noise.
Whereas the propeller noise is typically addressed by analysis of thickness and loading constituents, the influence of airfoil shape and the related airfoil self-noise cannot be neglected.
Airfoil self-noise is typically divided into five types [2]:
Turbulent boundary layer (TBL) trailing edge noise, where the noise is generated as turbulence passes over the trailing edge. This typically results in a broadband noise spectrum. Laminar boundary layer (LBL) vortex shedding (VS) instability noise, which occurs for low Reynolds number and low turbulence flows. Hydrodynamic instabilities generated on the laminar bubble develop vortices responsible for high levels of tonal noise, which is schematically depicted in Figure 1. Boundary layer separation (BLS), which generates noise due to shed turbulent vorticity. Increasing angle of attack shifts position of separation towards the leading edge, resulting in deep stall, which in turn generates low-frequency noise. Blunt trailing edge (BTE) noise is caused by VS behind the thick trailing edge. Tip vortex noise results from highly turbulent flow occurring near the tips of lifting blades or wings.
One of the first comprehensive studies on airfoil tonal noise was an experimental research that was conducted in the 1970s [3]. The ‘ladder-type’ variation of the peak tone frequency
In the light of these findings, LBL-VS is potentially important in the context of low Reynolds number propellers, which operate in low turbulence conditions and at moderate angles of attack [9]. Additionally, an appropriate computational model must be selected, which allows us to identify the existence of tonal noise on arbitrary airfoil, with the associated cost not ruling out the shape optimisation.
The computational analysis of airfoil noise, especially if broadband noise is considered, requires a high-fidelity computational fluid dynamics (CFD) model, capable of resolving scales of turbulence, such as the large eddy simulation (LES) or hybrid RANS-LES [10]. The tonal nature of LBL-VS noise requires appropriate modelling of laminar-TBL transition and formation of hydrodynamic instabilities, which propagate over the airfoil surface up to the trailing edge. Especially, late transition occurring at low turbulence conditions is important. The complexity of LBL-VS, further increased by the aeroacoustic feedback from the airfoil trailing edge, responsible for amplification of T–S waves, requires direct computational aeroacoustics (CAA) methods performed with direct numerical simulations (DNS) [11,12]. Whereas a direct application of CAA would resolve the acoustic field, the use of the relatively inexpensive unsteady Reynolds-averaged Navier–Stokes (RANS) model [13], in conjunction with the use of transition SST for turbulence closure, would suffice for modelling, with good accuracy, the existence of tonal noise as well as the levels of this noise. For the selected computational case of symmetric NACA airfoil, the method was shown whereby promising results, in agreement with the experiment described in the study of De Gennaro et al. [13], could be obtained. However, sufficient investigations have not been conducted regarding the performance associated with, as well as the limitations occurring in, a generic airfoil problem.
In the present work, the validity of an incompressible unsteady Reynolds-averaged Navier–Stokes (URANS)-based modelling approach for airfoil LBL-VS tonal noise was assessed. The generic, asymmetric S834 airfoil was investigated, and the results obtained were compared to the experiment [14]. As the exact reproduction of an experimental setup, i.e. modelling the finite-span airfoil, was not feasible, the two-dimensional analysis of an airfoil, i.e. assuming an airfoil of infinite span, was conducted. The focus was thus on the qualitative aspects of tonal noise rather than the quantitative comparison of noise levels. Specifically, the existence and frequency of tones were of primary interest.
The present article is structured as follows: in Section 2, the details of aerodynamic and acoustic models are described. In Section 3, aerodynamic and acoustic results are presented and discussed. In Section 4, the summary and conclusions are stated.
Navier–Stokes (N–S) equations are used to describe the viscous fluid motion, which in the present instance is represented by the following continuity equation:
So far as practical engineering applications are concerned, owing to the magnitude of the computational cost involved, derivation of direct numerical solutions to the N–S equations is not performed. The most popular simplification of N–S equations in industrial practice involves averaging; and among approaches involving averaging, the most popular, Reynolds averaging, gives RANS equations. RANS equations are of the form:
The selection of eddy-viscosity model is case-specific. In this work, the requirement to accurately model the laminar bubble and the associated transition to the TBL led to the selection of a transition SST model [15], a common choice in the literature [16]. Transition SST is a four-equations, correlation-based turbulence model, which solves transport equations for turbulence kinetic energy
In the present work, airfoil analysis at low Mach and moderate Reynolds number (
Case summary.
22.4, 32, 47.9 [m/s] | |
1.225 [kg/m3] | |
5×105 [–] | |
4.4 [°] | |
c | 0.2286 [m] |
TILE | 0.18%–0.99% |
The analyses were performed with a 2D model, which exploits the two-dimensional nature of LBL-VS and significantly reduces the computational cost. The domain extends over 20 chords upstream and almost 40 chords downstream from the airfoil. The problem was defined using common boundary conditions for incompressible flow, namely velocity-inlet with required freestream velocity and pressureoutlet in the downstream of domain (Fig. 2).
Due to the turbulence dissipation in the computational domain, its value is computed from the results of CFD model at the location directly above the leading edge of an airfoil. The turbulence at the inlet of the computational domain is controlled using the turbulence intensity and turbulent viscosity ratio.
Turbulence intensity is defined as:
The selected turbulence model (transition SST requires a mesh that allows resolution of the viscous sublayer, and specifically in such a way that y+ does not exceed unity. The mesh is particularly important in the region of laminar-turbulent transition (laminar bubble), where the hydrodynamic instabilities are formed. Moreover, the region downstream of the laminar bubble should be well-modelled, such that the VS directly influencing pressure oscillations on the airfoil surface is properly resolved. A structured C-mesh was used to achieve orthogonal mesh in the region surrounding the airfoil (Fig. 3). The mesh consists of 305,000 quadrilateral elements. The near wall mesh is presented in Figure 4, whereas the mesh near the trailing edge is presented in Figure 5. The airfoil trailing edge is blunt, with three elements across its thickness.
In this work, a commercial software Ansys Fluent 21 [18] was used compute the flowfield around the airfoil. The incompressible pressure-based solver was selected due to low freestream velocity (
The hybrid approach to modelling aeroacoustics is based on the separation of flow dynamics and acoustics. Firstly, the aerodynamic model is calculated. Then the noise sources are computed and propagated to the observer location using a selected method. The possible approaches include Curle analogy [19], Ffowcs Willimas–Hawkings analogy [20], Lighthill analogy [21,22], and others. The latter approach, schematically presented in Figure 6, is used in this work.
Lighthills analogy is formulated by rearranging N–S equations with acoustic density perturbation
In the present work, Lighthill analogy, with volume sources in the region surrounding the airfoil, was used. The analysis was performed using a commercial software,
To recreate the experimental environment, turbulence intensity not exceeding 1% had to be ensured at the airfoil [14]. However, the wind tunnel test section turbulence properties were not described in sufficient detail to allow more exact reproduction of the freestream conditions. For this reason, the influence of freestream turbulence, specifically turbulence intensity, on the computational model was investigated. In Figure 8, the example of turbulence dissipation is presented with turbulence intensity levels around the airfoil.
Several analyses with turbulence intensity at the leading edge of an airfoil ranging from 0.18% to 0.99% were conducted. Figure 9 presents the instantaneous pressure coefficient on the rear part of the airfoil for turbulence intensities of 0.18% and 0.55%. The hydrodynamic instabilities are formed at the position of laminar bubble, approximately at 60% chord and 87% chord for the suction and pressure sides, respectively. Hence, on the suction side the pressure oscillations reach their maximum amplitude and begin to decay. On the pressure side, the oscillations are only observed for the lower turbulence case.
A comparison of pressure contours for cases with different turbulence intensities is presented in Figure 10. Lower turbulence results in pressure oscillations at both pressure and suction sides, whereas for higher turbulence only suction side oscillations are present. It can be noted that the strength of vortices generated at suction side is higher for the lower turbulence (as assessed by density of the isopressure lines). Additionally, the vortices in the wake of an airfoil are only observed in the lower turbulence case.
The boundary layer in the vicinity of hydrodynamic instabilities formation is presented in Figure 11 for the low turbulence case (TI = 0.18%). The three snapshots were taken in equal time steps of 0.4
Similarly, in Figure 12, the vortices pattern is presented for the low turbulence case (TI = 0.18%), where the vortices are generated at both pressure and suction sides. The location and movement of these are visualised in the rear part of the airfoil using q-criterion. It can be observed that the vortex at the suction side, initially downstream of the pressure side vortex, is moving faster, and eventually passes the latter. Consequently, different frequencies of oscillations at both sides are expected to be observed during quantitative analysis.
In Figure 13, drag monitor frequency domain representation is shown. Two peaks are observed: the major at
A summary of pressure oscillations’ presence and strength is shown in Figure 14. The variation of peak amplitudes in spectral analysis of pressure signals, collected from four different points on the airfoil surface, is plotted against turbulence intensity. The general trend of decreasing amplitude with increasing turbulence intensity is observed. Vortices are only present on the pressure side for the lowest turbulence intensity case, while on the suction side they are observed up to the intensity of 0.65%. It should be noted that the amplitude of pressure oscillations on the suction side halves between the 87% chord and 96% chord locations, while on the pressure side the amplitudes are approximately doubled. In Table 2, a summary of peak tone frequencies and amplitudes for pressure signals collected at 87% chord and drag counts (whole airfoil) is presented. An important result is shown for the case with a turbulence intensity of 0.18%, where the lower pressure amplitude at 1,300 Hz results in the peak frequency for drag signal. This result agrees with the generally observed dominant role of pressure side instabilities in the generation of tonal noise. For higher turbulence levels, where the pressure side oscillations are very small, their frequency is the same as that of the suction side oscillations, which suggests that these oscillations are generated by flow pattern of the suction side.
Summary of peak tone frequencies and amplitudes for pressure signals collected at points at 96% chord and drag signal.
Turbulence Intensity | Suction side, 96% chord | Pressure side, 96% chord | Force (drag), airfoil | |||
---|---|---|---|---|---|---|
Frequency [Hz] | Amplitude [Pa] | Frequency [Hz] | Amplitude [Pa] | Frequency [Hz] | Amplitude (drag counts) | |
0.18% | 1,660 | 5.1 | 1,300 | 9.6 | 1,300 | 1.7 |
0.25% | 1,660 | 6.1 | 1,660 | 1.1 | 1,660 | 0.74 |
0.65% | 1,620 | 5.0 | 1,620 | 0.6 | 1,620 | 0.66 |
The final assessment of the URANS modelling is performed by comparing the frequency and amplitude of drag oscillations for S834 airfoil at different freestream speeds, and thus different Reynolds numbers; a summary of the comparison results is presented in Table 3. It is observed that the frequency of oscillations for different Reynolds numbers is well-matched with the peak tone frequency measured in the experiment. These results were obtained at the turbulence intensities of 0.44% and 0.69% for freestream velocity values of 22.4 m/s and 47.9 m/s, respectively. Hence, the intensities were in the regime, where the expected LBL-VS frequency matches the experiment tonal noise peak, as shown for the baseline case. The trend of increasing sound power level (PWL) with increasing amplitude of drag oscillations is also maintained; however, due to various turbulence levels between cases, further conclusions cannot be drawn.
Comparison of frequencies and amplitudes of drag oscillations for the baseline case (32 m/s) and additional cases at 22.4 m/s and 479 m/s, together with presentation of the experimental data pertaining to peak tone frequency (1/3 octave band centre) and PWL for comparison.
u∞ [m/s] | TI | fD [Hz] | Amplitude, drag counts | fpeak, EXP [Hz] | PWLEXP [dB] |
---|---|---|---|---|---|
22.4 | 0.44% | 1,120 | 0.18 | 1,000 | 69.36 |
32.0 | 0.65% | 1,620 | 0.74 | 1,600 | 71.32 |
47.9 | 0.69% | 2,560 | 1.88 | 2,508 | 71.44 |
The acoustic pressure field is presented in Figure 15 for the frequency of 1,600 Hz. The following conclusions can be drawn:
The noise source is located at the trailing edge of an airfoil. The hydrodynamic pressure oscillations appear at the laminar bubble and propagate downstream, over the trailing edge, and eventually form the wake. The wake does not contribute to the noise generation at the airfoil, and thus it is not required to participate in the calculation of aeroacoustic sources. The directivity patterns can be clearly observed with higher and lower amplitudes of acoustic pressure.
The quantitative comparison of noise levels obtained in the present work and in the experiment presented in the study of Oerlemans [14] is provided in Figure 16. The frequency of oscillations from the aerodynamic model was transferred to the noise tonal peak observed in the noise spectrum. The tone frequency from CAA matches the peak tone frequency measured in the experiment. The level of noise is, however, significantly higher than in the experiment. The difference can be attributed to the differences in the computational model setup, which assumes an ideal, infinite span airfoil, while the experiment was conducted on the finite wing. Noise levels’ differences at the remaining frequencies result from turbulence averaging in RANS. To model higher frequencies scale resolving CFD is needed, which was not in the scope of the present work.
In this work, the feasibility of using low-cost URANS simulations to model the airfoil LBL-VS tonal noise was investigated. In comparison with the LES or DNS approaches, which are typically used in conjunction with direct application of CAA, the simplified physics model captured with URANS constitutes one of the important limitations characterising this approach. The lack of aeroacoustic feedback typically responsible for triggering the hydrodynamic instabilities resulted in a lack of these for cases calculated with a turbulence intensity just below 1%. On the other hand, for lower turbulence intensity, the laminar bubble was sufficient to form the instabilities, which propagated downstream, forming a VS pattern on the suction side. The frequency of oscillations at moderate turbulence intensity matched that of the experiment. At the lowest turbulence intensities, the oscillations appeared at the pressure side, and with a different frequency compared to those at the suction side. The resulting force oscillations’ frequency was dominated by pressure side vortices, which represents a difference from the experimental results.
The main outcomes of this work are related to the sensitivity of analyses of LBL-VS to the freestream turbulence. The analyses showed a very significant sensitivity to freestream turbulence, with three different results within a 0%–1% turbulence intensity range: domination of pressure side vortices over suction side vortices (observed for TI = 0.18%), existence of suction side vortices only (observed for 0.25% ≤ TI ≤ 0.65%) and lack of vortices for higher turbulence. Even though some of the results match the experimental measurements very well, the usefulness of URANS in the capturing of airfoil tonal noise remains limited.