天然气重整器内部传递过程的数值模拟.docxVIP

· 121 · 第28卷第3期 天 然 气 工 业 维普资讯http://www.eqvip.cot 加工利用与安全环保 加工利用与安全环保 天然气重整器内部传递过程的数值模拟* 杨国刚1 岳丹婷1 袁金良2 (1.大连海事大学轮机学院 2.瑞典兰德大学热能工程系) 杨国刚等.天然气重整器内部传递过程的数值模拟，天然气工业，2008,28(3):121-123. 摘 要 天然气通过重整反应可以生成氢气，成为燃料电池氢原料的重要来源。目前紧凑型天然气重整器的 开发设计，主要集中在材料方面的研究，而内部机理方面的研究较少。对天然气重整器内部流动与传热现象进行 了模拟与分析。分析中考虑了化学反应的影响，采用了耦合的边界条件以及可变的热物性参数，建立了描述流通 管道和多孔催化剂层内部流动与传热现象的三维数学模型。采用 SIMPLE 算法对模型进行求解，得到了化学反应 速率、混合气浓度以及温度等的空间分布。结果表明，化学反应限制在厚度约为1 mm 的催化剂薄层内进行，沿主 流动方向，CH; 浓度由19.5%降到7%左右，H2 浓度由4.1%上升到13%左右，混合气温度上升约10℃。研究结 果对天然气重整器的开发、结构设计具有参考意义。 主题词 燃料电池 天然气 重整器 传热 传质 数值模拟 天然气(甲烷占97%以上)通过重整反应生成的 氢气，可成为燃料电池氢原料的重要来源。近年出 现的紧凑型天然气重整器由燃烧与重整两个平行管 道组成，具有热传导的距离以及气体扩散路径短的 优点，降低了重整器中的传热、传质阻力。因此，该 重整器不仅能量利用率较高，也使得体积大大缩 小[12]。目前紧凑型天然气重整器的内部机理方面 的研究较少[3-]。笔者将对单个重整管道的不同传 递过程进行三维模拟，不涉及燃烧管道，燃烧管道对 重整管道的影响将采用热边界条件来描述(相邻燃 烧管道固壁有热量流入，其他为绝热条件)。 一 、过 程 描 述 与 数 学 模 型 重整器结构见图1。取一半宽度管道进行研究， 坐标原点设定在复合管道的左下角处，坐标x、y、z 分别表示管道的长、高、宽，U、V、W 分别为对应3个 方向的速度分量 重整的化学过程用2个独立的平衡式来表示： CH?+H?O→CO+3H?(△h10ook=226 kJ/mol) ( 1) CO+H?O→→CO?+H?(△h100k=—35 kJ/mol) ( 2) 燃烧反应提供的热量(9) 图 1 紧凑型重整器结构示意图 式中：△h100ok表示反应温度为1000 K 时的反应热， 方程(1)称为水蒸气重整反应，为吸热反应；方程(2) 称为水一气变换反应，为放热反应。总反应吸热，热 量由燃烧管道内的催化燃烧反应提供。 控制方程包括连续方程、动量方程、能量方程等。 连续方程： ▽ (peftv)=0 ( 3) 动量方程：▽(pefw)=- ▽P+ ▽(μeff▽v)+Sai (4) 式中：paii为有效密度，v 为混合物的体积平均速度张 量；P 为压强；源项 Su 的计算公式[5.6]为： Sa =(μev/βat)-pahBV; |v | (5) *本文受到国家自然科学基金资助项目(编号、辽宁省2007年博士启动基金的资助。 作者简介：杨国刚，1972年生，副教授，工学博士；毕业于大连理工大学，主要从事气体流动与传热现象的数值模拟研究； 发表文章10余篇。地址：(116026)辽宁省大连市凌海路2号。电话：(0411E-mail:yanggg168@so- 加工利用与安全环保 加工利用与安全环保 天 然 气 工 业 维普资讯hitp// 2008年3月 式中： μd 、P 、B 分别为有效黏度、多孔层的有效渗 透率、公式常数。 能量方程表示为： ▽(pantCpeftvT)= ▽(kei ▽T)+Sr (6) Sr =R:△h;+R,△h 。 (7) 式中：Sr 为与重整反应有关的热源项；R 为水蒸气 重整反应与变换反应的摩尔速率；△h 为反应焓。 组分守恒方程采用以下的统一形式： ▽(pervp)=▽(penDert▽φ)+S, (8) 这里，φ为燃料气体的质量分数。 边界条件为： (1)在固壁底部 (y=0) 处 (φ为CH?、H?、CO、H?O、CO?) (2)在固壁顶部与侧面处 U=V=W=0,q=0,J,=0 (3)在平板中心位置(x=0) 处 采用控制容积积分法对数学模型进行离散，对流 项采用二次