Page 92 - 《真空与低温》2026年第2期

薛 璞等：正仲氢转化器双重介质渗流模型的 IMPTECR 算法 211

0.50 hydrogen conversion during liquefaction[J]. Chemical Engi-
实验数值
0.45 仿真数值
neering & Technology，2019，42（3）：669−679.
0.40
[6] WANG K， XU S， ZHENG J， et al. Experimental study on
0.35
catalytic conversion and flow characteristics of ortho-para
0.30
p-H 2 /− 0.25 hydrogen in tubular packed-bed converters at liquid nitrogen
γ B
0.20 temperature[J]. Cryogenics，2025，147：104041.
0.15 [7] WHITAKER S. The method of volume averaging： Vol. 13
0.10
[M]. Dordrecht：Springer Netherlands，1999.
0.05
[8] KUWAHARA F，SHIROTA M，NAKAYAMA A. A numer-
0
0 200 400 600 ical study of interfacial convective heat transfer coefficient in
u v /min −1
two-energy equation model for convection in porous media[J].
图 7 仿真数据与实验数据的对比 International Journal of Heat and Mass Transfer，2001，44（6）：
Fig. 7 Comparison of simulation data and experimental data 1153−1159.

[9] 代盟，刘志帆，张启勇. 基于双重介质模型的正仲氢转化器
3 结论
流动规律研究 [J]. 低温与超导，2023，51（4）：33−39.
针对正仲氢转化器双重介质渗流模型求解，采 [10] （美）哈利德·阿齐兹（AZIZ K），（加）安东尼·塞特瑞（SET-
用 IMPTECR 算法的模拟过程展现出良好的计算 TARI A）. 油藏数值模拟 [M]. 袁士义，王家禄，译. 北京：
收敛性，可稳定输出各物理量的连续分布特征；而 石油工业出版社，2004.

采用分离式求解方法时，受方程中组分比例时间导 [11] MONTEAGUDO J E P，FIROOZABADI A. Comparison of
数项的影响，数值稳定性显著下降，不仅出现明显 fully implicit and IMPES formulations for simulation of
的振荡现象，且误差随计算步长累积放大，最终导 water injection in fractured and unfractured media[J]. Inter-
致模拟结果与实际物理过程产生显著偏差。并基 national Jaurnal for Numerical Methods in Engineering，
于 Hutchinson 等发表的实验数据，通过多组工况下 2007，69（4）：698−728.
的模拟值与实测值对比分析，验证了本研究数值求 [12] 梁润红. 多孔介质中不可压多相流模拟的改进 IMPES 算
解模型的准确性与适用性。在未来的研究中，作者 法 [D]. 贵阳：贵州大学，2021.
将运用该算法对正仲氢转化器催化剂堆床的传热 [13] LEE S H，TENE M，DU S，et al. A conservative sequential
传质过程展开深入探究，以期为氢液化器设计提供 fully implicit method for compositional reservoir simulation
理论支持。 [J]. Journal of Computational Physics，2021，428：109961.
[14] 朱文昌. 多孔介质中单相多组分可压渗流的数值模拟研
参考文献：
究 [D]. 贵阳：贵州师范大学，2023.
[1] 付强，杨洸，金辉，等. 中国氢能产业链技术现状及发展趋 [15] 熊峰，姜清辉，陈胜云，等. 裂隙-孔隙双重介质 Darcy-
势 [J]. 油气与新能源，2024，36（4）：19−30. Forchheimer 耦合流动模拟方法及工程应用 [J]. 岩土工程
[2] 徐硕，余碧莹. 中国氢能技术发展现状与未来展望 [J]. 北 学报，2021，43（11）：2037−2045.
京理工大学学报（社会科学版），2021，23（6）：1−12. [16] 刘军文. 基于体积平均法对高渗多孔介质渗流规律研
[3] WILHELMSEN Ø，BERSTAD D，AASEN A，et al. Reduc- 究 [D]. 合肥：中国科学技术大学，2019.
ing the exergy destruction in the cryogenic heat exchangers [17] SINGLETON A H，LAPIN A，WENZEL L A. Rate model
of hydrogen liquefaction processes[J]. International Journal for ortho-parahydrogen reaction on a highly active catalyst
of Hydrogen Energy，2018，43（10）：5033−5047. [C]//Timmerhaus K D. Advances in cryogenic engineering.
[4] PARK J，LIM H，RHEE G H，et al. Catalyst filled heat ex- Boston，MA：Springer US，1995：409−427.
changer for hydrogen liquefaction[J]. International Journal of [18] BLIESNER R M. Parahydrogen-orthohydrogen conversion
Heat and Mass Transfer，2021，170：121007. for boil-off reduction from space stage fuel systems[D]. Wa-
[5] DONAUBAUER P J， CARDELLA U， DECKER L， et al. shington State University，2013.
Kinetics and heat exchanger design for catalytic ortho-para [19] LEACHMAN J W，JACOBSEN R T，PENONCELLO S G，