当氢气遇见镍壁：微尺度爆轰诞生记

几年前，我住在俄勒冈州波特兰市郊的一个小区。不止一次，当我以每小时30英里的悠闲速度行驶在本地道路上时，有人会以极高的速度从我身边呼啸而过。虽然他们可能并没有达到每小时100英里（就像我会激动地告诉同车伙伴的那样），但感觉上确实如此。

如今，我在一家从事燃烧分析的公司工作，而那段经历正是我自学DDT（deflagration to detonation transition：爆燃转爆轰）概念的契机。如果，在某个反乌托邦的宇宙中，我的现实与那位飙车者的现实合二为一，那辆新车将以30英里的时速平稳行驶，然后在不到一秒内突然加速到100英里，理论上就经历了一次DDT。

DDT被定义为一种这样的过程：一个缓慢移动的火焰（即我的车）迅速加速为一个超音速的爆轰波（即飙车者的车）。DDT发生前的微秒阶段被称为火焰加速（FA：Flame Acceleration），这两个现象通常是一起研究的。过去，FA和DDT的研究是在大规模场景中进行的，比如超新星爆炸、大型激波管或煤矿巷道。然而，随着排放法规的日益严格和对更紧凑能源系统需求的上升，人们也开始在更小的尺度上研究它们，比如微通道。

这些微通道设备具有更强的传热传质能力，且制造成本更低，广泛应用于电子冷却、生物系统和暖通空调设备等领域。然而，可燃气体混合物在通过微通道这种高度受限的通道时更容易发生爆轰，这些通道的尺寸大约只有一根头发丝的直径。研究微通道中的FA及其随后的DDT，有助于我们更好地理解引发爆轰的条件，并为高压系统提供更有效的控制策略。

现有的大量关于爆炸安全的文献都集中在热壁边界条件的影响上，这些条件通过影响传热损失、火焰稳定性和点火行为，在火焰传播中起着重要作用。另一个可能影响火焰传播和爆轰的因素是非均相化学反应，特别是催化壁上的表面反应。在微反应器中，反应性催化壁涂层可以改变并诱导壁面的化学交换，从而影响FA和DDT过程。催化壁为燃料/空气混合物提供了一个反应表面，这种非均相燃烧发生在催化剂表面，而不是气相中。

大多数催化燃烧的研究都集中在贵金属（如铂或铑）上的催化燃烧。相比之下，像镍这样的过渡金属仅被研究用于化学重整——一种改变烃类分子结构以生产其他化学组分的过程。在本研究中，明尼苏达大学双城分校的博士生Suryanarayan（Surya）Ramachandran与Suo Yang教授以及埃克森美孚科技公司的研究工程师合作，研究了在具有含镍催化壁的微通道中的氢气点火与火焰传播。在这种高度受限的微通道环境中，FA和DDT的问题变得更加突出。现在，我把发言权交给Surya，让他来介绍他的研究！

CONVERGE 案例研究  

污染、稀释与纯度   

在理想的氢燃烧系统中，燃料/氧化剂混合物应只含氢气以及空气中的氧气和氮气。然而在工业现场，部分燃烧产物（如水蒸气）会回流到燃料/氧化剂中，使混合物被 H₂O 严重“污染”。为了贴近真实工况，我们先用氢-氧-氮-水四组分混合气进行模拟，该算例命名为 C1，结果未出现爆轰。

既然课题组的目标是研究 DDT，C1 连爆轰都没发生就无从谈起。于是我开始思考：水蒸气并不参与放热反应，只是稀释剂，且比热容高，像“能量黑洞”一样吸走热量、降低火焰温度。干脆把水去掉！剔除 H₂O 后得到纯氢＋干燥空气的新工况，称为 C1d（氮气仍是稀释剂，但已无“污染”）。为做对比，又测试了纯 H₂/O₂体系，命名为 C1p。

数值方法   

要同时考虑气相（均相）与壁面（非均相）化学反应对 FA & DDT 的影响，我们选用 CONVERGE 进行数值分析。爆轰问题需要高分辨率网格和自适应加密（AMR）来捕捉火焰面，而 CONVERGE 的自动网格生成与均相/非均相表面化学反应功能的耦合恰好能满足我们的需求。

具体实现： 

- 用 CONVERGE 求解多组分反应 Navier–Stokes 方程，采用有限体积法（FVM）保证质量、动量、总能量及组分质量分数守恒。  

- 单元面速度采用混合格式（中心＋迎风）通量混合；压力-速度耦合采用用 PISO 算法，Rhie-Chow插值算法抑制网格的“棋盘”压力振荡。 

- 压力泊松方程采用 CONVERGE的 BiCGSTAB 线性求解器，直接解耦压力场。  

- 气相与表面反应均采用 SAGE 详细化学反应求解器，同步计算表面覆盖率与气相质量分数，实现壁面耦合反应。  

AMR 策略：  

- C1 火焰厚度 δf≈700 μm，未启用 AMR。  

- C1d δf≈27 μm，在基础网格上加 AMR，确保火焰面厚度内至少 6 个网格。 

- C1p δf≈20 μm，进一步加密至 16 个网格跨越火焰面，充分解析火焰结构。  

验证工作   

为增强结果可信度，我们对 CONVERGE 的气相与表面反应机理进行了系统验证。例如，将 CONVERGE 的气相 SAGE 详细化学反应与流场耦合结果与开源燃烧 CFD 代码 PeleC 对比，验证结果见图 1。

图 1. (A) CONVERGE 与 PeleC 结果对比：不同时刻的火焰前锋变化，以温度着色显示。 (B) CONVERGE 与 PeleC 得到的火焰顶端位置随时间变化曲线。

CONVERGE 的表面化学反应模块与Chen 等人³的经典论文进行了验证。该文针对预混 H₂/空气在催化微管内的气相-表面耦合反应做了模拟，主要关注火焰稳定而非 FA/DDT。CONVERGE 的结果与文献数据吻合良好¹。

结果   

在 C1 算例中，火焰既未加速，也未形成爆轰，而是以恒定火焰速度平稳传播。然而，一旦启用表面化学（即显式考虑壁面反应），传统抛物线形的火焰锋面便发生“倒置”，如图 2 所示：火焰优先贴着壁面传播，这正是催化表面化学作用的表现。

图2. 温度云图展示C1算例在气-固耦合化学作用下火焰面的演化过程

然而，当关闭表面化学反应后，火焰重新恢复为传统的抛物线形锋面，如图3所示。

图3. 温度云图展示C1算例仅考虑气相化学时火焰面的演化过程

在发现沿壁面有大量中间自由基 OH 和 O 生成后，我们得出结论：催化表面反应通过产生活性中间体，直接强化了气相燃烧，从而使火焰优先贴着壁面传播。换言之，燃料/氧化剂混合物在壁面发生表面反应并释放中间自由基，推动火焰沿催化壁面前进。此外，启用表面化学的 C1 算例温度分布普遍高于仅用气相化学的情况，这是因为表面反应额外释放了热量。

所有 C1 算例中火焰均未加速，归因于混合物中的稀释剂和污染组分降低了火焰速度，抑制了 FA/DDT⁴。因此，我们按同一流程对 C1d 混合物进行了模拟。剔除H2O后，火焰速度提高并出现明显加速，但仍未发展为爆轰。与受污染的 C1 不同，C1d 仅在“仅启用表面化学、关闭气相化学”时出现火焰倒置；一旦允许气相反应，锋面又恢复抛物线形。各 C1d 算例的火焰最终加速至约 0.1 Ma，却因驻留时间过短，壁面催化效应来不及与气相耦合，未见贴壁传播。C1d 出现了快速 FA，但所需 DDT 酝酿距离不足，未能转捩为爆轰。

相比之下，C1p算例在初期形成“郁金香”形锋面后迅速完成 DDT。火焰两分支优先沿壁面传播，最终汇合形成爆轰波，如图 4 所示。

图4. 温度云图展示C1p算例（仅气相化学）中DDT（爆燃转爆轰）过程的发展。

最后简要总结一下：Surya 团队携手埃克森美孚科技研究人员，利用 CONVERGE 对三种不同燃料混合物（H₂/O₂与 H₂/空气）在镍催化微通道内的火焰传播与加速过程进行了模拟。各混合物因气-固化学耦合程度不同，在火焰加速（FA）与爆燃转爆轰（DDT）上呈现出显著差异。更多细节请阅读 Surya 的论文！

总体而言，本研究首次在该领域内同时考虑了镍催化微通道中的气相与表面耦合反应，并以此评估 DDT，为今后针对工业场景、提升爆炸安全性的更精细化研究提供了新方向。

参考文献  

[1] Ramachandran, S., et al. “Flame Acceleration and Deflagration to Detonation Transition in a Microchannel with Catalytic Nickel Walls.”Physics of Fluids, 36(11), 2024, 116-143. https://doi.org/10.1063/5.0235540 

[2] Zhang, S., Zhao, X., and Bayyuk, S., “Generalized Formulations for the Rhie–Chow Interpolation.”Journal of Computational Physics, 258, 2014, 880–914. https://doi.org/10.1016/j.jcp.2013.11.006 

[3] Chen, G.-B., et al. “Effects of Catalytic Walls on Hydrogen/Air Combustion inside a Micro-Tube.”Applied Catalysis A: General, 332(1), 2007, 89–97 https://doi.org/10.1016/j.apcata.2007.08.011 

[4] Ramachandran, S., Srinivasan, N., Wang, Z., Behkish, A. and Yang, S., “A Numerical Investigation of Deflagration Propagation and Transition to Detonation in a Microchannel With Detailed Chemistry: Effects of Thermal Boundary Conditions and Vitiation.” Physics of Fluids, 35(7), 2023, https://doi.org/10.1063/5.0155645