Numerical Investigation of a Rectangular Jet Exhausting over a Flat Plate with Periodic Surface Deformations at the Trailing Edge

In an attempt to reduce jet noise, blended wing concept vehicle utilizes rectangular jet exhaust ports exiting from above the wing ahead of the trailing edge. In this study, we take another look at the rectangular exhaust port configuration with some notional modifications to the geometry of the trailing edge to determine if the emitted noise levels due to jet interactions can be reduced with respect to a baseline configuration. We consider various horizontal and vertical offsets of the jet exit with respect to a flat plate standing in for the aft wing surface. We then introduce a series of sinusoidal deformations to the trailing edge of the plate of varying amplitude and wave number. Our results show that the emitted sound levels due to the jet--surface interactions can be significantly altered by the proposed geometry modifications. While sound levels remained fairly consistent over many configurations, there were some that showed both increased and decreased sound levels in specific directions. We present results here for the simulated configurations which showed the greatest decrease in overall sound levels with respect to the baseline. These results provide strong indications that such geometry modifications can potentially be tailored to optimize for further reductions in sound levels.


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 largescale 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 threedimensional 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 k − ω 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.

Scalings
The governing equations considered here are expressed in terms of a generalized curvilinear coordinate system, subject to the following transformations: where ξ, η, and ζ are the body-aligned curvilinear coordinates corresponding to the streamwise, wall-normal, and spanwise directions, respectively, and x, y, and z are the (non-dimensional) global Cartesian coordinates of physical space. All dimensional spatial coordinates (x * , y * , z * ) are normalized by the reference length, D j , the height of the nozzle (See Figure 1) The fluid velocity values are scaled by the jet velocity, V j .
Pressure values are scaled by the dynamic pressure, ρ ∞ V 2 j , and temperature values by the freestream temperature, T ∞ . Reynolds number, Mach number, and Prandtl number are defined as: where µ ∞ , a ∞ , and k ∞ are the freestream dynamic viscosity, speed of sound, and thermal conductivity, respectively. C p is the specific heat at constant pressure. For all simulations, these values are initialized for air as an ideal gas.

Governing Equations
The filtered Navier-Stokes equations can be expressed in conservative form as: Here, the subscripts denote partial derivatives with respect to time and the curvilinear coordinates. Q is the vector of conserved quantities: ρ = ρ * /ρ ∞ is the non-dimensional density of the fluid, u = (u, v, w) is the non-dimensional velocity vector in physical space, and E is the total energy. The flux vectors, F, G and H, are given by: Here, the repeated indices indicate Einstein summation convention over i ∈ 1, 2, 3. U, V, W are the contravariant velocity components, τ is the shear stress tensor, and Θ is the heat flux, Again, repeated indices (that are not present on the left-hand side) indicate summation. The Jacobian of the curvilinear transformation from the physical space to computational space is denoted by J, and S is a vector of prescribed source terms. Pressure, temperature, and density are related by the ideal gas equation of state: The dynamic viscosity and thermal conductivity k are related to temperature using the Sutherland's equations in dimensionless form: For air at sea level, C 1 = 110.4 K, and C 2 = 194 K.
There are no explicit subgrid scale turbulence terms in Equation (4) (see Grinstein et al. [29]). Instead, the compressible Navier-Stokes equations are solved within the framework of an implicit large eddy simulation, where numerical filtering is applied to account for the missing sub-grid scale energy. The numerical solver uses high-order finite difference approximations for the spatial derivatives and explicit time marching. The time integration is performed using a second order Adams-Bashforth method [30]: Here, the constants β ν are chosen to give either the maximum order of accuracy [30] or the lowest dispersion and dissipation. L(Q) is the residual.
The spatial derivatives are discretized using the dispersion-relation-preserving schemes of Tam and Webb [31] or a high-resolution 9-point dispersion-relation-preserving optimized scheme of Bogey and Bailly [32]. To damp out the unwanted high wavenumber waves from the solution, high-order spatial filters, as developed by Kennedy and Carpenter [33], are used. No slip boundary conditions for velocity and adiabatic conditions for temperature are imposed at the solid surfaces. Sponge layers are imposed near the far-field boundaries in regions that are outside the flow domain of interest. The sponge layers combined with grid stretching act to damp-out unwanted waves returning from the farfield boundaries. For the LES simulations, the flow in proximity to the wall is modeled using the Werner-Wengler model (Werner and Wengle [34]), so some small turbulent flow structures may not captured adequately. At the inflow boundary, a constant mean flow in the core region is combined with a hyperbolic tangent function in the shear region, and disturbances in the form of a superposition of Fourier modes with random amplitudes, frequencies, and wavenumbers. These disturbances introduce three-dimensional flow variations into the jet flow.
The mean flows used to initialize the LES are obtained from RANS simulations, where a classical SST k −ω turbulence model (Menter [28]) is applied to account for missing fluctuations. This model was found to perform well for jet applications (see Mihaescu et al. [35]) In the present work, high-order, central-difference schemes are used to achieve increased resolution of the propagating disturbances. A shock-capturing technique suitable for simulations involving central differences in space is required to avoid unwanted oscillations that may propagate from discontinuities arising in supersonic flows. Shock capturing techniques are employed based on the general explicit filtering framework, a straightforward approach which introduces sufficient numerical viscosity in the area of the discontinuities, and negligible artificial viscosity in the rest of the domain. The technique of Bogey et al. [36] introduces selective filtering at each grid vertex to minimize numerical oscillations, and shock-capturing in the areas where discontinuities are present. This method has been proven to work efficiently for high-order accurate, nonlinear computations.

Ffowcs-Williams Hawkings Acoustic Analogy Method
The Ffowcs-Williams and Hawkings (FW-H) equation [37,38] is an inhomogeneous wave equation that can be derived by manipulating the continuity equation and the Navier-Stokes equations for a compressible fluid. The FW-H equation can be written as: where u i is the fluid velocity component in the x i direction, u n is the fluid velocity component normal to the surface at f = 0, v i is the surface velocity components in the x i direction, v n is the surface velocity component normal to the surface, δ(f ) is Dirac delta function, H(f )is Heaviside function, and p = p − p 0 is the magnitude of the sound pressure disturbance at the far field. The acoustic pressure in the farfield at a location, x, can be determined from: where Given a point, y, in the near-field acoustic source on surface S, r = | r| = | x − y| is the distance from the surface source location to the farfield position x; U n is the scalar product of U = ρ u/ρ 0 and the unit normal to the surface S at y; L r is the scalar product of L i = P ij n i + ρu i u n and the unit vector in the direction from y to x; P ij = −τ ij + pδ ij ; τ ij is the viscous stress; ρ is the instantaneous density; and ρ 0 and a 0 are the ambient density and sound speed, respectively.
In Equation (15), p Q ( x, t) is the contribution from the quadrupole, but it is neglected here. The integrals in (16) and (17) are calculated at the emission time τ = t − r/a, where a is the speed of sound and t is the reception time. The flow variables at τ must be interpolated from the flow data at multiple time iterations, t, due to its dependency on y, which varies over the surface S surrounding the jet.
For the numerical implementation of the FW-H method, a conical surface that includes the entire jet is considered. The primitive variables and viscous stresses are interpolated from the flow domain to the surface at every b time steps (the value of b depends on the desired time resolution) and written out to a large number of files for post-processing. Given a reception time t, the emission time τ is calculated for every small surface element ∆S on the conical surface. Within a specific surface element ∆S on the surface, an interpolation is necessary to obtain the flow data at the emission time τ , which is needed to calculate the integrals in (16) and (17). Next, we determine the time derivatives of U n and L r and evaluate the integrals in Equations (16) and (17). The acoustic pressure then is calculated from Equation (15).

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 (D j ). 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 0.8. The amplitude of velocity disturbances imposed at the nozzle exit are restricted to less than 3% of the jet velocity. Table 1: Position of the plate with respect to the nozzle (see also Figure 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 Figures 2-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 (d = 1.0D j , L = 5D j ), 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 (d = 0.65D j , L = 5D j ), 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 Figures 5 and 6, we compare cross-flow vertical profiles of mean velocity magnitude and turbulent kinetic energy among the five configurations, at four axial locations, x = 4D j ;, x = 10D j , x = 25D j , and x = 45D j ; 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 d = 1.0D j and L = 5D j , 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. Figures 5d and 6d indicate that the jet flow becomes more symmetric further downstream for all normalized profiles, except the case corresponding to d = 1.0D j , L = 5D j . (c) (d) Figure 6: Turbulent kinetic energy profiles along vertical direction, for different configurations:

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 x 0 = −1.2 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; σ x = 2 is a parameter controlling the 'thickness' of the tangent hyperbolic function; and L z 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 4D j 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.
is the rate-of-strain tensor, and Ω = 1/2[∇v − (∇v) T ] 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). 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 [41]). Results were analyzed by using the sound pressure level and overall sound pressure level plots. The Strouhal number, St e , was computed by the formula: where f is the frequency, D j is nozzle height, and V j is the jet velocity. Figure 9: Locations of the three Probes, and the Ffowcs-Williams integration surface. 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 (w z = 4, A = 0.14). 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 Figures 12  and 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 (d = 0.65D j , L = 5D j ). 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 (d = 0.75D j L = 4D j ). The average overall sound pressure level for the no deformation case was 121.367 dB. For the w z = 10 A = 0.04 treatment, the average was 120.467 dB, and, for the w z = 6 A = 0.08 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 (d = 0.65D j L = 5D j ), 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 w z = 10 A = 0.04 treatment, the average was 120.633 dB, and for the w z = 4 A = 0.14 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.

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.