航空学报 > 2022, Vol. 43 Issue (1): 626005-626005   doi: 10.7527/S1000-6893.2021.26005

激波/边界层干扰机理与控制专栏

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

沈鹏飞1, 刘朋欣1, 孙东1, 袁先旭1,2   

  1. 1. 空气动力学国家重点实验室, 绵阳 621000;
    2. 中国空气动力研究与发展中心 计算空气动力所, 绵阳 621000
  • 收稿日期:2021-06-23 修回日期:2021-12-13 出版日期:2022-01-15 发布日期:2021-09-22
  • 通讯作者: 袁先旭 E-mail:yuanxianxu@cardc.cn
  • 基金资助:
    国家重点研发计划(2019YFA0405300);国家自然科学基金(11802324)

Statistical characteristics of skin friction of shock wave/turbulent boundary layer interaction in hollow cylinder-flare configuration at Mach 6

SHEN Pengfei1, LIU Pengxin1, SUN Dong1, YUAN Xianxu1,2   

  1. 1. State Key Laboratory of Aerodynamics, Mianyang 621000, China;
    2. Computational Aerodynamics Institute, China Aerodynamics Research & Development Center, Mianyang 621000, China
  • Received:2021-06-23 Revised:2021-12-13 Online:2022-01-15 Published:2021-09-22
  • Supported by:
    National Key Research and Development Program of China (2019YFA0405300); National Natural Science Foundation of China (11802324)

摘要: 为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项在干扰后呈现明显增长的原因是该项在对数区及以上区域对时均摩阻有明显增益。

关键词: 激波/湍流边界层干扰, 直接数值模拟, 非定常流动, 摩阻分解, 统计分析

Abstract: To examine the statistical characteristics of skin friction of shock wave/turbulent boundary layer interaction, the direct numerical simulation method was used to obtain the exact flow field of the shock wave/turbulent boundary layer in a hollow cylinder-flare configuration at Mach 6. The decomposition formula of averaged skin friction was derived and studied together with the convection item, streamwise inhomogeneity item, molecular viscous item, curvature effect item, and turbulent kinetic energy dissipation item. The statistical characteristics of fluctuating skin friction and averaged skin friction were explored. Probability density functionresults indicate that the fluctuating skin friction deviates from the normal distribution in the interaction region with distinct intermittency. Power Spectrum Density results show that energy of skin friction plays the main role in middle frequency zone whose peal location is 0.14 before interaction while it does in higher frequency zone after interaction. After the decomposition of the averaged skin friction, the molecular viscous item and turbulent kinetic energy dissipation item mainly contribute to the averaged skin friction. Before the interaction region, the molecular viscous item plays the main role, while after the interaction region, the turbulent kinetic energy dissipation item dominates. After the interaction, the convection item becomes stronger than before the interaction, the curvature effect item remains the same, and the streamwise inhomogeneity item turns negative from positive before the interaction, which is caused by the change in the pressure gradient item. Moreover, the fact that different regions in the boundary layer indicate different flow features leads to different contributions to the averaged skin friction, particularly for the turbulent kinetic energy dissipation item.

Key words: shock wave/boundary layer interaction, direct numerical simulation, unsteady flow, decomposition of skin friction, statistical analysis

中图分类号: