提升管内气固双流体模型的计算模拟.pdfVIP

第1 卷 第3 期 过 程 工 程 学 报 Vol.1 No.3 2001 年 7 月 The Chinese Journal of Process Engineering July 2001 提升管内气固双流体模型的计算模拟 郑 雨 万晓涛 魏 飞 金 涌 (清华大学化工系 北京100084) 摘 要 对提升管内气固流动行为建立了k–e–k –e –Q 双流体模型. 利用计算软件CFX 模拟了内 p p 径0. 418 m 高18 m 的提升管内的冷态气固两相流动. 模拟结果与实验数据吻合较好. 关键词 提升管 气固流动 双流体模型 计算模拟 中图分类号 TQ021 文献标识码 A 文章编号 1009–606X(2001)03–0249–08 1 前 言 气固两相流的计算模型根据对颗粒相的处理不同可分为轨道模型和双流体模型. 双流体模型 基于拟流体假设 将颗粒处理为连续相. 颗粒动理学理论(Kinetic Theory of Granular Flow)类比稠 密气体的分子运动学理论而建立颗粒相的运动方程 认为颗粒相能量可能的耗散机制是由于颗粒 [1, 2] [3] 间的非弹性碰撞 并提出用“颗粒温度”的概念来反映颗粒相速度脉动 . 自Sinclair 首次将颗粒 动理学理论用于模拟提升管内的气固流动以来 这一理论结合多相流模型在气固流动的研究和模 拟的应用中得到了不断的发展[4–9] 但目前基本上还没有令人满意的气相湍动和颗粒相脉动的耦 合封闭. 本文在Cheng[9] 的k–e–k –Q 模型基础上进一步完善而建立了k–e–k –e –Q 气固湍动多相流模 p p p 型 利用计算软件CFX( 由英国AEA Technology 公司开发的流体单元模拟软件)模拟了内径0. 418 m 高 18 m 的提升管内的冷态气固流动 并将预测结果与实验数据进行了对照与分析 二者吻 合较好. 考虑到管流的对称性 计算域为轴对称二维截面的一半 采取 7020 网格 在 PC 机上 运算获得收敛结果需要5~8 h. 2 数学模型 2.1气固湍流模型 [10] [9] 周力行 提出k–e–k 模型 并在煤粉燃烧等稀疏颗粒湍流气粒体系中得到成功应用. Cheng p 将该模型与颗粒动理学理论结合 建立了 k–e–k –Q 模型 并成功地预测了循环流化床下行式反 p 应器内的气固流动行为. 本文在k–e–k –Q 模型的基础上 对固相的湍动模拟采用类比气相k–e模 p 型的推导方法 对固相守衡方程进行雷诺平均 建立了k–e–k –e –Q 模型 推导过程见文献[11]. p p 2.2柱坐标下的k–e–k –e –Q 模型方程组