Direct Numerical Simulation Research for Supersonic Shock Wave/Turbulent Boundary Layer Interactions

Abstract

激波/边界层干扰问题是可压缩湍流领域中一个基本的流动现象，具有十分重要的工程应用背景。干扰区内的流动非常复杂，其中包括边界层大尺度的分离与再附过程、激波低频振荡现象、大尺度涡结构的生成以及强压力脉动载荷等。本文采用直接数值模拟对来流马赫数Ma = 2.25，33.2&deg; 激波角的入射激波/平板湍流边界层干扰流动的分离特性、再附边界层结构以及摩擦阻力的生成演化机制进行了研究。此外，还对来流马赫数Ma = 2.9，24&deg; 偏转角的压缩-膨胀折角构型中的复杂流动进行了直接数值模拟，探讨了膨胀角高度变化对离心不稳定性以及干扰区的非定常运动特性的影响规律，还对壁面压力脉动场的空间分布结构进行了分析。主要的工作和研究成果如下： (1) 入射激波与平板湍流边界层干扰流动的分离泡结构具有三维空间分布特征，整体上呈现扁平的单峰造型。通过定量对比三个不同展向站位的分析结果，研究了展向三维结构对分离区非定常运动特性以及瞬态分离微团的几何形态的影响规律。研究发现，分离泡的非定常特性表征为低频的、大尺度的膨胀和收缩运动。壁面压力脉动的功率谱密度结果表明，展向三维结构对非定常运动的特征频率影响较小。通过经验模态分解方法对分离泡面积脉动信号进行了低频重构，从中提取出了膨胀和收缩过程，借此条件统计分析分离泡呼吸运动对瞬态分离微团几何形状的影响。此外，应用本征正交模态分解对流向速度场进行了低阶近似，分析发现分离泡呼吸运动与低阶模态密切相关，前十个低阶模态重构出的流场能够精确复现分离泡低频的、大尺度呼吸运动。再附边界层内湍流运动的统计特性研究发现，边界层外层的相干结构强度明显增加、具有较大的空间尺度。壁面剪切应力脉动与流向速度脉动的时-空相关性分析表明，涡结构倾角在再附边界层的恢复过程中逐渐减小。通过空间两点的振幅调制相关性，对边界层内、外层结构间的尺度干扰效应进行了分析，提出了一种可能的机制来解释激波干扰区下游的强调制作用。 (2) 考虑激波入射角分别为33.2&deg;&nbsp;和28&deg;&nbsp;的两种情况，对比研究激波干扰强度对壁面摩擦阻力生成、演化特性的影响规律。研究发现，摩擦阻力在流动产生分离时的恢复速率较为缓慢；而边界层在保持附着的情况下，在一个边界层厚度的流向范围内快速恢复到来流充分发展湍流边界层的水平。本研究将最新提出的平均摩擦阻力分解方法应用在激波干扰区内，定量地考察摩擦阻力在不同激波入射角情况下的演化规律。研究发现，与充分发展湍流边界层不同，再附边界层内湍流相关项的贡献显著增强，而压力梯度引起的空间发展项在很大程度上抑制了摩擦阻力的产生。此外，通过二维经验模态分解技术对湍流脉动实现尺度分解，定量评估不同展向空间尺度的湍流运动对摩擦阻力生成的贡献。结果表明，在分离区和再附区来自外层大尺度结构的贡献占主导地位，而内层小尺度结构的贡献相当有限。然而在弱干扰的情况下，外层大尺度结构的贡献显著减少，导致外层大尺度和内层小尺度运动的贡献相当。 (3) 对压缩-膨胀折角干扰构型中的激波/湍流边界层干扰问题进行了直接数值模拟研究。考虑了膨胀角高度为 / = 4.25, 1.22 两种情况。当膨胀角高度较高时，流动规律与传统的压缩折角构型一致；而对于高度较低的情况，激波干扰区受下游膨胀角的存在影响较大，分离区的长度显著减小。壁面压力脉动的功率谱密度结果表明，干扰区非定常运动的低频成分被抑制。极限流线确认了G&ouml;rtler涡在再附区附近存在。通过对边界层内流线的曲率半径和G&ouml;rtler 数分布分析发现，离心不稳定性在压缩角附近区域得到保留，但在膨胀角附近区域消失。膨胀折角下游的流动结构相对复杂，可以观察到附加的激波串结构和新的流向涡生成。此外，本研究还在膨胀角上、下游流动之间发现了信号之间的负反馈机制。</p