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:
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]:
where
s is the wavelet scale,
is a time shift,
is a constant that takes into account the mean value of
and
is the complex conjugate of the dilated and translated mother wavelet
. 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:
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
magnitude of the Screech tone that is persistent over time. A similar trend is followed for the Screech tone for the
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
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
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
, 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:
where
are the wavelet coefficients evaluated with Equation (
4) whilst the symbol
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].
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
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:
where
is the wavelet cross-spectrum of two consecutive signals in the axial direction, while
and
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
,
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:
where
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.