In this paper, aerobreakup in the stagnation region of high-Mach-number flow over a bluff body is studied with experiments and computations. Water drops of diameter 0.51–2.30 mm were acoustically levitated at sea level along the flight path of a rectangular $100\ {\rm mm} \times 150\ {\rm mm}$ rail-gun launched projectile. This enabled the study of aerobreakup at high Mach (3.03–5.12), post-shock Mach (1.5–1.9), Weber $(5 \times 10^4\unicode{x2013}4 \times 10^5)$ and Reynolds $(6 \times 10^4\unicode{x2013}3 \times 10^5)$ numbers. High-speed backlit shadowgraphy is used to record the flow structure. Computations are made for two cases, and it was found that the drop behaviour is not dominated by viscous or surface-tension effects and can be adequately captured by treating the gas as calorically perfect with the ratio of specific heats set to 1.3 to account for thermochemical effects. To assess drop surface stability at early breakup times, results from Newton's inclination method are used to determine the flow along the drop surface and input to a linear-stability analysis; from this, it was found that viscosity and surface tension can be neglected. Moreover, the acceleration term dominates the shear term at the stagnation point, a point accentuated as a drop flattens; this relation inverts closer to the drop equator. Linear-stability analysis was insufficient for modelling late-stage aerobreakup because the predicted wavelengths were too small and the expected aerobreakup times were non-physically short. To address this discrepancy, a nonlinear instability model with constant-rate growth is used that treats the accelerated drop surface as analogous to bubbles rising through a liquid; agreement with computations is good.

Developed in this study is a theoretical description of squeeze-film lubrication systems that involve the flexural oscillation of a thin plate near a parallel wall. Such systems were discovered in recent experiments to produce load-bearing attractive forces that are a thousandfold stronger than those generated by rigid oscillators, which typically favour repulsion. Analyses of squeeze-film gas flow driven by a presumed plate deformation reproduce the observed magnification of attractive load capacity, but exhibit serious discrepancies with crucial aspects of the experimental measurements – most importantly, the precise distribution of air pressure along the film. The discrepancies are resolved in this study by accounting for the presence of two-way-coupled fluid–structure interactions whereby the undulations of the plate, modelled here with use of the classical Kirchhoff–Love equation, are affected non-negligibly by the evolving pressure, described by a modified Reynolds lubrication equation that accounts for compressibility. The resulting problem of elastohydrodynamic lubrication is solved with use of perturbation methods that exploit the limit of small oscillation amplitudes. The analysis ultimately provides an explicit expression specifying the attractive load capacity of a squeeze-film system as a function of relevant operating parameters – including, in particular, the amplitude and frequency of the localized excitation force exerted on the plate. The rudimentary theory derived here may be readily generalized to guide the analysis and development of a wide variety of emerging engineering systems that exploit the vibration-induced squeeze-film effect – such as wall-climbing soft robots and contactless grippers.

In the astrophysics community it is common practice to model collisionless dust, entrained in a gas flow, as a pressureless fluid. However, a pressureless fluid is fundamentally different from a collisionless fluid – the latter of which generically possess a non-zero anisotropic pressure or stress tensor. In this paper we derive a fluid model for collisionless dust, entrained in a turbulent gas, starting from the equations describing the motion of individual dust grains. We adopt a covariant formulation of our model to allow for the geometry and coordinate systems prevalent in astrophysics, and provide a closure valid for the accretion disc context. We show that the continuum mechanics properties of a dust fluid corresponds to a higher-dimensional anisotropic Maxwell fluid, after the extra dimensions are averaged out, with a dynamically important rheological stress tensor. This higher-dimensional treatment has the advantage of keeping the dust velocity and velocity of the fluid seen, and their respective moments, on the same footing. This results in a simplification of the constitutive relation describing the evolution of the dust rheological stress.

Compound flows consist of two or more parallel compressible streams in a duct and their theoretical treatment has gained attention for the analysis and modelling of ejectors. Recent works have shown that these flows can experience choking upstream of the geometric throat. While it is well known that friction can push the sonic section downstream of the throat, no mechanism has been identified yet to explain its displacement in the opposite direction. This study extends the existing compound flow theory and proposes a one-dimensional (1-D) model, including friction between the streams and the duct walls. The model captures the upstream and downstream displacements of the sonic section. Through an analytical investigation of the singularity at the sonic section, it is demonstrated that friction between the streams is the primary driver of upstream displacement. The 1-D formulation is validated against axisymmetric Reynolds averaged Navier–Stokes simulations of a compound nozzle for various inlet pressures and geometries. The effect of friction is investigated using an inviscid simulation for the isentropic case and viscous simulations with both slip and no-slip conditions at the wall. The proposed extension accurately captures the displacement of the sonic section, offering a new tool for in-depth analysis and modelling of internal compound flows.

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

The gas dynamics of shock-induced gas filtration through densely packed granular columns with vastly varying shock intensity and the structural parameters are numerically investigated using a coupled Eulerian–Lagrangian approach. The results shed fundamental light on the thermal effects of the shock-induced gas filtration manifested by a distinctive self-heating hot gas layer traversing the medium. The characteristics of the thermal effects in terms of the thermal intensity and uniformity are found to vary with the shock Mach number, Ms, and the filtration coefficient of the granular media, Π. As the incident shock transitions from weak to strong, and (or) the filtration coefficient increases from O(10−5) to O(104), the heating mechanisms transition between three distinct heating modes. A phase diagram of heating modes is established on the parameter space (Ms, Π), which enables us to predict the characteristics of the thermal effect in different shock-induced gas filtrations. The thermal effects markedly accelerate the pressure diffusion due to the additional heat influx when the time scale of the former is smaller than or comparable to the latter. Based on the contour map displaying the coupling degree of the thermal effects and the pressure diffusion, we identify a decoupling criterion whereby the isothermal assumption holds if only the pressure diffusion is concerned. The thermal effects may well bring about considerable thermal shocks which pose a great threat to the integrity of the solid skeleton and further reduce the overall shock resistance performance of the porous media.

The dynamics of a shock-induced separation unit generated by a 20$^\circ$ sharp fin placed on a cylindrical surface in a Mach 2.5 flow was investigated. Specifically, the present work investigated the mechanisms that govern the mid-frequency range of separation shock unsteadiness in the fin shock wave–boundary layer interaction (SBLI) unit. Two-dimensional pressure fields were obtained over the cylinder surface spanning the entire fin SBLI unit using high-bandwidth pressure-sensitive paint at 40 kHz imaging rate that allowed probing the low- through mid-frequency ranges of the separation shock unsteadiness. The mean pressure field showed a progressive weakening of the separation shock with downstream distance, which is an artifact of the three-dimensional relief offered by the curved mounting surface. The root-mean-square (r.m.s.) pressure field exhibited a banded structure with elevated $p_{r.m.s.}$ levels beneath the intermittent region, separation vortex and adjacent to the fin root. The power spectral density (PSD) of the surface pressure fluctuations obtained beneath the intermittent region revealed that the separation shock oscillations exhibited the mid-frequency content over the majority of its length. Interestingly, neither the PSD nor the length of the intermittent region varied noticeably with downstream distance, revealing a constant separation shock foot velocity along the entire SBLI. The pressure fluctuation PSD beneath the separation vortex also exhibited the broadband peak at the mid-frequency range of the separation shock motions over the majority of its length within the measurement domain. By contrast, the region adjacent to the fin root exhibited pressure oscillations at a substantially lower frequency compared with the separation shock and the separation vortex. Two-point coherence and cross-correlation analysis provided unique insights into the critical sources and mechanisms that drive the separation shock unsteadiness. The separation vortex and separation shock dynamics were found to be driven by a combination of convecting perturbations that originated from the vicinity of the fin leading edge and the local interactions of the separated flow with the incoming boundary layer. The boundary layer locally strengthened or weakened the convecting pressure perturbations depending on the local momentum fluctuations within the boundary layer.

Shock refraction in a gas–liquid interface is ubiquitous in nature and engineering. This study investigates the shock refraction phenomena in air–water interfaces for various inclination angles. The interface inclination angles are achieved using a tiltable vertical shock tube. The time-resolved schlieren images are compared with numerical simulations performed using the BlastFoam solver in the OpenFOAM software. The stiffened gas equation of state is used to model water in the simulations. The shock polar analysis using modified shock relations for a stiffened gas is used to elucidate the refraction patterns. A regular refraction pattern with a reflected shock wave and a bound precursor refraction with a regular reflection are observed experimentally for the first time in an air–water system. Further, a new free precursor refraction pattern with a Mach reflection is observed. The transition criteria and the corresponding boundaries for each refraction pattern are demarcated in the ($M_S, \theta _w^c$)-plane. The refraction sequence and the range for various incident shock strength regimes are also identified.

Standing acoustic waves have been known to generate Eulerian time-mean ‘streaming’ flows at least since the seminal investigation of Lord Rayleigh in the 1880s. Nevertheless, a recent body of numerical and experimental evidence has shown that inhomogeneities in the ambient density distribution lead to much faster flows than arise in classical Rayleigh streaming. The emergence of these unusually strong flows creates new opportunities to enhance heat transfer in systems in which convective cooling cannot otherwise be easily achieved. To assess this possibility, a theoretical study of acoustic streaming in an ideal gas confined in a rectangular channel with top and bottom walls maintained at fixed but differing temperatures is performed. A two time scale system of equations is utilized to efficiently capture the coupling between the fast acoustic waves and the slowly evolving streaming flow, enabling strongly nonlinear regimes to be accessed. A large suite of numerical simulations is carried out to probe the streaming dynamics, to highlight the critical role played by baroclinically generated wave vorticity and to quantify the additional heat flux induced by the standing acoustic wave. Proper treatment of the two-way coupling between the waves and mean flow is found to be essential for convergence to a self-consistent steady state, and the variation of the resulting acoustically enhanced steady-state heat flux with both the amplitude of the acoustic wave and the $O(1)$ aspect ratio of the channel is documented. For certain parameters, heat fluxes almost two orders of magnitude larger than those realizable by conduction alone can be attained.

Cavity evolution in granular media is crucial in explosion-driven gas–particle flows, particularly in many engineering applications. In this study, a theoretical model was first proposed to describe the cavity evolution in granular media by extending the classical Rayleigh–Plesset model. A closed equation set comprising the radius, pressure and gas leak-off velocity equations was built by considering the gas expansion and non-Darcy gas-penetration effects. Both centrally symmetric and non-centrally symmetric cases of gas injection into granular media were investigated. Especially for modelling the non-symmetric scenario, the radius and gas leak-off velocity equations were proposed in each radial direction with angle $\theta$, and then the pressure equation was built up based on the integral gas leak-off along the cavity outline. Through non-dimensionalizing the theoretical equations, four key dimensionless numbers $\varPi_1,\ \varPi_4$ were obtained to characterize the response time of cavity expansion and the intensity of non-Darcy effects for both cases. This allowed us to determine a scaling law of effective cavity radius $R_{eff}^*=\sqrt {2\varPi _2/(7{\rm \pi} )}t^{*1/2}$ and the critical time $t_{cr}^*=\sqrt {12.5/\varPi _1}$ for two-dimensional cavity evolution. Additionally, the necessity of incorporating non-Darcy effects was ascertained under conditions of $\varPi _4>400$. The findings demonstrate that the proposed theoretical equations effectively predict the cavity evolution results under various operational conditions ($0.7<\varPi _1<7\times 10^2, 3<\varPi _4<1.1\times 10^3$), as validated by refined Euler–Lagrange numerical simulations.

Shock waves are of great interest in many fields of science and engineering, but the mechanisms of their formation, maintenance and dissipation are still not well understood. While all transport processes existing in a shock wave contribute to its compression and irreversibility, they are not of equal importance. To figure out the roles of viscosity and heat conduction in shock transition, the existence of smooth shock solutions and the counter-intuitive entropy overshoot phenomenon (the specific entropy is not monotonically increasing and exhibits a peak inside the shock front) are theoretically and numerically investigated, with emphasis on the effects of viscosity and heat conduction. Instead of higher-order hydrodynamics, the Navier–Stokes formalism is employed for its stability and simplicity. Supplemented with nonlinear thermodynamically consistent constitutive relations, the Navier–Stokes equations are adequate to demonstrate the general nature of shock profiles. It is found that heat conduction cannot sustain strong shocks without the presence of viscosity, while viscosity can maintain smooth shock transition at all strengths, regardless of heat conduction. Hence, the critical role in shock compression is played by viscosity rather than heat conduction. Nevertheless, the dispensability of heat conduction would not compromise its essential role in the emergence of an entropy peak. It is the entropy flux resulting from heat conduction that neutralises the positive entropy production and thus prevents the decreasing entropy from violating the second law of thermodynamics. This mechanism of entropy overshoot has not been addressed previously in the literature and is revealed using the entropy balance equation.

Two-dimensional compressible flows in radial equilibrium are investigated in the ideal dilute-gas regime and the non-ideal single-phase regime close to the liquid–vapour saturation curve and the critical point. Radial equilibrium flows along constant-curvature streamlines are considered. All properties are therefore independent of the tangential streamwise coordinate. A differential relation for the Mach number dependency on the radius is derived for both ideal and non-ideal conditions. For ideal flows, the differential relation is integrated analytically. Assuming a constant specific heat ratio $\gamma$, the Mach number is a monotonically decreasing function of the radius of curvature for ideal flows, with $\gamma$ being the only fluid-dependent parameter. In non-ideal conditions, the Mach number profile also depends on the total thermodynamic conditions of the fluid. For high molecular complexity fluids, such as toluene or hexamethyldisiloxane, a non-monotone Mach number profile is admissible in single-phase supersonic conditions. For Bethe–Zel'dovich–Thompson fluids, non-monotone behaviour is observed in subsonic conditions. Numerical simulations of subsonic and supersonic turning flows are carried out using the streamline curvature method and the computational fluid dynamics software SU2, respectively, both confirming the flow evolution from uniform flow conditions to the radial equilibrium profile predicted by the theory.

The choking condition of a one-dimensional steady-state diabatic flow, with wall friction, of a perfect gas through a constant cross-section pipe has been analysed by applying a recent analytical solution. Such a choking condition can be achieved for an initially subsonic flow and a supersonic one, even when the heat flux goes from the fluid to the wall. Entropy–enthalpy diagrams have been analysed, and a universal choking condition for compressible diabatic flows with wall friction has been determined. The analytical solution of a supersonic flow, in which a normal shock occurs, has been obtained for a diabatic flow with wall friction and compared with the results of a numerical model.

Green hydrogen is expected to have a key role in the transition towards a carbon-neutral society, particularly in the transportation sector. Exploration of novel solutions for the direct injection of hydrogen in internal combustion engines (ICEs) is a main topic for both academia and industry. Here, the authors explore the gas-dynamic and mixing features of hydrogen under-expanded jets exiting from non-axisymmetric cross-sections, with the aim to provide guidelines for designing novel generations of ICE hydrogen injectors. Triangular and star shapes have been compared with elliptical and rectangular sections with different aspect ratios. Differences in the shock wave systems are reported, and explanations of the gas-dynamic mechanisms developed by each section are proposed. Non-uniformity in the radial expansion of the jet boundary has been noticed in all of the non-axisymmetric sections, leading to an axis switching of the jet in some cases. Differences in the radial mixing and axial jet penetration have been reported too, and a robust correlation with the vorticity distribution along the jet boundary has been observed.

Indirect noise is a significant contributor to aircraft engine noise, which needs to be minimized in the design of aircraft engines. Indirect noise is caused by the acceleration of flow inhomogeneities through a nozzle. High-fidelity simulations showed that some flow inhomogeneities can be chemically reacting when they leave the combustor and enter the nozzle (Giusti et al., Trans. ASME J. Engng Gas Turbines Power, vol. 141, issue 1, 2019). The state-of-the-art models, however, are limited to chemically non-reacting (frozen) flows. In this work, first, we propose a low-order model to predict indirect noise in nozzle flows with reacting inhomogeneities. Second, we identify the physical sources of sound, which generate indirect noise via two physical mechanisms: (i) chemical reaction generates compositional perturbations, thereby adding to compositional noise; and (ii) exothermic reaction generates entropy perturbations. Third, we numerically compute the nozzle transfer functions for different frequency ranges (Helmholtz numbers) and reaction rates (Damköhler numbers) in subsonic flows with hydrogen and methane inhomogeneities. Fourth, we extend the model to supersonic flows. We find that hydrogen inhomogeneities have a larger impact on indirect noise than methane inhomogeneities. Both the Damköhler number and the Helmholtz number markedly influence the phase and magnitude of the transmitted and reflected waves, which affect sound generation and thermoacoustic stability. This work provides a physics-based low-order model which can open new opportunities for predicting noise emissions and instabilities in aeronautical gas turbines with multi-physics flows.

We present a theoretical and numerical stability analysis for a piston-driven planar shock against two-dimensional perturbations. The results agree with the well-established theory for isolated planar shocks: in the range of $h_c< h<1+2 {\mathcal {M}}_2$, where $h$ is the D'yakov–Kontorovich (DK) parameter related to the slope of the Rankine–Hugoniot curve, $h_c$ is its critical value corresponding to the onset of the spontaneous acoustic emission (SAE) and ${\mathcal {M}}_2$ is the downstream Mach number, non-decaying oscillations of shock-front ripples occur. The effect of the piston is manifested in the presence of additional frequencies occurring by the reflection of the sonic waves on the piston surface that can reach the shock. An unstable behaviour of the shock perturbation is found to be possible when there is an external excitation source affecting the shock, whose frequency coincides with the self-induced oscillation frequency in the SAE regime, thereby being limited to the range $h_c< h<1+2 {\mathcal {M}}_2$. An unstable evolution of the shock is also observed for planar shocks restricted to one-dimensional perturbations within the range $1< h<1+2 {\mathcal {M}}_2$. Both numerical integration of the Euler equations via the method of characteristics and theoretical analysis via Laplace transform are employed to cross-validate the results.

We propose a physical model to predict indirect noise generated by the acceleration of compositional inhomogeneities in nozzles with viscous dissipation (non-isentropic nozzles). First, we derive the quasi-one-dimensional equations from the conservation laws of multicomponent flows. Second, we validate the proposed model with the experimental data available in the literature for binary mixtures of four gases. Third, we calculate the nozzle transfer functions for different Helmholtz numbers and friction factors, in both subsonic and supersonic flows with/without shock waves. We show that friction and dissipation have a significant effect on the generation of indirect noise, for which the physical mechanism is identified and explained. Fourth, we find a semi-analytical solution with path integrals, which provides an asymptotic expansion with respect to the Helmholtz number. Fifth, we introduce the compositional-noise scaling factor, which is applied to quickly estimate compositional noise from the knowledge of only one single-component gas transfer function. The approximation error is less than $1\,\%$. The proposed low-order model provides accurate estimates of the transfer functions and physical insight into indirect noise for multicomponent gases. This opens up new possibilities to accurately predict, and understand, sound generation in gas turbines.

In this paper, the reflection of curved shock waves over a symmetry plane in planar supersonic flow is studied. This includes stable Mach reflection (MR) and the regular reflection (RR) to MR transition process. Curved shock theory (CST) is applied to derive the high-order parameters in front of and behind the shock wave. The method of curved shock characteristics is used to establish an analytical model to predict the wave configurations. The shock structures provided by the proposed model agree well with the numerical results. Flow structures, such as the height of the Mach stem and the shape of the shock wave and slip line, are studied by applying the analytical model. Isentropic waves generated from a curved wall are found to significantly influence the flow patterns. It appears that the compression waves obstruct the formation of the sonic throat and increase the Mach-stem height. The expansion waves have the opposite effect. The evolution mechanism of the Mach stem is found in conjunction with the RR-to-MR transition process. The CST is extended to a moving frame and used to model the transition. The time history of the moving triple point illustrates the effects of the incident shock angle and isentropic waves on the transition process.