The interaction between a weakly turbulent free stream and a hypersonic shock wave is investigated theoretically by using linear interaction analysis (LIA). The formulation is developed in the limit in which the thickness of the thermochemical nonequilibrium region downstream of the shock, where relaxation toward vibrational and chemical equilibrium occurs, is assumed to be much smaller than the characteristic size of the shock wrinkles caused by turbulence. Modified Rankine–Hugoniot jump conditions that account for dissociation and vibrational excitation are derived and employed in a Fourier analysis of a shock interacting with three-dimensional isotropic vortical disturbances. This provides the modal structure of the post-shock gas arising from the interaction, along with integral formulas for the amplification of enstrophy, concentration variance, turbulent kinetic energy (TKE), and turbulence intensity across the shock. In addition to confirming known endothermic effects of dissociation and vibrational excitation in decreasing the mean post-shock temperature and velocity, these LIA results indicate that the enstrophy, anisotropy, intensity, and TKE of the fluctuations are much more amplified through the shock than in the thermochemically frozen case. In addition, the turbulent Reynolds number is amplified across the shock at hypersonic Mach numbers in the presence of dissociation and vibrational excitation, as opposed to the attenuation observed in the thermochemically frozen case. These results suggest that turbulence may persist and get augmented across hypersonic shock waves despite the high post-shock temperatures.
We propose a linear asymptotic theory to describe the propagation of planar detonation fronts through heterogeneous mixtures of reactive gases consisting of random fluctuations in the fuel mass fraction. The analysis starts with the derivation of the transfer functions that relate the upstream fuel mass fraction inhomogeneities with the burnt-gas perturbations via normal mode analysis. These results are then used in a Fourier analysis of a detonation wave interacting with two- and three-dimensional isotropic heterogeneous fields. This yields integral formulae for the turbulent kinetic energy, sonic energy and averaged vorticity and entropy production rates. Second-order corrections for the turbulent Rankine–Hugoniot conditions are also obtained, along with analytical expressions for the deviation of the detonation velocity with respect to that of the equivalent homogeneous mixture. Upstream inhomogeneities are found to speed up the detonation front in the vast majority of scenarios studied, with a velocity amplification factor that depends on the properties of the fuel–air mixture, particularly on the variation of the density and the heat release with the fuel mass fraction.
Hypersonic flight involves a variety of complex flow phenomena that directly impact the aerothermodynamic loading of high-speed vehicles. The turbulence encountered during a typical flight trajectory influences and interacts with the shock waves on and around the surface of a vehicle and its propulsion system, affecting both aerodynamic and power plant performance. These interactions can be studied by isolating a turbulent flow convected through a normal shock, commonly referred to as the canonical shock-turbulence interaction (STI) problem. Scale-resolving computational fluid dynamics (CFD) and linear interaction analysis (LIA) have been crucial in studying this problem and formulating scaling laws that explain the observed behavior. In this work, an extensive review of the theoretical (LIA) and numerical (CFD) work on the canonical STI is presented. The majority of the work conducted to date has focused on calorically perfect gases with constant heat capacities. However, in hypersonic flows, chemical and thermal non-equilibrium effects may alter the nature of the interaction. As a result, relevant LIA and CFD studies addressing high-enthalpy phenomena are also succinctly discussed.
The dynamic response of structural elements subjected to blast loading is a problem of growing interest in the field of defense and security. In this work, a novel computational tool for the rapid evaluation of the effects of explosions, hereafter referred to as SimEx, is presented and discussed. The classical correlations for the reference chemical (1 kg of TNT) and nuclear ( kg of TNT) explosions, both spherical and hemispherical, are used together with the blast wave scaling laws and the International Standard Atmosphere (ISA) to compute the dynamic response of Single-Degree-of-Freedom (SDOF) systems subject to blast loading. The underlying simplifications in the analysis of the structural response follow the directives established by UFC 3-340-02 and the Protective Design Center Technical Reports of the US Army Corps of Engineers. This offers useful estimates with a low computational cost that enable in particular the computation of damage diagrams in the Charge Weight–Standoff distance (CW–S) space for the rapid screening of component (or building) damage levels. SimEx is a computer application based on Matlab and developed by the Fluid Mechanics Research Group at University Carlos III of Madrid (UC3M). It has been successfully used for both teaching and research purposes in the Degree in Security Engineering, taught to the future Guardia Civil officers at the Spanish University Center of the Civil Guard (CUGC). This dual use has allowed the development of the application well beyond its initial objective, testing on one hand the implemented capacities by undergraduate cadets with the end-user profile, and implementing new functionalities and utilities by Masters and PhD students. With this experience, the application has been continuously growing since its initial inception in 2014 both at a visual and a functional level, including new effects in the propagation of the blast waves, such as clearing and confinement, and incorporating new calculation assistants, such as those for the thermochemical analysis of explosive mixtures; crater formation; fragment mass distributions, ejection speeds and ballistic trajectories; and the statistical evaluation of damage to people due to overpressure, body projection, and fragment injuries.
View Video Presentation: https://doi.org/10.2514/6.2023-0075.vid The interaction between a weakly turbulent free stream of air and a hypersonic shock wave is investigated theoretically by using linear interaction analysis (LIA). The perturbation-free jump conditions across the shock are computed using Combustion Toolbox, an in-house thermochemical code capable of capturing high-temperature phenomena such as dissociation, ionization, and recombination in multi-species mixtures, which are found to be dominant effects in hypersonic shocks in air. The formulation is developed in the limit in which the thickness of the thermochemical nonequilibrium region is assumed to be much smaller than the characteristic size of the shock wrinkles caused by turbulence. Similarly to our previous work [Physics of Fluids, 33(8), 086111 (2021)] which only accounted for vibrational and dissociation effects of single-species diatomic gases, the LIA results presented here for air indicate that the enstrophy, anisotropy, intensity, and turbulent kinetic energy (TKE) of the fluctuations are more amplified through the shock than in the thermochemical frozen case. Moreover, the turbulent Reynolds number is also amplified across the shock at hypersonic Mach numbers in the presence of dissociation and vibrational excitation, as opposed to the attenuation observed in the thermochemical frozen case. Multi-species effects reshape the TKE curve by rendering two maxima that fit fairly well within the O2 and N2 dissociation processes.
19 December 2025 in Computer Physics Communications