基于扩散近似法与双通量法耦合的多孔介质内燃烧辐射换热计算方法.pdf

第 30卷 第 2期 计 算 物 理 Vo1．30．No．2 2013年 3月 CHINESE JOURNAL OF COMPUTATIONALPHYSICS Mar．．2013 文章编号 ：1001—246X(2013)02-0203-06 基于扩散近似法和双通量法耦合的多孔介质 内 燃烧辐射换热计算方法 李德波 ， 宋正昶 (1．广东 电网公司电力科学研究院，广州 510060 2．中国矿业大学 电力学院 ，徐 州 221116) 摘 要 ：将 Rosseland辐射模型和双通量法结合 ，对瓦斯气体在泡沫陶瓷中预混燃烧辐射换热计算方法进行研究 数值模拟 的结果与试验结果吻合 比较好 ，证明提出的辐射换热计算方法是有效的． 关键词 ：辐射换热计算方法 ；Rosseland扩散近似法 ；双通量法 ；瓦斯预混燃烧 ；多孔介质 中图分类号 ：TK124 文献标识码 ：A O 引言 多孔介质中的预混合燃烧是一个包含化学反应及导热、对流和热辐射三种换热方式互相耦合的复杂换 热过程 ¨ ．从理论上研究多孔介质预混合燃烧 ，建立燃烧的数学模型是众多研究人员努力的方 向．大部分 研究者都在实验研究的基础上 ，通过适当简化 ，建立燃烧数学模型，再通过数值计算 的方法，得到和实验基本 一 致的结果 ，验证燃烧模型的有效性． Kotani。。等人用小直径 的陶瓷管束作为研究对象，建立了数学模型，在模型中提出对气体和固体采用独 立的能量方程 ，这是他们对多孔介质预}昆合燃烧理论研究的重要贡献．但是模型中没有考虑辐射作用 ，而辐 射换热在多孔介质预混合燃烧 中对火焰位置 、火焰温度 、燃烧速度等重要参数有很大的影响．Churchill 的 数学模型考虑 了壁面之间的辐射作用 、壁面内的导热作用 、壁面与气体之间的对流作用以及层流到湍流的过 渡状态热交换的强化作用，并认为气体的辐射作用与固体相 比可以忽略． 对辐射传热的处理上 ，研究者大多采用适用于光学厚介质的Rosseland辐射模型，理论上讲 ，在多孔介质 燃烧器出口，光学厚介质 的假设是不满足的，因此 Rosseland辐射模型显然是不适用 的；有的采用双通量辐射 模型，在处理多孔介质燃烧器 内部计算区域时显得过于复杂．本文将 Rosseland辐射模型和双通量法结合 ，克 服各 自的缺点，对瓦斯气体在泡沫陶瓷 中预混燃烧辐射换热计算方法进行研究． 1 数学模型 对于泡沫陶瓷的燃烧过程 ，气 、固相的传热能力有 明显差异 ，气体燃料燃烧释放 出来 的热量不可能立刻 完全传递给 固体基质 ，因此两相之间存在局部温差 ，即二者处于局部非热平衡状态 ，应分别建立能量输运方 程 ，并通过两相之间的表面对流换热系数将这两个方程耦合起来 ，这一点与普通的预混燃烧不同．本文研究 引入多孔介质 的体积孔 隙率 以描述 固相介质 的存在对燃烧 的影 响．基于上面 的假设和分析 ，参考 Williams 和周力行 ¨。。的论著 ，根据多组分反应流体的基本定律和基本方程，建立控制方程 ： 连续性方程 ： 0 ， (1) a 式 中，p为气体混合物的密度 ；“为气体流速 ；s为多孔介质孔隙率 ； 为轴 向距离． 投稿 日期 ：2012—04—23；修回 日期 ：2012—07—25 基金资助 ：国家 自然科学基金重点项 目资助项 目 作者简介 ：李德波 (1983一)，男 ，土家族 ，湖北宜昌，工程师 ，工学博士 ，主要从事煤粉燃烧污染物控制 ，百万超临界燃煤机组调试 、试验，煤 粉燃烧高级数值模拟，大规模并行数值计算方法和程序开发等研究，E-mail：ldbyx@126．eom 计 算 物 理 第 3O卷