高超音速飞行器的激波特性
高超音速飞行器（Mach > 5）在大气中飞行时，机体前缘、机翼前缘以及喷管出口等部位会产生强激波。激波是一种薄层的不可逆压缩波，跨越激波后流体的压强、密度、温度和流速都会出现突变。由于激波厚度仅为分子自由程的数量级，宏观上可视为一个几何面，常用 Rankine–Hugoniot 关系描述跨激波的状态变化。
激波的形成机制
- 超声速流动的特征：当局部马赫数 (M>1) 时，信息只能沿流向传播，扰动无法逆流。
- 几何诱导：机体前缘、凸起或喷管收缩等导致流线收敛，使局部马赫数进一步升高，最终在某一面上形成激波。
- 热力学不连续：激波面上满足质量、动量、能量守恒，但由于不可逆过程，熵会增加，从而产生温度升高和气体离解。
常见激波类型
| 类型 | 形状 | 典型出现位置 | 主要特征 |
|---|---|---|---|
| 正激波 | 平面或略呈弧形 | 机体前缘、喷管收缩端 | 压强、温度显著升高，马赫数下降至次声或低超声速 |
| 斜激波 | 斜面 | 机翼前缘、机体侧面 | 跨激波后仍保持超声速，马赫数下降幅度取决于激波倾角 |
| 螺旋激波 | 螺旋形 | 机体尾迹、喷管膨胀段 | 伴随强旋转流动，常导致局部热负荷不均匀 |
| 冲击波-膨胀波组合 | 交替出现 | 可变截面喷管 | 形成系列压缩-膨胀过程，影响推力和排气温度 |
激波对气动热负荷的影响
- 温度升高：跨正激波后温度可提升至 3000 K 以上，导致气体电离、化学反应加剧。
- 热流分布：激波位置决定热流热点，常在机体前缘形成局部最高热流密度。
- 材料选型：需要选用高温耐热材料（如碳-碳复合材料）并配合主动冷却技术。
示例：假设机体前缘马赫数为 8，入口温度 220 K，跨正激波后温度可用下式估算
[
\frac{T_2}{T_1}= \frac{(2\gamma M_1^2-(\gamma-1))((\gamma-1)M_1^2+2)}{(\gamma+1)^2M_1^2}
]
取 (\gamma=1.4)，则 (T_2\approx 2800) K。
数值仿真方法
网格划分
- 在激波可能出现的区域采用 局部细化（(y^+) < 1），保证激波捕捉精度。
- 使用 非结构化混合网格（六面体+棱柱体）兼顾几何复杂度和计算效率。
求解器选择
- 显式有限体积（如 ANSYS Fluent）适用于瞬态激波捕获。
- 隐式耦合求解（如 NASA FUN3D）在稳态高马赫数流动中更为稳健。
湍流模型
- 对于高温化学非平衡流动，推荐 SST‑k‑ω 与 化学反应模型 联合使用。
- 若需捕捉激波细节，可采用 大涡模拟（LES） 或 混合 RANS‑LES（DES）。
# 示例：Fluent 中的激波捕获设置
mesh refine region -name "ShockZone" -size 0.001
solver set scheme second-order-upwind
physics set model compressible-viscous
turbulence model sst-k-omega
chemistry enable
solve initialize
solve iterate 5000
实验验证与工程应用
- 冲击管实验：利用高压气体瞬时释放产生可控激波，测量压强、温度和光谱，验证数值模型的化学反应率。
- 风洞试验：在 高超音速风洞（Mach 6–10）中装配热流计和红外摄像，获取机体表面热流分布。
- 飞行试验：如 X‑37B、HTV‑2 项目，通过机载传感器记录激波位置和热负荷，实现模型的全尺度验证。
小结
高超音速飞行器的激波特性是决定气动热负荷、结构材料和制导控制的核心因素。通过 激波形成机理、激波类型、热负荷影响 的系统分析，结合 高分辨率数值仿真 与 实验验证，可以实现对激波行为的精准预测，为高超音速飞行器的设计提供可靠的技术支撑。未来的研究方向包括更精细的化学非平衡模型、基于机器学习的激波捕捉算法以及新型高温材料的协同优化。