1. Introduction
Methods to reduce the noise generated by aircraft propulsion systems have long been a topic of research within the aerospace community. There have been many challenges faced by designers attempting to solve these types of problems due to the complexity of the flow phenomena as well as the costs involved with testing potential solutions. In some cases, noise reductions have only come at the expense of aircraft performance in other metrics (e.g., weight, cost).
Jet–surface interaction noise is usually associated with an increase in low frequency noise as a result of the interaction between a turbulent jet and a flat surface that is parallel to the jet axis. Jet noise is generally considered to be the result of two sources: scrubbing noise and trailing edge scattering noise (Brown and Wernet [
1]). Scrubbing is generated by pressure fluctuations in the turbulent boundary layer impinging on a flat surface. Trailing edge scattering, which is the dominant noise source at low frequency (Brown [
2], Podboy [
3]), is produced when the hydrodynamic component of the upstream turbulence is scattered into acoustic waves after interacting with the trailing edge. The total installation noise can propagate upstream or downstream of the trailing edge.
High-speed jet exhausts interacting with nearby surfaces can generate flow distortions that significantly increase the overall noise emitted. In one such example, an aircraft engine placed on top of the wing in an attempt to shield the jet noise that would propagate to the ground while also improving aircraft performance actually resulted in an overall increase in sound levels. Another example of note is the tremendous noise generated by the jet exhausts of military aircraft interacting with the various surfaces of an aircraft carrier deck during take-off and landing operations.
Considerable research has been conducted on axisymmetric jets (Bridges [
4], Brown and Bridges [
5]), but, since the 1980s, a number of studies have shifted their focus to non-circular jet engine exit designs as potential passive flow control devices (Gutmark and Grinstein [
6]). Geometric corrugations applied at the nozzle exit plane (e.g., tabbed studied by Gao et al. [
7], cruciform by El Hassan et al. [
8], and chevrons by Violato et al. [
9]) have been observed to break down the large-scale structures present in the jet exhaust flow. This subsequently reduces the mixing between the jet and ambient-fluid. It was noted that these effects were most pronounced for larger wetted perimeters of the jet exit (Shakouchi and Iryama [
10]) and an increased number of corners/sides (Gutmark and Grinstein [
6]). The effect of external expansion ramps on supersonic flow emitted from high-aspect-ratio rectangular nozzles was analyzed by Malla et al. [
11] to determine if the ramps could be used to reduce noise. Research conducted by Mancinelli et al. [
12] and Proenca et al. [
13] show further results regarding the interaction of rectangular jets with flat plates. Other aircraft propulsion concepts have utilized acoustically treated engine parts (inlet ducts, exhaust ducts, and inner walls), but this type of treatment tends to have a negative overall impact due to additional cost and weight.
While static chevron configurations at the jet exit have been shown to provide a significant reduction in noise without noticeable loss in thrust, some researchers have been looking into potential acoustic benefits that could be gained by exciting the chevrons with piezoelectric actuators. Mechanical perturbations of the flow by the actuated chevrons are believed to increase the production of small-scale disturbances while diminishing the large-scale turbulent structures thought to be responsible for the dominant portion of jet mixing noise. Mohan et al. [
14] introduced piezoelectric actuators to chevrons in an attempt to modify the growth rate of the mixing layer. Butler and Calkins [
15] examined four different types of nozzles that exist: round, faceted, faceted with static chevrons, and faceted with active chevrons. Their results demonstrated a 2 to 4 dB reduction in noise with the static chevrons, and an additional 2 dB reduction with the actuated chevron configurations.
Zaman [
16] studied the integration of microjets near the engine exhaust port and found that a clear noise reduction was observed as the microjet pressure increased. The results showed that turbulent mixing noise reduction, as monitored by the overall sound pressure level at a shallow angle, correlated with the ratio of the microjet to the primary jet driving pressures (normalized by the ratio of the corresponding diameters). Semlitsch et al. [
17] used implicit large eddy simulation to analyze the relationship between screech tone frequency and fluidic injection pressure of the microjets.
Rego et al. [
18] utilized a Lattice–Boltzmann Method (LBM) to carry out numerical simulations of a flat plate placed in an irrotational hydrodynamic field of a jet with a round nozzle. Three cases were investigated with input flow characteristics determined from NASA wind tunnel experiments (see Brown [
19]). Heated and cooled jets were analyzed at different Mach numbers, and the far-field noise was computed using the Ffowcs-Williams and Hawkings equation. Far-field spectral results showed a large noise increase for low to mid frequencies due to upstream hydrodynamic waves in the jet that propagate along the plate, ’scrubbing’ over the surface, until it reaches the trailing edge of the plate. (This is referred to as the ’gust solution’ in Afsar et al. [
20]). As these waves interact with the trailing edge and the downstream turbulent mixing layer, they are scattered in the far-field as noise. A noticeable change in the far-field noise was observed when the trailing edge location was changed relative to the fixed nozzle. As the length of the plate was increased, the trailing edge was positioned in a region dominated by large-scale structures and a greater degree of low-frequency amplification was observed. When the plate was moved closer to the jet in the radial direction, the sound levels at lower frequencies increased.
Behrouzi and McGuirk [
21] conducted experiments at the Loughborough University high-pressure nozzle test facility on a rectangular nozzle exhausting over a rectangular plate acting as an aft-deck. (A configuration similar to the one used in the current study.) The focus of the experiment was to study the effect of the rectangular plate on flow field development and to determine if any acoustic benefits could be obtained from the shielding effect that the aft-deck might produce. Schlieren imaging, pneumatic Pitot probes, and nonintrusive Laser Doppler Anemometry measurements were used to capture detailed information on the flow structures. The data revealed that the aft-deck created an asymmetry in the entrainment characteristics of the shear layers, inducing a net transverse pressure force on the plume by changing the inviscid shock cell structure. The presence of the aft-deck altered the net pressure force exerted on the flow which led to a transverse deflection of the jet. The aft-deck configuration produced a dramatic effect on the plume development, extending the potential core length slightly and reducing turbulence levels in the plume near field.
Bridges [
22] further investigated this configuration by studying whether beveled edges on the rectangular nozzle would make a difference in the sound generated when compared to an aft-deck being added to the propulsion configuration (see also Zaman [
23] and Zaman et al. [
24]). Nozzles were fabricated with different aspect ratios of 2:1, 4:1, and 8:1. For each aspect ratio, three nozzle designs were tried: one basic nozzle with no extension and two beveled variants with bevel lengths of 1.3 and 2.7 inches. A rectangular nozzle with an aft-deck surface was also considered in which the rectangular nozzles had plates fitted to them such that when the plate surface was even with the inner lip of the nozzle, the surface and nozzle were effectively one piece. Five aft-deck lengths were tested: 1.3, 2.7, 4, 8, and 12 inches long. When the beveled nozzle was analyzed, it was found that the noise levels decreased in both azimuthal planes as the aspect ratio was increased, with much more dramatic reductions in the nozzle’s major axis plane. These reductions in noise were lost when one lip of the nozzle was extended to make an aft deck. Low and high frequencies were examined with the aft deck extensions. For low frequencies, the sound levels were roughly the same on both sides of the plate. However, there were noticeable differences in the sound levels for higher frequencies on each side of the plate. This is due to the effect of the noise being shielded and reflected. The shielding effect seemed to increase as the length of the aft deck increased, but this was not the case for reflection. In addition, it is important to note that, for the largest plate lengths, the low-frequency amplification above the baseline jet noise was not dependent on the aspect ratio.
Seiner and Manning [
25] considered the interaction between a supersonic jet and a flat plate from a rectangular nozzle. They showed that the distance between the nozzle exit and flat surface was an important parameter which can impact the screech noise. Ibrahim et al. [
26] studied the effect of turbulence characteristics of jet flows on the radiated jet noise. Berland et al. [
27] used compressible large eddy simulations to investigate the generation of screech tones from an under-expanded jet with a three-dimensional planar geometry.
In the present work, we aim to extend these results by examining the effect of sinusoidal surface deformations added to the trailing edge of the aft-deck plate. To this end, a suite of large eddy simulations have been performed, each targeting different geometrical configurations and flow conditions. Sound levels in the farfield were then computed to determine whether these deformations could potentially reduce the total noise emitted by the installed configuration. In these simulations, we consider a high aspect ratio rectangular jet. To simplify the nozzle geometry, we used a spanwise slice of width four times the height of the nozzle with periodic boundary conditions applied in the spanwise direction. Flow disturbances were imposed at the inflow boundaries to introduce three-dimensional flow structures into the jet flow. A high-order accurate solver was used to discretize the unsteady, compressible, conservative form of the filtered Navier–Stokes equations. The mean flow was determined using Reynolds Averaged Navier–Stokes (RANS) equations via a SST
turbulence model (Menter [
28]). The farfield noise radiation was evaluated using Ffowcs-Williams and a Hawkings acoustic analogy method.
We provide further details of the numerical formulation in
Section 2, and present our results in
Section 3.
3. Results and Discussion
Results are reported and discussed in this section for a high aspect ratio rectangular jet exhausting over a flat surface located underneath the jet flow. Because the aspect ratio of the nozzle exit is high, we simplify the problem by considering only a portion from the nozzle along the spanwise direction. This simplification is based on the assumption that the flow is statistically two-dimensional (an assumption also considered by [
24]). Thus, the spanwise length of the flow domain is set to be four times the height of the nozzle, and periodic boundary conditions are imposed in the lateral direction.
The layout of the configuration is illustrated schematically in
Figure 1a. The flat plate is oriented parallel to the jet axis and located beneath the jet exit. The vertical separation of the plate from the jet axis (
d) and the distance from the jet exit to the trailing edge of the plate (
L) were both varied in this study. These distances are specified below and in
Table 1 in terms of multiples of the nozzle height (
). The vertical separation considered were 0.65, 0.75, and 1 times the nozzle height, while the horizontal distances to the trailing edge considered were 4, 5, and 6 times the nozzle height. A 3D view of the geometry of the nozzle and the plate with deformed trailing edge are shown in
Figure 1. The deformations in the spanwise direction are prescribed as a sine function while the amplitude is modulated by a hyperbolic tangent function in the streamwise direction. Two combinations of wavenumbers and amplitudes are shown in
Figure 1b,c.
The Reynolds number, based on the jet velocity and the height of the nozzle, is 500,000, and the acoustic Mach number is . The amplitude of velocity disturbances imposed at the nozzle exit are restricted to less than 3% of the jet velocity.
3.1. RANS Results
The mean flow and turbulent kinetic energy (TKE) used to initialize the LES are obtained from separate two-dimensional RANS simulations using a
turbulence model. The RANS results shown in
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6 are reported first to demonstrate the effect of the plate location on the jet mean flow and TKE. Plots of mean velocity magnitude (left) and turbulent kinetic energy (right) are shown in
Figure 2 for two selected cases, corresponding to the largest and the smallest offset distance of the plate from the jet axis. These color plots indicate that there is a visible effect from the plate on both the mean flow and the turbulent kinetic energy. For the first configuration (
,
), which corresponds to the largest distance
d from the center of the nozzle, the jet is significantly deviated downward by the plate (resembling the ’Coanda’ effect). This deviation is also present in other cases, but they do not appear to be as significant as it is for this case. The smallest deviation from the jet axis occurs for the configuration (
,
), which corresponds to the smallest offset distance,
d.
From the plots of TKE shown on the right side of
Figure 2, there appears to be an attenuation of the TKE in the bottom shear layer of the jet as a result of the interaction with the plate boundary layer. This is potentially due to the viscous dissipation in the boundary layer that is developing on the surface of the plate. Further downstream of the plate trailing edge, the magnitude of the TKE appears to recover somewhat as the two shear layers show similar intensities. This is further revealed by the profiles of TKE across the jet shown in
Figure 6.
Center-line mean velocity distributions are plotted in
Figure 3 for several tested configurations along with the results obtained from the baseline free jet (no plate). The presence of the plate seems to slow down the flow inside the potential core as seen in the zoomed-in plot shown in
Figure 3b. While all curves seem to follow the same trend in the downstream region, there is a slight acceleration for the cases that include the plate. This does seem to depend on the distance of the jet axis to the plate,
d; the smaller this distance, the higher the acceleration. TKE distributions along the jet center line that are plotted in
Figure 4 reveal a decrease of TKE in the downstream region for the cases that involve the plate. This decrease is more significant for the configurations with the smallest
d and smallest
L (yellow and red curves in
Figure 4).
In
Figure 5 and
Figure 6, we compare cross-flow vertical profiles of mean velocity magnitude and turbulent kinetic energy among the five configurations, at four axial locations,
;,
,
, and
; the results from the free jet case are also included. They all show that the largest deviation of the jet from the original position occurs for
and
, corresponding to the largest distance of the plate from the jet (this was also observed in the contour plots of the mean velocity and TKE).
Figure 6 shows that the reduction of the TKE in the lower shear layer is more significant for the configuration with the smallest distance from the trailing edge.
Figure 5d and
Figure 6d indicate that the jet flow becomes more symmetric further downstream for all normalized profiles, except the case corresponding to
,
.
3.2. LES Results
A suite of LES runs were carried out corresponding to the geometrical configurations given in
Table 1, and the trailing edge deformations characterized by the function
where
is the streamwise location on the plate where the deformation begins to gradually increase from zero to the maximum amplitude at the trailing edge location;
is a parameter controlling the ’thickness’ of the tangent hyperbolic function; and
is the spanwise width of the plate. The amplitude,
A, and wavenumber,
k, are varied as given in
Table 2. (See also
Figure 1 for a visual representation of the trailing edge deformations.) As mentioned previously, the spanwise length of the flow domain is four times the height of the nozzle, and periodic boundary conditions are imposed in the lateral direction with the assumption that the aspect ratio of the rectangular nozzle is very high. This lateral dimension of
was selected to ensure that all relevant turbulent flow structures are captured accurately in the spanwise direction. The mesh consists of approximately 10 millions grid points, clustered mostly in proximity to the jet region. Grid stretching is used in the farfield along with an imposed sponge layer condition at the outflow boundary to ensure that the fluctuations are gradually dissipated while they leave the flow domain.
Figure 7 shows iso-surfaces of Q-criterion colored by the streamwise velocity component. (
, where
is the rate-of-strain tensor, and
is the vorticity tensor.) Acoustic waves radiating from the jet are plotted in gray contours. (The white lines are caused by gaps between blocks from the MPI decomposition). The attenuation of acoustic waves propagating from the jet underneath the plate can be observed. This is a result of the shielding provided by the plate causing the waves to not be as intense in this region.
The mean centerline velocity distributions are compared against RANS results in
Figure 8. While there are some slight differences between the two in the potential core region, overall the agreement is good, showing similar decay in the downstream flow. (The LES profiles show some fluctuations from insufficient time spans used for the time averaging).
Acoustic spectra have been calculated at three probe locations relatively close to the jet (in the nearfield). These probes are located at large angles with respect to the jet axis since it was observed that jet surface interaction noise dominates at large polar angles (see [
39] for a study of AR8 nozzle and [
20] where various rectangular nozzle jets were considered). The Ffowcs-Williams Hawkins surface is shown in
Figure 9. The vertical side in the downstream was placed at 30 nozzle diameters from the nozzle exit plane. (Varying this distance did not affect the results significantly.) Two probes are located above the jet while the third is located underneath the jet (see
Figure 9). The coordinates of the probes are listed in
Table 3. A fourth probe was placed in the far field at 100 equivalent diameters and 50 deg angle with respect to the jet axis to compare the numerical results with experimental results collected by Bridges et al. [
40] (see also Bridges [
22]). Results were analyzed by using the sound pressure level and overall sound pressure level plots. The Strouhal number,
, was computed by the formula:
where
f is the frequency,
is nozzle height, and
is the jet velocity.
Figure 10 shows the comparison between our acoustic spectrum and the spectrum obtained from measurements (Bridges et al. [
40]). Bridges et al. [
40] conducted an extensive set of measurements to acquire acoustic data on single-flow convergent rectangular nozzles of aspect ratios 2:1, 4:1, and 8:1, at Mach numbers 0.7 and 0.9. Our comparison is for the largest aspect ratio 8:1 and Mach number 0.9. The agreement in
Figure 10 is fairly good in the low frequency range of the spectrum. For the high frequency range the spectrum curve obtained from numerical simulations is slightly lower, which is most likely due to the mesh resolution of the LES not being sufficient to fully resolve the smallest flow scales that are responsible for the highest frequency noise.
In general, the effect of trailing edge deformation on the radiated noise was small. We noticed both an increase and a decrease in noise, depending on the combination of amplitudes and wavenumbers characterizing the deformations. A visible reduction in noise was achieved for treatment 5 given in
Table 2 (
,
). The acoustic spectra at the three probe points are plotted for this case in
Figure 11a–c.
The overall sound pressure levels (OASPL) observed at the probe points are plotted in
Figure 12 and
Figure 13. In
Figure 12, the OASPL’s are plotted for the five configurations that do not involve trailing edge deformations. For the probes that are located above the jet, there is a significant increase in the radiated noise for the case corresponding to the smallest distance from the nozzle center (
,
). This is to be expected for this case since the end of the potential core is the closest to the trailing edge of the plate. For the probe located under the jet, there is an increase in the radiated noise for the case corresponding to the shortest plate length. This is an indication of a reduced shielding effect.
Figure 13 displays the OASPL results at all three probe locations for each of the tested configurations and prescribed trailing edge deformations. Each plot displays the results for all of the treatments given in
Table 2 at one of the fixed plate configurations given in
Table 1. When the data are broken down by each plate placement location, there were two apparent plate configurations that yielded a noticeable reduction in sound levels at the probe locations. We will look further into these configurations to determine if the trailing edge deformations provide an additional positive effect.
There appeared to be an overall reduction in noise for most of the treatments in case 1 (
). The average overall sound pressure level for the no deformation case was 121.367 dB. For the
treatment, the average was 120.467 dB, and, for the
case, the average was 120.4 dB. The results for the Overall Sound Pressure Level data can be seen in
Table 4. Although there was a reduction in noise for all of the deformations at this plate placement location, the two deformations mentioned above demonstrated the greatest reduction in noise.
For case 5 (
), there was also a reduction in noise for all of the deformations applied to the trailing edge of the plate. The average overall sound pressure level for the probes for the no deformation case was 122.067 dB. For the
treatment, the average was 120.633 dB, and for the
treatment the average was 121.133 dB. The Overall Sound Pressure Level results for this plate placement location can be seen in
Table 5.
The OASPL results from these configurations yielded the best results for overall noise reduction among those tested. It would be worthwhile to consider other deformations that are close to these treatments as well as other nearby plate placements to determine if further noise reduction could be achieved.
4. Conclusions
In this work, jet flow over various surface deformations at the trailing edge of a plate installed under a high aspect ratio rectangular jet were simulated to determine the effect of the jet-surface interaction on the noise radiated to the farfield. RANS calculations and a suite of LES simulations were performed with a high-order accurate flow solver. Numerous configurations and flow conditions were investigated, and the emitted noise from the jet-surface interactions calculated using Ffowcs-Williams and the Hawkings acoustic analogy method.
The RANS results, consisting of centerline velocity distributions, showed that the presence of the plate appeared to slow down the flow inside the potential core. For some of the tested cases, a slight acceleration of the flow further in the downstream was observed. The amount of acceleration appears to depend on the offset distance of the plate from the center of the jet, d. TKE distributions along the jet centerline revealed a decrease in TKE magnitude in the downstream region, which was expected. This decrease was more significant for the configuration with the smallest d and smallest L. Cross-flow profiles of mean velocity magnitude and turbulent kinetic energy at different axial locations showed a clear deviation of the jet from the axis of symmetry near the trailing edge. The TKE profiles indicated that the jet is slightly skewed around the trailing edge of the plate in a manner resembling the Coanda effect. Despite this deviation, the cross-flow profile becomes more symmetric further downstream. A reduction of the TKE in the lower shear was observed for all cases with the plate under the jet. This is most significant for the configurations with the smallest distance to the trailing edge.
For the LES analysis, various results consisting of iso-surface of Q-criterion, centerline velocity, acoustic spectra, and overall sound pressure levels were presented and discussed. Overall, good agreement was found between the centerline velocities from RANS and LES, with some small discrepancies in the potential core region. Acoustic spectra plotted for the configurations without trailing edge deformations confirmed a reduction of the noise under the plate from the shielding effect. For most cases involving the trailing edge deformations, there did not appear to be a significant change to overall sound levels. However, some cases did demonstrate a noticeable increase or decrease in emitted noise. In the best cases, the overall sound pressure levels were decreased from 1 to 5 dB.
The results presented in this work provide evidence that further refinement of the configurations described herein may yield additional noise reductions. Future research aimed at analyzing the interactions of the jet flow with different types of surface deformations (e.g., non-periodic and/or streamwise-oriented) may uncover a greater degree of noise reduction, beyond what our results have shown.