# Effects of freestream turbulence on cavity tone and sound source.

1. IntroductionFlows over open cavities such as the sunroofs of automobiles and landing gear configurations of airplanes often generate self-sustained oscillations and intense tonal sound. Therefore, to reduce the tonal sound from cavity is one of the most important issues for development of rapid transport vehicles such as airplanes or automobiles.

A mechanism for oscillations with fluid-acoustic interactions in compressible cavity flows was proposed by Rossiter [1]. The vortex impinging at the downstream edge of the cavity generates acoustic waves. These acoustic waves induce vortex shedding from the upstream edge of the cavity again. Sarohia [2] measured velocity fluctuations around cavities for laminar flows over shallow cavities (depth-to-length ratio of cavities D/L < 0.35 and freestream Mach number M < 0.07). He found that instability in the shear layer of the cavities was amplified by self-sustained oscillations. East [3] measured the acoustic pressure of deep cavities in the turbulent boundary layer (1.0 < D/L < 8.5 and M < 0.2). It was clarified that the self-sustained oscillations were amplified by coupling between instability in the shear layer and the acoustic mode in the direction of cavity depth. A cavity flow without acoustic resonance does not generate self-sustained oscillations or intense tonal sound in cavity flows at low Mach numbers in the turbulent boundary layer. Recently, direct numerical simulations were performed to investigate the three-dimensional stability of cavity flows [4].

The control of self-sustained oscillations in cavity flows and sound generated from cavities have also been investigated. Huang and Zhang [5] used plasma actuators implemented at the leading edge ahead of cavities in a recent study on the control of cavity flows and obtained significantly reduced levels of tonal sound. The flow induced by actuators traversed in the streamwise direction makes the shear layer thickness nonuniform in the spanwise direction.

Lusk et al. [6] investigated the attenuation effects of steady normal mass injection from the upstream of the leading edge of cavities on oscillations in supersonic cavity flows. The blowing flows from multisegmented slots spanning the cavity width broke up large-scale structures in cavity flows and reduced the level ofpressure fluctuations in the cavities. These results indicate that disturbance in the incoming boundary layer is important to reduce the cavity tone.

The disturbances in the boundary layer are also possibly caused by freestream turbulence, freestream vortices, and acoustic waves, as reviewed by Saric et al. [7]. Brandt et al. [8] studied the effects of the integral scale of freestream turbulence on bypass transitions of the boundary layer by using direct numerical simulations. They found that transitions occurred earlier for larger-scale freestream turbulence. The condition of the incoming boundary layer affects the cavity tone, and this could cause the past measured scatter of sound pressure spectra of cavity tone particularly at a low Mach number [9]. However, the effects of freestream turbulence on cavity flows and generated sound have not yet been clarified.

Also, it has been considered to be important to clarify the structures of the sound source to understand the mechanism for reducing the cavity tone. Larsson et al. [10] directly simulated two-dimensional flows and acoustic fields and estimated the sound source based on the Curle's acoustic analogy [11] by using fluctuating surface pressure. They found that the sound source of the cavity tone was maximal at the downstream edge of the cavity. Ask and Davidson [12] discussed the cancellation of the sources at the cavity bottom and downstream wall. However, the effects of the intensity of the Lighthill's acoustic source [13] and spanwise coherence of the source on the sound have not been investigated.

The present study directly simulated flow and acoustic fields for cavity flows with various intensities of freestream turbulence. The sound source of the cavity tone based on Lighthill's acoustic analogy [13] was estimated to clarify the reduction mechanism for the cavity tone by freestream turbulence. The effects of the weakening of the intensity of the sound source and the lowering of spanwise coherence of the sound source on the cavity tone were separately investigated.

2. Flow Configurations

Figure 1 outlines the configurations for flows around a two-dimensional cavity. The origin of the coordinate system is located at the upstream edge of the cavity. The streamwise direction is the x-axis, the vertical direction is the y-axis, and the spanwise direction intersecting these two axes is the z-axis. The cavity length, L, is 20 mm. The freestream velocity, [U.sub.[infinity]], is 30 m/s.

Table 1 summarizes the computational and experimental conditions. The depth-to-length ratio of the cavity is D/L = 0.5 and 2.5. In addition, acoustic resonance in the direction of cavity depth occurs in the flow over the cavity with D/L = 2.5, as explained in Section 5.2. The freestream Mach number, M [equivalent to] [U.sub.[infinity]]/c, is 0.09, where c denotes the speed of sound. The Reynolds number based on the cavity length, [Re.sub.L] [equivalent to] [U.sub.[infinity]]L/v, is 4.0 x [10.sup.4], where [nu] denotes the kinetic viscosity and that based on the cavity depth, [Re.sub.D] [equivalent to] [U.sub.[infinity]]D/[nu], is 2.0 x [10.sup.4] for D/L = 0.5 and 9.9 x [10.sup.4] for D/L = 2.5.

The velocity profile in the flow over the flat plate without a cavity at the streamwise position corresponding to the location of the upstream edge of the cavity with and without freestream turbulence is in good agreement with that of the laminar Blasius boundary layer. The ratio of the momentum thickness in the upstream boundary layer to the cavity length was slightly changed by freestream turbulence (turbulent intensity of Tu = 0.0-2.3%) and varied in the range of [theta]/L = 0.0071-0.0074 in the computation. This was approximately the same as that in the experiment ([theta]/L = 0.0072), where the freestream turbulence is 0.6%.

The features of freestream turbulence in the present paper are described by the turbulent intensity, Tu, and the integral scale, [L.sub.e]. The turbulent intensity, Tu, is calculated by using the rms values of the streamwise velocity fluctuations, [u.sub.rms]/[U.sub.[infinity]], in y/L = 0.2-0.5, where even the height of this lower limit is three times greater than the boundary layer thickness. The integral scale of the freestream turbulence, [L.sub.e], is calculated by using the autocorrelation function values of the streamwise velocity, u, in the same range of height. The intensity and integral scale are approximately constant in this range of height.

The computations were performed for turbulent intensity Tu = 0.0%, 0.9%, 1.5%, 1.8%, and 2.3%, where the integral scale was [L.sub.e]/L = 0.095 for Tu = 0.9% and [L.sub.e]/L = 0.111 for Tu = 2.3%. In addition, the ratio of the integral scale to the momentum thickness of the boundary layer, [L.sub.e]/[theta], was 13.3 for Tu = 0.9% and 15.5 for Tu = 2.3%. The turbulent intensity, Tu, was 0.6% in the experiment, and the integral scale, [L.sub.e]/L, was 0.165.

It has been confirmed that no laminar-turbulent transitions in the boundary layer occurred with this freestream turbulence in both computations and experiments. Also, the preliminary computations for the flow over the cavity with D/L = 0.5 for freestream turbulence of [L.sub.e]/L = 0.078 and 0.137 with turbulent intensity Tu = 1.5% indicated that the effects of difference of the integral scale on the reduced cavity tone were within 2 dB.

3. Experimental Methods

The experiments [14] were carried out using a suction-type, low-noise wind tunnel with a rectangular test section with a cross-section of 150 mm x 75 mm, as outlined in Figure 2. The intensity of freestream turbulence was less than 0.6% and nonuniformity of the freestream was less than [+ or -] 0.1% at the freestream velocity of 30 m/s. The background noise level was 58 dB (A) at the same velocity.

The test section was terminated in a spanwise manner by the end walls composed of acrylic plate and porous plate to visualize the flow within the cavity and suppress acoustic resonance in the spanwise direction. The velocity along the center of the cavity (x/L = 0.5) was measured using a hotwire anemometer. The acoustic pressure in the far field (x/L = 6.75, y/L = 23.5) was measured with a nondirectional 1/2 inch microphone.

4. Numerical Methods

4.1. Governing Equations and Finite Difference Formulations. To clarify the fluid-acoustic interactions in the cavity flows, flow and acoustic fields were simulated simultaneously by directly solving the three-dimensional compressible Navier-Stokes equations in the conserved form.

The spatial derivatives were evaluated using a sixth-order-accurate compact finite difference scheme [15] (fourth-order-accurate scheme at the boundaries). The time integration was performed by the third-order-accurate Runge-Kutta method.

In order to reduce the computational cost, large-eddy simulations (LES) were performed in the present study. A tenth-order-accurate spatial filter dissipated the energy that should be transferred to subgrid scale vortices. This filter also removed numerical instabilities due to the compact finite difference scheme [16]. The details on these computational methods have been presented in Yokoyama and Kato [17].

4.2. Computational Grids. Figure 3 shows the computational grid for the flow over the cavity with D/L = 2.5. The spanwise length of the computational domain is W/L = 1.0 and the spanwise grid resolution is [DELTA]z/L = 0.0125. It has been preliminarily confirmed that the computational results with W/L = 1.0 approximately agree with those with a wider computational domain (W/L = 2.0). The computational grid for the flow over the cavity with D/L = 0.5 was the same as that for the range of y/L [greater than or equal to] -0.5 for the flow over the cavity with D/L = 2.5. The total grid points are 18 million.

Figure 4 outlines the computational domain that is divided into vortical, acoustic, and buffer regions. The vortical region is from the upstream edge to the downstream edge of the cavity (0.0 [less than or equal to] x/L [less than or equal to] 1.0 and -0.1 [less than or equal to] y/L [less than or equal to] 0.1), where the shear layer is spread. The grid resolutions are [DELTA]x/L = [DELTA]y/L [less than or equal to] 0.01 and sufficiently fine to capture the large-scale vortical structures and active fine scale vortices in the cavity flow. The predicted velocity spectra will be discussed later in Section 5.1.

The region surrounded by the dashed line (-12.5 [less than or equal to] x/L [less than or equal to] 8.5 and y/L [less than or equal to] 23.5) in Figure 4 is the acoustic region. The acoustic region has an extent that is three and two times larger than the acoustic wavelength of the fundamental frequency of the cavity tone with D/L = 0.5 and 2.5, respectively. More than 10 points are used per acoustic wavelength that is two and four times greater than the fundamental frequency of the cavity tone with D/L = 0.5 and 2.5, respectively, in the grid resolutions in the acoustic region.

In the buffer region, the grid is sufficiently stretched to dissipate acoustic and voritcal disturbances near the artificial nonreflecting boundary conditions.

As discussed in Section 5, the comparison of the predicted flow and sound with those measured clarifies that it is possible to capture acoustic waves radiating from the cavity flow by using this computational grid and these numerical methods.

4.3. Boundary Conditions. Figure 4 also shows the boundary conditions. Nonreflecting boundaries [18-20] are used at the inflow and outflow boundaries. On the lower boundary in the upstream of the cavity, the boundary condition is changed from a slip wall to a nonslip and adiabatic wall. The position of this change was determined in order to set the boundary layer thickness to be that measured. Periodic boundary conditions are used in the spanwise direction.

A uniform steady flow is asymptotically imposed in the inflow buffer region (x/L [less than or equal to] -12.5) for the cavity flow without freestream turbulence (Tu = 0.0%). In the case with freestream turbulence, a homogeneous turbulence field is imposed on the uniform steady flow in the inflow buffer region.

This homogeneous turbulence field was separately predicted as box turbulence, where the computation starts from initial flow fields based on the von Karman spectrum [21]. These methods were described in detail in the previous paper [17]. The turbulent intensity and scale were controlled by using the parameters of the von Karman spectrum. Figure 5 plots the energy spectrum of homogeneous box turbulence. The Taylor scale of homogeneous turbulence is almost 20 mm.

4.4. Prediction of Acoustic Far Fields. The porous Ffowcs Williams and Hawkings (FW-H) method [22-24] was used to predict the acoustic pressure at the measurement point (x/L = 6.75 and y/L = 23.5) because predicting acoustic pressure in the far field with direct simulation greatly expends computational resources and time. The FW-H formula can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where r = [x.sub.2] - [x.sub.1], r = [absolute value of r], [??] = r/r, and [??] represent the outward unit normal vector of surface S. The subscript "ret" denotes consideration of retarded time [tau] = t - r/c. The pressure, p, was sampled at the semicylindrical surface of the radius, r = 2.5L, from the center of the cavity (x/L = 0.5 and y = 0), as outlined in Figure 6.

The spanwise computational domain, W = L, was smaller than that of the experiment, [W.sub.e] = 7.5L, to reduce computational resources. To take the effects of this difference into consideration, the sound pressure levels, SPL(f), were estimated by using the equivalent coherent length, [L.sub.c](f), following (2) [25-27]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

where [SPL.sub.s](f) is the sound pressure level predicted with the above porous FW-H method and the second and third terms on the right hand side are the correction terms. This length was determined by the coherence function values of the vertical velocity along the center of the cavity (x/L = 0.5, y = 0) in the spanwise direction, [[gamma].sup.2](f, [DELTA]z) = 1/3.

The coherence function was computed by the same way as that in the reference by Van Der Velden et al. [28] as indicated in (3):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

where [[PHI].sub.11] and [[PHI].sub.22] denote autopower spectra and [[PHI].sub.12] denote the cross-power spectrum between two points along a given dimensional line [DELTA]z = [z.sub.2] - [z.sub.1]. The equivalent coherent length, [L.sub.c](f), also has been confirmed to approximately agree with the distance for the phase difference in Lighthill's stress tensor at the same position of [phi] = [pi]/4 (45 degrees).

When the coherence is larger than 1/3 in the entire computational region, the equivalent coherent length was computed by linear extrapolation using the values in the range of [DELTA]z/L = 0.25-0.5, where the gradient of the coherence is approximately constant.

5. Validation of Computational Accuracy

5.1. Flow Field. The predicted flow and acoustic fields for the cavity flow without freestream turbulence (Tu = 0.0%) were compared with those measured because the measured turbulent intensity in the boundary layer was negligibly small.

Figure 7 shows the predicted and measured profiles of mean velocity and turbulent intensity at the center of the cavity (x/L = 0.5) for the cavity flow with D/L = 0.5, where the contributions of vertical velocity v to the velocity measured by a hot-wire anemometer were considered as [u.sub.h] = [([u.sup.2] + [(0.5v).sup.2]).sup.0.5] in the same way as our previous research [27]. It should be noted that the markers in Figure 7 in the predicted results represent only data to sample in the computational duration time, where the grid resolutions are finer. It is clarified that the predicted results agree well with those measured.

Figure 8 shows the predicted and measured power spectra of the velocity, [u.sub.h]/[U.sub.[infinity]], at the center of the cavity and the height for the maximal turbulent intensity (x/L = 0.5 and y/L = 0.02) for the cavity flow with D/L = 0.5. The frequency resolution of all spectral analyses is [DELTA]St [equivalent to] [DELTA]fL/[U.sub.0] = 0.104. The predicted spectrum was averaged 20 times, where the samples were overlapped by 50% and the duration time of the samples for spectral analysis was 100L/[U.sub.0]. Also, the experimental spectrum was averaged 9400 times with respect to time.

As shown in Figure 8, there is the main peak of [St.sub.f] = 1.56 (fundamental frequency) in both predicted and measured results, where the cavity tone radiates at this frequency. Moreover, the predicted power at the fundamental frequency is in good agreement with that measured.

To clarify the effects of the width of the computational domain, W, on the predicted flow fields, the coherence of the normal velocity at x/L = 0.5, y = 0 was computed for the flow fields predicted with different domains of W/L = 1.0 and 2.0. Figure 9 shows the predicted coherence at the fundamental frequency. It is presented that the distributions of the coherence predicted with the original width are in good agreement with those with the double width. This means that the independency of the computational results on the width of the domain supports the adoption of the periodic conditions.

5.2. Sound Pressure Level. Figures 10(a) and 10(b) show the predicted and measured sound pressure spectra at the points of the far field (x/L = 6.75 and y/L = 23.5) for D/L = 0.5 and 2.5, respectively. The measured level was affected by the background noise of the wind tunnel in the range of St < 0.5. The bars in the figures present how the results with the same averaging number of 20 as that of simulation vary around the measured value with the long averaging number of 9400.

The predicted spectra in both cavity flows with D/L = 0.5 and 2.5 are in good agreement with those measured. Also, the fundamental frequency for the cavity flow with D/L = 2.5 and [St.sub.f] = 0.94 is corresponding to the Rossiter mode of n = 2 [1], which is the streamwise number of vortices in the shear layer, while that for D/L = 0.5 and [St.sub.f] = 1.56 is corresponding to n = 3. The present peak also agrees with the resonant frequency of one-quarter wavelength mode in the direction of cavity depth, [St.sub.res] = 0.84 and 0.98, which are estimated by using the semiempirical formula by East [3] and by general open correction of a circular closed pipe of 8(D/2)/3[pi] [29]. The slight difference with the resonant frequency estimated by East [3] is possibly due to the conditions of the incoming boundary layer, where the cavity was placed in a turbulent boundary layer in the experiments by East.

To confirm the occurrence of the acoustic resonance in the cavity flow with D/L = 2.5, the phase of the pressure fluctuations in the cavity was investigated. Figure 11 shows the phase differences in the pressure fluctuations, [C.sub.p] = p/0.5 [[rho].sub.[infinity]][U.sub.[infinity].sup.2], at the fundamental frequency vertically along x/L = 0.5. The reference point for these phase differences is at the bottom of the cavity (y/L = -2.5). The phase of pressure fluctuations in the cavity is constant in the cavity (y/L < -0.1); that is, the standing waves due to acoustic resonance are generated in the cavity. These results indicate that acoustic resonance can be captured in the present computation.

To clarify the dependency of the predicted results on computational grid, the computations with two different meshes were also performed for cavity flow of D/L = 0.5 without freestream turbulence (Tu = 0.0%). One mesh consists of finer mesh of double grid points in the normal direction in the cavity, where the normal stretch of grid resolution is lower than that of original mesh. The other one is twice as wide as the present mesh as mentioned in Section 5.1.

Figure 12 shows the predicted sound pressure spectra with these meshes. As shown in this figure, the tonal sound is predicted at the same fundamental frequency of [St.sub.f] = 1.56 and the difference of level is within 5 dB. Therefore, it has been concluded that the dependency of the predicted results on the mesh is small.

It has been concluded along with the discussions of flow field in Section 5.1 that the present computations could adequately capture the flow and acoustic fields.

6. Results and Discussion

6.1. Shear Layer in Cavity Flow. Figure 13 shows the predicted mean velocity and turbulent intensity at the center of the cavity (x/L = 0.5) for the cavity flow with D/L = 0.5 with and without freestream turbulence. As shown in Figure 13(a), it is clarified that the shear layer is thickening with freestream turbulence.

The momentum thickness [[delta].sub.m] at x/L = 0.5 was computed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where [y.sub.min]/L and [y.sub.max]/L were set to be -0.2 and 0.2. This range is the same as that in reference [30]. The momentum thickness was clarified to become thicker from [[delta].sub.m]/L = 0.019 (Tu = 0.0%) to 0.027 (Tu = 2.3%) by the freestream turbulence. Figure 13(b) also shows that the turbulent intensity becomes more intense by the freestream turbulence.

Figure 14 shows the predicted profiles for the cavity flow with D/L = 2.5. As shown in this figure, the shear layer becomes thicker, where the momentum thickness was changed from [[delta].sub.m]/L = 0.041 (Tu = 0.0%) to 0.049 (Tu = 2.3%) by the freestream turbulence. Figure 14 shows that the peak of the turbulent intensity becomes slightly weaker by the freestream turbulence. This is because the velocity fluctuations at the fundamental frequency become weaker. The effects of the freestream turbulence on velocity spectra are discussed in the next section. Also, as shown in Figure 14, the region for intense turbulent intensity is spread in the normal direction by the freestream turbulence of Tu = 2.3%. As a result, the mixing in the shear layer is enhanced and the shear layer becomes thicker.

6.2. Vortical Structures. Figure 15 shows instantaneous flow fields around the cavity for Tu = 0.0% and 2.3% for the cavity flow without acoustic resonance (D/L = 0.5). Figure 16 shows those for the cavity flow with acoustic resonance (D/L = 2.5).

Large-scale vortices are found for the cavity flow without freestream turbulence (Tu = 0.0%) for both cavity flows with and without acoustic resonance. These vortices are shed from the upstream edge of the cavity at the fundamental frequency and contribute to the periodic acoustic radiation from the cavity.

The large-scale vortical structures, which are coherent in the spanwise direction, are dominant for Tu = 0.0% particularly in the cavity flow with acoustic resonance. Also, streamwise rib structures and fine scale eddies in the shear layer can be observed for Tu = 0.0% for the cavity flow without acoustic resonance. As the turbulent intensity is increased, the fine scale vortices become active regardless of acoustic resonance.

Figures 17(a) and 17(b) show the power spectra of the vertical velocity, v/[U.sub.[infinity]], at the center of the cavity (x/L = 0.5 and y = 0), where large-scale vortices are observed without freestream turbulence, for the cavity flow without acoustic resonance and that with acoustic resonance, respectively.

The power spectra for freestream turbulence Tu = 0.9% are almost the same as that for the cavity flow without freestream turbulence (Tu = 0.0%) in both cases with and without acoustic resonance. The level of velocity fluctuations at the fundamental frequency becomes significantly lower for Tu = 2.3%. This level corresponds to the power of velocity fluctuations due to the large-scale vortical structures, which became the sound source of the cavity tone.

These spectra are also compared with the curve for -5/3 power. All the slopes of the power spectra of velocity fluctuations with and without freestream turbulence in the range of St = 1-5 are along the curve of -5/3 power for the cavity flow without acoustic resonance. Also, the level in that range is approximately independent of freestream turbulence.

As shown in Figure 17(b), for the cavity flow with acoustic resonance, the slope of the spectrum becomes close to the curve for -5/3 power by adding freestream turbulence. Figures 18(a) and 18(b) show the power spectra near the downstream edge of the cavity (x/L = 0.9 and y = 0) for the cavity flow without acoustic resonance and that with acoustic resonance, respectively. Figure 18(b) shows that the slope of the power is clearly along the curve of-5/3 power in the range of St = 1-5 particularly for Tu = 2.3% for the cavity flow with acoustic resonance. Also, the power at high frequencies of St [greater than or equal to] 3 becomes more intense with the freestream turbulence of Tu = 2.3%.

The above-mentioned results indicate that turbulent mixing in the shear layer is enhanced by freestream turbulence, and the turbulent transition in the shear layer occurs due to freestream turbulence.

6.3. Effects of Freestream Turbulence on Cavity Tone. Figure 19 shows sound pressure spectra in the acoustic far field (x/L = 6.75 and y/L = 23.5) for cavity flows with and without acoustic resonance. The differences in the tonal sound pressure levels at the fundamental frequency, [DELTA]SPL, with reference to that for the cavity flow without freestream turbulence are also presented.

The [DELTA]SPL is -2.4 dB and -15.1 dB for Tu = 0.9 and 2.3%, respectively, for the cavity flow without acoustic resonance. The [DELTA]SPL is -0.3 dB and -20.0 dB for Tu = 0.9 and 2.3%, respectively, for the cavity flow with acoustic resonance.

The effects of freestream turbulence on the cavity tone are less than 2.5 dB for both cavity flows with and without acoustic resonance for freestream turbulence of Tu = 0.9%. Tonal sound generated from the cavity flow with acoustic resonance is greatly affected by freestream turbulence for Tu = 2.3%.

The weakening of the intensity and lowering of the spanwise coherence of the sound source are considered to be responsible for the reduced cavity tone. The following sections discuss the effects of the freestream turbulence on the intensity and spanwise coherence of the sound source.

6.4. Intensity of Sound Source. Aerodynamic sound is related to the intensity of fluctuating pressures/velocities and their structures (coherent motion) according to acoustic analogies [11, 13, 22]. In addition, the velocity fluctuations near the downstream edge of the cavity contribute to the sound pressure level of the cavity tone since the cavity tone is generated by impinging large-scale vortical structures at the downstream edge of the cavity [17].

The Lighthill's equations [13] are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

[T.sub.ij] [approximately equal to] [rho][u.sub.i][u.sub.j], (6)

where tensor [T.sub.ij] is Lighthill's stress tensor. This tensor is approximately equal to the Reynolds stress term at the low Mach number [31] as shown in (6). In the present section, the intensity of the sound source was estimated by using the Reynolds stress term of Lighthill's stress tensor [13] at the fundamental frequency near the downstream edge (0.75 [less than or equal to] x/L [less than or equal to] 1.0). The range of 0.75 [less than or equal to] x/L [less than or equal to] 1.0 is approximately equal to the scale of the large-scale vortical structures.

Figures 20 and 21 show the intensity of the abovementioned sound source (the Reynolds stress term of Lighthill's stress tensor) at the fundamental frequency, around the cavity. Also, Figure 22 shows the power level of the sound source at the fundamental frequency, [P.sub.f], along y = 0 from the upstream edge to the downstream edge. As shown in the figure, the acoustic source is weakened by the freestream turbulence.

Moreover, the sound reduction level of the cavity tone due to the weakening of the intensity of the sound source, [DELTA][SPL.sub.i], was calculated by the integration of the power of source, [P.sub.f], in the range of 0.75 [less than or equal to] x/L [less than or equal to] 1.0. The equation of [DELTA][SPL.sub.i] is as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

where the subscript "0" means the condition of the cavity flow without freestream turbulence (Tu = 0.0%).

Table 2 shows the values of [DELTA][SPL.sub.i] in each condition. The value of [DELTA][SPL.sub.i] is -8.5 dB and -8.2 dB for Tu = 2.3% with and without acoustic resonance, respectively, which represent greater effects than those for Tu = 0.9% (-0.3 dB and 0.8 dB). Also, this means that the cavity tone is greatly reduced by weakening of the sound source for Tu = 2.3%.

6.5. Spanwise Coherence of Sound Source. Two-dimensionality in the spanwise direction of large-scale vortical structures as shown in Figures 15 and 16 is quantitatively discussed. In the present section, the effects of freestream turbulence on the spanwise coherence of the sound source are discussed by computing the coherence of the vertical velocity, v/[U.sub.[infinity]], at the fundamental frequency at the center of the cavity (x/L = 0.5andy = 0), which is shown in Figure 23.

As shown in Figure 23, the coherence is decreased by the freestream turbulence for both cases with and without acoustic resonance. The equivalent coherent length, [L.sub.c]/L, which was discussed in Section 4.4, is 0.22, 0.12, and 0.04 for the flow with freestream turbulence of Tu = 0.0,0.9, and 2.3%, respectively, in the cavity flow without acoustic resonance, whereas it is 6.40, 5.76, and 0.05 for Tu = 0.0, 0.9, and 2.3% for the cavity flow with acoustic resonance.

Equivalent coherent length for the cavity flow with acoustic resonance without freestream turbulence, [L.sub.c]/L = 6.40, is greater than that for the flow without acoustic resonance ([L.sub.c]/L = 0.22). This means that the two-dimensionality of the large-scale vortices related with cavity tone becomes higher due to the acoustic resonance. As shown in Figures 15 and 16, three-dimensional fine scale vortices are less active for cavity flow with acoustic resonance.

In Section 4.4, the method for correction of SPL for the difference between the computational and experimental spanwise extent is discussed. With a similar method, the effects of lowering spanwise coherence of the sound source on the cavity tone, [DELTA][SPL.sub.c], are estimated as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

where [[lambda].sub.v] is the streamwise scale of vortices related to cavity tone and estimated to be approximately L/(4n) by the distance of the streamwise phase difference of velocity fluctuations of [DELTA][phi] < [pi]/2 in the shear layer.

Table 2 also shows the value of [DELTA][SPL.sub.c] for each case. It is -0.6 dB and -4.2 dB for Tu = 0.9 and 2.3% in the cavity flow without acoustic resonance, while it is -0.5 dB and -17.1 dB for Tu = 0.9 and 2.3% for the cavity flow with acoustic resonance. These results present that the effects of the freestream turbulence in the cavity flow with acoustic resonance were greater than those in the flow without acoustic resonance. This is possibly because the turbulent transition occurs in the free shear layer of the flow with acoustic resonance while the shear layer is turbulent even without freestream turbulence in the flow without acoustic resonance.

6.6. Contributions of Intensity and Spanwise Coherence of Sound Source. Both the intensity and spanwise coherence of the sound source affect variation of the sound pressure level by freestream turbulence, as was previously mentioned. The contributions of these two factors are discussed in the present section. Table 2 and Figure 24 help to explain the contributions of various factors to the cavity tone. The total difference in the cavity tone due to contributions by these factors, [DELTA][SPL.sub.total], is calculated as

[DELTA][SPL.sub.total] = [DELTA][SPL.sub.i] + [DELTA][SPL.sub.c]. (9)

Figure 24 shows that the reduction level estimated with all factors, [DELTA][SPL.sub.total], is comparable to the predicted one from the direct simulation, [DELTA]SPL, for each case. Also, we found qualitative agreement such that the effects of freestream turbulence on the cavity tone were greater in the cavity flow with acoustic resonance than that without acoustic resonance. This is because the reduced sound pressure level due to the reduction of the coherence is greater for the cavity flow with acoustic resonance, where the turbulent transition occurs. Also, this agreement shows that the present method for estimating the effects of the intensity and coherence of the sound source on the cavity tone is reasonable.

7. Conclusions

Direct simulations of flow and acoustic field were carried out for cavity flows with and without acoustic resonance (D/L = 0.5 and 2.5) under various freestream turbulent conditions to clarify the effects of freestream turbulence on the cavity tone. The freestream Mach number was M = 0.09 and the Reynolds number based on the cavity length was [Re.sub.L] = 4.0 x [10.sup.4]. The incoming boundary layer was laminar with momentum thickness [theta]/L = 0.0071-0.0074. The effects of turbulent intensities of Tu = 0.0-2.3% on flow and acoustic fields were investigated.

The sound source was estimated by using the flow field based on Lighthill's acoustic analogy. Moreover, the reduced cavity tone was decomposed into the effects of the intensity and spanwise coherence of the sound source to enable the reduction mechanism of the cavity tone to be understood. The total reduction level by summing these contributions by the intensity and coherence was clarified to be comparable to that predicted by direct simulations. Also, these analytical results can explain well the effects of the acoustic resonance and freestream turbulence on the reduced level, which were predicted from direct simulations. Thus, this indicates that the method for estimating the effects of the intensity and coherence of the sound source on the cavity tone in the present study is reasonable.

The effects of freestream turbulence on the cavity tone in a cavity flow with acoustic resonance, -20.0 dB, were greater than those in a cavity flowwithout acoustic resonance, -15.1 dB, for a cavity flow with freestream turbulence of Tu = 2.3%. This was caused by the difference in the flow field for the cavity flow without freestream turbulence Tu = 0.0%. That is, the large-scale vortical structures in the cavity flow with acoustic resonance were extremely coherent in the spanwise direction due to acoustic resonance, whereas they were less coherent even for Tu = 0.0% for the cavity flow without acoustic resonance. As a result of this difference in coherence, the cavity flow with acoustic resonance was affected more by freestream turbulence.

The analyses of the intensity and spanwise coherence of the sound source clarified that the effects of the intensity and spanwise coherence of the sound source, on the other hand, were almost the same as those without acoustic resonance. The effects of spanwise coherence of the sound source on the cavity tone were greater in the cavity flow with acoustic resonance than in that without resonance.

http://dx.doi.org/10.1155/2016/7347106

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The present study was supported by JSPS KAKENHI Grant nos. 24760134 and 26820044 and through the Next-Generation Supercomputer Strategy Program made available by the Ministry of Education, Culture, Sports, Science, and Technology of Japan (MEXT).

References

[1] J. E. Rossiter, "Wind-tunnel experiments on the flow over rectangular cavities at subsonic and transonic speeds," Aeronautical Research Council Report 3438, 1964.

[2] V. Sarohia, "Experimental investigation of oscillations in flows over shallow cavities," AIAA Journal, vol. 15, no. 7, pp. 984-991, 1977.

[3] L. F. East, "Aerodynamically induced resonance in rectangular cavities," Journal ofSound and Vibration, vol. 3, no. 3, pp. 277-287, 1966.

[4] G. A. Bres and T. Colonius, "Three-dimensional instabilities in compressible flow over open cavities," Journal of Fluid Mechanics, vol. 599, pp. 309-339, 2008.

[5] X. Huang and X. Zhang, "Streamwise and spanwise plasma actuators for flow-induced cavity noise control," Physics of Fluids, vol. 20, no. 3, Article ID 037101, 2008.

[6] T. Lusk, L. Cattafesta, and L. Ukeiley, "Leading edge slot blowing on an open cavity in supersonic flow," Experiments in Fluids, vol. 53, no. 1, pp. 187-199, 2012.

[7] W. S. Saric, H. L. Reed, and E. J. Kerschen, "Boundary-layer receptivity to freestream disturbances," Annual Review of Fluid Mechanics, vol. 34, no. 1, pp. 291-319, 2002.

[8] L. Brandt, P. Schlatter, and D. S. Henningson, "Transition in boundary layers subject to free-stream turbulence," Journal of Fluid Mechanics, vol. 517, pp. 167-198, 2004.

[9] T. Colonius and S. K. Lele, "Computational aeroacoustics: progress on nonlinear problems of sound generation," Progress in Aerospace Sciences, vol. 40, no. 6, pp. 345-416, 2004.

[10] J. Larsson, L. Davidson, M. Olsson, and L.-E. Eriksson, "Aeroacoustic investigation of an open cavity at low mach number," AIAA Journal, vol. 42, no. 12, pp. 2462-2473, 2004.

[11] N. Curle, "The influence of solid boundaries upon aerodynamic sound," Proceedings of the Royal Society. London. Series A. Mathematical, Physical and Engineering Sciences, vol. 231, pp. 505-514, 1955.

[12] J. Ask and L. Davidson, "Sound generation and radiation of an open two-dimensional cavity," AIAA Journal, vol. 47, no. 6, pp. 1337-1349, 2009.

[13] M. J. Lighthill, "On sound generated aerodynamically. I. General theory," Proceedings of the Royal Society. London. Series A. Mathematical, Physical and Engineering Sciences, vol. 211, pp. 564-597, 1952.

[14] K. Terao, H. Yokoyama, Y. Ogoe, and A. Iida, "Proposition of new formula for frequency prediction based on generation mechanism of aerodynamic sound in cavity flows," Transactions of the Japan Society of Mechanical Engineers B, vol. 77, no. 779, pp. 1522-1532, 2011 (Japanese).

[15] S. K. Lele, "Compact finite difference schemes with spectral-like resolution," Journal of Computational Physics, vol. 103, no. 1, pp. 16-42, 1992.

[16] K. Matsuura and C. Kato, "Large-eddy simulation of compressible transitional flows in a low-pressure turbine cascade," AIAA Journal, vol. 45, no. 2, pp. 442-457, 2007

[17] H. Yokoyama and C. Kato, "Fluid-acoustic interactions in self-sustained oscillations in turbulent cavity flows, I. Fluid-dynamic oscillations," Physics of Fluids, vol. 21, no. 10, Article ID 105103, 2009.

[18] K. W. Thompson, "Time dependent boundary conditions for hyperbolic systems," Journal of Computational Physics, vol. 68, no. 1, pp. 1-24, 1987

[19] T. J. Poinsot and S. K. Lele, "Boundary conditions for direct simulations of compressible viscous flows," Journal of Computational Physics, vol. 101, no. 1, pp. 104-129, 1992.

[20] J. W. Kim and D. J. Lee, "Generalized characteristic boundary conditions for computational aeroacoustics," AIAA Journal, vol. 38, no. 11, pp. 2040-2049, 2000.

[21] J. O. Hinze, Turbulence, McGraw-Hill, New York, NY, USA, 1975.

[22] J. E. F. Williams and D. L. Hawkings, "Sound generation by turbulence and surfaces in arbitrary motion," Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering, vol. 264, no. 1151, pp. 321-342, 1969.

[23] A. S. Lyrintzis, "Surface integral methods in computational aeroacoustics--from the (CFD) near-field to the (Acoustic) far-field," International Journal of Aeroacoustics, vol. 2, no. 2, pp. 95-128, 2003.

[24] M. Shur, P. Spalart, and M. Strelets, "Noise prediction for increasingly complex jets--part I: methods and tests," International Journal of Aeroacoustics, vol. 4, no. 3, pp. 213-246, 2005.

[25] C. Kato, A. Iida, Y. Takano, H. Fujita, and M. Ikegawa, "Numerical prediction of aerodynamic noise radiated from low Mach number turbulent wake," in Proceedings of the 31st Aerospace Sciences Meeting and Exhibit, Reno, Nev, USA, 1993.

[26] K. W. Chang, J. H. Seo, Y. J. Moon, and M. Roger, "Prediction of flat plate self-noise," in Proceedings of the 12th AIAA/CEAS Aeroacoustics Conference, pp. 1451-1464, May 2006.

[27] H. Yokoyama, K. Kitamiya, and A. Iida, "Flows around a cascade of flat plates with acoustic resonance," Physics of Fluids, vol. 25, no. 10, Article ID 106104, 2013.

[28] W. C. P. Van DerVelden, A. H. Van Zuijlen, A. T. De Jong, and H. Bijl, "On the estimation of Spanwise pressure coherence of a turbulent boundary layer over a flat plate," in Proceedings of 11th World Congress on Computational Mechanics (WCCM '14), pp. 5710-5721, July 2014.

[29] J. W. S. Rayleigh, The Theory of Sound, Dover, Mineola, NY, USA, 1945.

[30] L. Larcheveque, P. Sagaut, I. Mary, O. Labbe, and P. Comte, "Large-eddy simulation of a compressible flow past a deep cavity," Physics of Fluids, vol. 15, no. 1, pp. 193-210, 2003.

[31] M. E. Goldstein, Aeroacoustics, McGraw-Hill, New York, NY, USA, 1976.

Hiroshi Yokoyama, Hiroshi Odawara, and Akiyoshi Iida

Department of Mechanical Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku, Aichi 441-8580, Japan

Correspondence should be addressed to Hiroshi Yokoyama; h-yokoyama@me.tut.ac.jp

Received 6 July 2016; Accepted 10 November 2016

Academic Editor: William W. Liou

Caption: Figure 1: Configurations for flows around two-dimensional cavity. (a) D/L = 0.5. (b) D/L = 2.5.

Caption: Figure 2: Schematics for experimental setup.

Caption: Figure 3: Computational grids for cavity flow with D/L = 2.5 (every second grid line is shown).

Caption: Figure 4: Computational domain and boundary conditions.

Caption: Figure 5: Energy spectra of homogeneous box turbulence.

Caption: Figure 6: Schematics for prediction of acoustic far field by porous FW-H.

Caption: Figure 7: Predicted and measured profiles of mean velocity (a) and turbulent intensity (b) for the cavity flow with D/L = 0.5.

Caption: Figure 8: Predicted and measured power spectra of velocity fluctuations (x/L = 0.5 and y/L = 0.02) for cavity flow with D/L = 0.5.

Caption: Figure 9: Predicted coherence of velocity, v/[U.sub.[infinity]], at fundamental frequency in spanwise direction (x/L = 0.5 and y = 0) for cavity flow (D/L = 0.5) with original and wider computational domains of W/L = 1.0 and 2.0.

Caption: Figure 10: Comparison of predicted and measured sound pressure spectra (x/L = 6.75 and y/L = 23.5). (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 11: Predicted phase differences of pressure fluctuations, [C.sub.p] [equivalent to] p/0.5 [[rho].sub.[infinity]][U.sub.[infinity].sup.2], at fundamental frequency vertically along x/L = 0.5 for cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 12: Comparison of predicted sound pressure spectra (x/L = 6.75 and y/L = 23.5) with original, finer, and wider meshes for cavity flow without freestream turbulence (D/L = 0.5).

Caption: Figure 13: Predicted profiles of mean velocity (a) and turbulent intensity (b) for the cavity flow with D/L = 0.5.

Caption: Figure 14: Predicted profiles of mean velocity (a) and turbulent intensity (b) for the cavity flow with D/L = 2.5.

Caption: Figure 15: Instantaneous isosurfaces of second invariant of velocity tensor gradient (q/[([U.sub.[infinity]]/L).sup.2] = 20) for cavity flows without acoustic resonance (D/L = 0.5). (a) Cavity flow without freestream turbulence (Tu = 0.0%). (b) Cavity flow with freestream turbulence, Tu = 2.3%.

Caption: Figure 16: Instantaneous isosurfaces of second invariant of velocity tensor gradient (q/[([U.sub.[infinity]]/L).sup.2] = 20) for cavity flows with acoustic resonance (D/L = 2.5). (a) Cavity flow without freestream turbulence (Tu = 0.0%). (b) Cavity flow with freestream turbulence, Tu = 2.3%.

Caption: Figure 17: Power spectra of velocity v/[U.sub.[infinity]] (x/L = 0.5 and y = 0). (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 18: Power spectra of velocity v/[U.sub.[infinity]] (x/L = 0.9 and y = 0). (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 19: Sound pressure spectra (x/L = 6.75 and y/L = 23.5) and difference in tonal sound pressure levels due to effects of freestream turbulence. (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 20: Contours of real part of power of Lighthill's stress tensor at fundamental frequency, [[partial derivative].sup.2][T.sub.ij]/[partial derivative][x.sub.i][partial derivative][x.sub.j]/([[rho].sub.[infinity]][U.sub.[infinity].sup.2]/[L.sup.2]), around cavity for cavity flows without acoustic resonance (D/L = 0.5). (a) Cavity flow without freestream turbulence (Tu = 0.0%). (b) Cavity flow with freestream turbulence, Tu = 2.3%.

Caption: Figure 21: Contours of real part of power of Lighthill's stress tensor at fundamental frequency, [[partial derivative].sup.2][T.sub.ij]/[partial derivative][x.sub.i][partial derivative][x.sub.j]/([[rho].sub.[infinity]][U.sub.[infinity].sup.2]/[L.sup.2]), around cavity for cavity flows with acoustic resonance (D/L = 2.5). (a) Cavity flow without freestream turbulence (Tu = 0.0%). (b) Cavity flow with freestream turbulence, Tu = 2.3%.

Caption: Figure 22: Power level of Lighthill's stress tensor at the fundamental frequency, [P.sub.f], along y = 0. (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 23: Coherence of velocity v/[U.sub.[infinity]] at fundamental frequency in spanwise direction (x/L = 0.5 and y = 0). (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Caption: Figure 24: Effects of intensity and spanwise coherence of sound source on variation of cavity tone. (a) Cavity flow without acoustic resonance (D/L = 0.5). (b) Cavity flow with acoustic resonance (D/L = 2.5).

Table 1: Computational and experimental conditions. D/L Tu [%] [L.sub.e]/L [theta]/L M Computation 0.5 0.0 -- 0.0071 0.09 0.9 0.095 0.0072 1.5 0.078 0.0074 1.8 0.090 0.0077 2.3 0.111 0.0078 Computation 2.5 0.0 -- 0.0071 0.09 0.9 0.095 0.0072 1.5 0.078 0.0074 1.8 0.090 0.0077 2.3 0.111 0.0078 Experiment 0.5 0.6 0.165 0.0072 0.09 2.5 D/L [Re.sub.L] [Re.sub.D] Computation 0.5 4.0 x [10.sup.4] 2.0 x [10.sup.4] Computation 2.5 4.0 x [10.sup.4] 9.9 x [10.sup.4] Experiment 0.5 4.0 x [10.sup.4] 2.0 x [10.sup.4] 2.5 9.9 x [10.sup.4] Table 2: Contributions of intensity and spanwise coherence of sound source to the sound pressure level. D/L Tu [%] [L.sub.c]/L [DELTA]SPL [dB] [DELTA][SPL.sub.i] [dB] 0.5 0 0.22 -- -- 0.9 0.12 -2.4 -0.3 2.3 0.04 -15.1 -8.5 2.5 0 6.40 -- -- 0.9 5.76 -0.3 0.8 2.3 0.05 -20.0 -8.2 D/L Tu [%] [DELTA][SPL.sub.c] [dB] [DELTA][SPL.sub.total] [dB] 0.5 0 -- -- 0.9 -2.7 -3.0 2.3 -4.2 -12.7 2.5 0 -- -- 0.9 -0.5 0.3 2.3 -17.1 -25.3

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Yokoyama, Hiroshi; Odawara, Hiroshi; Iida, Akiyoshi |

Publication: | International Journal of Aerospace Engineering |

Date: | Jan 1, 2016 |

Words: | 8145 |

Previous Article: | Numerical and experimental investigation of computed tomography of chemiluminescence for hydrogen-air premixed laminar flames. |

Next Article: | Erratum to "aeroelastic analysis of wings in the transonic regime: planform's influence on the dynamic instability". |