喷嘴流
喷嘴流的理论基础
概述
老师，拉瓦尔喷嘴是用收缩-膨胀形状来产生超音速的，对吧？为什么这种形状能超过音速？
这是个好问题。可压缩流体力学的基本定理，断面积变化与速度变化的关系为
表示。亚音速（M<1）时缩小面积会加速，超音速（M>1）时扩大面积会加速。因此通过达到音速的喉部（最小截面积）可以过渡到超音速。
喉部恰好达到M=1。这就是壅塞条件吗？
完全正确。壅塞（choking）是指喉部马赫数达到1且质量流量最大的状态。这个临界质量流量由以下公式确定：
无论背压降低多少，通过喉部的质量流量都不会超过这个值。
截面积-马赫数关系
从喷嘴截面积怎样求出马赫数？
在等熵流假设下，截面积比与马赫数的关系为
对于给定的 $A/A^*$，存在亚音速解和超音速解两个解。哪个解会实现取决于背压条件。
这是非线性方程，无法解析求解吧？
可以用牛顿-拉夫逊法轻松求解。将初始值设在亚音速侧（M<1）或超音速侧（M>1）中的任一方，就会收敛到对应的解。结合等熵关系
可得温度比和压力比。对于空气（$\gamma=1.4$），当M=2时，$T_0/T=1.8$，$p_0/p=7.824$。
喷嘴内的激波
背压与设计值不同时会发生什么？
背压过高时，在膨胀部会产生垂直激波。激波前后的压力比由Rankine-Hugoniot关系确定
激波后流动变为亚音速，然后减速至喷嘴出口。继续降低背压会使激波向喷嘴出口移动，最终形成喷嘴外的膨胀波-压缩波模式。
火箭发动机排气中看到的钻石图案就是这个吗？
正是。过膨胀（出口压 < 背压）时产生斜激波，不足膨胀（出口压 > 背压）时产生Prandtl-Meyer膨胀波。那个马赫盘的位置由背压比确定。
拉瓦尔喷嘴的发明——收缩膨胀"节流"的发现拯救了蒸汽涡轮机
收缩膨胀喷嘴（de Laval nozzle）由瑞典工程师Gustav de Laval在1880年代设计。当时他在改进蒸汽涡轮机效率，通过实验发现要将蒸汽加速到超音速，需要在喉部先达到音速，然后再膨胀。虽然当时可压缩气体理论尚不完善，但他通过反复试验最终找到了这个形状。现代火箭发动机和喷气发动机的喷嘴都遵循这一原理。这是理论前于实验发现"答案"的典范之一。
喷嘴流的数值计算方法
准一维喷嘴的数值解法
准一维喷嘴流是CFD入门的最佳选择，对吧？
完全同意。这是求解一维Euler方程加截面积变化项的问题，可以实验可压缩CFD的所有基本方案。保存形式为
右侧的 $p \, dA/dx$ 作为源项发挥作用。
一般用什么格式来求解？
从教育角度看，MacCormack法（预测-修正两阶段显式格式）是经典做法，这也是Anderson教科书的标准问题。实际应用中用Roe法或AUSM计算通量，用MUSCL重构实现二阶精度。
二维三维喷嘴的网格策略
二维或三维喷嘴分析如何生成网格？
轴对称喷嘴用二维轴对称网格就足够了。建议用结构化网格，喉部附近要加密。
| 区域 | 网格策略 | 理由 |
|---|---|---|
| 收缩部 | 壁面法向至少20层 | 边界层分辨率 |
| 喉部附近 | 轴向高密度 | 精确捕捉音速条件 |
| 膨胀部 | 激波位置自适应加密 | 清晰捕捉激波 |
| 中心轴 | 避免奇异性（楔形单元等） | 轴对称计算数值稳定 |
三维还有推力偏转喷嘴（TVC）之类的，对吧？
推力偏转喷嘴和多喷嘴阵列必须用三维计算。喷嘴壁面曲率大的部分要用棱柱层（inflation layer）分辨边界层，中心用多面体或六面体填充。
边界条件的设定
喷嘴计算的边界条件怎样设置？
这是喷嘴CFD的核心部分。
- 入口：固定全温 $T_0$ 和全压 $p_0$（燃烧室条件）。流向为轴向
- 出口：超音速流出时所有变量外推（信息不向上游传播）。存在亚音速域时指定背压 $p_b$
- 壁面：绝热壁（adiabatic）或等温壁。无滑移条件
- 对称轴：轴对称条件或对称面条件
Fluent中用"Pressure Inlet"设置 $p_0, T_0$，用"Pressure Outlet"设置背压。超音速出口时背压值实际上被忽略。
亚音速-超音速过渡的数值处理
喉部的M=1过渡在数值上会造成问题吗？
很好的观察。Euler方程的双曲性在M=1处改变，某些数值格式在过渡点附近可能出问题。需要Roe法的熵修正（entropy fix）或Harten-Hyman修正。AUSM系列本身就对这个问题有很好的设计。
压力基础求解器不能解超音速吗？
最新的压力基础耦合求解器（如Fluent的coupled scheme）支持全速度范围，但对含激波的喷嘴流，密度基础求解器更稳定、精度更高。特别是要精确预测喷嘴内激波的位置和强度，非得用密度基础不可。
喷嘴数值解法中"密度基础还是压力基础"的困扰由来
对拉瓦尔喷嘴的超音速区域进行数值计算时，初学者经常困于"密度基础求解器与压力基础求解器选哪个"的问题。当亚音速到超音速连续求解喷嘴流时，由于马赫数跨越1，压力基础求解器的Presto!和邻近修正会失效。一般而言，对"超音速、可压缩性占优"的问题，密度基础更稳定，尤其采用时间前进（显式密度基础）时激波捕捉特别清晰。但壁面附近的Y+管理不当也会导致边界层计算崩溃，因此要从网格和求解方法的整体考虑。
喷嘴流的实际应用
火箭喷嘴的设计分析
火箭喷嘴设计CFD是怎样的流程？
看看实际的流程吧。
1. 等熵设计：用MOC（特性线法）确定最短长度喷嘴或钟形喷嘴的壁面轮廓
2. 冷态CFD：先用理想气体验证流场。检查激波或喉部附近的流动模式
3. 高温气体CFD：考虑燃烧气体的实在气体效应（离解·再结合）的计算
4. 壁面冷却：再生冷却通道或膜冷却的传热分析
5. 结构耦合：热负荷和压力负荷的结构应力分析
用MOC（特性线法）确定喷嘴形状。
指定设计点的马赫数（例如M=3.0）和喉部半径，追踪特性线的交点，就能得到能在出口实现均匀平行流的喷嘴壁面形状。Anderson《现代可压缩流》有详细的算法说明。
背压变动的参数化分析
改变背压时喷嘴内流动如何变化？
根据背压比 $p_b/p_0$，会出现以下工作阶段。
| $p_b/p_0$ | 状态 | 流动特征 |
|---|---|---|
| $= 1$ | 无流动 | 等压 |
| $>$ p*/p_0 | 全域亚音速 | 仅有文丘里效应 |
| $=$ p*/p_0 | 喉部M=1 | 壅塞开始 |
| 中间值 | 喷嘴内垂直激波 | 激波位置依赖背压 |
| 设计值 | 完全膨胀 | 出口达设计马赫数 |
| 低于设计值 | 过膨胀/不足膨胀 | 喷嘴外的波模式 |
其中 $p^*/p_0 = (2/(\gamma+1))^{\gamma/(\gamma-1)}$，对空气约为0.528。
CFD参数化扫描背压，激波在喷嘴内移动的样子应该很清楚吧。
完全同意。Fluent可以创建设计点（Design Points），用背压作参数进行批量计算，然后用ParaView制作动画，从教育和可视化角度都非常有效。
三维效应与非对称性
实际喷嘴有时不是轴对称的吧？
过膨胀喷嘴会出现"流动分离"——激波从壁面剥离。这个分离如果非对称化，就产生横向力（side load），这是喷嘴结构设计的大问题。SSME（航天飞机主发动机）的设计中也是个重大课题。
非对称剥离的计算需要三维吗？
验证要点
喷嘴CFD结果怎样验证？
要确认以下几点。
- 质量流量：与理论值 $\dot{m}_{max}$ 偏差在1%以内
- 喉部压力：确认 $p^*/p_0 \approx 0.528$（空气）
- 出口马赫数：与截面积比计算的理论值一致
- 壁面压力分布：与实验数据或等熵解对比
- 推力：用推力方程 $F = \dot{m} V_e + (p_e - p_a) A_e$ 验证
用"手工计算"优化喷嘴形状的年代
当代喷嘴轮廓设计往往结合CFD和优化算法，但在1950-60年代，设计人员用"特性线法（Method of Characteristics）"手工计算来确定喷嘴形状。在方格纸上一条一条绘制特性线（马赫线），逐步推进以求出实现等熵膨胀的壁面轮廓。NASA初期火箭发动机喷嘴就这样设计的。现在教科书仍然介绍这个方法，深入理解它就能直观看出"CFD为什么要这样优化形状"的道理。
喷嘴流的软件比较
喷嘴流分析工具
喷嘴流分析用什么软件比较合适？
喷嘴流是可压缩CFD的典型问题，主要商用和开源工具都有对应功能。但功能完整度有差异。
| 工具 | 密度基础 | 实在气体 | 化学反应 | MOC设计 | 适用场景 |
|---|---|---|---|---|---|
| Ansys Fluent | Density-Based | NIST RGP | Species Transport | 无 | 通用喷嘴CFD |
| Ansys CFX | Coupled | 实在气体EOS | 反应流 | 无 | 涡轮机械喷嘴 |
| STAR-CCM+ | Coupled Flow | 表格/EOS | 详细化学 | 无 | 燃烧喷嘴 |
| OpenFOAM (rhoCentralFoam) | 密度基础 | 有限 | reactingFoam | 无 | 研究·教育 |
| NPSS (NASA) | 0D/1D | 各种 | CEA耦联 | 有 | 发动机循环分析 |
| NASA CEA | 0D | 多种 | 化学平衡 | 仅输出 | 燃烧气体组成计算 |
NASA CEA是什么？
Chemical Equilibrium with Applications。输入燃料、氧化剂及其混合比，就能计算化学平衡状态的气体组成和热力性质。火箭喷嘴的入口条件设定时必须用，Web界面免费公开。
Ansys Fluent中的喷嘴设置
用Fluent做喷嘴计算的具体设置请告诉我。
基本设置如下。
- 求解器：Density-Based, Steady
- 空间离散化：2nd Order Upwind（初期收敛困难时从1st Order开始）
- 通量类型：Roe-FDS（清晰捕捉激波）
- 气体模型：Ideal Gas（冷态）/ Real Gas Property Table（高温气体）
- 乱流：SST k-omega模型（注重壁面近处精度）
- 入口BC：Pressure-Inlet（指定 $p_0, T_0$）
- 出口BC：Pressure-Outlet（指定背压。超音速出口设Supersonic/Initial Gauge Pressure）
超音速出口时背压值没有意义，对吧？
完全同意。整个出口面都是超音速时，所有变量都从内部外推，背压被忽略。但过渡计算中部分亚音速或背压会生效的情况需要注意。
OpenFOAM中的喷嘴计算
OpenFOAM用什么求解器？
rhoCentralFoam是最适合的。这是基于KNP（Kurganov-Noelle-Petrova）格式的密度基础求解器，对含激波的喷嘴流稳定计算。设置例子如下。
thermophysicalProperties:hePsiThermo+perfectGasfvSchemes:localEulerddt（定常）、vanLeerdiv（二阶TVD）fvSolution: smoothSolver + GaussSeidel0/: p, T, U的初值和边界条件
边界条件：
- inlet:
totalPressure+totalTemperature - outlet:
waveTransmissive（非反射边界） - walls:
zeroGradient（绝热壁）
waveTransmissive边界是什么？
基于特性线的非反射边界条件，能阻止波从计算域出口反射。对喷嘴出口避免人为反射波很重要。
各种工具验证喷嘴流时最初检验的基准算例
对喷嘴流计算工具的性能评估，CFD业界多年来都用"一维等熵喷嘴流"作基准。因为存在解析解，可以直接比较计算值和严格解。改变压力比使喷嘴内产生垂直激波的条件也有解析解，可以精确到kPa级验证激波位置和强度。FLUENT或OpenFOAM等商用代码开发团队，每次新增格式时都要回归测试该喷嘴算例。导入工具时首先验证喷嘴流是工程师的标准做法。
喷嘴流的前沿研究
双铃形喷嘴与自适应喷嘴
喷嘴设计的最新趋势是什么？
随着上升高度变化大气压，喷嘴膨胀比动态优化的"高度补偿喷嘴"受关注。主要概念有三种。
- 双铃形喷嘴：内外两层铃形结构。低高度时在内铃形末端流动分离，高高度时膨胀到外铃形
- 塞式喷嘴（航空尖峰喷嘴）：中心塞周围膨胀。外侧自动与大气压适应
- 伸展-偏转喷嘴（E-D喷嘴）：省略壁面一部分以获得高度补偿效果
SpaceX的Raptor发动机也优化了喷嘴吧？
两相流喷嘴
固体火箭的喷嘴里还有固体粒子混在里面吧？
完全正确。固体推进剂燃烧气体含有Al₂O₃粒子（直径1-10μm），导致比推力降低和喷嘴侵蚀。这种两相流喷嘴计算采用欧拉-拉格朗日方法。
- 连续相（气体）：Euler方程
- 分散相（粒子）：Lagrangian追踪。粒子-气体间阻力、热传递耦合
Fluent的DPM（离散相模型）或STAR-CCM+的Lagrangian Multiphase均可处理。粒子通过喉部时的速度滞后（velocity lag）导致推力损失。
推力偏转控制（TVC）的CFD
推力偏转CFD怎样预测？
推力偏转方式有多种，CFD处理方法也各不相同。
| TVC方式 | CFD方法 | 难题 |
|---|---|---|
| 万向节 | 网格变形/重叠网格 | 可动部件网格 |
| 二次流体喷射（SITVC） | 横喷嘴的喷流干涉 | 激波-边界层干涉 |
| 柔性喷嘴 | FSI耦合 | 结构变形耦合 |
| 活塞式喷嘴 | 可动活塞网格变形 | 喉部面积变化 |
二次流体喷射是研究最活跃的方向。在主流超音速膨胀流中横向喷入亚音速或超音速喷流，诱发激波产生横向力。精确预测激波-边界层干涉需要DES级以上乱流模型。
形状最优化
CFD和喷嘴形状最优化结合的研究吗？
伴随法（adjoint method）基于梯度的灵敏度分析效果非常强。以推力或比推力作目标函数，能高效计算喷嘴壁面形状参数（贝塞尔曲线控制点等）的梯度。
- Ansys Fluent：内置adjoint求解器。与壁面形状网格变形联动
- SU2：开源CFD。伴随法设计最优化实力强。由斯坦福大学开发
- DAFoam：基于OpenFOAM的离散伴随最优化框架
SU2在航空宇宙很受欢迎吧。
完全同意。SU2的Euler/RANS伴随法很成熟，翼型和喷嘴形状最优化论文很多。免费使用，大学研究也很理想。
"航空尖峰喷嘴"——下一代火箭重新关注的理由
普通钟形喷嘴只在特定高度（设计高度）达到最优效率。随大气压变化，就会过膨胀或不足膨胀，推力下降。航空尖峰喷嘴（aerospike nozzle）中心是尖峰形，周围燃气流动。理论上全高度自动实现最优膨胀的"高度补偿功能"。NASA在1970年代已研究但冷却设计复杂没实用化。现在3D打印能制造复杂冷却流路，民间航天公司又在开发。CFD形状最优化在这个领域重放光彩。
喷嘴流的故障排除
常见故障与对应方法
喷嘴CFD初学者容易遇到什么问题，请告诉我。
来看喷嘴特有的故障模式。
1. 喉部不壅塞
症状：喉部马赫数达不到1
原因和对策：
- 背压过高 → 降低背压到 $p_b/p_0 < 0.528$（空气）
- 入口用Mass-flow-inlet而非Pressure-Inlet → 指定超过壅塞流量会导致求解器自动提高全压。应该用全温全压设置入口
- 网格太粗 → 加细喉部附近的轴向网格
2. 激波振荡不收敛
残差不下降，激波反复摇摆…
原因：喷嘴内激波对背压很敏感，数值振荡会摇晃激波位置。
对策：
- 降低CFL数（Fluent设0.5-1.0左右）
- 先1阶精度收敛，再切换2阶精度
- 用Fluent的FMG（全重网格）初始化
- 残差完全定常化不了，则切换URANS（可能物理上非定常）
3. 质量流量与理论值不符
CFD的质量流量和理论值差几个百分点，正常吗？
2-3%的偏差可以用壁面边界层的有效截面积减少（排移厚度效应）解释。但5%以上就要检查。
- 喉部截面积定义对否（确认CAD喉部径）
- 气体常数 $R$ 和比热比 $\gamma$ 输入有无错误
- 二维轴对称计算时楔角和轴条件对否
- 入口全温全压设置对否
4. 壁面压力呈锯齿状
壁面压力分布图呈锯齿波形…
原因：网格质量问题。特别是壁面法向方向单元配列不均匀。
对策：
- 沿壁面切向方向的单元尺寸变化率（growth rate）控制在1.2以下
- 喉部曲率大处加密网格
- 非结构网格的壁面加棱柱层
- 后处理时局部平均滤波（不改物理但改善可视化）
5. 过膨胀喷嘴的壁面分离
过膨胀条件下喷嘴壁面分离预测不出来…
RANS用来精确预测过膨胀喷嘴的壁面分离有困难。下面的经验关联式可用来估计分离位置后与CFD比较。
Schmucker-Summerfield关联式：
这里 $p_{sep}$ 是分离点壁面压力，$p_a$ 是周围大气压。SST k-omega模型中壁面压力恢复到该值的位置可估计为分离点。
RANS有局限时怎样？
喷嘴流"堵住"还是"膨胀"——过膨胀与不足膨胀的判法
火箭喷嘴CFD中遇到"改变背压流量还是不变"的现象，多数是进入"壅塞"状态。喷嘴喉部达到M1时，下游信息无法向上游传播，背压再怎样改也流量不变——这是正常工作。关键是分辨出口处是过膨胀还是不足膨胀。出口压力低于背压就是过膨胀，会在外侧产生斜激波；高于背压是不足膨胀，膨胀扇展开。CFD可视化看出口附近压力分布就能判出是哪种，这是排故的铁律。