Compressible jets impinging on a perpendicular surface can produce high-intensity, discrete-frequency tones. The character of these tones is a function of nozzle shape, jet Mach number, impingement-plate geometry, and the distance between nozzle and plate. Though it has long been recognised that these tones are associated with a resonance cycle, the exact mechanism by which they are generated has remained a topic of some debate. In this work, we present evidence for a number of distinct tone-generation mechanisms, reconciling some of the different findings of prior authors. We demonstrate that the upstream-propagating waves that close resonance can be confined within the jet, or external to it. These waves can be either weak and relatively linear, or strong and nonlinear from their inception. The waves can undergo coalescence or merging, and in some configurations, pairs of waves rather than singletons appear. We discuss both historical and new evidence for multiple distinct processes by which upstream-propagating waves are produced: direct vortex sound, shock leakage, wall-jet-boundary fluctuations, and wall-jet shocklets. We link these various mechanisms to the disparate collection of upstream-propagating waves observed in the data. We also demonstrate that multiple mechanisms can be provoked by a single vortex, providing an explanation as to why sometimes pairs of waves or merging waves are observed. Through this body of work, we demonstrate that rather than being in opposition, the various pieces of past research on this topic were simply identifying different mechanisms that can support resonance.

In the present work, we experimentally investigate the transverse injection of elliptic liquid jets into a supersonic cross-flow ($M_\infty$ = 2.5). The primary focus is to understand the effect of injection orifice aspect ratio ($\textit{AR}$ = spanwise/streamwise dimension), on the liquid jet breakup mechanism, the flow field around the liquid jet and the resulting droplet sizes formed downstream, for three $\textit{AR}$ cases ($\textit{AR}$ = 0.3, 1, 3.3). We find that the $\textit{AR}$ = 0.3 case has large unsteadiness in the spray core due to relatively large wavelength Rayleigh–Taylor (RT) waves formed on the liquid jet surface. However, the primary jet breakup occurs through Kelvin–Helmholtz (KH) instabilities formed on the large lateral surfaces, as in coaxial liquid jet breakup. This leads to a higher Sauter mean diameter (SMD) of the droplets in the spray core with a wider range of droplet sizes compared with the circular case ($\textit{AR}$ = 1.0). However, in the case of $\textit{AR}$ = 3.3, the RT waves are more intense and of smaller wavelength due to the large drag on the liquid jet, which results in a direct catastrophic breakup of the liquid jet by the RT waves. This results in a relatively steady liquid jet and shock structure with the formation of a fine spray and smaller droplets in the spray core than for the $\textit{AR}=1.0$ case. The study shows the importance of the orifice $\textit{AR}$ on the flow around, and the spray downstream of, the liquid jet injection into supersonic cross-flow.

Supersonic free jets are extensively employed across a range of applications, especially in high-tech industries such as semiconductor processing and aerospace propulsion. Due to the difficulties involved in flow measurement, previous research on supersonic free jets has primarily focused on investigating near-field shockwave structures, with quantitative experimental analysis of the far-field zone being relatively scarce. However, physical understanding of the far-field flow, particularly post-shockwave energy dissipation, holds significant importance for the application and utilisation of these jets in vacuum environments. Therefore, this study aims to provide a robust experimental foundation for a rarefied supersonic free jet through the analysis of the flow field in both the near- and far-field zones. Nanometre-sized tracer particles and molecules were utilised to measure the rarefied supersonic jet flow field using particle image velocimetry and acetone molecular tagging velocimetry, respectively. The experiments revealed that in rarefied conditions, the supersonic jet exhibits a one-barrel shockwave structure in the near field, and after passing the Mach disk, a long annular viscous layer develops downstream. Experimental data on the jet velocity profile and width demonstrated a transition to a laminar flow regime in the far-field zone. This transition aligns with the theoretically inferred flow regimes based on the complex Reynolds number. The velocity profile and potential core length of the laminar flow regime could be modelled using a bi-modal distribution, which represents the summation of symmetric Gaussian distributions.

The investigation of shock/blast wave diffraction over various objects has garnered significant attention in recent decades on account of the catastrophic changes that these waves inflict on the environment. Equally important flow phenomena can occur when the moving expansion waves diffract over bodies, which has been hardly investigated. To investigate the effect of expansion wave diffraction over different bodies, we conducted shock tube experiments and numerical simulations to visualise the intricate wave interactions that occur during this process. The current investigation focuses on the phenomenon of expansion wave diffraction across three distinct diffracting configurations, namely the bluff, wedge and ogive bodies. The diffraction phenomenon is subsequently investigated under varying expansion wave strengths through the control of the initial diaphragm rupture pressure ratios. The shock waves generated by the expansion wave diffraction in the driver side of the shock tube, which was initially identified in numerical simulations by Mahomed & Skews (2014 J. Fluid Mech., vol. 757, pp. 649–664), have been visualised in the experiments. Interesting flow features, such as unsteady shock generation, transition, and symmetric/asymmetric vortex breakdown, have been observed in these expansion flows. An in-depth analysis of such intricate flow features resulting from expansion wave diffraction is performed and characterised in the current study.

Görtler vortices induced by concave curvature in supersonic turbulent flows are investigated using resolvent analysis and large-eddy simulations at Mach 2.95 and Reynolds number $ Re_{\delta }=63\,500$ based on the boundary-layer thickness $ \delta$. Resolvent analysis reveals that the most amplified coherent structures manifest as streamwise counter-rotating vortices with optimal spanwise wavelength $ 2.4\delta$ at cut-off frequency $f\delta /{u}_{\infty } =0.036$, where $ {u}_{\infty }$ is the freestream velocity. The leading spectral proper orthogonal decomposition modes with spanwise wavelength approximately $ 2\delta$ align well with the predicted coherent structures from resolvent analysis at $f\delta /{u}_{\infty } =0.036$. These predicted and extracted coherent structures are identified as Görtler vortices, driven by the Görtler instability. The preferential spanwise scale of the Görtler vortices is further examined under varying geometric and freestream parameters. The optimal spanwise wavelength is insensitive to the total turning angle beyond a critical value, but sensitive to the concave curvature $ K$ at the same turning angle. A limit spanwise wavelength $ 1.96\delta$, corresponding to an infinite concave curvature as $ K\rightarrow \infty$, is identified and validated. Increasing the freestream Mach number or decreasing the ratio of wall temperature to freestream temperature reduces the optimal wavelength normalised by $ \delta$, while variations in freestream Reynolds number have negligible impact. Additionally, a modified definition of the turbulent Görtler number $ G_{T}$ based on the peak eddy viscosity in boundary layers is proposed and employed to assess the occurrence of Görtler instability.

The supersonic wake of a circular cylinder in Mach 3 flow was studied through spectral proper orthogonal decomposition (SPOD) of high-speed focussing schlieren datasets. A wavenumber decomposition of the SPOD eigenvectors was found to be an effective tool for isolating imaging artefacts from the flow features, resulting in a clearer interpretation of the SPOD modes. The cylinder wake consists of both symmetric and antisymmetric instabilities, with the former being the dominant type. The free shear layers that form after the flow separates from the cylinder surface radiate strong Mach waves that interact with the recompression shocks to release significant disturbances into the wake. The wake shows a bimodal vortex shedding behaviour with a purely hydrodynamic instability mode around a Strouhal number of 0.2 and an aeroacoustic instability mode around Strouhal number of 0.42. The hydrodynamic mode, which is presumably the same as the incompressible case, is weaker and decays rapidly as the wake accelerates due to increasing compressibility. The aeroacoustic mode is the dominant shedding mode and persists farther into the wake because of an indirect energy input received through free-stream acoustic waves. A simple aeroacoustic feedback model based on an interaction between downstream propagating shear-layer instabilities and upstream propagating acoustic waves within the recirculation region is shown to accurately predict the shedding frequency. Based on this model, the vortex shedding in supersonic flows over a circular cylinder occurs at a universal Strouhal number (based on approach free-stream velocity and feedback path length) of approximately 0.3.

In the present study, we investigate the modulation effects of particles on compressible turbulent boundary layers at a Mach number of 6, employing high-fidelity direct numerical simulations based on the Eulerian–Lagrangian point-particle approach. Our findings reveal that the mean and fluctuating velocities in particle-laden flows exhibit similarities to incompressible flows under compressibility transformations and semi-local viscous scaling. With increasing particle mass loading, the reduction in Reynolds shear stress and the increase in particle feedback force constitute competing effects, leading to a non-monotonic variation in skin friction, particularly in turbulence over cold walls. Furthermore, dilatational motions near the wall, manifested as travelling-wave structures, persist under the influence of particles. However, these structures are significantly weakened due to the suppression of solenoidal bursting events and the negative work exerted by the particle feedback force. These findings align with the insight of Yu et al. (J. Fluid. Mech., vol. 984, 2024, A44), who demonstrated that dilatational motions are generated by the vortices associated with intense bursting events, rather than acting as evolving perturbations beneath velocity streaks. The attenuation of travelling-wave structures at higher particle mass loadings also contributes to the reduction in the intensities of wall shear stress and heat flux fluctuations, as well as the probability of extreme events. These results highlight the potential of particle-laden flows to mitigate aerodynamic forces and thermal loads in high-speed vehicles.

This article delves into the dynamics of inviscid annular supersonic jets, akin to those exiting converging–diverging nozzles in over-expanded regimes. It focuses on the first azimuthal Fourier mode of flow fluctuations and examines their behaviour with varying mixing layer parameters and expansion regimes. The study reveals that two unstable Kelvin–Helmholtz waves exist in all cases, with the outer-layer wave being more unstable due to differences in the velocity gradient. The inner-layer wave is more sensitive to changes in base flow and extends beyond the jet, potentially contributing to nozzle resonances. The article also investigates upstream propagating guided-jet modes, which are found to be robust and not highly sensitive to changes in base flow, which makes them essential for understanding jet dynamics. A simplified model is used to obtain ideal base flows but with realistic shape in order to study the effects of varying nozzle pressure ratios on the dynamics of the waves supported by the jet.

This work is a numerical study of a transitional shock wave boundary layer interaction (SWBLI). The main goal is to improve our understanding of the well-known low-frequency SWBLI unsteadiness and especially the suspected role of triadic interactions in the underlying physical mechanism. To this end, a direct numerical simulation is performed using a high-order finite-volume scheme equipped with a suitable shock capturing procedure. The resulting database is then extensively post-processed in order to extract the main dynamical features of the interaction zone dynamics (involved characteristic frequencies, characteristics of the vortical structures, etc.). The dynamical organisation and space–time evolution of the flow at dominant frequencies are then further characterised by mean of an spectral proper orthogonal decomposition analysis. In order to study the role of triadic interactions occurring in the interaction region, a bispectral mode decomposition analysis is applied to the database. It allows us to extract the significant triadic interactions, their location and the resulting physical spatial modes. Strong triadic interactions are detected in the downstream part of the separation bubble whose role on the low-frequency unsteadiness is characterised. All the results of the various analyses are then discussed and integrated to formulate a possible mechanism fuelling low-frequency SWBLI unsteadiness.

The study examines supersonic square jets in a twin nozzle configuration with the aim of identifying and characterising emergent instability modes during overexpanded operation. Unlike screeching rectangular jets that undergo strong fluctuations normal to the wider jet dimension, the equilateral nature of the exit geometry in square nozzles leads to multiple instability states dictated by shock–turbulence interactions and nozzle operating conditions. Furthermore, strong coupling modes between the jets were identified that led to either phase locked or out of phase interactions of the inner shear layers. Results from experimental studies were examined using spatial and temporal decomposition techniques based on spectral methods to identify the resultants from triadic shock–turbulence interactions. The primary instability mode across both operating conditions were driven by optimal interactions while the harmonics were found to be associated with the suboptimal shock–turbulence interactions.

Although active flow control based on deep reinforcement learning (DRL) has been demonstrated extensively in numerical environments, practical implementation of real-time DRL control in experiments remains challenging, largely because of the critical time requirement imposed on data acquisition and neural-network computation. In this study, a high-speed field-programmable gate array (FPGA) -based experimental DRL (FeDRL) control framework is developed, capable of achieving a control frequency of 1–10 kHz, two orders higher than that of the existing CPU-based framework (10 Hz). The feasibility of the FeDRL framework is tested in a rather challenging case of supersonic backward-facing step flow at Mach 2, with an array of plasma synthetic jets and a hot-wire acting as the actuator and sensor, respectively. The closed-loop control law is represented by a radial basis function network and optimised by a classical value-based algorithm (i.e. deep Q-network). Results show that, with only ten seconds of training, the agent is able to find a satisfying control law that increases the mixing in the shear layer by 21.2 %. Such a high training efficiency has never been reported in previous experiments (typical time cost: hours).

We investigate the effects of thermal boundary conditions and Mach number on turbulence close to walls. In particular, we study the near-wall asymptotic behaviour for adiabatic and pseudo-adiabatic walls, and compare to the asymptotic behaviour recently found near isothermal cold walls (Baranwal et al. 2022. J. Fluid Mech. 933, A28). This is done by analysing a new large database of highly-resolved direct numerical simulations of turbulent channels with different wall thermal conditions and centreline Mach numbers. We observe that the asymptotic power-law behaviour of Reynolds stresses as well as heat fluxes does change with both centreline Mach number and thermal condition at the wall. Power-law exponents transition from their analytical expansion for solenoidal fields to those for non-solenoidal field as the Mach number is increased, though this transition is found to be dependent on the thermal boundary conditions. The correlation coefficients between velocity and temperature are also found to be affected by these factors. Consistent with recent proposals on universal behaviour of compressible turbulence, we find that dilatation at the wall is the key scaling parameter for these power-law exponents, providing a universal functional law that can provide a basis for general models of near-wall behaviour.

The reflection of multiple incident shock waves that converge to a single point on the reflecting surface is studied in this paper. The number of the incident shocks, denoted $K$, is arbitrary. The interaction between the reflected shock of one incident shock and the other incident shocks may produce various possible configurations, such as type-I, type-II and type-IV shock interferences. The number of possible reflection configurations is shown to be an exponential function of ($K-1$) with base 2. The possibility of pre-, middle- and post-Mach reflections, which means Mach reflection occurs for the first, middle and last incident shock, is revealed through numerical simulation for $K=3$. For the particular case where the incident shocks are produced by equal variation of wedge surface deflection, the conventional von Neumann condition and detachment condition for the $k\mathrm{th}$ incident shock to have Mach reflection are derived. It is shown that the von Neumann condition for regular reflection is lowered and the detachment condition for Mach reflection is elevated as $k$ increases. The shock reflection patterns for $ K=1,2,\ldots ,10$ are obtained by numerical simulations. We observe a shock interaction train structure, where we have pre-Mach reflection followed by ($K-1$) type-I or type-II shock interferences. We also observe that the Mach stem height decreases with $K$ well above the von Neumann condition and becomes non-monotonic near the von Neumann condition.

A combination of physics-based and data-driven post-processing techniques is proposed to extract acoustic-related shear-layer perturbation responses directly from spatio-temporally resolved schlieren video. The physics-based component uses momentum potential theory to extract the irrotational (acoustic and thermal) component from density gradients embedded in schlieren pixel intensities. For the unheated shear layer, the method filters acoustic structures and tones not evident in the raw data. The filtered data are then subjected to an efficient data-driven dynamic mode decomposition reduced-order model, which provides the forced acoustic perturbation response for broad parameter ranges. A shear layer comprising Mach 2.461 and 0.175 streams, corresponding to a convective Mach number 0.88 and containing shocks, is adopted for illustration. The overall perturbation response is first obtained using an impulse forcing in the wall-normal direction of the splitter plate, extending into subsonic and supersonic streams. Subsequently, impulse and harmonic forcings are independently applied in a pixel-by-pixel manner for a precise receptivity study. The acoustic response shows a convective wavepacket and acoustic burst from the splitter plate. The interaction with the shock and associated wave dispersion emits a second, slower, acoustic wave. Harmonic forcing indicates higher frequency-dependent sensitivity in the supersonic stream, with the most sensitive location near the outer boundary-layer region, which elicits an order of magnitude larger acoustic response compared with disturbances in the subsonic stream. Some receptive forcing regions do not generate significant acoustic waves, which may guide excitation with low noise impact.

We investigate the stability of a compressible boundary layer over an impedance wall for both constant impedances and a frequency-dependent porous wall model. For an exponential mean flow profile, the solution of the Pridmore-Brown equation, i.e. the linearised Euler equations for compressible shear flows, is expressed exactly in confluent Heun functions and, with the boundary condition of acoustic wall impedance, reduced to a single algebraic eigenvalue equation. This, in turn, is solved asymptotically and numerically and provides the complete inviscid eigenvalue spectrum without spurious modes. The key finding is that impedance walls not only have a desirable stabilising effect on inviscid disturbances, but also induce new instabilities. The type of the destabilised mode and therefore also the direction of propagation of the modes with maximum growth rate as well as the destabilised wavenumbers depend significantly on the porous wall properties, in particular on the porous wall layer thickness. For small porous layer thicknesses, the impedance-induced instability is observed as a second mode instability, where we find above a critical porosity growth rates exceeding those present in the rigid-wall case.

Supersonic impinging tones have been attracting significant interest because high-intensity discrete-frequency tones pose substantial risks to structural safety in applications such as rocket launch and recovery, and space vehicle attitude adjustment. However, various issues remain to be addressed regarding the jet oscillation and tone generation mechanism. In this study, a numerical simulation of the supersonic impinging jet with a nozzle pressure ratio of 4.03 and an impingement distance of 2.08 times the nozzle exit diameter is conducted. The results show good consistency with the reference data by other researchers. A phase-locked averaging analysis of 2960 flow field snapshots is employed to investigate jet structure oscillation dynamics and the tone generation mechanism. The phase-locked averaged images reveal that the pressure variation induced by Kelvin–Helmholtz vortices as they pass through the reflected shock results in the periodic motions of the reflected shock and Mach disk. The periodically oscillating Mach disk generates high-pressure fluid masses driving recirculation bubbles through a cyclic ‘compression–generation–merging’ oscillation. The streamline oscillation and sound-ray analyses reveal there are two distinct tone source regions: the impinging zone and the wall jet region. Consequently, it is proposed that vortex collapse in conjunction with wall jet oscillations coexist to generate the tone. According to the directivity, the tone emitted from the wall jet source region is believed to contribute to the feedback loop. These findings collectively contribute to an improved understanding of the jet plume oscillation and tone generation mechanisms of the supersonic impinging jet.

Wave propagation in channels with area changes is a topic of significant practical interest that involves a rich set of coupled physics. While the acoustic wave problem has been studied extensively, the shock propagation problem has received less attention. In addition to its practical significance, this problem also introduces deep fundamental issues associated with how energy in propagating large-amplitude disturbances is redistributed upon interaction with inhomogeneities. This paper presents a study of shock scattering and entropy and vorticity coupling for shock wave propagation through discrete area changes. It compares results from computational fluid dynamics to those of one-dimensional quasi-steady calculations. The solution space is naturally divided into five ‘regimes’ based upon the incident shock strength and area ratio. This paper also presents perturbation methods to quantify the dimensionless scaling of physical effects associated with wave reflection/transmission and energy transfer to other disturbances. Finally, it presents an analysis of the ‘energetics’ of the interaction, quantifying how energy that initially resides in dilatational disturbances and propagates at the shock speed is redistributed into finite-amplitude reflected and transmitted waves as well as convecting vortical and entropy disturbances.

In this paper, the reflection of shock waves with downstream expansion fan interference in two-dimensional, inviscid flow is investigated, including steady Mach reflection (MR) and the unsteady transition process from regular reflection (RR) to MR. A threshold for the configuration based on non-dimensional wedge length is proposed. The analytical model for the steady MR and RR$\rightarrow$MR transition process is established based on the classical shock and expansion wave relations, whose prediction agrees well with results obtained through inviscid numerical simulation. It is found that the expansion fan interference significantly influences the steady flow patterns, especially the height of the Mach stem and the shape of the slip line. The interaction accelerates the formation of the sonic throat, stabilizing the flow structure rapidly, and results in generally small Mach stem heights. The exposure of the triple point to the expansion fan eliminates the inflection point on the slip line, whose slope increases smoothly. The interaction further affects the time evolution of the Mach stem during the multiple-interaction stage of the RR$\rightarrow$MR transition process. It appears that the modifications come from the curvature of the incident shock brought by the wave interference. During the multiple-interaction stage, the triple point moves upstream along the curved incident shock, where the incident shock angle changes according to the curvature, resulting in the variation of the evolution velocity.

The interaction between planar incident shocks and cylindrical boundary layers is prevalent in missiles equipped with inverted inlets, which typically leads to substantial three-dimensional flow separation and the formation of vortical flow. This study utilizes wind-tunnel experiments and theoretical analysis to elucidate the shock structure, surface topology and pressure distributions induced by a planar shock with finite width impinging on a cylinder wall at Mach 2.0. In the central region, a refraction phenomenon occurs as the transmitted shock bends within the boundary layer, generating a series of compression waves that coalesce into a shock, forming a ‘shock triangle’ structure. As the incident shock propagates backward along both sides, it gradually evolves into a Mach stem, where the transmitted shock refracts the expansion wave. The incident shock interacts with the boundary layer, resulting in the formation of a highly swept separation region that yields a pair of counter-rotating horseshoe-like vortices above the separation lines. These vortices facilitate the accumulation of low-energy fluid on both sides. Although the interaction of the symmetry plane aligns with free-interaction-theory, the separation shock angle away from the centre significantly deviates from the predicted value owing to the accumulation of low-energy fluids. The primary separation line and pressure distribution jointly exhibit an elliptical similarity on the cylindrical surface. Furthermore, the potential unsteady behaviour is assessed, and the Strouhal number of the low-frequency oscillation is found to be 0.0094, which is insufficient to trigger significant alterations in the flow field structure.