导航

ACTA AERONAUTICAET ASTRONAUTICA SINICA ›› 2022, Vol. 43 ›› Issue (7): 125437-125437.doi: 10.7527/S1000-6893.2021.25437

• Fluid Mechanics and Flight Mechanics • Previous Articles     Next Articles

Direct numerical simulation of separation bubble in shock wave/turbulent boundary layer interaction

TONG Fulin1,2,3, DONG Siwei1,3, DUAN Junyi2,4, LI Xinliang2,4   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    2. State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China;
    3. Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    4. School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China
  • Received:2021-03-02 Revised:2021-04-03 Published:2021-04-08
  • Supported by:
    National Natural Science Foundation of China (11972356, 91852203); National Key Research and Development Program of China (2019YFA0405300)

Abstract: Characteristics of separation bubbles in the interaction of a supersonic turbulent boundary layer at Mach number 2.25 with an impinging shock wave of 33.2° are investigated by means of Direct Numerical Simulation (DNS). After verifying the reliability of the numerical results, fundamental mechanisms associated with separation bubbles, including unsteadiness, separation micro-clusters geometries features and coherent structures, at three different spanwise locations are quantitatively compared to analyze the influence of the three-dimensionality in the spanwise direction. It is found that the separation bubble is highly three-dimensional, with the streamwise extent significantly larger than the wall-normal height and spanwise width. In the spanwise direction, the bubble height is generally large in the middle and small on both sides, exhibiting a single flat peak behavior. The pre-multiplied power spectrum density of the fluctuating separation bubble area suggests that the separation bubble unsteadiness is characterized by large-scale low-frequency contraction and dilation, which is less affected by the spanwise three-dimensionality. The bubbles on both sides lag slightly behind that in the middle. Conditional analysis based on Empirical Mode Decomposition (EMD) is performed to analyze the influence of the bubble dilation and contraction on geometries features of the separation micro-clusters. The statistical results indicate no essential changes in both motions, where the probability peak of the aspect ratio appears around 0.1, and the area and the normal height of the micro-clusters approximately satisfy the quadratic distribution. In addition, the Proper Orthogonal Decomposition (POD) analysis of the fluctuating streamwise velocity indicates that the unsteady motion of the separation bubble is strongly related to the low-order modes, whereas the contribution from the high-order modes is rather small. With the first ten low-order modes, the low-frequency dilation and contraction process of separation bubbles is accurately reconstructed.

Key words: shock wave/turbulent boundary layer interaction, separation bubbles, low-frequency dilation and contraction, empirical mode decomposition, proper orthogonal decompositionhttp

CLC Number: