水雾抑制氢–空气爆轰的稳态传播区划与爆轰极限

《International Journal of Hydrogen Energy》：Water spray mitigation of hydrogen-air detonations: Steady-state regimes and detonability limits

　　

氢–空气爆轰因其超声速传播和极端压力与温度而构成重大安全危险。水雾喷射过去已被提出作为一种有效的缓解策略。本研究采用一维稳态多连续介质（multi-continua）模型，结合详细化学动力学以及两相质量、动量和能量传递，考察了氢–空气爆轰与单分散水雾之间的相互作用。模型中纳入了非量热完全气体（non-calorically perfect gas）效应和液滴蒸发。结果表明，在简化模型假设下，爆轰行为主要由水质量载荷（water mass loading）和液滴尺寸控制。小液滴快速蒸发，使得爆轰速度可通过考虑潜热的全局能量平衡加以预测，且Chapman–Jouguet（CJ）条件在热化学平衡处得到满足。在此区划内增加水载荷可导致爆轰完全猝灭。对于较大液滴，则出现病理爆轰（pathological detonations），此时蒸发和放热与驱动机制解耦，传播速度由内部爆轰结构而非全局能量平衡决定。极大液滴蒸发过慢，无法显著影响爆轰，传播速度趋近于干燥氢–空气混合物的值。这些发现阐明了氢安全应用中有效爆轰缓解的参数区划，并使研究人员能够在存在液态水分散相的条件下重新绘制氢–空气混合物的理论爆轰极限。

**论文解读**

**一、研究背景与问题提出**

氢气被公认为未来可持续能源系统中的关键组分，然而其独特的燃烧特性——包括宽广的可燃范围（flammability range）和低点火能量——带来了严重的安全挑战，其中爆轰风险尤为突出。爆轰以超声速传播并伴随极端压力和温度，一旦发生将造成灾难性后果。因此，理解氢爆轰的传播极限及其抑制条件是危险缓解研究的核心议题。

在爆轰混合物中引入稀释剂可显著影响爆轰特性，包括波前传播速度和Mach数。水雾（water mist）因其能够吸收大量热量（显热和蒸发潜热）、具有高比表面积且成本低廉，被视为有效的爆轰缓解手段。然而，水雾与爆轰的相互作用涉及多个耦合物理过程：液滴首先升温至沸点，随后经历吸热的蒸发相变，显著降低激波后气体温度；水蒸气本身作为稀释剂降低混合物可燃性；水蒸气还与中间自由基（如原子氢H）发生化学作用，有助于抑制燃烧反应；气液两相之间的动量和能量传递（包括传热和曳力）减弱冲击波强度。这些效应的综合结果决定了爆轰的最终状态——从稳定传播到完全熄灭。此前研究表明，水雾的有效性强烈依赖于水质量载荷比（water mass loading ratio，即单位质量气体对应的液体质量）和初始液滴尺寸。然而，对于不同液滴尺寸和水载荷组合下爆轰行为的系统性区划，以及相应的爆轰极限理论，尚缺乏深入认识。这正是本研究的出发点。

**二、研究内容与意义**

研究人员采用一维稳态多连续介质模型，将爆轰描述为一层无限薄的非反应激波后跟随两相反应流。气相采用Zel’dovich–von Neumann–D?ring（ZND）模型描述，液相则假定为稀薄、单分散的水雾，采用Euler方法处理。模型完整考虑了液滴蒸发、气动曳力和传热引起的两相质量、动量和能量交换，并采用San Diego详细化学动力学机制（包含8种组分和21个基元反应）描述H₂/O₂燃烧过程。通过数值积分稳态爆轰剖面，系统考察了不同水质量载荷比和激波后液滴半径下爆轰的内部结构。

该研究的核心意义在于：首次在统一的参数空间中识别出四种截然不同的爆轰传播区划——全局平衡区划（Regime I）、无效区划（Regime II）、猝灭区划（Regime III）和病理爆轰区划（Regime IV），并据此重新绘制了含水雾条件下氢–空气混合物的理论爆轰极限图。这一成果为氢安全领域中水雾防护系统的设计提供了关键参数依据，尤其是明确了有效缓解所需的液滴尺寸和水载荷范围。

**三、关键技术方法**

本研究采用一维稳态多连续介质模型，以Euler方法描述稀薄单分散水雾相，气相则基于ZND模型框架。激波后状态由Rankine–Hugoniot（RH）跃变条件确定，假定组分冻结且两相交换可忽略。数值求解采用二分法寻找满足CJ自持条件的传播Mach数特征值（eigenvalue），积分使用SciPy软件包的Radau隐式Runge–Kutta方法。全局平衡CJ Mach数通过热化学平衡求解器（Combustion Toolbox）独立计算作为对照。化学动力学采用San Diego机制，涵盖H₂/O₂系统的21个基元反应，用于获取诱导长度（induction length）、传播速度和自由基行为。液滴蒸发模型基于Stefan流和Clausius–Clapeyron关系，曳力模型采用基于Reynolds数（Re）的阻力系数关联式。爆轰极限分析则基于Belles准则及其修正形式，通过比较激波后温度与由链分支反应与终止反应竞争决定的交叉温度（crossover temperature）来判定。

**四、主要研究结果**

**（1）全局平衡分析与无穷小液滴极限**

当液滴足够小以致调节和蒸发时间远小于特征反应时间时，可在激波后薄层内完成蒸发。此时全局守恒方程可积分得到修正的反应RH方程，CJ传播Mach数可通过热化学平衡求解器或简化解析表达式计算。结果表明，增加水质量载荷持续降低自持传播Mach数。在无穷小液滴极限下，修正Belles准则同时计入水蒸气的化学效应（增强链终止反应）和蒸发冷却效应，使抑制爆轰所需的水质量载荷较经典准则降低近一个数量级。

**（2）两相ZND结构与病理爆轰**

对于α_u=0.2、a_N=0.2μm的工况，液滴在到达声速点前完全蒸发，传播Mach数与全局平衡解一致（M_CJ=4.39），属于规则爆轰。当液滴增大至a_N=1μm（相同水载荷），诱导长度与液滴调节长度可比，内部结构发生质变：声速点在液滴完全蒸发之前出现，蒸发吸热发生在与导激波声学隔离的区域，导致实际传播速度（M_CJ=4.52）高于全局平衡预测——即病理爆轰。对于a_N?20μm的大液滴，惯性过大，蒸发和调节时间远长于化学时间，液相几乎不产生影响，传播速度趋近于纯气相爆轰。

**（3）传播Mach数与诱导长度的参数研究**

对于a_N≤1μm的液滴，传播Mach数几乎不受液滴半径影响，与无穷小液滴极限一致。超过此阈值后，M_CJ随半径增大渐进趋近气相爆轰极限。诱导长度对液滴存在的敏感度远高于传播速度：对于a_N≤0.1μm的小液滴，诱导长度可在毫米至米量级变化，关联于爆轰猝灭。

**（4）四区划图和爆轰极限**

在水质量载荷比-液滴半径参数空间中识别出四种区划：区划I（低水载荷、小液滴）——传播速度可由全局平衡预测，接近爆轰极限处诱导长度对水载荷极度敏感；区划II（低水载荷、中等液滴）——蒸发发生在燃烧区和声速点下游，爆轰几乎不受水雾影响；区划III（中高水载荷、小液滴）——蒸发吸热导致猝灭；区划IV（中等水载荷、中等液滴）——病理爆轰，声速点出现在非平衡区有限距离处，M_CJ在log(a_N)–log(α_u)空间中呈准线性标度，临界液滴半径a*_N略低于1μm。爆轰极限图表明，液滴越小、水载荷越高，可爆轰区域越窄。

**五、讨论与结论**

本研究揭示了液滴尺寸决定了分散相的响应时间尺度（调节、加热、蒸发），与基于Stokes数的无量纲分析一致；而水质量载荷比则调控分散相的集体效应。病理爆轰区划的发现说明，仅凭全局能量平衡不足以预测某些参数组合下的传播速度，必须解析内部结构。此外，气动破碎效应可能使大液滴局部恢复小液滴区划，导致缓解有效性对液滴尺寸的非单调依赖——存在一个中间尺寸范围使水雾缓解效果最差。

研究结论指出：本工作建立了一个基于尺度的理论框架，可指导氢–空气–液态水混合物的爆轰极限探索。通过采用修正Belles准则，证明了计入蒸发能量预算和水蒸气与中间自由基的化学相互作用，可使抑制所需的水质量载荷降低近一个数量级。这些发现阐明了有效缓解的参数条件，为氢能应用中水雾安全系统的设计提供了关键依据。未来研究将拓展至多分散水雾和气动破碎的详细建模，以及多维非稳态效应的引入——横向波和三维胞格结构有望引入额外的猝灭机制并改变液滴–流动相互作用。