射流悬浮床气固两相流动的数值模拟.pdfVIP

第 30 卷 第 2 期 北 京 化 工 大 学 学 报 Vol. 30 , No. 2 2003 年 JOURNAL OF BEIJ ING UNIVERSITY OF CHEMICAL TECHNOLOGY 2003 射流悬浮床气固两相流动的数值模拟 杨　烨　王翠苹　李定凯　吕子安 (清华大学煤的清洁燃烧技术国家重点实验室 ,北京　100084) ε 摘　要 : 从多相流多流体模型、气相湍流 k 模型与颗粒湍流 kp 模型出发 ,对底部向上射流悬浮床内气固两相流动 进行数值模拟 ,得到了不同入口气速下各截面的颗粒相速度场、颗粒相脉动速度场的分布 ,结果与实验定性一致。计 算结果与实验值对比的结果表明 ,数值模拟对进一步的研究有指导意义。数值结果与实测值间存在的差异则有待于 进行直接三维流动模拟、进一步考察更为合理的边界条件及应用更先进的模型来解决。 关键词 : 底部向上射流 ;气固两相流 ;多流体模型; kp 模型 中图分类号: O359 ; TK229. 6 + 6 ε 中可舍去时间项。在采用 k 模型描述气相湍流特 引　言 性 ,并在模型方程中加入固相的影响因素后可以得 在应用循环床进行脱硫的实践中 ,为了降低系 到气相时均控制方程组为 统阻力 ,将传统的布风板改为大孔射流板 ,直接将含 ( ) 1 气相连续方程 灰的烟气射入反应塔内。由于反应塔内固体颗粒较 ( ρ) ( ) 细 ,尺寸比较均匀 ,可使床内的颗粒分布上下比较均 xj vj = 0 1 匀 ,形成悬浮床。入口结构改进后的效果及进一步 气相动量方程 优化的方向需要对反应塔内气固两相流动特性进行 p vi vj (ρ ) μ 研究后才能确定。目前塔内的气固两相流动特性已 xj vj vi = - x x + xj e xj + x i - f p i [1 ] ( ) 经有实验研究的结果 可供参考。为了克服实验 2 研究条件的局限性 ,更深入地了解条件变化对反应 k 方程 μ 塔内气固两相流动特性的影响 ,有必要进一步对射 ( ρ ) e k ρε ( )