一种新型仙人掌形状反应器用于氨裂解：结合数值模拟与响应面优化方法

《ENERGY CONVERSION AND MANAGEMENT》：A novel cactus-shaped reactor design for ammonia cracking: Using an integrated numerical modelling and response surface optimisation method

　　

通过催化分解从氨（NH₃）中实现最大氢气（H₂）产量且能耗最小仍是该领域的主要挑战。在本数值研究中，研究人员设计并优化了一种新型反应器。该反应器几何形状受巨型仙人掌（Mighty Saguaro cactus）启发。这种分形形状结构增加了传热表面积并防止反应器内出现冷点。设计参数包括几何特征如分支角度和直径，以及操作条件如氨流量、进口温度和反应器壁温。采用实验设计（DOE）方法生成设计点，并通过基于有限体积法（FVM）的流场、传热和传质模拟评估每个点的反应器性能。目标函数定义为最大氢气（H₂）摩尔分数、最低操作温度和最小压降。通过将数值结果与响应面方法（RSM）耦合，预测了最优设计点，并通过数值模拟验证了预测的反应器在建议最优条件下的性能准确性。最优点的敏感性分析表明，氢气（H₂）产量对进口速度（?83.73%）和反应器壁温（9.57%）最为敏感。在最优仙人掌形状反应器中，转化效率达到88.12%，而在相同操作条件下，传统反应器仅达到74%效率。尽管仙人掌形状反应器增加了克服压降所需的泵功率，但最优条件下氢气（H₂）产量的相应提高使这种权衡可接受且合理。总体而言，当前优化的新型反应器使用22.36 kWh（包括加热负荷和泵功率）每天可生产2.93 kg氢气（H₂）。

现代社会对可持续能源的需求日益增加，氢气（H₂）作为无碳能源载体备受关注，但其低体积能量密度使储存和运输面临挑战。氨（NH₃）因高氢含量和易液化特性被视为理想氢载体。氨裂解生成氢气和氮气（N₂）是强吸热反应（ΔH=46 kJ/mol），需在400–500°C下通过催化剂进行。传统固定床反应器中，催化剂导热率低（如商用钌基催化剂ACTA 10010导热率0.1 W/m·K），导致径向温度梯度大，形成冷点，限制转化效率。工程上采用多根小直径反应器并联虽能缓解问题，但增加系统复杂性和成本。已有研究通过加装翅片或分形结构强化传热，但针对吸热氨分解反应器的仿生设计仍较少。因此，研究人员受巨型仙人掌（Mighty Saguaro cactus）几何形状启发，设计了一种新型反应器，旨在增加传热面积、改善温度均匀性，提高氨裂解效率并降低能耗。该研究发表于《ENERGY CONVERSION AND MANAGEMENT》。

研究人员采用三维数值模拟，基于有限体积法（FVM）求解流场、传热和传质控制方程。催化剂床层视为均匀多孔介质，使用Ergun方程计算动量损失，有效导热系数按体积平均求得。氨分解反应采用一阶幂律动力学（Tamaru模型）。通过实验设计（DOE）和拉丁超立方采样（LHS）生成设计点，参数范围包括分支直径（d_b）0.013–0.0254 m、分支角度（θ）35°–60°、氨流量1.29–10.66 Nm3/h、进口温度673–773 K、壁温673–773 K。目标函数为最大氢气摩尔分数、最小压降和最低操作温度。数值结果与响应面方法（RSM）耦合，建立二次多项式代理模型进行多目标优化，并通过CFD验证最优解。

主要关键技术方法包括：有限体积法（FVM）求解连续性、动量、能量和物种输运方程；多孔介质模型采用Ergun方程（孔隙率0.4，颗粒直径3 mm，导热率0.1 W/m·K）；氨分解采用一阶幂律动力学（A=823 s^-1，E_a=2.85×10⁷ J/kmol）；实验设计（DOE）结合拉丁超立方采样（LHS）生成25个设计点；响应面方法（RSM）建立代理模型进行多目标优化；CFD验证最优解。

研究结果如下：

**4.1 方差分析（ANOVA）**：评估响应面模型拟合质量。氢气（H₂）摩尔分数模型R2=0.994，调整R2=0.974，表明预测与模拟高度一致。压降模型R2=0.963，操作温度模型R2=0.940，反映其强非线性。基于编码系数的全局敏感性分析显示：氢气产量对进口流量敏感性为-83.76%，对壁温敏感性为9.57%；压降对进口流量敏感性为64.94%；操作温度对壁温、进口温度和分支直径敏感性分别为46.37%、28.96%和17.74%。分支角度对各目标函数影响不显著。

**4.2 优化**：通过RSM预测最优设计点：分支直径0.0133 m，分支角度35°，氨流量1.29 Nm3/h，进口温度673 K，壁温673 K。CFD验证结果：氢气摩尔分数预测误差3.52%，压降误差13.7%，操作温度误差1.95%，均在可接受范围。最优条件下转化效率88.12%，压降3262 Pa，操作温度489.14 K。

**4.3 流场**：氨（NH₃）从底部进入，在分支处因温度升高分解加快，氢气（H₂）和氮气（N₂）摩尔分数升高。与传统管式反应器相比，仙人掌形状反应器中冷点区域显著减少，分支处局部转化率更高。

**4.4 轴向演变**：沿中心线温度分布显示，管式反应器温度单调下降至约0.52倍进口温度；仙人掌反应器分支处温度恢复并保持较高（出口约0.93倍）。密度分布相应变化。氢气摩尔分数沿轴向在分支处急剧上升，中心线平均氢摩尔分数达0.649，而管式中心线仅0.389。但整体转化效率提升12.8个百分点（从74%到88.12%），低于中心线提升幅度，因分支通道减小导致流速增加、停留时间降低。

**4.5 转化效率、能量效率和泵功率**：最优流量（1.29 Nm3/h）下，仙人掌反应器转化效率88.12%，能量效率74.1%，均高于管式反应器（74%和62.1%）。随流量增加，两种反应器效率均下降，高流量（6.45 Nm3/h）时仙人掌反应器能量效率略低于管式，因流速加快缩短停留时间。泵功率分析表明，最优条件下仙人掌反应器泵功率2.393 W，远低于高流量时的458 W，低流量操作能平衡转化效率和泵耗。

总结讨论部分：研究结果表明，仙人掌形状反应器通过分支结构有效增强传热，减少冷点，在最优条件下显著提升氨裂解转化效率。主要结论包括：（a）氨进口流量对氢气产量和压降影响最显著，氢气产量对进口流量敏感性为-83.76%，对壁温敏感性为9.57%；（b）分支角度对目标函数影响不显著；（c）操作温度主要受壁温、进口温度和分支直径影响；（d）最优点转化效率88.12%，压降3262 Pa，操作温度489.14 K，预测误差可接受；（e）与传统管式反应器相比，优化后效率从74%提升至88.12%，但泵功率增加；（f）随流量增加，转化效率和能量效率下降，高流量下仙人掌反应器能量效率甚至低于管式。该反应器在低流量操作窗口优势显著。尽管研究基于数值模拟，但实验验证可行，主要挑战包括分支结构的高精度制造和壁温控制。