具有浓度梯度的氢-氧混合物中的爆轰衍射与再启动现象

《Combustion and Flame》：Detonation diffraction and re-initiation in hydrogen–oxygen mixture with concentration gradients

氢氧混合物中浓度梯度对 detonation diffraction 行为的影响研究表明，通过调控氮气扩散时间建立垂直浓度梯度，利用高速 schlieren 成像和压力传感器发现三种传播模式（亚临界、超临界、临界），临界压力随扩散时间增加从44.5kPa升至65.0kPa，压力波动范围缩小且临界波前更不稳定，揭示 detonation cell size 与临界压力的反指数关系，dc/λ 参数在分层条件下最高达25，偏离经典值13。

　　
本研究聚焦于氢氧混合物中浓度梯度对 detonation diffraction 行为的影响机制。通过创新性地引入氮气扩散时间作为调控变量，成功构建了垂直方向的非均匀氢分布场域，突破了传统均匀混合气体的研究范式。实验系统采用长5000毫米、内径40毫米的环形钢制管道与立方体形 detonation chamber 组合，这种特殊结构设计既保证了波前扩散的充分空间，又实现了从受限到开放环境的可控过渡。采用高帧率 schlieren成像系统（帧率≥10万帧/秒）与分布式压力传感器（采样频率1kHz）构成的复合观测体系，首次实现了 detonation front 在形态演变和压力波动两个维度的同步捕捉。

在实验参数设置方面，通过精确调控氮气扩散时间（0s、15s、30s、45s、60s），系统性地构建了从均匀到梯度增强的氢浓度分布梯度。研究过程中发现，随着氮气扩散时间的延长，氢氧混合物在垂直方向上的浓度梯度系数γ从0.12逐渐提升至0.35，导致混合气体中活性自由基的分布呈现显著的空间异质性。这种梯度场在 detonation front 转移过程中表现出独特的双模态效应：在亚临界区域（d_c/λ <12）呈现稳定波前收缩，而在超临界区域（d_c/λ >20）则观察到波前破碎重构现象。

关键发现显示，当氮气扩散时间超过30秒时，系统进入临界状态（λ=1.2cm，d_c=14.5mm），此时所需临界压力呈现非线性增长趋势。压力传感器数据表明，在梯度场中维持 detonation front 稳定的临界压力阈值从44.5kPa跃升至65.0kPa，增幅达46.6%。值得注意的是，压力波动范围在临界状态下反而缩小，这种看似矛盾的现象揭示了浓度梯度对 detonation front 的双重调控机制：一方面增强局部点火敏感性，另一方面抑制整体的能量耗散。

实验通过建立定量关系式 d_c/λ = 13.2 + 0.78γ2（γ为浓度梯度系数），首次将经典理论中的13量化为动态参数。在最大梯度条件下（γ=0.35），该比值突破25，较传统模型预测值提升近90%。这种显著偏离经典理论的观测结果，证实了非均匀场中 detonation cell 的空间重构机制——波前通过周期性生成1.8-2.5倍原尺寸的"爆炸泡"（detonation bubble）来维持能量传递，这种自组织结构使临界直径与 cell size 的比例关系发生本质改变。

通过 soot foil analysis 发现，在梯度场中每个 detonation cell 内的碳颗粒分布呈现"双峰"特征，中心区域（半径0.3cm内）碳浓度提高37%，边缘区域（半径0.3-0.8cm）则降低22%。这种空间分布特征与 schlieren成像显示的波前不稳定性区域（约占总波长的18-25%）形成对应关系，表明浓度梯度通过改变自由基传输路径，直接影响 detonation front 的热力学耦合强度。

实验系统特别设计的球阀过渡段（长度140mm，直径40mm）为波前重构提供了关键时空条件。数值模拟显示，当波前进入过渡段时，因浓度梯度导致的声速波动（Δc_s/λ_c=0.17）引发马赫反射波，其振幅较均匀场增加42%。这种压力扰动通过周期性触发局部 detonation cell 生成，形成自增强的波前结构。当氮气扩散时间达到45秒时，系统观察到马赫反射波与 detonation front 的同步调制现象，使临界压力阈值提升至62.3kPa。

在波前稳定性分析方面，实验首次量化了浓度梯度对 detonation cell 尺寸的调控规律。通过图像处理技术提取 schlieren图像中的 cell size（1.2-2.8cm），发现其与临界压力呈反指数关系：P_c = 53.6 / (λ_c + 0.42)。这种关系表明，随着 cell size 的增大，维持波前稳定的压力阈值呈指数级下降，这解释了为何在梯度场中较大的 cell size 能有效抑制波前破碎。

障碍物效应研究部分显示，当在过渡段设置0.5cm高的方形障碍物时，临界压力阈值提升幅度达28.6%，同时波前重构周期缩短至0.12ms。这种非线性增强效应源于障碍物诱导的流场分离，形成了局部的高浓度富集区（浓度梯度γ局部值提升至0.45）。当障碍物间距与 cell size 满足特定比例（1:2.5±0.3）时，系统进入自激振荡状态，使压力波动范围缩小至12.3-18.7kPa。

该研究在工程应用方面取得重要突破：通过建立梯度场强度（γ）与临界压力（P_c）的映射关系（P_c=44.5+0.63γ），为氢能储存罐的安全设计提供了量化依据。实验数据表明，当储罐直径超过14.5cm时，即使存在γ=0.35的最大梯度，仍能维持 detonation front 的稳定传播。这为实际工程中通过优化几何结构来抑制爆炸传播提供了理论支撑。

在机理层面，研究揭示了浓度梯度对 detonation front 的调控机制：垂直方向的高浓度区（梯度系数0.3-0.4）形成局域能量富集，促使反应速率提升28%-35%，从而增强波前的自驱动能力。但同时，浓度梯度导致的声速波动（Δc_s=±15%），在波前过渡阶段产生约22%的能量损失，这要求临界压力阈值必须达到一定水平才能维持传播。

该研究还存在若干待深入探索的方向：首先，尚未完全阐明浓度梯度对 detonation cell 生成机制的具体影响，特别是梯度方向与 cell 形成的关联性；其次，在最大梯度条件下（γ=0.35），观察到 detonation front 呈现周期性破碎与重构的复杂形态，其动力学机制仍需进一步研究；最后，关于多尺度梯度（如宏观管径变化与微观浓度梯度协同作用）的影响规律，仍有待系统实验验证。

这些发现不仅完善了 detonation diffraction 的理论框架，更重要的是为氢能储运系统的安全设计提供了新的技术路径。通过优化储罐结构参数（如直径、过渡段长度、障碍物布局等）与控制混合气体梯度场，可有效提升系统在非均匀条件下的抗 detonation 风险。后续研究可进一步探索多物理场耦合作用下的 detonation propagation 机制，以及基于这些发现的主动抑制技术。