Next Article in Journal
A Multi-Node Lumped Parameter Model Including Gravity and Real Gas Effects for Steady and Transient Analysis of Heat Pipes
Next Article in Special Issue
Evaluation of Turbulent Jet Characteristic Scales Using Joint Statistical Moments and an Adaptive Time-Frequency Analysis
Previous Article in Journal
The Functional Relationship of Sediment Transport under Various Simulated Rainfall Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Wavelet-Based Time-Frequency Analysis on the Supersonic Jet Noise Features with Chevrons

by
Stefano Meloni
1,* and
Hasan Kamliya Jawahar
2
1
Department of Engineering, University of Roma Tre, 00146 Rome, Italy
2
Department of Aerospace Engineering, University of Bristol, Bristol BS8 1TR, UK
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(3), 108; https://doi.org/10.3390/fluids7030108
Submission received: 11 February 2022 / Revised: 4 March 2022 / Accepted: 11 March 2022 / Published: 16 March 2022

Abstract

:
A detailed investigation of the statistical properties of the near-field pressure fluctuations induced by an under-expanded jet, by varying the nozzle exit shapes has been presented. Experiments using different convergent Chevron nozzles were carried out in the anechoic chamber at the University of Bristol to assess the importance of the Chevron shape on the near pressure field emitted by a single stream under-expanded jet. Measurements were carried out through an axial microphone array traversed radially to various positions for jet in an under-expanded condition at Mach number M = 1.3. The intermittent behavior is investigated considering the standard statistical indicators and a wavelet-based conditional approach, including the phase angle. The intermittent degree of various features related to different scales, such as Screech tones and broadband shock associate noise were estimated. A series of recently developed wavelet-based tools were assessed as a viable approach to investigate the noise emitted by under-expanded jets.

1. Introduction

One of the key aircraft noise sources that should be taken into account during the design of modern aircraft is the jet noise which dominates the take-off and the cruise phase [1]. The generation of sound by a jet exiting into a fluid medium is of great interest for several aeronautical applications and it has long been recognized by several previous studies [2,3]. A large body of literature has clarified that the dominant noise source in the subsonic case is the turbulent mixing nose, which is generated by the large-scale turbulence structures/instability waves of the jet flow [4,5]. During takeoff, when a high level of thrust is required, jet exiting flow from the engine’s nozzle is under-expanded. This could also happen during the cruise phase due to the combined effect of low external static pressure and required thrust to maintain the flight Mach number in maneuvers or changing the lane [6,7]. An under-expanded jet plume is characterized by a shock cell train, which creates a series of compression and expansion into the flow, generating a further shock-associated noise [7,8,9]. This form of noise generated by non-ideally expanded supersonic jet comprises of two components: Screech tones and Broadband Shock Associated Noise (BBSAN). Screech tones are discrete tones that originate from an acoustic feedback loop between the shock cell train and the nozzle lip [10]. Screech is unusual among resonance phenomena, in that the resonance is entirely contained within the flow itself [11].
BBSAN, on the other hand, is generated from a weak interaction between downstream propagating large scale turbulent structures and the quasi-periodic shock cells in the jet plume. Contrary to Screech, BBSAN is present in both under and over-expanded jets. One of the characteristics of BBSAN (unlike Screech) is that the peak frequency varies as a function of observer position, a phenomenon observed in experimental data. The inclusion of these noise components in the analysis is essential to solve the problem and thus provide tolerable noise levels in the cabin to improve passenger comfort.
Aircraft manufacturers developed various technologies to reduce jet noise, and one of the most common jet technologies is the use of Chevrons at the nozzle exit [12,13,14,15]. They reduce the noise by reducing the velocity gradients in the jet shear layer and in the supersonic case they disrupt the feedback loop responsible of the Screech tone [16]. Bridges and Brown [17] showed that the increase in Chevron number and penetration achieves considerable noise reduction at lower frequencies for convergent jet.
The velocity gradient variation in the jet shear layer, induced by the Chevron presence, has been expected to modify the stochastic nature of the noise source in terms of intermitency of the aforementioned noise components. The role of intermittent events in the near-field of a jet is fundamental to evaluate the noise emitted in far-field and its role in the noise generation mechanism has been recognized in the subsonic case by several works [4,5,18,19].
In this paper, for the first time statistical analysis of the intermittency events related to the under-expanded acoustic signatures using a wavelet-based approach has been carried out. The wavelet technique is an efficient tool for extracting and analyzing the time evolution of the frequency-localized acoustic signatures. Different wavelet-based indicators, directly related to intermittency, like the local intermittency measure (LIM) and its square (LIM2), which is correlated to the flatness factor [4], was used to analyze the modification of the time evolution of the acoustic signatures detected in the baseline configurations. The time-frequency evolution of the phase-angle has been included in the analysis using a multivariate wavelet-based approach, being its modifications related to the modification of the noise directivity
To carry out this investigation, we used a database of near-field pressure signals acquired from the Bristol Jet Aeroacoustic Research Facility (BJARF) at the University of Bristol. A total of three different nozzles having different shapes of Chevrons were considered, one baseline and two chevron nozzles. The choice of nozzles were based on the availability of flow field, which is well characterized in the literature by several works (see [20,21]). According, to previous works [16,17] we consider the best nozzle configuration for noise SMC006, which has a slight penalty on the effective nozzle exhaust diameter directly correlated to the thrust and SMC002, which provided a lower noise reduction incrementing the effective nozzle exhaust diameter. The analyses were performed at M = 1.3 using a near-field microphone array positioned at h/D = 3.5 with 10 microphones covering a distance starting at the nozzle exit x/D = 0 up to x/D = 18 (D = 16.8 mm). The array position was chosen to avoid high hydrodynamic effects due to the flow grazing over the microphones. The study mainly focused on the zone close to the nozzle exit, as it is mostly dominated by Screech and BBSAN.
The paper is organized as follows: the experimental setup is reported in Section 2 and the results are illustrated in Section 3. Conclusions are presented in Section 4.

2. Experimental Setup

The experiments were conducted at the BJARF at the University of Bristol. The flow in the BJARF facility is conditioned and silenced using three different in-line silencers to create a quiet flow. Two silencers were placed right after the control valve outside the anechoic chamber and has a diameter of 0.3 m and a height of 1.5 m each. The third large silencer, which also acts as a plenum, was placed inside the anechoic chamber and has a diameter of 0.457 m and a height of 1.9 m. The silencers were equipped with perforated tubes for the flow with the remaining area packed with glass wool. The anechoic chamber where the tests were carried out has dimensions of 7.9 m in length, 5.0 m in width and 4.6 m in height, including the surrounding acoustic walls [16].
The tests were carried out using a round convergent nozzle (SMC 000) and two different Chevron nozzles (SMC 002 and SMC 006) with a different number of lobes, 4 and 6, respectively, as shown in Figure 1. Details about the nozzle geometries have been reported in Table 1. The Chevron nozzles were chosen from the detailed study carried out by Bridges and Brown [17].
The tested nozzles were 3:1 down-scaled version of the nozzles used by Bridges and Brown [17], which corresponds to an exit diameter of D = 16.933 mm for the SMC 000 round convergent nozzle. The tests were carried out for a wide range of supersonic flows ranging from M = 1.1 up to 1.3 . Near-field unsteady pressure measurements were carried out using 1 / 4 -inch G.R.A.S 40PL microphone that has a corrected flat frequency response at frequencies from 100 Hz to 20 kHz with a dynamic range of 150 dB. The near-filed measurements were carried with a linear array of 10 microphones placed at a distance of 2 D away from each other, moved to 12 different heights ( h ) radially away from the jet centerline using a traverse. The closest height was chosen to be h = 3.5 D to avoid grazing of flow on the microphone. The data were acquired using a National Instrument PXle-4499 for t = 12 s at a sampling frequency of f = 2 17 Hz. To clearly show the experiment, a photo of the experimental setup is reported in Figure 2.

3. Results and Discussions

3.1. Near-Field Spectral Levels

The results for the Sound Pressure Level (SPL) for several streamwise locations ranging from x/D = 0–18 at radial location h/D = 3.5 for M = 1.3 are presented in Figure 3. The SPL was evaluated using the following equation:
S P L = 10 log 10 PSD Δ f ref P ref 2 ,
where PSD denotes the power spectral density evaluated using Welch’s method, Δ f ref is the frequency bandwidth and P ref is the reference pressure in air (equal to 20 μ Pa).
The results for the baseline configuration SMC000 clearly show the presence of Screech tone, after which the BBSAN spectral hump is accurately captured in Figure 3a. The results for the Chevron configurations in Figure 3b,c do not show the presence of Screech, however the BBSAN hump is present. According to the literature, [16,17] and as reported in Figure 3 the Chevrons effectively eliminate the Screech tone, which is known to be generated due to the feedback loop between the shock cells and the nozzle lip. The bump related to BBSAN results shifted to the higher frequencies for the Chevron configuration, which might be most likely due to the variation of the shear layer velocity gradient. For both the Baseline and Chevron configurations, the locations closer to the nozzle, the spectra show a flatter trend, however, at downstream locations, the spectra show a low-frequency spectral hump due to the large scale structures that are often found at the mixing region. Overall, the results show the capability of Chevron to eliminate the feedback loop, thus the Screech tone. Moreover, the results show that the number of Chevron lobes also plays a significant role in the noise reduction mechanism as the SMC 006 with six lobes clearly shows improved noise reduction compared to SMC 002 with just four lobes.

3.2. Global Intermittency Analysis

One of the main aims of the present work is to provide a qualitative description of the time evolution of the near-field pressure signatures modified by the presence of the Chevrons. To have a global picture of intermittency, the third and four order statistical moments, named skewness (s) and kurtosis (k), were taken into account, which are defined in the following equations:
s = E [ p μ ] 3 σ p 3 ,
k = E [ p μ ] 4 σ p 4 ,
where μ is the mean of the signal p and E[ ] is the expected value. As shown in Figure 4a, the evolution of the skewness factor in the stream-wise direction varies consistently between the baseline and Chevron configurations. Both negative and positive skewness can be observed for the baseline nozzle in the jet potential core region, dominated by positive pressure events within the potential core region. Predominantly, the presence of Chevrons moves the skewness to zero in the jet potential core region, and this is ascribed to the disruption of the Screech feedback loop. At the region downstream of the potential core, the near-field pressure signals have negative skewness, which is often observed in the ideally expanded case (subsonic in our case) where the Screech tone does not appear [22,23]. Negative values of the third-order statistical moments were identified in the fully turbulent jet zone at the downstream locations, which can be attributed to the development of a fully turbulent jet flow. The amplitude of these negative values increases with the Chevron presence and seems amplified by the higher number of lobes. This could be related to the increase in the velocity gradient in the jet shear layer. However, further studies about the turbulent characteristics of the nozzle are needed to clarify this assumption.
The kurtosis trend has been reported in Figure 4b for different axial locations and different nozzle configurations. The kurtosis values were identified to be lower than three solely in the baseline configuration (SMC000) for axial positions of the microphones within the jet potential core. This could be related to the Screech presence, which is a probably persistent tonal component, further details on this will be provided with the wavelet analysis. It is important to note that the kurtosis close to 3 and the skewness close to zero indicate that the presence of Chevrons generates a Gaussian distribution of the fluctuating pressure events, suggesting a low presence of the hydrodynamic contribution.
At locations farther downstream of the nozzle exit higher kurtosis is observed, especially in the fully turbulent jet zone increasing with the number of lobes. These kurtosis values suggest that Chevron nozzles disrupt the feedback loop related to the Screech in the first few axial locations increasing the jet flow development.
To further improve our understanding on the nature of the time signature, an analysis of the signal stochastic behavior has been carried out on a few representative cases by using the probability density functions (PDF) reported in Figure 5a,b, considering all the presented nozzle exit configuration and two different locations of the microphones in the axial direction. The pressure variable is expressed in reduced form, i.e., normalized to have zero mean value and unitary standard deviation [22].
A quasi-bimodal PDF has been observed in Figure 5a for the SMC000 configuration, which is due to the presence of the Screech tone. As expected, PDF shape is modified by the use of Chevron nozzles. These results also follow a Gaussian form, but with slightly larger PDF tails, which could be related to the disappearance of the Screech tone and to the quite flow developing due to the higher velocity gradient in the jet shear layer. To confirm this, at downstream location, Figure 5b, we observed a large number of pressure events in the PDF tails, especially for the SMC006 configuration, which is characterized by a higher number of Chevron lobes.

3.3. Single-Point Wavelet Analysis

The wavelet transform is a very proficient tool when it comes to analyzing intermittent or time-dependent features. The wavelet transform of the signal p(t) is obtained by the following expression [24,25]:
w ( s , t ) = C ψ 1 2 s 1 2 p ( τ ) ψ * ψ * t τ s ,
where s is the wavelet scale, τ is a time shift, C ψ 1 2 is a constant that takes into account the mean value of ψ ( t ) and ψ * t τ s is the complex conjugate of the dilated and translated mother wavelet ψ ( t ) . In this analysis, we have applied the continuous wavelet transform (CWT) using the Morlet mother wavelet. The wavelet scalograms for two streamwise locations for all the three tested cases are presented in Figure 6. The first location was chosen to be within the jet potential core at x/D = 2 (see left column in Figure 6) and the second location was chosen in the turbulence mixing region at x/D = 18 (see the right column in Figure 6).
The wavelet scalogram results for the baseline configuration shows a time-persistent high energetic signature for a range of frequencies between 7 kHz and 10 kHz at x/D = 2 in Figure 6a. According to the spectral analysis seen in the previous sections, this feature is related to the Screech tone. At a higher frequency range, the same microphone location within the potential core is dominated by intermittent signatures related to the BSSAN. Moving to the characteristics of the downstream microphone at x/D = 18, Figure 6b, the intensity of the Screech signature has reduced, and energetic events intermittent in time appear at the lower frequencies because of the jet development and larger flow structures at the downstream location. It is interesting to note that although the Screech tone was observed in the spectral levels in Figure 3a at downstream locations, wavelet scalogram results show a substantial change in its intermittent characteristics compared to the upstream location. The scalogram results for the Chevron configurations, at location x/D = 2 in Figure 6c,e, show the absence of Screech signature and seems to increase the energy of the intermittent features associated with the BBSAN. According to the spectra presented in Figure 3b,c the presence of Chevrons seems to disrupt the high energetic time signatures detected in Figure 6b at higher axial locations where the jet is developed (i.e., x/D = 18) (see Figure 6d,f). This could be the cause of the reduction of the jet mixing noise and seems amplified by the increase of the Chevron lobes.
To further understand the nature of the Screech tone and the multiple tone generation mechanism seen in the current study, the power spectral density of the time signal and the absolute values of the wavelet coefficients were calculated. The time evolution of the Screech tone frequencies at the two different axial locations were selected using the results presented in Figure 3a, considering the Screech tone and its first harmonic for the baseline configuration within the potential core region and downstream location, respectively, x/D = 2 and x/D = 18.
The results presented in Figure 7 essentially show the modulation frequency of the selected tone. The results for the first fundamental tone within the potential core shows strong frequency modulation at f = 7.7 kHz for x/D = 2. However, at x/D = 18, the spectral energy of the modulation is absent.
In order to present a better comparison of the fluctuating characteristics of the unsteady near-field pressure components, the statistical representation of the wavelet coefficient moduli is presented in terms of the arithmetic mean and standard deviation, according to [26] have been evaluated as follows:
μ ( | w x | ) = Σ i = 1 N | w x | N ,
σ ( | w x | ) = Σ i = 1 N ( | w x | i μ ( | w x | ) ) 2 N ,
The mean μ and standard deviation σ of the wavelet coefficient moduli over time for each frequency were calculated using Equations (5) and (6). The results for the same for the three nozzle configurations for axial location x/D = 2 and 5 is shown in Figure 8.
At first glance, Figure 8a at location x/D = 2 shows a peak in the Screech related frequency with the nozzle in the baseline configuration, which confirms the high levels of μ ( | w x | ) magnitude of the Screech tone that is persistent over time. A similar trend is followed for the Screech tone for the σ ( | w x | ) results. When considering the chevron configuration, the broadband humps possess a higher magnitude compared to the baseline, with SMC006 showing high levels of magnitude at high frequencies. However, in the case of σ ( | w x | ) for the chevron configurations, the results for the SMC002 show high levels of dispersion of pressure from the mean compared to the SMC006 and SMC000. At the downstream location, x/D = 18, the μ ( | w x | ) shows no signs of the Screech tone for all the three configurations; however, high levels of fluctuations could be observed at low and mid-frequency range for the baseline configuration compared to the chevrons. This could be attributed to the high levels of turbulence at the downstream location. A similar trend is followed for σ ( | w x | ) , but with higher levels of dispersion of the fluctuations for the SMC002 compared to the SMC006. Overall, the results show low levels of fluctuation intensity and dispersion from the mean for the chevron configurations compared to the baseline.
To remove the dependence on the local feature energy, the so-called Local Intermittency Measure (LIM) [27] that represents a normalized version of the wavelet scalogram was used. Its formal definition is the following:
L I M ( s , t ) = w 2 ( s , t ) < w 2 ( s , t ) > t ,
where w 2 ( s , t ) are the wavelet coefficients evaluated with Equation (4) whilst the symbol < > t indicates time average. Intermittent features were identified by a LIM higher than one, while persistent features have a LIM lower or equal to this threshold. A series of LIM contour maps are presented in Figure 9, for both the baseline and Chevron configurations reported in the previous sections. At the near-field axial location x/D = 2 presented in Figure 9a a wide zone dominated by Screech having LIM = 1 in the baseline configuration can be observed. This zone spans from f = 7 kHz up to f = 10 kHz. As expected, this zone disappears in Figure 9c,e because Chevrons disrupt the feedback loop. Interestingly, the use of Chevron nozzles also seems to increase the degree of intermittency at the lower frequencies, by increasing the number of events with LIM higher than 1. At downstream location x/D = 18 LIM values results in higher levels due to the jet development, moreover, the presence of Screech is undetectable for the baseline configuration Figure 9b. In the case of the Chevron nozzles increased LIM value at low frequencies can be observed, this could be closely connected to a rise in the degree of intermittency and probably related to a more significant jet development, Figure 9d,f.
Following previous studies [4] highly intermittent events could be identified using the square of the LIM, named LIM2 as reported in the following equation, which is the definition suggested by [28].
L I M 2 ( s , t ) = w 4 ( s , t ) < w 4 ( s , t ) > t ,
LIM2 represents a convenient tool to extract those features contributing to the deviation from Gaussianity of the wavelet coefficients and can be interpreted as a time scale-dependent measure of the flatness factor of the signal. Indeed, LIM2 is equal to 3 if the PDF is Gaussian, and consequently, the condition LIM2 > 3 identifies only those events contributing to the deviation of the PDF from the Gaussian distribution of the wavelet coefficients. Otherwise, if LIM2 is lower than 3, it identifies only those events that induce a bi-modal distribution of the wavelet coefficients. This is an efficient tool to highlight the statistical content of a signal when intermittency or bi-modal contributions are confined at specific frequencies.
In the presented case, the Screech feature is characterized in the baseline configuration (see Figure 10a) by LIM2, which varies between 1 and 2 with no event above 3. This showed as expected that the Screech signatures has a bi-modal distribution in time, which is generated by the resonant phenomenon of the Screech. However, Chevron nozzles create a series of intermittent events in this frequency region, contributing to the deviation from a Gaussian distribution, increasing the global kurtosis of the signal (see Figure 10c,e). Considering the other frequencies, various effects can be observed: the higher frequency region seems slightly influenced by the Chevrons, while frequency at around 10 3 increased the number of events with larger LIM2 for the Chevron configuration. On the other hand, at x/D = 18, an increase of the LIM2 in all the three nozzle configurations (see Figure 10b,d,f) were observed in the higher frequency region while in the lower frequencies, LIM2 values decreased, with the larger part of it close to zero. These results are very evident in the baseline configuration. In contrast, some higher LIM2 values are detected in both Chevron configurations, which are probably responsible for the higher kurtosis observed in this location (see Figure 2). According to [4], it is important to underline that this behavior was entirely missed by the classical global statistical indicator reported in the first part of the present work.

3.4. Bi-Variate Wavelet Analysis

The multivariate wavelet analysis has been performed in the frequency domain using the wavelet coherence evaluated as follows:
γ ( 2 f , t ) = ψ | C x * ( f , t ) C y ( f , t ) | 2 ψ | C x ( f , t ) | 2 · ψ | C y ( f , t ) | 2
where | C x * ( f , t ) C y ( f , t ) | is the wavelet cross-spectrum of two consecutive signals in the axial direction, while C x ( f , t ) and C y ( f , t ) denote the continuous wavelet transforms of x and y (the two signals) at frequency f and time position t and ψ indicates the smoothing factor. Without the smooth function ψ , γ ( f , t ) will be equal to one everywhere. Hence, the wavelet coherence is a normalized scalogram [29,30], which describes the common power of two signals. The advantage in using a cross-wavelet transform with respect to computing a direct coherence is in the locality of the wavelet transform and the different resolutions achievable at the different scales. A series of wavelet coherence contour maps are reported in Figure 11 they are evaluated using signals from consecutive microphones located at different axial locations from the nozzle exit.
At first look high levels of coherence can be observed for the baseline configuration, especially at the Screech frequency in Figure 11a when considering two consecutive microphones that are located close to the nozzle exit at x/D = 2–4. The resulting coherence signature is observed to be persistent in time. As expected, when the nozzle configuration is changed to Chevron, the Screech characteristics fade and the coherence value for the above mentioned signature results reduced Figure 11c. These effects remain very similar when the Chevron penetration angle is varied, as shown in Figure 11e. For higher frequencies (>10 kHz) related to the BBSAN, the time step between two consecutive high coherence values results increased with the use of Chevrons. The variation of the Chevrons lobes modifies the low-frequency coherence, reducing it, especially at the mid-frequency range (see Figure 11e,f). This could be ascribed to the variation of the jet velocity gradient, but this needs to be clarified with further investigations with aerodynamic measurements. The analysis has also been repeated for two axial microphones located in the well-developed jet region (x/D = 16–18). There, Screech signature results are less evident with lower coherence levels (see Figure 11b). High coherence levels in low frequency can be observed in this region due to the large scale structures. For the baseline nozzle, high coherence is found at a frequency as low as 100 Hz. Interestingly, for Chevron nozzles, the coherence levels increase to higher frequencies compared to the baseline configuration. At the downstream location, the effect of the Chevrons results was indistinct on the coherence due to the well-developed jet flow, Figure 11d,f.
The following analysis has been carried out to understand the phase relationship between the two microphones at the frequency related to Screech and BBSAN for the tested configurations. The complex argument can be interpreted as the local relative phase between the two considered pressure signals, in this way, the wavelet cross-spectrum allows the evaluation of the phase angle between two consecutive signals in both time and frequency domains.
The phase angle θ has been computed using the following equation:
θ ( f , t ) = a r c t a n i m a g ( w c s ( f , t ) ) r e a l ( w c s ( f , t ) )
where w c s ( f , t ) are the wavelet cross spectrum coefficients.
In contrast to an average phase shift provided by the Fourier analysis, the wavelet phase provides a measure of phase shift in the time domain for each frequency localized feature. Figure 12 shows the phase angle, normalized by π , between two consecutive microphones in both the time and frequency domain. A phase angle close to − π has been detected in the frequency zone dominated by the Screech in the baseline configuration, Figure 12a. This could be due to a persistent upstream event related to the Screech presence with no phase dependency in time. The presence of Chevrons portrays intermittent phase angles as seen in Figure 12c,e. At higher axial locations in the downstream region, the phase angles are characterized by a positive signature at the higher frequencies and a negative signature at the lower frequencies Figure 12b. As for the coherence magnitude at this axial location, the phase angle seems only slightly influenced by the Chevron presence Figure 12d,f. On the other hand, the increasing of the Chevron numbers seems to increase the dominance of the negative phase angle at the lower frequencies, especially close to the nozzle exit. The physical understanding of these results is not trivial and needs further analysis. Regardless of that, in this paper there is the presence of highly intermittent events in both magnitude and phase that should be accounted for in the model to predict both near-field and far-field has been highlighted.

4. Conclusions

An in-depth time-frequency analysis using wavelet transform was carried out to deduce and report the stochastic behavior of the jet flow features in the under-expanded supersonic regime with and without Chevrons. Experiments were conducted at the Bristol Jet Aeroacoustic Research Facility at the University of Bristol using a microphone array to attain pressure fluctuations at the vicinity of the jet. Three different nozzle configurations were considered, a baseline case which a classical circular nozzle and two different Chevron nozzles with a different number of Chevron lobes. Data were acquired with the convergent nozzles in under-expanded conditions at M = 1.3.
The conventional spectral level analysis shows the presence of the characteristic Screech tone and its presence in all the presented axial locations. The use of Chevron nozzles eliminated the Screech tone due to the disruption of the feedback loops that drive the tone. The BBSAN was also captured well for the baseline configuration, and it was found to be still persistent in the Chevron configuration. The global statistical analysis highlighted that the Chevron presence modified the statistical distribution of the pressure events close to the nozzle exit by varying the PDF from a bi-modal trend to a Gaussian distribution. Additionally, larger PDF tails were observed at downstream location x/D = 18, ascribed to the increase of the flow development due to the Chevron lobes.
Consecutive analyses were focused on the time evolution of the features responsible for these statistical changes by using the continuous wavelet transform. The resulting wavelet scalogram was dominated by increased wavelet coefficients at the Screech frequency close to the nozzle exit, and the use of Chevrons diminished these signatures. The LIM contours underline that this feature is very persistent in time for the round jet and that Chevrons increase the intermittent events at these frequencies. At x/D = 18 from the nozzle exit, the main effect of the Chevrons were the disruption of the low-frequency energetic signatures detected in the baseline configuration. This effect was amplified with increased Chevron lobes, which could be linked to the reduction of the jet mixing noise. An increase in intermittency was also observed in the LIM contour maps in this frequency range. These observations were also reconfirmed by the LIM2 analysis, which measures the deviation from a Gaussianity distribution of the PDF events related to frequency localized features. Screech was yet again observed to have a bi-modal distribution in time, while the presence of Chevrons disrupts these signatures and generates a series of intermittent events contributing to a deviation from a Gaussian distribution, increasing the signal kurtosis.
Bi-variate wavelet analysis was included to show the local coherence of two consecutive signals in the near-field region within the potential core and in the jet mixing region. The baseline nozzle showed high levels of time-persistent coherence in the Screech related frequency close to the nozzle exit. The Screech coherence was absent for the Chevron nozzles, other time-frequency signatures were similar for both the tested Chevron nozzles. On the other hand, at the mid-frequency range, the increased number of Chevron lobes reduces the coherence magnitude compared to the round jet. Finally, Screech is less evident in the developed jet region, and the effect of Chevrons results are indistinct. The coherence analysis has also been used to evaluate the local relative phase angle, measuring the detected features’ phase shift in both time and frequency domains. The most important findings were the observation of a phase angle close to π at the Screech frequency for the baseline configuration, which became intermittent with the Chevron presence alternating between positive and negative phase angles. No particular effects were detected at higher axial locations.
Overall, an extensive analysis of the intermittent statistics of the supersonic jet noise features has been reported for the first time. A series of efficient tools that can be applied in a more extensive analysis, oriented at providing all the information needed to model the supersonic jet noise sources has been showcased.

Author Contributions

Conceptualization, S.M. and H.K.J.; Methodology, S.M.; Experiments H.K.J.; Data Analysis, S.M.; Writing—Original Draft Preparation, S.M. and H.K.J.; Writing—Review & Editing, S.M. and H.K.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the authors.

Acknowledgments

Hasan K. J. would like to acknowledge the Engineering and Physical Sciences Research Council (EPSRC) for supporting this research (Grant No. EP/S000917/1). The authors would like to acknowledge Mahdi Azarpeyvand for allocating time for the tests at the aeroacoustic facility in University of Bristol.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meloni, S.; Proença, A.R.; Lawrence, J.L.; Camussi, R. An experimental investigation into model-scale installed jet–pylon–wing noise. J. Fluid Mech. 2021, 929, A4. [Google Scholar] [CrossRef]
  2. Lighthill, M.J. On sound generated aerodynamically I. General theory. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1952, 211, 564–587. [Google Scholar]
  3. Lilley, G. On the noise from air jets. AGARD CP 1974, 131, 13.1–13.2. [Google Scholar]
  4. Camussi, R.; Bogey, C. Intermittent statistics of the 0-mode pressure fluctuations in the near field of Mach 0.9 circular jets at low and high Reynolds numbers. Theor. Comput. Fluid Dyn. 2021, 35, 229–247. [Google Scholar] [CrossRef]
  5. Micci, G.L.; Camussi, R.; Meloni, S.; Bogey, C. Intermittency and Stochastic Modeling of Low- and High-Reynolds-Number Compressible Jets. AIAA J. 2022, 60, 1983–1990. [Google Scholar] [CrossRef]
  6. Huber, J.; Fleury, V.; Bulté, J.; Laurendeau, E.; Sylla, A.A. Understanding and Reduction of Cruise Jet Noise at Aircraft Level. Int. J. Aeroacoust. 2014, 13, 61–84. [Google Scholar] [CrossRef]
  7. Camussi, R.; Ahmad, M.K.; Meloni, S.; de Paola, E.; Di Marco, A. Experimental analysis of an under-expanded jet interacting with a tangential flat plate: Flow visualizations and wall pressure statistics. Exp. Therm. Fluid Sci. 2022, 130, 110474. [Google Scholar] [CrossRef]
  8. Clem, M.; Zaman, K.; Fagan, A. Background Oriented Schlieren Applied to Study Shock Spacing in a Screeching Circular Jet. In Proceedings of the 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Nashville, TN, USA, 9–12 January 2012. [Google Scholar]
  9. Clem, M.; Brown, C.; Fagan, A. Background Oriented Schlieren Implementation in a Jet-Surface Interaction Test. In Proceedings of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Grapevine, TX, USA, 7–10 January 2013. [Google Scholar]
  10. Tam, C. Supersonic jet noise. Ann. Rev. Fluid Mech. 1995, 27, 17–43. [Google Scholar] [CrossRef]
  11. Edgington-Mitchell, D. Aeroacoustic resonance and self-excitation in screeching and impinging supersonic jets—A review. Int. J. Aeroacoust. 2019, 18, 118–188. [Google Scholar] [CrossRef]
  12. Humphrey, N.J.; Edgington-Mitchell, D. The effect of low lobe count chevron nozzles on supersonic jet screech. Int. J. Aeroacoust. 2016, 15, 294–311. [Google Scholar] [CrossRef]
  13. Nikam, S.; Sharma, S. Effect of chevron nozzle penetration on aero-acoustic characteristics of jet at M = 0.8. Fluid Dyn. Res. 2017, 49, 065506. [Google Scholar] [CrossRef] [Green Version]
  14. Massey, S.; Elmiligui, A.; Hunter, C.; Thomas, R.; Pao, P.; Mengle, V. Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics. In Proceedings of the 12th AIAA/CEAS Aeroacoustics Conference (27th AIAA Aeroacoustics Conference), Cambridge, MA, USA, 8–10 May 2006. [Google Scholar]
  15. Schlinker, R.; Simonich, J.; Shannon, D.; Reba, R.; Colonius, T.; Gudmundsson, K.; Ladeinde, F. Supersonic jet noise from round and chevron nozzles: Experimental studies. In Proceedings of the 15th AIAA/CEAS Aeroacoustics Conference (30th AIAA Aeroacoustics Conference), Miami, FL, USA, 11–13 May 2009; p. 3257. [Google Scholar]
  16. Jawahar, H.K.; Meloni, S.; Camussi, R.; Azarpeyvand, M. Experimental Investigation on the Jet Noise Sources for Chevron Nozzles in Under-expanded Condition. In Proceedings of the AIAA AVIATION 2021 FORUM, Virtual Event, 2–6 August 2021. [Google Scholar] [CrossRef]
  17. Bridges, J.; Brown, C. Parametric testing of chevrons on single flow hot jets. In Proceedings of the 10th AIAA/CEAS Aeroacoustics Conference, Manchester, UK, 10–12 May 2004; p. 2824. [Google Scholar]
  18. Camussi, R.; Mancinelli, M.; Di Marco, A. Intermittency and stochastic modeling of hydrodynamic pressure fluctuations in the near field of compressible jets. Int. J. Heat Fluid Flow 2017, 68, 180–188. [Google Scholar] [CrossRef]
  19. Juvé, D.; Sunyach, M.; Comte-Bellot, G. Intermittency of the noise emission in subsonic cold jets. J. Sound Vib. 1980, 71, 319–332. [Google Scholar] [CrossRef]
  20. Tide, P.; Srinivasan, K. Novel chevron nozzle concepts for jet noise reduction. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2009, 223, 51–67. [Google Scholar] [CrossRef]
  21. Tide, P.; Srinivasan, K. Effect of chevron count and penetration on the acoustic characteristics of chevron nozzles. Appl. Acoust. 2010, 71, 201–220. [Google Scholar] [CrossRef]
  22. Meloni, S.; Di Marco, A.; Mancinelli, M.; Camussi, R. Wall pressure fluctuations induced by a compressible jet flow over a flat plate at different Mach numbers. Exp. Fluids 2019, 60, 48–60. [Google Scholar] [CrossRef]
  23. Meloni, S.; Di Marco, A.; Mancinelli, M.; Camussi, R. Experimental investigation of jet-induced wall pressure fluctuations over a tangential flat plate at two Reynolds numbers. Sci. Rep. 2020, 10, 9140. [Google Scholar] [CrossRef]
  24. Meloni, S.; Lawrence, J.L.; Proença, A.R.; Self, R.H.; Camussi, R. Wall pressure fluctuations induced by a single stream jet over a semi-finite plate. Int. J. Aeroacoust. 2020, 19, 240–253. [Google Scholar] [CrossRef]
  25. Camussi, R.; Meloni, S. On the Application of Wavelet Transform in Jet Aeroacoustics. Fluids 2021, 6, 299. [Google Scholar] [CrossRef]
  26. Kamliya Jawahar, H.; Meloni, S.; Camussi, R.; Azarpeyvand, M. Intermittent and stochastic characteristics of slat tones. Phys. Fluids 2021, 33, 025120. [Google Scholar] [CrossRef]
  27. Farge, M. Wavelet transforms and their applications to turbulence. Annu. Rev. Fluid Mech. 1992, 24, 395–458. [Google Scholar] [CrossRef]
  28. Meneveau, C. Analysis of turbulence in the orthonormal wavelet representation. J. Fluid Mech. 1991, 232, 469–520. [Google Scholar] [CrossRef]
  29. Camussi, R.; Robert, G.; Jacob, M.C. Cross-wavelet analysis of wall pressure fluctuations beneath incompressible turbulent boundary layers. J. Fluid Mech. 2008, 617, 11–30. [Google Scholar] [CrossRef] [Green Version]
  30. Torrence, C.; Compo, G.P. A Practical Guide to Wavelet Analysis. Bull. Am. Meteorol. Soc. 1998, 79, 61–78. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the various nozzle configurations used in the present study: (a) SMC000 or baseline nozzle; (b) SMC002; and (c) SMC006.
Figure 1. Schematic of the various nozzle configurations used in the present study: (a) SMC000 or baseline nozzle; (b) SMC002; and (c) SMC006.
Fluids 07 00108 g001
Figure 2. A photo of the experimental setup with the baseline nozzle configuration.
Figure 2. A photo of the experimental setup with the baseline nozzle configuration.
Fluids 07 00108 g002
Figure 3. Near-field SPL spectra at M = 1.3 h/D = 3.5 and various axial locations. (a) SMC000 nozzle; (b) SMC002 nozzle; (c) SMC006 nozzle.
Figure 3. Near-field SPL spectra at M = 1.3 h/D = 3.5 and various axial locations. (a) SMC000 nozzle; (b) SMC002 nozzle; (c) SMC006 nozzle.
Fluids 07 00108 g003
Figure 4. Axial evolution at M = 1.3 and h/D = 3.5 of: (a) Skewness; (b) Kurtosis.
Figure 4. Axial evolution at M = 1.3 and h/D = 3.5 of: (a) Skewness; (b) Kurtosis.
Fluids 07 00108 g004
Figure 5. PDF comparison at M = 1.3 and h/D = 3.5: (a) At x/D = 2; (b) At x/D = 18.
Figure 5. PDF comparison at M = 1.3 and h/D = 3.5: (a) At x/D = 2; (b) At x/D = 18.
Fluids 07 00108 g005
Figure 6. Wavelet scalograms at M = 1.3, h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18 respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18 respectively.
Figure 6. Wavelet scalograms at M = 1.3, h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18 respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18 respectively.
Fluids 07 00108 g006
Figure 7. Fourier spectra of the wavelet coefficients absolute values related to the Screech tone in the baseline configuration. Results are at the same Mach number of the previous plots.
Figure 7. Fourier spectra of the wavelet coefficients absolute values related to the Screech tone in the baseline configuration. Results are at the same Mach number of the previous plots.
Fluids 07 00108 g007
Figure 8. Stochastic analysis for the various frequencies of interest. (a) μ ( | w x | ) at x/D = 2; (b) μ ( | w x | ) at x/D = 2; (c) σ ( | w x | ) at x/D = 2; (d) σ ( | w x | ) at x/D = 18.
Figure 8. Stochastic analysis for the various frequencies of interest. (a) μ ( | w x | ) at x/D = 2; (b) μ ( | w x | ) at x/D = 2; (c) σ ( | w x | ) at x/D = 2; (d) σ ( | w x | ) at x/D = 18.
Fluids 07 00108 g008
Figure 9. LIM contour maps at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18, respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18, respectively.
Figure 9. LIM contour maps at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18, respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18, respectively.
Fluids 07 00108 g009
Figure 10. LIM2 contour maps at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18, respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18, respectively.
Figure 10. LIM2 contour maps at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle at x/D = 2 and x/D = 18, respectively; (c,d) SMC002 nozzle at x/D = 2 and x/D = 18, respectively; (e,f) SMC006 nozzle at x/D = 2 and x/D = 18, respectively.
Fluids 07 00108 g010
Figure 11. Wavelet Coherence contour maps between two consecutive microphones at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18 respectively; (c,d) SMC002 nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (e,f) SMC006 nozzle having the reference microphone at x/D = 2 and x/D = 18, respectively.
Figure 11. Wavelet Coherence contour maps between two consecutive microphones at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18 respectively; (c,d) SMC002 nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (e,f) SMC006 nozzle having the reference microphone at x/D = 2 and x/D = 18, respectively.
Fluids 07 00108 g011
Figure 12. Wavelet phase angle contour maps between two consecutive microphones at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (c,d) SMC002 nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (e,f) SMC006 nozzle having the reference microphone at x/D = 2 and x/D = 18, respectively; the reported contours are normalized by π .
Figure 12. Wavelet phase angle contour maps between two consecutive microphones at M = 1.3 and h/D = 3.5: (a,b) Baseline nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (c,d) SMC002 nozzle having the reference microphone at x/D = 2–4 and x/D = 16–18, respectively; (e,f) SMC006 nozzle having the reference microphone at x/D = 2 and x/D = 18, respectively; the reported contours are normalized by π .
Fluids 07 00108 g012
Table 1. Parameters used for the tested chevron nozzles.
Table 1. Parameters used for the tested chevron nozzles.
Nozzle IDNLength (mm)Angle ( )Penetration (mm)D e (mm) Γ
SMC 0000 16.93330.089
SMC 002410.666750.465017.86670.089
SMC 00667.533318.21.175015.90000.292
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Meloni, S.; Jawahar, H.K. A Wavelet-Based Time-Frequency Analysis on the Supersonic Jet Noise Features with Chevrons. Fluids 2022, 7, 108. https://doi.org/10.3390/fluids7030108

AMA Style

Meloni S, Jawahar HK. A Wavelet-Based Time-Frequency Analysis on the Supersonic Jet Noise Features with Chevrons. Fluids. 2022; 7(3):108. https://doi.org/10.3390/fluids7030108

Chicago/Turabian Style

Meloni, Stefano, and Hasan Kamliya Jawahar. 2022. "A Wavelet-Based Time-Frequency Analysis on the Supersonic Jet Noise Features with Chevrons" Fluids 7, no. 3: 108. https://doi.org/10.3390/fluids7030108

APA Style

Meloni, S., & Jawahar, H. K. (2022). A Wavelet-Based Time-Frequency Analysis on the Supersonic Jet Noise Features with Chevrons. Fluids, 7(3), 108. https://doi.org/10.3390/fluids7030108

Article Metrics

Back to TopTop