2.2 C-J 假设与 ZND 结构：解剖爆轰波
2.2 C-J 假设与 ZND 结构：解剖爆轰波
C-J 理论用一行条件——波后流动速度恰等于当地声速——锁定了爆轰速度的下限，让爆速成为可预测的物性参数；ZND 模型则进一步给出爆轰波的内部结构：引领激波压缩未反应炸药、其后跟随一个化学反应区。本节用勘测数据检验这两个模型，顺带引出炸药状态方程。
上一节说明了爆轰为什么能自持，这一节回答"自持的波内部长什么样、速度由什么决定"。它是第 5 章数值模拟爆轰时喂给程序的物理内核。
C-J 条件：一行式子锁定爆速
把第 1 章的三大守恒跨过爆轰波阵面写出来（激波间断关系），再补上炸药的反应释热和产物状态方程，理论上能解出波后状态。但守恒关系只给出一条曲线——雨贡纽曲线，曲线上每一点都对应一种"可能的爆轰"，速度并不唯一。
查普曼与朱格特的贡献是补上了缺失的选择规则：稳定的爆轰对应雨贡纽曲线上产物流速恰等于当地声速的那一点（C-J 点）。物理直觉是：比它慢的爆轰，波后压缩波能追上前沿给波"加油"使其加速；比它快的，波后稀疏波会追上前沿"抽油"使其减速。只有 C-J 点上扰动恰好与波平行传播，速度自锁。
用 Python 从产物状态方程出发算一次 C-J 爆速（简化演示，JWL 状态方程）：
# TNT 的 C-J 参数与 JWL 状态方程快速核对 # JWL: p = A*exp(-R1*V) + B*exp(-R2*V) + C/V**omega （V 为相对体积） A, B, C = 3.7377e11, 3.747e9, 2.734e9 # 帕 R1, R2, omega = 4.15, 0.9, 0.35 rho0 = 1.63 # 装药密度 g/cm3 D_cj = 6.93 # 实测 C-J 爆速 km/s p_cj = 21.0 # C-J 压力 GPa（经验值） V_cj = 1 - 1.7*0 # 演示：C-J 相对体积约 0.76（手册值） V = 0.76 p = A*(1 - R1*V/R1)*0 # 占位说明：完整求 C-J 点需解方程组，见下方注释 # C-J 点数值解的思路（伪代码）： # 1) 在 V 上扫描，用 JWL 算 p(V)、e(V) # 2) 代入雨贡纽关系（含爆热 Q）检验能量闭合 # 3) 找流速=声速的 V，对应 D 与 p_cj print(f"TNT 装药密度 {rho0} g/cm3，实测爆速 {D_cj} km/s，C-J 压力 {p_cj} GPa") print("手册规律：p_cj ≈ rho0 * D**2 / 4（数量级核对用）") print(f" 核对：1.63e3 kg/m3 * (6.93e3 m/s)**2 / 4 = {1630*6.93e3**2/4/1e9:.1f} GPa")
输出显示数量级核对式算出约 19.5 吉帕，与手册值 21 吉帕相差约 7%——对一条经验式子来说相当好。勘测者未必亲手解雨贡纽方程，但一定要会用这条数量级核对式：现场拿到爆速估计（例如从残留坑尺寸或电探针测量），乘密度平方除四，就能立刻判断产物压力是否合理，筛掉不可能的证词。
ZND 结构：波内的一小段燃烧室
C-J 理论把爆轰波当成一个"面"，但真实爆轰波有厚度。ZND 模型画出了这个厚度里的三明治结构：
- 引领激波：厚约一个分子自由程，把未反应炸药绝热压缩到冯·诺依曼状态（压力可达 C-J 压力的两倍），温度瞬时冲到两千开以上；
- 化学反应区：厚约 0.1 到 1 毫米，被压缩的炸药在这里完成诱导与放热，压力从冯·诺依曼尖峰松弛到 C-J 状态；
- 泰勒稀疏波：反应区之后，产物膨胀驱动周围介质。
图：ZND 一维结构。注意压力的"先峰后谷"：最高压力出现在完全反应之前——这是很多工程计算用 C-J 压力会偏保守的原因。
图：ZND 压力剖面。尖峰出现在反应完成之前，多数工程计算取 C-J 压力作保守近似即可。
ZND 是理想化的一维模型，真实凝相炸药的反应区是三维的：前沿激波并非平直，而是由一系列马赫杆与入射波组成的胞格结构，波后压力在胞格的不同位置周期起伏。这就是"爆轰胞格"，第 3 章讲马赫反射后你会更清楚胞格里的波系几何。勘测上，爆轰胞格可以在烟熏铝板上留下痕迹——胞格尺寸是估算炸药临界起爆直径的关键数据。
状态方程：让模型"落地"的最后一环
无论 C-J 还是 ZND，都要喂一个产物状态方程才能算数。工程上最常用的是 JWL（Jones-Wilkins-Lee）方程，它以相对体积的指数函数拟合爆轰产物的等熵膨胀线，覆盖从 C-J 点到后期膨胀的全过程；研究气相爆轰则常用 BKW 或更简单的多方气体方程。选哪个状态方程，直接决定第 5 章数值模拟里冲击波出流场的压力曲线。勘测者不必推导 JWL，但要记住它的适用边界：标定于特定装药密度，密度偏离超过 5% 时参数需要重新标定，套用旧参数算出的远场超压可能偏差两成。
⚠️ 常见坑：把冯·诺依曼尖峰压力当驱动压力用。对近距离接触爆炸（第 7 章聚能射流、第 8 章爆破），局部载荷确实由尖峰主导；但炸药自由场冲击波参数用 C-J 压力即可，因为尖峰只维持亚毫米尺度。混用两者是新手模拟结果对不上的高发原因。
三点五、从参数表到勘测：C-J 数据的现场用法
爆轰参数表（装药密度、爆速、C-J 压力、C-J 温度）是勘测者的常备工具，三类用法值得记牢。其一，密度自洽检查：手册 JWL 参数标定于特定装药密度，现场取样密度偏差超过百分之五，爆速估算应按密度与爆速的近似线性关系修正，每偏离零点一克每立方厘米，爆速变化量级为每秒百米。其二，殉爆能力粗估：起爆冲量与 C-J 压力及脉冲宽度正相关，用隔板试验数据类比现场刺激强度时，两者的装药密度必须先对齐。其三，产物做功能力核对：水下爆炸试验的气泡能与冲击波能之和，除以装药爆热得到能量利用率，这一栏对不上，优先怀疑反应完全度而非测量。
本节要点回顾
- C-J 条件：产物流速等于当地声速，速度自锁，使爆速成为可预测的炸药物性；
- 数量级核对式：C-J 压力约等于密度乘爆速平方除四，勘测现场的快速验算器；
- ZND 三层结构：引领激波、反应区、泰勒稀疏波；最高压力出现在反应完成之前；
- 状态方程是落地件：JWL 参数与装药密度绑定，越界使用误差显著。
下一节处理最贴近工程事故的问题：起爆感度与爆燃转爆轰——一次意外，往往就发生在这两个概念的交叉处。