极超音速流
极超音速流的理论基础
概述
老师，极超音速流是指马赫数5以上的流体吗？普通超音速有什么不同呢？
决定性的区别在于冲击波后方的温度达到数千K，气体的化学性质发生变化。空气发生离解和电离，理想气体的假设完全破裂。再入舱和超燃冲压发动机的设计中无法回避这个领域。
温度太高导致空气分子被破坏，这太危险了…
是的。N₂和O₂离解成原子，更高温时会电离。因此需要考虑热化学非平衡的控制方程。
控制方程
具体会出现什么方程？
首先是垂直冲击波的基本关系，Rankine-Hugoniot关系式。对于理想气体，冲击波前后的压力比为
温度比为
当M=20时温度比大约是多少？
用理想气体计算的话T₂/T₁会是几百倍，但实际上离解反应会吸收能量，温度上升会被抑制。这正是实际气体效应。牛顿修正压力系数在极超音速中也很重要
是钝头物体滞点压力的很好近似。此外滞点焓
是TPS（热防护系统）设计的出发点。在马赫25的轨道速度下，$h_0 \approx 30$ MJ/kg。
30 MJ/kg是个巨大的能量…
热化学非平衡模型
空气离解时用什么模型？
需要添加5种（N₂, O₂, NO, N, O）或11种（包括离子种和电子）的化学种输运方程。每个化学种 $Y_s$ 的输运方程为
其中 $\dot{\omega}_s$ 是化学反应源项，用Arrhenius反应速率常数
计算。Park (1990) 的2温度模型中，将平动-旋转温度 $T_{tr}$ 和振动-电子温度 $T_{ve}$ 分开处理。
用两个温度吗！我想象冲击波后面振动温度滞后上升。
完全正确。振动松弛时间用Millikan-White相关式估计。这个非平衡区域对冲击波的驻止距离和壁面加热率有很大影响，所以CFD中适当的建模是必不可少的。
壁面加热与耐热设计
再入舱的表面加热程度有多大？
滞点加热率可用Fay-Riddell公式粗略估算。
鼻尖半径越大加热率越低，所以再入体采用钝头形状。Apollo型舱在峰值加热率时可达几MW/m²。
所以再入体才是圆形的。尖头的话加热会集中。
是的。这就是极超音速气动设计的基本思想，也称为钝体悖论。
航天飞机再入——马赫25的空气变成等离子体
极超音速流（马赫5以上）最大的难点是"空气不再是理想气体"。航天飞机大气圈再入时，飞行器前方的空气在相当于马赫25的冲击波中瞬间加热到6000K以上。在这个温度下，氮气分子（N₂）离解成N原子，进一步电离形成等离子体。要在CFD中精确计算，需要在通常的流动方程基础上加入化学反应方程——"反应气体模型"。哥伦比亚号事故的热防护分析中，这种化学反应CFD在事故原因调查中发挥了重要作用。
极超音速流的数值计算方法
数值格式的选择
极超音速流的CFD能用普通的压缩性求解器计算吗？
基本的有限体积法（FVM）框架是一样的，但通量计算格式的选择至关重要。要稳定捕获强冲击波，标准做法是用Roe法或AUSM+系列的格式配合限制器（Van Leer, Barth-Jespersen等）。
AUSM+是什么缩写？
Advection Upstream Splitting Method。将质量通量和压力通量分离评估的方法，能有效降低极超音速域的carbuncle现象。特别是AUSM+-up格式在全速度范围内稳定性很高。
空间离散化与冲击波捕获
在网格上处理冲击波很重要吧？
完全同意。在冲击波捕获法（shock capturing）中，数值粘性将不连续面平滑到数个网格单元宽度。为提高精度，用MUSCL重建达到2阶精度，同时用TVD（Total Variation Diminishing）限制器抑制振荡。
这里 $\phi(r)$ 是限制器函数。MinMod限制器最富有扩散性但稳定，Superbee很陡峭但有超调的风险。
极超音速下限制器的选择会很敏感吗？
非常敏感。在M>10的强弓形冲击波中，Roe法单独使用容易产生carbuncle不稳定。H-CUSP法或Rotated Roe格式等考虑多维效应的方法很有效。
时间积分与显式方法的约束
时间积分用显式还是隐式方法？
求定常解用隐式法（LU-SGS、DPLR型point-implicit等）更高效。CFL数可取100以上，收敛很快。而非定常现象（舱体振荡等）需用2阶隐式法（BDF2）或显式二阶Runge-Kutta。不过显式法的CFL限制为
极超音速下音速 $a$ 因高温增大，时间步长会非常小。
化学反应源的刚性问题
我听说包含化学反应的计算很难收敛。
化学反应的时间尺度比流体时间尺度小好几个数量级，导致源项"stiff"。解决办法是用point-implicit法隐式处理化学种源项的雅可比矩阵，或用operator splitting法分离流体和化学。NASA的DPLR码和US3D码标准使用point-implicit法。
DPLR是什么？
Data Parallel Line Relaxation。NASA Ames开发的极超音速CFD码，内置Park热化学模型。是再入体解析的行业标准之一。此外还有LAURA（NASA Langley）和US3D（明尼苏达大学）比较著名。
商业求解器不支持吗？
ANSYS Fluent也支持实际气体和有限速度化学反应，但极超音速专用码的验证实绩更多。STAR-CCM+有Reacting Flow模型，但5种和11种空气模型的实现可能需要用户侧定制。
极超音速CFD"网格就是生命"的真实原因
极超音速CFD中最关键的是网格设计。冲击波厚度达到平均自由程（分子间距离）的量级，接近连续体假设的极限。同时边界层厚度比超音速情况薄得多。必须同时解析多个薄薄的层。实际中使用"冲击波拟合（shock fitting）"技术——将冲击波位置显式地纳入网格——能大幅减少冲击波后方的数值扩散。结构格子与非结构格子的混合方法正成为现场标准。
极超音速流的实际应用
分析工作流
极超音速CFD分析具体怎么进行？
典型工作流如下。
1. 形状定义：准备再入体或波浪形飞行器的CAD形状。鼻尖圆角半径必须设为有限值
2. 网格生成：用结构网格覆盖包含弓形冲击波的区域。壁面第一层y⁺<1
3. 边界条件：自由流（M, T∞, p∞）、壁面（等温壁或辐射平衡壁温）、出口（超音速流出）
4. 物理模型：5种空气或11种空气、Park 2温度模型、催化壁条件的选择
5. 求解：用CFL逐步增加（初始CFL=0.1到某个值）开始计算
6. 后处理：检验壁面加热率、压力分布、冲击波形状
CFL逐步增加是什么意思？
即使隐式法，突然用很大的CFL数也会因为初期强非线性而发散。所以从CFL=0.1开始，随着残差下降逐步增加到100或1000。Fluent有自动调整CFL数的选项。
网格设计要点
网格最需要注意什么？
有几个重要要点。
| 区域 | 推荐单元数 | 原因 |
|---|---|---|
| 冲击波层（墙面到冲击波） | 50-80层 | 冲击波形状与加热率精度 |
| 壁面边界层 | 30-50层（y⁺<1） | 壁面热流量的准确预测 |
| 冲击波对齐网格 | shock-aligned为理想 | 数值扩散最小化 |
| 鼻尖周围 | 高密度集中 | 精确捕获最大加热率 |
怎样做冲击波对齐网格？
事先用牛顿法或MOC（特性曲线法）估算冲击波形状，将网格外边界与该形状对齐。NASA的GRIDGEN2D或Pointwise的椭圆型网格生成功能可用。非结构网格的情况下，可用AMR（自适应网格细化）检测冲击波并自动细化。
壁面边界条件的设置
"催化壁"是什么意思？
极超音速中壁面附近离解的原子是否在壁面上重组，对加热率有很大影响。完全催化壁（FCW）下所有原子在壁面上重组并释放能量，加热率最大。非催化壁（NCW）下不重组，加热率最小。实际的TPS材料介于二者之间，用催化效率 $\gamma_{cat}$ 作为参数。
设计应该用哪一种？
保守的设计用FCW。SiC表面或RCC（强化碳复合材）的实际催化效率在NASA TP中有公布数据。CFD验证时最好两种都计算，确认灵敏度。
验证基准
计算结果用什么数据验证？
代表性基准如下。
- FIRE II飞行实验：1965年再入舱。壁面加热率飞行数据公开
- RAM-C II：电子密度飞行实测数据。用于电离等离子体的验证
- 双圆锥/双楔：CUBRC（Calspan）冲击波隧道实验数据
- HIFiRE飞行实验：极超音速边界层转变数据
必须与实验数据对比。仅凭计算不能保证可信度。
完全同意。极超音速CFD物理模型的不确定性很大，所以V&V（Verification & Validation）特别重要。
TPS（热防护系统）设计时CFD工程师的困境
极超音速舱的设计中最严峻的是"怎样确定热防护系统（TPS）的厚度"。用CFD计算的热流量乘以安全系数决定材料厚度，但实际气体模型的不确定性±20～30%，导致安全系数异常之大。结果是要装载沉重的隔热材，挤压有效载荷。实际方法是"细致对标有飞行历史的再入体实测数据与CFD，掌握自己代码的系统偏差，从而缩小不确定性"。虽然是地道的工作，却能实现TPS的轻量化。
极超音速流的软件比较
极超音速支持工具比较
极超音速流解析用什么软件？
极超音速涉及化学反应、高温气体，通常的CFD求解器不够用。对比专用码和通用码。
| 工具 | 种类 | 化学反应模型 | 热非平衡 | 催化壁 | 备注 |
|---|---|---|---|---|---|
| DPLR (NASA Ames) | 专用 | 5/11/13种空气 | Park 2温度 | 支持 | 再入体行业标准 |
| LAURA (NASA Langley) | 专用 | 多种支持 | 多温度 | 支持 | Orion/Mars Entry设计用 |
| US3D (U. Minnesota) | 专用 | 任意化学系 | 多温度 | 支持 | 非结构网格支持 |
| ANSYS Fluent | 通用 | 有限速度反应 | UDF扩展 | UDF | 物种运输模型支持 |
| STAR-CCM+ | 通用 | 化学反应 | 有限 | 定制 | 反应流模块 |
| OpenFOAM (hy2Foam) | 开源 | 5/11种空气 | Park 2温度 | 支持 | 斯特拉思克莱德大学开发 |
| Eilmer (U. Queensland) | 开源 | 多种 | 多温度 | 支持 | Python+Lua界面 |
NASA专用码很强。但获取很难吧？
DPLR和LAURA是NASA内部码，不公开。但US3D有大学授权，hy2Foam在GitHub上公开。用商业求解器做极超音速的话，多数用Fluent的实际气体模型+UDF实现Park模型。
ANSYS Fluent的设置示例
Fluent做极超音速的注意点？
主要要点整理如下。
- 求解器：必须用Density-Based求解器（Pressure-Based不适用）
- 通量格式：Roe-FDS或AUSM
- 气体模型：物种运输中定义N₂, O₂, NO, N, O五种
- 反应机制：有限速率化学，输入各反应的Arrhenius系数
- 运输系数：动理论或Sutherland公式（高温区推荐Blottner曲线拟合）
- CFL控制：初期0.5 → 逐步增加直至收敛
Fluent能处理热非平衡吗？
标准功能只有单温度模型。2温度模型需在UDF中添加能量方程。这要花不少开发时间，有时hy2Foam会更便捷。
hy2Foam（OpenFOAM）的应用
hy2Foam能实用吗？
由斯特拉思克莱德大学的Vincent Casseau博士开发，基于OpenFOAM，已实现Park 2温度模型和多种空气模型。有FIRE II和RAM-C的验证论文。用于研究足够实用，源代码也是学习热化学模型实现的好教材。
开源就能做极超音速，太吸引人了。
但要理解没有图形界面，全是文本配置。并行计算可扩展性和网格规模实绩不如商业码。
选择极超音速解析软件时要确认"热化学反应数据库"
选商用工具做极超音速CFD时，GUI的美观度不重要——重要是"热化学反应数据库包含的物种数量和精度"。马赫10以上时氮气、氧气多阶段离解，产生一氧化氮、电子等。你的求解器只有5种空气模型，还是支持11种（N₂, O₂, N, O, NO, NO⁺, N₂⁺, O₂⁺, N⁺, O⁺, e⁻），这直接影响热荷重计算精度。特别是TPS设计中，这个选择关系到安全系数，不容妥协。
极超音速流的前沿研究
极超音速边界层转变
极超音速研究现在最热的课题是什么？
最大课题之一是边界层转变。层流到湍流的转变位置影响壁面加热率3～8倍，直接关系到TPS重量。Mack第2模式不稳定性在极超音速占主导，这是声波在边界层内捕获的不稳定。
Mack第2模式是什么？
马赫4以上边界层中出现的高频不稳定模式。与墙面附近相对绝对压力梯度导致的第1模式不同，第2模式因为边界层相对超音速区域的声波反射和放大而增长。典型频率几百kHz，波长约为边界层厚度的两倍。
CFD直接预测很难吧？
RANS无法直接预测转变，要用eN方法或DNS/LES。用NASA的eN码（LASTRAC）或PSE（抛物线型稳定方程）做线性稳定分析，N因子达到临界值（通常N=5～10）的位置为转变点。最近也有用DNS完整计算非线性阶段的研究。
烧蚀耦合分析
TPS材料融化过程也用CFD计算吗？
是的。再入体用PICA（酚醛浸渍碳烧蚀剂）或SiC等烧蚀材。材料热分解气体从表面喷出（吹气效应），推开边界层，降低加热率。这种耦合分析需要CFD和材料响应码（FIAT、TITAN等）的耦联。
表面形状变化，网格也要更新。
对。用ALE（任意拉格朗日欧拉）法或overset grid法处理表面后退。NASA的CHAR（炭化烧蚀响应）码与DPLR耦联过。
机器学习的应用
极超音速也用AI吗？
研究很活跃。主要应用有三块。
- 代理模型：从形状参数（鼻尖半径、外倾角等）瞬时预测加热率的神经网络
- 湍流模型修正：用DNS数据训练的机器学习模型修正RANS的Reynolds应力张量
- 边界层转变预测：用实验数据库学习N因子和转变位置的分类器
DNS计算量太大，用机器学习代替很合理。
但外推性（超出训练数据范围的泛化）是课题。极超音速实验数据有限，所以物理约束神经网络（PINN）来自支配方程的方法受到关注。
极超音速流的问题处理
常见问题与解决方案
极超音速CFD计算上不去，原因是什么？
极超音速特有的问题模式整理如下。
1. Carbuncle不稳定
听说过carbuncle。弓形冲击波崩溃的现象吧？
Roe法计算钝头物体时，滞点附近的弓形冲击波非物理地变成凹凸形。格子对称线上容易发生。
对策：
- 切换到AUSM+或HLLC格式（比Roe法耐性强）
- 如用Roe法，添加H修正（entropy fix）
- 在冲击波法向方向加微小人工粘性
- 避开格子对称线（引入微弱的非对称性）
2. 负温度负密度的出现
温度或密度变成负数，什么情况？
强膨胀波或冲击波干涉附近，从守恒量转换为原始变量时产生非物理值。特别有化学反应时，质量分数会出现<0或>1。
对策：
- 降低CFL数（从0.1以下开始）
- 先用1阶精度得到初值，收敛后切换到2阶
- 添加质量分数剪裁处理（$\max(0, \min(1, Y_s))$）
- 使用保正限制器（Positivity-preserving limiter）
3. 壁面加热率的网格依赖性
改网格加热率就变了…
这是极超音速CFD最常见的问题。壁面加热率正比于温度梯度 $\partial T / \partial n |_{wall}$，强烈依赖壁面附近的网格分辨率。
对策：
- 壁面法向至少50～80层网格
- 第一层高度 $y^+ < 1$（极超音速边界层很薄，可能达微米级）
- 3水准以上网格收敛性验证（Richardson外推最佳）
- 冲击波与壁面间网格均匀分布
4. 化学反应发散
加入化学反应计算就发散，怎么办？
通常是化学源项的刚性引起。
对策：
- 用operator splitting分离流体和化学，化学物种用stiff ODE求解器（CVODE等）
- Point-implicit法隐式处理源项雅可比矩阵
- 先计算冻结流（frozen flow）收敛，然后逐步启动反应
- 设温度下限（100K以下反应速率式可能失效）
5. 冲击波驻止距离不匹配
冲击波位置与实验不符？
冲击波驻止距离对气体模型很敏感。理想气体模型驻止距离偏大，化学平衡模型偏小，常常非平衡模型最接近实验。
检查清单：
- 气体模型（理想/平衡/非平衡）选择恰当？
- 自由流条件（M, T, p）正确？
- 反应速率常数数据源合理？（Park 1990 vs. Gupta 1990等）
- 网格充分解析冲击波？（冲击波宽度内至少5个单元）
极超音速是物理模型的选择直接决定结果，问题处理也要回归物理。
完全同意。不仅数值方法调试，还要问"这个条件下真的需要非平衡效应吗""壁面催化条件合理吗"这样的物理判断。这就是极超音速CFD难的地方，也是有意思的地方。
极超音速计算的"carbuncle现象"——冲击波变成怪物的数值bug
极超音速CFD的臭名昭著的故障是"carbuncle（硬结肿）现象"。强弓形冲击波前缘附近，现实中不存在的非物理突起凭空长出，用Roe格式或Godunov格式常发生。名字来自"carbuncle（脓肿）"——长得跟这个一样可怕的bug。对付办法包括切换到HLLE系有适度数值粘性的格式，或用冲击波检测器做局部格式切换。极超音速代码的开发和使用要领，这个症状是必须早早铭记的。
有帮助
更多
错误