Informed by LES data and resolvent analysis of the mean flow, we examine the structure of turbulence in jets in the subsonic, transonic, and supersonic regimes. Spectral (frequency-space) proper orthogonal decomposition is used to extract energy spectra and decompose the flow into energy-ranked coherent structures. The educed structures are generally well predicted by the resolvent analysis. Over a range of low frequencies and the first few azimuthal mode numbers, these jets exhibit a low-rank response characterized by Kelvin-Helmholtz (KH) type wavepackets associated with the annular shear layer up to the end of the potential core and that are excited by forcing in the very-near-nozzle shear layer. These modes too the have been experimentally observed before and predicted by quasi-parallel stability theory and other approximations-they comprise a considerable portion of the total turbulent energy. At still lower frequencies, particularly for the axisymmetric mode, and again at high frequencies for all azimuthal wavenumbers, the response is not low rank, but consists of a family of similarly amplified modes. These modes, which are primarily active downstream of the potential core, are associated with the Orr mechanism. They occur also as sub-dominant modes in the range of frequencies dominated by the KH response. Our global analysis helps tie together previous observations based on local spatial stability theory, and explains why quasiparallel predictions were successful at some frequencies and azimuthal wavenumbers, but failed at others.

The purpose of this paper is to characterize and model waves that are observed within the potential core of subsonic jets and relate them to previously observed tones in the near-nozzle region. The waves are detected in data from a large-eddy simulation of a Mach 0.9 isothermal jet and modelled using parallel and weakly non-parallel linear modal analysis of the Euler equations linearized about the turbulent mean flow, as well as simplified models based on a cylindrical vortex sheet and the acoustic modes of a cylindrical soft duct. In addition to the Kelvin–Helmholtz instability waves, three types of waves with negative phase velocities are identified in the potential core: upstream- and downstream-propagating duct-like acoustic modes that experience the shear layer as a pressure-release surface and are therefore radially confined to the potential core, and upstream-propagating acoustic modes that represent a weak coupling between the jet core and the free stream. The slow streamwise contraction of the potential core imposes a frequency-dependent end condition on the waves that is modelled as the turning points of a weakly non-parallel approximation of the waves. These turning points provide a mechanism by which the upstream- and downstream-travelling waves can interact and exchange energy through reflection and transmission processes. Paired with a second end condition provided by the nozzle, this leads to the possibility of resonance in limited frequency bands that are bound by two saddle points in the complex wavenumber plane. The predicted frequencies closely match the observed tones detected outside of the jet. The vortex-sheet model is then used to systematically explore the Mach number and temperature ratio dependence of the phenomenon. For isothermal jets, the model suggests that resonance is likely to occur in a narrow range of Mach number,$0.82<M<1$.

To investigate the effects of the nozzle-exit conditions on jet flow and sound fields, large-eddy simulations of an isothermal Mach 0.9 jet issued from a convergent-straight nozzle are performed at a diameter-based Reynolds number of $1\times 10^{6}$. The simulations feature near-wall adaptive mesh refinement, synthetic turbulence and wall modelling inside the nozzle. This leads to fully turbulent nozzle-exit boundary layers and results in significant improvements for the flow field and sound predictions compared with those obtained from the typical approach based on laminar flow in the nozzle. The far-field pressure spectra for the turbulent jet match companion experimental measurements, which use a boundary-layer trip to ensure a turbulent nozzle-exit boundary layer to within 0.5 dB for all relevant angles and frequencies. By contrast, the initially laminar jet results in greater high-frequency noise. For both initially laminar and turbulent jets, decomposition of the radiated noise into azimuthal Fourier modes is performed, and the results show similar azimuthal characteristics for the two jets. The axisymmetric mode is the dominant source of sound at the peak radiation angles and frequencies. The first three azimuthal modes recover more than 97 % of the total acoustic energy at these angles and more than 65 % (i.e. error less than 2 dB) for all angles. For the main azimuthal modes, linear stability analysis of the near-nozzle mean-velocity profiles is conducted in both jets. The analysis suggests that the differences in radiated noise between the initially laminar and turbulent jets are related to the differences in growth rate of the Kelvin–Helmholtz mode in the near-nozzle region.

Direct numerical simulations are performed to investigate the three-dimensional stability of compressible flow over open cavities. A linear stability analysis is conducted to search for three-dimensional global instabilities of the two-dimensional mean flow for cavities that are homogeneous in the spanwise direction. The presence of such instabilities is reported for a range of flow conditions and cavity aspect ratios. For cavities of aspect ratio (length to depth) of 2 and 4, the three-dimensional mode has a spanwise wavelength of approximately one cavity depth and oscillates with a frequency about one order of magnitude lower than two-dimensional Rossiter (flow/acoustics) instabilities. A steady mode of smaller spanwise wavelength is also identified for square cavities. The linear results indicate that the instability is hydrodynamic (rather than acoustic) in nature and arises from a generic centrifugal instability mechanism associated with the mean recirculating vortical flow in the downstream part of the cavity. These three-dimensional instabilities are related to centrifugal instabilities previously reported in flows over backward-facing steps, lid-driven cavity flows and Couette flows. Results from three-dimensional simulations of the nonlinear compressible Navier–Stokes equations are also reported. The formation of oscillating (and, in some cases, steady) spanwise structures is observed inside the cavity. The spanwise wavelength and oscillation frequency of these structures agree with the linear analysis predictions. When present, the shear-layer (Rossiter) oscillations experience a low-frequency modulation that arises from nonlinear interactions with the three-dimensional mode. The results are consistent with observations of low-frequency modulations and spanwise structures in previous experimental and numerical studies on open cavity flows.

Gudmundsson and Colonius (J. Fluid Mech., vol. 689, 2011, pp. 97-128) have recently shown that the average evolution of low-frequency, low-azimuthal modal large-scale structures in the near field of subsonic jets are remarkably well predicted as linear instability waves of the turbulent mean flow using parabolized stability equations. In this work, we extend this modelling technique to an isothermal and a moderately heated Mach 1.5 jet for which the mean flow fields are obtained from a high-fidelity large-eddy simulation database. The latter affords a rigourous and extensive validation of the model, which had only been pursued earlier with more limited experimental data. A filter based on proper orthogonal decomposition is applied to the data to extract the most energetic coherent components. These components display a distinct wavepacket character, and agree fairly well with the parabolized stability equations model predictions in terms of near-field pressure and flow velocity. We next apply a Kirchhoff surface acoustic propagation technique to the near-field pressure model and obtain an encouraging match for far-field noise levels in the peak aft direction. The results suggest that linear wavepackets in the turbulence are responsible for the loudest portion of the supersonic jet acoustic field.

Acoustic analogy methods are used as post-processing tools to predict aerodynamically generated sound from numerical solutions of unsteady flow. The Ffowcs WilliamsHawkings (FW-H) equation and related formulations, such as Farassat's Formulations 1 and 1A, are among the commonly used analogies because of their relative low computation cost and their robustness. These formulations assume the propagation of sound waves in a medium at rest. The present paper describes a surface integral formulation based on the convective wave equation, which takes into account the presence of a mean flow. The formulation was derived to be easy to implement as a numerical postprocessing tool for computational fluid dynamics codes. The new formulation constitutes one possible extension of Farassat's Formulation 1 and 1A based on the convective form of the FW-H equation.

No abstract

Coherent features of a turbulent Mach 0.9, Reynolds number$10^{6}$jet are educed from a high-fidelity large eddy simulation. Besides the well-known Kelvin–Helmholtz instabilities of the shear layer, a new class of trapped acoustic waves is identified in the potential core. A global linear stability analysis based on the turbulent mean flow is conducted. The trapped acoustic waves form branches of discrete eigenvalues in the global spectrum, and the corresponding global modes accurately match the educed structures. Discrete trapped acoustic modes occur in a hierarchy determined by their radial and axial order. A local dispersion relation is constructed from the global modes and found to agree favourably with an empirical dispersion relation educed from the simulation data. The product between direct and adjoint modes is then used to isolate the trapped waves. Under certain conditions, resonance in the form of a beating occurs between trapped acoustic waves of positive and negative group velocities. This resonance explains why the trapped modes are prominently observed in the simulation and as tones in previous experimental studies. In the past, these tones were attributed to external factors. Here, we show that they are an intrinsic feature of high-subsonic jets that can be unambiguously identified by a global linear stability analysis.