Page 86 - 《真空与低温》2026年第2期
P. 86

第  32 卷    第  2 期                          真空与低温
                 2026 年 3 月                           Vacuum and Cryogenics                                205



                      正   仲  氢   转   化   器  双   重   介   质  渗   流   模   型   的   IMPTECR         算   法



                                                              2
                                             薛 璞    1,2 ,张启勇 ,刘世贸     1,2 ,刘志帆  2*
                    (1. 中国科学技术大学,合肥 230026;2. 中国科学院合肥物质科学研究院,合肥 230031)


                     摘要:为研究多孔介质中正仲氢催化转化过程的传热传质特性,基于双重介质模型分别构建颗粒孔隙与堆床
                  空隙的质量守恒方程和能量守恒方程,将关键物理效应转化为控制方程组中的源项进行表征。针对该过程压力
                  场、温度场、组分比例场不同的特征时间,仲氢(或正氢)组分比例场变化更为剧烈,控制方程中的非定常项易造
                  成数值弥散或振荡。论文引入 IMPES 算法相似的隐式-显式分步求解思想,提出隐式压力及温度-显式组分比例
                  (IMPTECR)算法,核心目的在于分离非定常项得到形式上类似于稳态问题的压力场和温度场控制方程,采用“压
                  力场隐式求解-温度场隐式求解-组分比例场显式求解”的分步求解策略,再结合控制容积法对控制方程进行离散
                  求解。数值算例表明,IMPTECR 算法在绝热与等温工况下均表现出良好的收敛性和稳定性,可避免数值振荡与
                  误差累积。基于      Hutchinson 等报道的正仲氢催化转化实验数据,对所建立的数值求解模型进行了验证,通过对比
                  模型预测值与实验实测值,证实了该数值求解模型的准确性。研究为正仲氢转化器的数值模拟及优化设计提供
                  了可靠算法。
                     关键词:正仲氢转化;双重介质模型;数值模拟;多孔介质流动
                     中图分类号:TB61                       文献标志码:A       文章编号:1006-7086(2026)02-0205-08
                     DOI:10.12446/j.issn.1006-7086.2026.02.012


                     IMPTECR Algorithm for the Dual-porosity Seepage Model of Ortho-parahydrogen Converter


                                               1,2            2          1,2        2*
                                        XUE Pu ,ZHANG Qiyong ,LIU Shimao ,LIU Zhifan
                                 (1. University of Science and Technology of China, Hefei 230026, China;
                           2. Hefei Institutes of Physical Science,Chinese Academy of Sciences, Hefei 230031, China)


                     Abstract:To explore the heat and mass transfer laws of ortho-para hydrogen catalytic conversion in porous media, this
                  study establishes mass and energy conservation equations for particle pores and bed voids separately based on a dual-porosi-
                  ty medium model. Key physical effects in the conversion process, such as catalytic reaction kinetics and interphase heat/mass
                  transfer,  are characterized as source terms in the governing equations,  enabling accurate description of complex physical
                  processes in the multi-scale pore structure. Aiming at the core problem in the multiphase seepage system of ortho-para hydro-
                  gen converters—significant time-scale differences among the pressure field, temperature field, and component ratio field,
                  with the ortho/para hydrogen ratio field changing sharply and its steep gradient regions easily causing numerical
                  dispersion and oscillations—this study introduces the implicit-explicit stepwise solution idea of the IMPES algorithm
                  and proposes the Implicit Pressure and Temperature-Explicit Component Ratio(IMPTECR) algorithm. Through equivalent
                  algebraic transformation, the algorithm eliminates the time derivative term of the component ratio, and adopts a stepwise
                  strategy of "implicit solution for pressure field-implicit solution for temperature field-explicit solution for component ratio
                  field" to effectively decouple strongly coupled unknowns. Subsequently, the control volume method is used for discretiza-
                  tion and numerical solution of the governing equations. Numerical results show that the IMPTECR algorithm exhibits excel-
                  lent convergence and numerical stability under both adiabatic and isothermal conditions, which can effectively suppress nu-
                  merical oscillations and error accumulation. Validated against the ortho-para hydrogen catalytic conversion experimental da-
                  ta reported by Hutchinson et al., the model predictions of conversion rate and bed temperature distribution are in good agree-


              收稿日期:2026−01−04
              基金项目:聚变堆主机关键系统综合研究设施项目(CRAFT)(2018-000052-73-01001228)
              作者简介:薛璞,硕士研究生。E-mail:xuepu@mail.ustc.edu.cn
              通信作者:刘志帆,副研究员。E-mail:zhifan.liu@ipp.ac.cn
   81   82   83   84   85   86   87   88   89   90   91