This study investigates the correlation between the fluctuating wall heat flux, and the distribution and transport of Reynolds shear stress and turbulent heat flux in compressible boundary layers at Mach number 5.86 and friction Reynolds number 420, with a relatively weaker and a stronger wall cooling imposed. As illustrated from the probability density functions of the wall-heat-flux perturbations, with increasing wall cooling, the extreme wall heat flux is intensified and tends to be more negatively skewed. To examine the role of the extreme events in the transport of the momentum and heat, conditional analysis of the extreme positive and negative wall-heat-flux-perturbation events is conducted. In most regions of the boundary layer, the positive events are predominantly associated with an increase in Reynolds shear stress and a decrease in turbulent heat flux. Joint probability density functions of velocity and wall-heat-flux perturbations in the near-wall region indicate that the extreme positive events tend to be more correlated with ejections, which is particularly evident in the stronger wall-cooling case. To further shed light on the underlying mechanisms of the connections between wall heat flux and transport budgets, a transport equation for turbulent heat flux is derived, in a similar manner to that for Reynolds shear stress. The energy balance is inspected, with conditional analysis applied to budget terms and mean flow properties so as to quantify the correlation between wall-heat-flux fluctuations and energy evolution.

In convergent geometry, the effect of convergence and compression on the Rayleigh–Taylor instability (RTI) and Richtmyer–Meshkov instability (RMI) modifies the growth rate and behaviour of the instabilities. In order to better understand how compression/expansion caused by axial strain rates (i.e. strain rates normal to the interface) change the instability dynamics, axial strain rates are applied to RMI in planar geometry, isolated from the effects of convergence. Potential flow theory for the linear regime shows the growth rate of the instability is modified to include the background velocity difference of the instability's width. Resolved two-dimensional simulations of single-mode RMI showed the potential flow model is accurate whilst the amplitude is small compared with the wavelength. The application of strain rate to an RMI-induced mixing layer was investigated using three-dimensional implicit large eddy simulations (ILES) of the quarter-scale $\theta$-group case by Thornber et al. (Phys. Fluids, vol. 29, 2017, 105107). Whilst the background strain rate contributed to the mixing layer's growth, it was to a smaller extent than expected. The shear production of axial turbulent kinetic energy from the strain rate modified the rate of bulk entrainment, affecting the mixing layer's growth and mixedness, such that the strained simulations no longer attained the same self-similar state. The capability of the buoyancy-drag model by Youngs & Thornber (Physica D, vol. 410, 2020, 132517) to predict the mixing layer width was investigated, using a model calibrated to the unstrained case. New terms were introduced into the buoyancy-drag model, which correspond to the shear production of turbulent kinetic energy.

The influence of outer large-scale motions (LSMs) on near-wall structures in compressible turbulent channel flows is investigated. To separate the compressibility effects, velocity fluctuations are decomposed into solenoidal and dilatational components using the Helmholtz decomposition method. Solenoidal velocity fluctuations manifest as near-wall streaks and outer large-scale structures. The spanwise drifting of near-wall solenoidal streaks is found to be driven by the outer LSMs, while LSMs have a trivial influence on the spanwise density of solenoidal streaks, consistent with the outer LSM impacts found in incompressible flows (Zhou et al., J. Fluid Mech., vol. 940, 2022, p. A23). Dilatational motions are characterized by the near-wall small-scale travelling-wave packets and the large-scale parts in the outer region. The streamwise advection velocity of the near-wall structures remains at $16 \sim 18u_{\tau }$, hardly influenced by Mach numbers, Reynolds numbers and wall temperatures. The spanwise drifting of near-wall dilatational structures, quantified by the particle image velocimetry method, follows a mechanism distinct from solenoidal streaks. This drifting velocity is notably larger than those of the solenoidal streaks, and the influence of outer LSMs is not the primary trigger for this drifting.

We systematically study the dissipative anomaly in compressible magnetohydrodynamic (MHD) turbulence using direct numerical simulations, and show that the total dissipation remains finite as viscosity diminishes. The dimensionless dissipation rate $\mathcal {C}_{\varepsilon }$ fits well with the model $\mathcal {C}_{\varepsilon } = \mathcal {C}_{\varepsilon,\infty } + \mathcal {D}/R_L^-$ for all levels of flow compressibility considered here, where $R_L^-$ is the generalized large-scale Reynolds number. The asymptotic value $\mathcal {C}_{\varepsilon,\infty }$ describes the total energy transfer flux, and decreases with increase of the flow compressibility, indicating non-universality of the dimensionless dissipation rate in compressible MHD turbulence. After introducing an empirically modified dissipation rate, the data from compressible cases collapse to a form similar to the incompressible MHD case depending only on the modified Reynolds number.

In the present study, we perform direct numerical simulations of compressible turbulent boundary layers at free stream Mach numbers $2\unicode{x2013}6$ laden with dilute phase of spherical particles to investigate the Mach number effects on particle transport and dynamics. Most of the phenomena observed and well-recognized for inertia particles in incompressible wall-bounded turbulent flows – such as near-wall preferential accumulation and clustering beneath low-speed streaks, flatter mean velocity profiles, and trend variation of the particle velocity fluctuations – are identified in the compressible turbulent boundary layer as well. However, we find that the compressibility effects are significant for large inertia particles. As the Mach number increases, the near-wall accumulation and the small-scale clustering are alleviated, which is probably caused by the variations of the fluid density and viscosity that are crucial to particle dynamics. This can be affected by the fact that the forces acting on the particles with viscous Stokes number greater than 500 are modulated by the comparatively high particle Mach numbers in the near-wall region. This is also the reason for the abatement of the streamwise particle velocity fluctuation intensities with the Mach numbers.

The turbulent boundary layer (TBL) is a widely existing flow phenomenon in nature and engineering applications. Its strong mixing effect can achieve more sufficient material mixing, heat transport, etc. The understanding of the entrainment process and mechanism of irrotational fluids entering the turbulent region can be promoted by studying the geometric and dynamic characteristics of turbulent${/}$non-turbulent interfaces (TNTI). In compressible flow, it is unclear whether the properties of TNTI will change and whether the entrainment will show different features due to the influence of compressibility. Based on the direct numerical simulation results of supersonic compressible plate TBLs with Mach number of 2.9, the geometric and dynamic characteristics of TNTI are investigated in this paper. The interface is identified by the enstrophy method, and the height, thickness, fractal dimension, enstrophy transportation and entrainment characteristics of the interface are investigated. It is found that for the enstrophy transportation in a TBL, the contribution of compressibility-related terms accounts for approximately 13.4 % of the total enstrophy transportation, which tends to transfer the enstrophy of turbulence near the interface to both directions vertical to the interface. This promotes the expansion of the turbulent region towards the non-turbulent region, and the mean height, thickness and entrainment velocity are increased by approximately 3.7 %, 7.0 % and 8.5 %, respectively, while the fractal dimension is basically unaffected. Different from the incompressible flow, the contribution of the compressibility-related terms to the entrainment velocity is independent of the local curvature, and the intense entrainment process is more likely to occur on a highly curved concave surface.

Analysing the dynamics of velocity gradients is useful for understanding various nonlinear turbulence processes. This work focuses on how the vibrational non-equilibrium of the constituent molecules in a gaseous medium affects the dynamics of the velocity gradient and the pressure Hessian tensors. We first derive the exact evolution equation of the pressure-Hessian tensor in the presence of the vibrational non-equilibrium process. Subsequently, we perform several direct numerical simulations of compressible isotropic turbulence, including the vibrational relaxation process therein. Using flow fields extracted from these simulations, we conduct several parametric studies over a range of the Damköhler number (ratio of the relevant fluid time scale to that of the mean vibrational relaxation process) and the initial ratio of the vibrational temperature to the mean local temperature. We find that a variation in the initial Damköhler number does influence the evolution of the pressure-Hessian and the velocity gradient tensors. As the vibrational relaxation process becomes more rapid (an increase in the value of the initial Damköhler number), it causes a decrease in the strength of the pressure-Hessian tensor and simultaneous suppression of dilatational fluctuations in the flow field. On the other hand, a variation in the initial value of the ratio of the vibrational temperature to the local temperature does not seem to affect the pressure-Hessian or the velocity gradient tensor. These findings are expected to aid in the development of closure models for the pressure-Hessian tensor in compressible flows under vibrational non-equilibrium conditions.

This work investigates the compressible turbulence induced by Richtmyer–Meshkov (RM) instability using high-resolution Navier–Stokes simulations. Special attention is paid to the characteristics of RM turbulence including the mixing width growth, the turbulent kinetic energy (TKE) decay, the mixing degree, inhomogeneity and anisotropy. Three distinct initial perturbation spectra are designed at the interface to reveal the initial condition imprint on the RM turbulence. Results show that cases with initial large-scale perturbations present a stronger imprint on statistical characteristics and also a quicker growth of mixing width, whereas cases with small-scale perturbations present a faster TKE decay, greater mixing level, higher isotropy and homogeneity. A thorough analysis on the inter-scale energy transfer in RM turbulence is also presented with the coarse-graining approach that exposes the two subfilter-scale (SFS) energy fluxes (i.e. deformation work and baropycnal work). A strong correlation between the nonlinear model of baropycnal work (Fluids, 4(2), 2019) and the simulation results is confirmed for the first time, demonstrating its potential in modelling the RM turbulence. Two primary mechanisms of baropycnal work (the straining and baroclinic generation processes) are explored with this nonlinear model. The evolutions of two SFS energy fluxes exhibit distinct behaviours at various filter scales, in different flow regions and under various flow motions (strain and rotation). It is found that all three cases share the common inter-scale energy transfer dynamics, which is important for modelling the RM turbulence.

Direct numerical simulations of hypersonic turbulent Couette flows are performed for top-wall Mach numbers of 6, 7 and 8, inspired by non-reactive high-enthalpy wind tunnel free-stream conditions, with the goal of analysing the physical processes driving the sub-filter-scale (SFS) stresses to inform development of large-eddy simulation techniques for hypersonic wall-bounded flows. Semi-local scaling laws collapse mean profiles and second-order turbulent statistics well, in spite of the strong wall-normal gradients of temperature and density. On the other hand, the SFS shear stresses exhibit an unexpected profile characterized by a region of pronounced shear stress deficit, which becomes more pronounced for higher Mach numbers. Instantaneous visualizations suggest that the SFS shear stress deficit is induced by the counter-gradient resolved momentum transport driven by residual velocity motions at the interface between high-density low-speed streaks being ejected away from the wall, and low-density high-speed ones replacing the displaced fluid, qualifying this as a compressibility effect. It is shown that the SFS shear stresses are primarily driven by second-order interactions between residual velocities, in spite of their triply nonlinear nature. This, in turn, motivated a statistical quadrant analysis revealing the presence of a SFS shear stress deficit, that is, SFS processes driving momentum transport of the resolved field towards the top wall. Additionally, higher-Reynolds-number simulations reveal that there is an upper limit of spatial filter widths to resolve large-scale structures, and such deficit is observable at any Reynolds number when a reasonable spatial filter is applied.

This paper extends the work of Tamano & Morinishi (J. Fluid Mech., vol. 548, 2006, pp. 361–373) by simulating supersonic turbulent channel flow with asymmetric thermal walls using a larger computational domain and a finer mesh. Direct numerical simulation is carried out for four cases with different thermal wall boundaries at the top wall at fixed $Ma=1.5$, $Re=6000$ and $Pr=0.72$, while the bottom wall is maintained at a constant temperature of $T_L$ equal to the reference temperature. These cases are referred to as the adiabatic case TAd, where the top wall is adiabatic; the pseudo-adiabatic case T32, where the top wall is isothermal with temperature $T_{w,t}=T_A$; the sub-adiabatic case T25, with $T_{w,t}=0.77T_A$; and the super-adiabatic case T40, with $T_{w,t}=1.24T_A$. Here, $T_A=3.234$ is the mean temperature at the adiabatic wall in the TAd case. The objective of this study is to compare and contrast the TAd case with its corresponding T32 case, and to investigate the effect of the wall temperature difference between the two isothermal walls. Comparisons of the basic turbulent statistics, the heat transfer between the Favre-averaged mean-flow kinetic energy, the Favre-averaged turbulent kinetic energy and the Favre-averaged mean internal energy, as well as the wall heat transfer properties, indicate that the TAd case and its corresponding T32 case are generally equivalent. The only discernible difference is in the region very close to the top wall for the temperature-fluctuation-related quantities. The analysis reveals that the asymmetry of the thermal walls causes asymmetry in the flow and thermal fields. In addition, the transfer of the heat generated by the pressure dilatation and the viscous stress is facilitated by the turbulent heat flux term and the mean molecular heat flux term.

An improved resolvent analysis is proposed in the regime of compressible turbulent boundary layers. To better model nonlinear processes in the input, the resolvent framework is augmented by adding eddy viscosity. To this end, we propose two eddy-viscosity models: a modified Cess eddy-viscosity model coupling the compressibility transformation and outer-layer correction, and a new eddy-viscosity model based on an empirical relationship and mixing-length theory. Both are incorporated into the resolvent operator to examine the performance of the eddy-viscosity-improved resolvent-based reduced-order modelling. Results of the augmented resolvent analysis are compared qualitatively and quantitatively with the first leading mode of spectral proper orthogonal decomposition, by checking the profiles and cross-spectral densities of velocities, density and temperature in two hypersonic turbulent boundary layers under different wall conditions. Higher accuracy of the turbulence prediction is achieved by adding the proposed eddy-viscosity models, particularly for the energetic cycle in the outer-layer region where strong nonlinear energy transfer exists.

Direct numerical simulations of temporally developing mixing layers have been performed to investigate the effects of compressibility on statistics and structures near the interfaces of high- and low-speed large-scale structures (LSSs), covering a range of convective Mach numbers from $M_c=0.2$ to $1.8$ and Taylor Reynolds numbers up to 290. The interfaces of LSSs are directly defined by the isosurface of zero fluctuating streamwise velocity. The characteristic velocity jump at the interfaces grows rapidly in the transition stage and then decreases until reaching a gradual self-similar stage. The evolution of velocity jump is evidently delayed as the convective Mach number increases. The interface layer is formed by the shearing motion of neighbouring LSSs. A small vortical motion characterized by the Kolmogorov scale is induced in the interface layer by shear-dominated outer regions. The strengths of outer shearing motion and central vortical motion are greater at the leading edge, but smaller at the trailing edge, which is also reflected in a larger velocity jump at the leading edge. As the convective Mach number increases, the small-scale vortical structure is obviously suppressed by compressibility. At high convective Mach number $M_c=1.8$, the compressive shear-dominated outer regions are linked with a sheet-like structure passing through the centre of the expansion region corresponding to a small-scale vortical structure. The compressibility and shearing strength near the interface are highly dependent on the interface orientation.

Very-large-scale motions are commonly observed in moderate- and high-Reynolds-number wall turbulence, constituting a considerable portion of the Reynolds stress and skin friction. This study aims to investigate the behaviour of these motions in high-speed and high-Reynolds-number turbulent boundary layers at varying Mach numbers. With the aid of high-precision numerical simulations, numerical experiments and theoretical analysis, it is demonstrated that the very-large-scale motions are weakened in high-Mach-number turbulence at the same friction Reynolds numbers, leading to the reduction in turbulent kinetic energy in the outer region. Conversely, the lower wall temperature enhances the very-large-scale motions but shortens the scale separation between the structures in the near-wall and outer regions.

The paper studies the statistics of total and static temperature (total $h_t$ and static $h$ enthalpy) in compressible turbulent plane channel flow using direct numerical simulations (DNS) data covering the range of centreline Mach numbers $0.3\lessapprox \bar {M}_{{{CL}}_x} \lessapprox 2.5$ and Huang–Coleman–Bradshaw friction Reynolds numbers $100\lessapprox Re_{\tau ^\star }\lessapprox 1000$. For this class of very-cold-wall flows, the DNS data for correlation coefficients and joint probability density functions (p.d.f.s) show that $h_t'$ is invariably very strongly correlated with the streamwise velocity fluctuation $u'$, in contrast to static temperature (static enthalpy $h'$) whose correlation with $u'$ weakens rapidly with increasing wall distance. We study various correlations and joint p.d.f.s of $h_t'$ and $h'$ in relation to the fluctuating velocity field, including the turbulent Prandtl number $Pr_{{T}}$, and discuss the predictions of Reynolds analogy. The scaling of the mean enthalpy and the fluctuating enthalpy variance and fluxes with respect to inner and outer velocity scales is investigated. The complex behaviour and scaling of different terms in the transport equations for the enthalpy variance and fluxes are discussed.

Linear models, based on stochastically forced linearized equations, are deployed for spectral linear stochastic estimation (SLSE) of the velocity and temperature fluctuations in compressible turbulent channel flows with a bulk Mach number of 1.5. Through comparison with the direct numerical simulation (DNS) data, an eddy-viscosity-enhanced model (eLNS) outperforms the one not enhanced (LNS) in computing the coherence and amplitude ratio of streamwise velocity at different wall-normal heights, but they both largely deviate from DNS regarding the temperature prediction. For further investigation, the eigenspectra and pseudospectra of the linear operators are scrutinized. The eddy viscosity is shown to stabilize the eigenmodes and decrease the non-normality of the vortical modes. Consequently, the relative importance of acoustic and entropy modes increases, and they can contribute 20 % to 55 % of the response growth, which is not supported by DNS. Hence, it is an intrinsic defect of the eLNS model introduced by turbulence modelling. After a procedure of cospectrum decomposition, the contributions of acoustic and entropy components are filtered out. The resulting SLSE quantities for velocity, temperature and their coupling are basically agreeable with DNS, demonstrating that the coherent temperature fluctuation is dominated by advection and other vortical motions, instead of the compressibility effects. Moreover, a parameter study of Reynolds and Mach numbers (from 0.3 to 4) is conducted. The semi-local units are shown to well collapse the velocity SLSE quantities to the incompressible case for streamwise-elongated structures of high coherence.

In the present study, we investigate the compressibility effects in supersonic and hypersonic turbulent boundary layers under the influence of wall disturbances by exploiting direct numerical simulation databases at Mach numbers up to 6. Such wall disturbances enforce extra Reynolds shear stress on the wall and induce mean streamline curvature in rough wall turbulence that leads to the intensification of turbulent motions in the outer region. The turbulent and fluctuating Mach numbers, the density and the velocity divergence fluctuation intensities suggest that the compressibility effects are enhanced by the increment of the free-stream Mach number and the implementation of the wall disturbances. The differences between the Reynolds and Favre average due to the density fluctuations constitute approximately $9\,\%$ of the mean velocity close to the wall and $30\,\%$ of the Reynolds stress near the edge of the boundary layer, indicating their non-negligibility in turbulent modelling strategies. The comparatively strong compressive events behaving as eddy shocklets are observed at the free-stream Mach number of $6$ only in the cases with wall disturbances. By further splitting the velocity into the solenoidal and dilatational components with the Helmholtz decomposition, we found that the dilatational motions are organized as travelling wave packets in the wall-parallel planes close to the wall and as forward inclined structures in the form of radiated waves in the vertical planes. Despite their increased magnitudes and higher portion in the Reynolds normal and shear stresses, the dilatational motions show no tendency of contributing significantly to the skin friction and the production of turbulent kinetic energy due to their mitigation by the cross-correlation between the solenoidal and dilatational velocity components.

Buoyancy-driven turbulent convection leads to a fully compressible flow with a prominent top-down asymmetry of first- and second-order statistics when the adiabatic equilibrium profiles of temperature, density and pressure change very strongly across the convection layer. The growth of this asymmetry and the formation of an increasingly thicker stabilized sublayer with a slightly negative mean convective heat flux $J_c(z)$ at the top of the convection zone is reported here by a series of highly resolved three-dimensional direct numerical simulations beyond the Oberbeck–Boussinesq and anelastic limits for dimensionless dissipation numbers, $0.1 \le D\le 0.8$, at fixed Rayleigh number $Ra=10^6$ and superadiabaticity $\epsilon =0.1$. The highly stratified compressible convection regime appears for $D > D_{crit}\approx 0.65$, when density fluctuations collapse to those of pressure; it is characterized by an up to nearly 50 % reduced global turbulent heat transfer and a sparse network of focused thin and sheet-like thermal plumes falling through the top sublayer deep into the bulk.

Direct numerical simulation of a compressible round jet is carried out at Mach number of 0.9 and Reynolds number of 3600 and the data are used to perform velocity gradient tensor (VGT) analysis for different regions of the spatially developing jet. For the developed portion of the jet, the classical teardrop shape is observed for the joint probability density function (p.d.f.) of Q and R (second and third invariants of the VGT). In the region just after the potential core, between $X = 10$ and 15 $r_0$ ($r_0$ is the jet inlet radius), an inclination towards the third quadrant is observed in the Q–R joint p.d.f. which represents the presence of tube-like structures. It is also shown that this inclination in the turbulent/non-turbulent (T/NT) boundary and interface towards the third quadrant is mainly a contribution of points that lie in regions with negative dilatation. Regions with weak expansion also show this third quadrant inclination to some extent. Points that lie in regions with relatively higher positive dilatation show no such inclination towards the third quadrant but are inclined towards the fourth quadrant which indicates the presence of sheet-like structures. Similarly for the domain segment $X = 15$ to 20 $r_0$, it is observed that points that lie in the regions with positive dilatation have a joint p.d.f. with an inclination towards fourth quadrant, which suggests the presence of sheet-like structures at the T/NT boundary and interface. Points that lie in regions with negative dilatation show the appearance of a third quadrant lobe.

We present pore-resolved compressible direct numerical simulations of turbulent flows grazing over perforated plates, that closely resemble the acoustic liners found on aircraft engines. Our direct numerical simulations explore a large parameter space including the effects of porosity, thickness and viscous-scaled diameter of the perforated plates, at friction Reynolds numbers $\textit {Re}_\tau = 500\unicode{x2013}2000$, which allows us to develop a robust theory for estimating the added drag induced by acoustic liners. We find that acoustic liners can be regarded as porous surfaces with a wall-normal permeability and that the relevant length scale characterizing their added drag is the inverse of the wall-normal Forchheimer coefficient. Unlike other types of porous surfaces featuring Darcian velocities inside the pores, the flow inside the orifices of acoustic liners is fully turbulent, with a magnitude of the wall-normal velocity fluctuations comparable to the peak in the near-wall cycle. We provide clear evidence of a fully rough regime for acoustic liners, also confirmed by the increasing relevance of pressure drag. Once the fully rough asymptote is reached, canonical acoustic liners provide an added drag comparable to that of sand-grain roughness with viscous-scaled height matching the inverse of the viscous-scaled Forchheimer permeability of the plate.

Compressible turbulent plane Couette flows are studied via direct numerical simulation for wall Reynolds numbers up to $Re_w=10\ 000$ and wall Mach numbers up to $M_w=5$. Various turbulence statistics are compared with their incompressible counterparts at comparable semilocal Reynolds numbers $Re^*_{\tau,c}$. The skin friction coefficient $C_f$, which decreases with $Re_w$, only weakly depends on $M_w$. On the other hand, the thermodynamic properties (mean temperature, density and others) strongly vary with $M_w$. Under proper scaling transformations, the mean velocity profiles for the compressible and incompressible cases collapse well and show a logarithmic region with the Kárman constant $\kappa =0.41$. Compared with wall units, the semilocal units yield a better collapse for the profiles of the Reynolds stresses. While the wall-normal and spanwise Reynolds stress components slightly decrease in the near-wall region, the inner peak of the streamwise component notably increases with increasing $M_w$ – indicating that flow becomes more anisotropic when compressible. In addition, the near-wall turbulence production decreases as $M_w$ increases – due to rapid wall-normal changes of viscosity caused by viscous heating. The streamwise and spanwise energy spectra show that the length scale of near-wall coherent structures does not vary with $M_w$ in semilocal units. Consistent with those in incompressible flows, the superstructures (the large-scale streamwise rollers) with a typical spanwise scale of $\lambda _z/h\approx 1.5{\rm \pi}$ become stronger with increasing $Re_w$. For the highest $Re_w$ studied, they contribute about $40\,\%$ of the Reynolds shear stress at the channel centre. Interestingly, flow visualization and correlation analysis show that the streamwise coherence of these structures degrades with increasing $M_w$. In addition, at comparable $Re^*_{\tau,c}$, the amplitude modulation of these structures on the near-wall small scales is quite similar between incompressible and compressible cases – but much stronger than that in plane Poiseuille flows.