We report the characteristics of wall shear stress (WSS) and wall heat flux (WHF) from direct numerical simulation (DNS) of a spatially developing zero-pressure-gradient supersonic turbulent boundary layer at a free-stream Mach number M∞ = 2.25 and a Reynolds number Reτ = 769 with a cold-wall thermal condition (a ratio of wall temperature to recovery temperature Tw/Tr = 0.75). A comparative analysis is performed on statistical data, including fluctuation intensity, probability density function, frequency spectra, and space–time correlation. The root mean square fluctuations of the WHF exhibit a logarithmic dependence on Reτ similar to that for the WSS, the main difference being a larger constant. Unlike the WSS, the probability density function of the WHF does not follow a lognormal distribution. The results suggest that the WHF contains more energy in the higher frequencies and propagates downstream faster than the WSS. A detailed conditional analysis comparing the flow structures responsible for extreme positive and negative fluctuation events of the WSS and WHF is performed for the first time, to the best of our knowledge. The conditioned results for the WSS exhibit closer structural similarities with the incompressible DNS analysis documented by Pan and Kwon [“Extremely high wall-shear stress events in a turbulent boundary layer,” J. Phys.: Conf. Ser. 1001, 012004 (2018)] and Guerrero et al. [“Extreme wall shear stress events in turbulent pipe flows: Spatial characteristics of coherent motions,” J. Fluid Mech. 904, A18 (2020)]. Importantly, the conditionally averaged flow fields of the WHF exhibit a different mechanism, where the extreme positive and negative events are generated by a characteristic two-layer structure of temperature fluctuations under the action of a strong Q4 event or a pair of strong oblique vortices. Nevertheless, we use the bi-dimensional empirical decomposition method to split the fluctuating velocity and temperature structures into four different modes with specific spanwise length scales, and we quantify their influence on the mean WSS and WHF generation. It is shown that the mean WSS is mainly related to small-scale structures in the near-wall region, whereas the mean WHF is associated with the combined action of near-wall small-scale structures and large-scale structures in the logarithmic and outer regions.

No abstract

No abstract

Direct numerical simulations of a supersonic turbulent boundary layer on a flat plate interacting with an impinging shock wave are carried out with two different incident shock angles at Mach 2.25. The effect of the interaction strength on the recovery process in the downstream region is systematically studied, including the turbulence evolution, the statistical and structural properties of wall pressure fluctuations, and the generation of mean skin friction and wall heat flux. The variations of the Reynolds stress components, the anisotropy tensor, and the turbulent kinetic energy budget in the two flow cases highlight a slow reversal tendency and an increasingly pronounced importance of the outer-layer large-scale structures in the relaxation region of the strong interaction. We find that the effect of increasing the interaction strength on the fluctuating wall pressure is reflected by a decrease in the characteristic frequencies, an increase in the spatial extent, and a decrease in the convection velocity. We decompose the mean skin friction and wall heat flux into different physically informed contributions and reveal that the mean wall heat flux generation is the same regardless of the interaction strength; in contrast, the generation mechanism of mean skin friction is found to be fundamentally changed. A novel scale-decomposition method is used to quantify the effect of the increased interaction strength on the leading components, and it is demonstrated that the energetic outer-layer large-scale structures are the dominant contributor in the recovery process as the interaction strength is increased.

Direct numerical simulations are carried out to identify the effects of shock impingement on the behavior of bump flow at free-stream Mach number of 2.25. Two cosine-shaped bump cases, with and without an impinging oblique shock at an angle of 33.2{degree sign}, are compared. The shock impingement exhibits a remarkable influence on the pattern of the shock system and on the size of the separation region. Spectral analysis finds that low-frequency unsteadiness is significantly enhanced by the impingement interaction, and the proper orthogonal decomposition highlights the low-frequency breathing motion of the separation bubble, which is accurately reconstructed using only the first ten low-order modes. Downstream of the bump, the Reynolds stress components and the turbulence kinetic energy both exhibit a general amplification, with the peaks reoccurring at outer wall-normal locations. Turbulent kinetic energy budget analysis shows the greatly increased production in the outer layer which is balanced by turbulent transport and dissipation. Anisotropy-invariant map analysis identifies enhanced isotropic turbulence in the vicinity of the bump, which is qualitatively modified into a two-component axisymmetric state around the reattachment point. In addition, the mean skin friction decomposition suggests that the shock impingement has little influence on the predominant contribution of turbulence kinetic energy production, apart from the spatial growth dominance at the bump summit in the absence of the impinging shock. Interestingly, a scale-decomposed analysis quantitatively demonstrates that the contributions of small-scale structures are attenuated, but those of large-scale ones are relatively increased, with a contribution of more than 80% with shock impingement.