高超音速零压力梯度边界层中的破裂机制与热传递超调现象研究
作者及发表信息
本研究由斯坦福大学航空航天系的Kenneth J. Franko和Sanjiva K. Lele合作完成,发表于2013年的*Journal of Fluid Mechanics*(J. Fluid Mech.)第730卷,页码491-532。论文标题为《Breakdown Mechanisms and Heat Transfer Overshoot in Hypersonic Zero Pressure Gradient Boundary Layers》,聚焦高超音速流动中边界层从层流到湍流的转捩(transition)机制及其对热传递的影响。
学术背景
高超音速流动(hypersonic flow,马赫数≥5)的边界层转捩机制与亚音速和超音速流动存在本质差异。传统转捩预测依赖线性扰动增长理论,但实验观察到转捩过程中壁面热传递和摩擦系数会先超越湍流值(即“超调”现象),随后才趋于稳定。这一现象可能导致局部高热流,威胁飞行器结构安全。然而,超调机制尚不明确,且实验受风洞噪声干扰难以精确捕捉转捩过程。因此,作者通过直接数值模拟(DNS, Direct Numerical Simulation)对比三种转捩机制(第一模态斜波破裂、第二模态斜波破裂、第二模态基频共振),旨在揭示非线性破裂过程与热传递超调的关系。
研究流程与方法
1. 模拟设置与数值方法
- 研究对象:马赫数6的平板边界层,采用等温壁条件(壁温与边界层外缘温度比Tw/Te=6.5)。
- 扰动引入:通过壁面吹吸(blowing/suction)激发三种转捩机制:
- 第一模态斜波破裂(First Mode Oblique Breakdown):强迫频率ω=0.9、展向波数β=0.3的两对斜波。
- 第二模态基频共振(Second Mode Fundamental Resonance):强迫二维第二模态(ω=2.4, β=0),并叠加低幅三维扰动。
- 第二模态斜波破裂(Second Mode Oblique Breakdown):强迫斜向第二模态(ω=2.4, β=0.1)。
- 数值算法:求解可压缩Navier-Stokes方程,采用六阶紧致格式空间离散、二阶隐式时间积分,并引入局部人工耗散(LAD, Localized Artificial Diffusivity)捕捉激波。
主要结果
1. 转捩机制对比
- 第一模态斜波破裂:转捩长度最短(约300δ*in),热传递超调显著(峰值超湍流值20%)。斜波相互作用生成强流向涡(模态(0,2)),通过二次不稳定性(secondary instability)快速破裂为湍流。
- 第二模态机制:转捩长度延长(>500δ*in),超调现象较弱。二维第二模态饱和后通过基频共振激发三维扰动,但流向涡强度较低,导致延迟破裂。
热传递超调成因
壁面冷却影响
结论与价值
1. 科学意义
- 揭示了高超音速边界层转捩中热传递超调的直接机制——流向涡增强的动量与热输运,弥补了传统线性稳定性理论的不足。
- 证实第一模态斜波破裂是超调现象的主因,为实验观测(如HIFiRE-1飞行器数据)提供了理论解释。
研究亮点
1. 方法创新:首次通过可控扰动DNS对比三种转捩机制,明确了非线性阶段涡动力学与热传递的关联。
2. 发现新颖:提出流向涡强度与展向尺度是超调的关键参数,挑战了传统间歇性(intermittency)模型的解释。
3. 数据详实:提供高分辨率湍流谱(图11)、模态能量演化(图12-19)及壁面热流分布(图25),支持结论的可重复性。
其他价值
- 研究指出,雷诺类比(Reynolds analogy)在转捩区失效(图26),需发展更精确的热流预测模型。
- 对壁面冷却效应的分析(§3.3)为高超声速飞行器主动热管理提供了理论依据。