爆轰（爆轰波）
爆轰（爆轰波）的理论基础
概述
教授，爆轰与爆炸有什么区别？
燃烧传播有两种形式。爆燃（缺陷燃烧）是亚音速火焰传播，普通燃烧器火焰属于此类。另一方面爆轰（爆轰波）是冲击波与燃烧波结合形成的超音速传播现象。传播速度在甲烷/空气中约1800 m/s，在氢气/空气中约2000 m/s。
冲击波和燃烧一起传播，是这样吗？
是的。先行冲击波对未燃混合气进行绝热压缩以升高温度，高温下化学反应迅速进行，反应能量维持冲击波。这种自我维持机制是爆轰的本质。
Chapman-Jouguet理论
Chapman-Jouguet（CJ）理论是什么？
CJ理论是热力学地求出爆轰波传播速度的理论。在冲击波前后的Rankine-Hugoniot关系式中加入化学反应放热量 $q$，应用CJ条件（爆轰波后流速恰好等于局部音速）。
用简化形式，CJ爆轰速度大致可以写成：
这里 $\gamma$ 是比热比，$q$ 是单位质量的放热量 [J/kg]。
还有CJ马赫数的概念吧？
CJ爆轰的马赫数由以下公式给出：
氢气/空气（当量比1.0）的情况下 $M_{CJ} \approx 5.0$，甲烷/空气为 $M_{CJ} \approx 5.2$。
ZND结构
爆轰波的内部结构是什么样的？
ZND（Zel'dovich-von Neumann-Doering）模型中，爆轰波具有三层结构。
1. 冲击波面（von Neumann spike）：未反应气体受到冲击压缩。压力达到CJ值的约2倍
2. 诱导带（Induction zone）：化学反应进行前的延迟区间。对应著火延迟时间
3. 反应带：迅速进行化学反应，达到CJ状态
诱导带越短爆轰越稳定，是这样吗？
是的。诱导带长度 $\Delta_i$ 直接关系到细胞大小 $\lambda$，有经验公式 $\lambda \approx (10-30)\Delta_i$。当这个细胞大小相对于爆轰管直径足够小（管直径 > 数个$\lambda$）时，稳定爆轰传播才能维持。
爆轰理论是冲击波力学和化学反应速率论的融合呢。
完全正确。用CFD处理需要既能准确捕捉冲击波的数值方法，又能处理高温高压下化学反应速率。
爆轰——爆轰波以「音速5～10倍」传播的原因
燃烧有「爆燃（缺陷燃烧：音速以下火焰传播）」和「爆轰（爆轰波：冲击波与燃烧一体化的超音速传播）」两种。爆轰波由Chapman-Jouguet（CJ）理论（1899～1905年）描述，燃烧气体在音速下的流动对应「CJ面」，其速度是特性值。氢-空气混合气约2 km/s，天然气-空气约1.8 km/s的爆轰速度。CAE中将欧拉方程与化学反应模型结合，用高解像度数值方案分析爆轰波传播、细胞结构、转移（DDT：爆燃到爆轰的转移）。
爆轰（爆轰波）的数值计算方法
数值方法的详细
用CFD求解爆轰需要什么数值方法？
爆轰的数值分析与普通RANS燃烧分析大不相同。需要准确捕捉冲击波的高解像度方案，同时兼顾高温高压下的详细化学反应。
空间离散化
捕捉冲击波适合用什么方案？
冲击波捕捉常用高阶精度TVD（Total Variation Diminishing）方案或WENO（Weighted Essentially Non-Oscillatory）方案。
| 方案 | 精度 | 特点 | 应用 |
|---|---|---|---|
| Roe + Minmod | 2阶 | 稳定但数值扩散大 | 初期研究 |
| HLLC | 2阶 | 分析接触不连续面 | 通用 |
| WENO-5 | 5阶 | 高精度但计算成本大 | DNS/高精度计算 |
| MUSCL-Hancock | 2阶 | 成本-精度平衡良好 | 实用爆轰计算 |
WENO方案最理想，但实务上用MUSCL就足够了吗？
正是。MUSCL+HLLC黎曼求解器的组合是实务的折中方案。但网格大小要保证为爆轰细胞宽度 $\lambda$ 的1/20以下。氢/空气（当量比1.0, 1 atm）时 $\lambda \approx 10$ mm，所以网格宽度需要0.5 mm以下。
时间积分
时间积分怎么处理？
爆轰波传播是微秒级现象，通常用显式时间积分。但化学反应部分是Stiff的，所以用算子分裂（Strang splitting）分离流体输送和化学反应。
我给出典型设定值。
| 参数 | 推荐值 | 备注 |
|---|---|---|
| CFL数 | 0.3-0.5 | 冲击波捕捉需保守选择 |
| 化学反应求解器 | CVODE (BDF) | Stiff系统必需 |
| 最小时间步 | $10^{-9}$ s | von Neumann尖峰分析 |
| 网格宽度 | $\lambda/20$ 以下 | 细胞结构分析最低条件 |
AMR（适应格子细分化）
均匀的细网格会导致计算成本巨大吧？
这时AMR就大显身手了。只在爆轰波面附近网格细分，未燃和已燃区域保持粗网格。CONVERGE内置自动AMR，OpenFOAM也可用 dynamicRefineFvMesh 实现。AMR可将记忆体和计算时间减少1-2个数量级。
AMR的细分基准用什么？
通常用温度梯度 $|\nabla T|$ 或压力梯度 $|\nabla p|$ 作为细分传感器。OH质量分率梯度也有效用于追踪反应带。
爆轰的数值计算是冲击波捕捉方案 + Stiff化学反应 + AMR的三位一体呢。
完全同意。少了其中任何一个都无法进行实用爆轰计算。
爆轰计算需要「超音速方案」的原因——CFL数0.3的壁障
爆轰数值计算最初面对的困难是离散化方案。爆轰波以超音速（马赫5～10等效）传播，普通燃烧CFD用的2阶精度中心差分方案会立即发散。解决办法是必须使用WENO或Roe法等冲击波捕捉方案，CFL数也需限制在0.3以下。这意味著时间步比普通燃烧计算小一个数量级，计算时间增加10倍以上。「爆轰是CAE计算中特别难处理的一类」的原因在此。
爆轰（爆轰波）的实务应用
实践指南
请教我实际进行爆轰分析的步骤。
爆轰分析的计算工作量比普通燃烧CFD大得多。让我整理一下实用流程。
分析流程
1. 0D/1D预备计算 -- 用Cantera或SD Toolbox计算CJ条件（传播速度、压力比）。从ZND结构估出诱导带长度和细胞大小
2. 网格设计 -- 细胞大小$\lambda$的1/20以下的分析度。使用AMR时设定最大细分级别
3. 爆轰开始条件 -- 设置Driver section（高压高温区）或用热点直接著火
4. 时间积分设置 -- CFL = 0.3-0.5，总计算时间直到爆轰波到达管端
5. 后处理 -- x-t图、细胞结构可视化、压力温度分布
RDE（回转爆轰引擎）的设定示例
最近受关注的RDE如何分析？
RDE中环形燃烧室的爆轰波沿周向传播。从2D展开模型（将环形展开为平面）开始计算，再扩展到3D比较有效率。
给出RDE分析的典型设定。
| 项目 | 设定值 | 备注 |
|---|---|---|
| 燃料 | H2/空气，当量比1.0 | CH4著火延迟太长，RDE不稳定 |
| 环形直径 | 100-200 mm | 研究规模 |
| 网格 | 0.2-0.5 mm | 推荐使用AMR |
| 反应机构 | Li et al. 9种/21反应 | H2详细机构已足够紧凑 |
| 入口条件 | 总压指定 + 混合气组成 | 推荐choked inlet防止逆流 |
| 计算时间 | 数毫秒（5-10转） | 确认定常旋转状态 |
H2较多是因为著火延迟短吧。
是的。甲烷著火延迟是氢的100倍以上，难以维持回转爆轰波。实用化方面正在检讨甲烷/氢混合、乙烯（C2H4）等。
常见失败与对策
| 症状 | 原因 | 对策 |
|---|---|---|
| 爆轰不传播 | 网格太粗 | 细分到$\lambda/20$以下 |
| 传播速度与CJ值不符 | 数值扩散大 | 改用高阶方案(MUSCL/WENO) |
| 看不到细胞结构 | 2D计算分析度不足 | 用0.1 mm以下网格重新计算 |
| RDE波消失 | 入口混气供给跟不上 | 增加入口面积或总压 |
爆轰分析是计算成本的战争呢。
是的。用2D AMR掌握概况，再对必要部分进行3D详细计算是现实做法。
爆轰管爆炸事故——「爆燃到爆轰的转移（DDT）」的恐怖
化学工厂管道中氢氧混合气泄漏时，最初是温和燃烧（爆燃），但火焰在管内传播时会加速并转化为爆轰（爆轰波），这称为DDT（Deflagration-to-Detonation Transition），是配管设计人员最担心的现象之一。爆轰转化的瞬间，压力会跳升到爆燃的10～20倍。爆轰CFD实践中，为了进行安全评估，必须准备足够细腻的网格来正确模拟DDT，并准确追踪从点火开始的时间演变。
爆轰（爆轰波）的软体比较
商用工具比较
能进行爆轰分析的CFD工具有哪些？
爆轰需要冲击波捕捉和详细化学反应两个条件，所以对应工具比较有限。
| 工具 | 冲击波方案 | AMR | 详细化学反应 | 适性 |
|---|---|---|---|---|
| CONVERGE | PISO + ALE | 自动AMR | SAGE | 最优 |
| Ansys Fluent | Density-Based | 无（UDF外部） | Stiff Chemistry | 可能但需要技巧 |
| OpenFOAM (rhoCentralFoam) | KNP/KT | dynamicRefine | ode | 灵活，需编码 |
| STAR-CCM+ | Coupled Flow | 无 | DARS | 有限 |
| 专用代码 (AMROC, PeleC) | Godunov类 | Block-AMR | 详细/缩约 | 最适研究 |
CONVERGE对爆轰最优出乎意料，印象中它是内燃机的工具呢。
CONVERGE的自动AMR和网格切割技术对爆轰波面追踪非常相容。省去了大量网格生成工作。实际上RDE和脉冲爆轰引擎（PDE）研究中CONVERGE应用增加。
Ansys Fluent中的注意点
用Fluent进行爆轰计算时的注意点？
OpenFOAM中的实现
OpenFOAM用哪个求解器？
rhoCentralFoam 是可压流的标准求解器，但需自行实现与化学反应耦合。社区版有 rhoReactingCentralFoam。用 dynamicRefineFvMesh 启用AMR，以温度、压力梯度为细分基准是定石。
选定指南
最后选哪个比较好？
专用代码的存在不知道。爆轰特殊要求的工具选择很重要呢。
是的。直接套用普通燃烧CFD工具会因数值扩散导致爆轰细胞结构消失。要根据工具特性选择。
爆轰CAE的「客户」是防卫、太空、安全评估——商用工具支持滞后的原因
爆轰模拟市场很小。客户主要限于防卫机构、太空机构、工厂安全评估，所以通用商用CFD工具（Fluent、STAR-CCM+）的爆轰专用功能有限。ANSYS的Fluent中可选WENO方案，但爆轰设置示例在官方手册中很少。实际上爆轰计算常用OpenFOAM的XiFoam或blastFoam（开源），或SANDIA、Lawrence Livermore自主开发代码。商用工具「对爆轰投资多少」取决于市场规模——CAE世界也由经济原理驱动。
爆轰（爆轰波）的前沿研究
前沿话题与研究动向
爆轰最前沿研究有哪些？
主要有三个方向。(1)回转爆轰引擎（RDE）实用化，(2)爆轰驱动的压力增益燃烧，(3)倾斜爆轰波引擎（ODWE）。
RDE（回转爆轰引擎）
为什么RDE受关注？
RDE的爆轰波在环形燃烧室连续回转，可提供定常推力。相对传统爆燃燃烧器，熵增加较小，理论上热效率提升5-15%。
| 机构 | 研究成果 | 燃料 |
|---|---|---|
| JAXA | 世界首次太空RDE实证（S-520火箭） | 乙烯/O2 |
| 美空军研究所 (AFRL) | 与GE共同将RDE用于喷气引擎燃烧器 | JP-8/空气 |
| 名古屋大学 | 连续回转爆轰稳定化方法 | H2/空气 |
| KAUST | LES+详细化学反应的RDE数值分析 | H2/空气 |
JAXA在太空进行实证实验了呢。
2021年S-520-31号火箭实现了世界首次太空RDE作动。CFD对设计做出了巨大贡献。
倾斜爆轰波引擎（ODWE）
ODWE是什么？
这是为极超音速飞行体（马赫5-10）设想的推进装置，在楔形燃烧室内驻定倾斜爆轰波。是SCRAMJET之后的概念。驻波不需RDE的回转，但需针对马赫数优化楔角度。
DNS（直接数值模拟）
进行爆轰DNS吗？
近年Peta-FLOPS级HPC使得2D爆轰细胞结构可用DNS分析。使用PeleC或S3D（Sandia）。3D DNS还需Exascale级计算资源，现由美DOE的ECP（Exascale Computing Project）推进中。
机器学习的应用
爆轰研究也用机器学习吗？
爆轰研究直接关系航空推进最前沿呢。
是的。特别是RDE作为实用化技术各国竞相研究。CFD设计支持的重要性今后还会增加。
回转爆轰引擎——爆炸「持续回转」的逆转思维
连续运用爆轰的回转爆轰引擎（RDE）从2000年代起研究加速。普通燃烧是「定压燃烧」，爆轰接近「定容燃烧」，故热效率理论上高10～15%。JAXA与各国防卫机构投入资金的原因。但实用化困难重重，在环形室内维持超音速爆轰波稳定运行难度极高。现在前沿研究试图用机器学习对爆轰波动态行为进行即时控制。「控制爆炸」这一矛盾命题正成为该领域正面攻关的主题。
爆轰（爆轰波）的故障排除
故障排除
请教爆轰计算常见的故障。
爆轰计算比普通燃烧CFD敏感得多。整理主要故障。
1. 爆轰波不开始
可能原因：
- Driver section压力温度不足
- 网格太粗，冲击波与反应带无法分析
- 反应机构未能正确再现该条件著火延迟
对策：
- 设定Driver条件为 $p > 30\,p_1$, $T > 2500$ K
- 先从1D验证开始条件再进行2D/3D扩展
- 用ZND计算确认所需分析度（诱导带 $\Delta_i$ 至少20点）
2. 传播速度与CJ值偏差
传播速度与CJ理论值不符的情况怎么办？
速度低的情况：数值扩散过大，反应带能量散逸。改用高阶方案（MUSCL、WENO），细分网格。
速度高的情况：过驱动状态，初始Driver影响仍存。让爆轰波传播足够距离后再测（管直径10倍以上）。
| 偏差 | 允许范围 | 要因 |
|---|---|---|
| $D/D_{CJ} < 0.95$ | 需改善 | 数值扩散、网格不足 |
| $0.95 < D/D_{CJ} < 1.02$ | 良好 | -- |
| $D/D_{CJ} > 1.05$ | 需改善 | 过驱动、过渡状态 |
3. RDE特有的问题
RDE计算特有的问题有吗？
4. 计算成本巨大
对策：
- 活用AMR（波面附近高分析，其他粗网格）
- 先2D掌握概况再进行3D详细计算
- H2情况用Li et al.的9种/21反应就足够（无需缩约）
- 利用GPU并列化求解器（PeleC支援GPU）
爆轰计算是分析度和成本的权衡很严峻呢。
是的。所以用1D ZND分析估算必要分析度，再用AMR降低成本的步骤很重要。直接上3D均匀网格的话，计算可能数月不完成。
「爆轰波消失了」——爆轰CFD重启困难的原因
爆轰计算故障排除的特有问题是「计算中途爆轰波消失，降级为爆燃（低速燃烧）」的现象。这种「quenching（消炎失速）」在实验中也发生，但计算中多因数值黏性导致爆轰波人为衰弱。典型是网格太粗无法分析爆轰细胞结构（detonation cell），经验法则「需要细胞宽度1/5以下的网格」。但实际爆轰细胞宽度只有数毫米至数厘米，实机规模全域细化会成本不现实，形成矛盾。
相关话题
价值
更详细
错误