航空学报 > 2022, Vol. 43 Issue (7): 125437-125437   doi: 10.7527/S1000-6893.2021.25437

激波/湍流边界层干扰分离泡直接数值模拟

童福林1,2,3, 董思卫1,3, 段俊亦2,4, 李新亮2,4   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 绵阳 621000;
    2. 中国科学院 力学研究所 高温气体动力学国家重点实验室, 北京 100190;
    3. 中国空气动力研究与发展中心 计算空气动力研究所, 绵阳 621000;
    4. 中国科学院大学 工程科学学院, 北京 100049
  • 收稿日期:2021-03-02 修回日期:2021-04-03 发布日期:2021-04-08
  • 通讯作者: 董思卫,E-mail:dswayb@126.com E-mail:dswayb@126.com
  • 基金资助:
    国家自然科学基金(11972356,91852203);国家重点研发计划(2019YFA0405300)

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)

摘要: 采用直接数值模拟(DNS)方法对来流马赫数为2.25、33.2°激波角的入射激波/平板湍流边界层干扰分离泡进行了数值研究。在验证了计算结果可靠性的基础上,通过分析比较3个不同展向站位分离泡的非定常运动特性、分离微团几何特征和相干结构等,定量考察了三维展向结构差异的影响规律。研究发现,分离泡存在复杂的三维结构,其流向长度明显大于法向高度和展向宽度,整体上沿展向呈现中间高两边低的扁平型单峰结构。分离泡面积脉动预乘功率谱结果表明,分离泡的非定常运动表征为大尺度低频膨胀/收缩过程,其展向三维结构对峰值频率的影响较小,且分离泡两侧略滞后于中间。采用经验模态分解(EMD)方法对分离泡低频膨胀/收缩过程进行了条件统计分析。统计结果表明,膨胀和收缩运动对分离微团几何特征没有实质影响,各展向站位分离微团高度/长度比值的概率峰值出现在0.1附近,同时分离微团面积和法向高度近似满足二次方分布。此外,流向速度脉动场的本征正交分解(POD)分析指出,分离泡的非定常运动与低阶模态密切相关,而高阶模态的贡献相对较小。采用前10个低阶模态可以准确重构出分离泡的低频膨胀/收缩过程。

关键词: 激波/湍流边界层干扰, 分离泡, 低频膨胀/收缩, 经验模态分解, 本征正交分解

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

中图分类号: