清华大学煤的清洁燃烧技术国家重点实验室评估材料.pptVIP

* 湍流气固两相流体力学 ——理论与数值模拟 祁海鹰 清华大学煤的清洁燃烧技术国家重点实验室 北京 2003.3.5 主题 气固两相湍流模型 流化床曳力理论模型 数值模拟 概述 气 液 固 气液两相流动 ——锅炉汽水循环 气固两相流动 ——沙尘暴 固液两相流动 ——泥石流 流体力学发展 ? …… ? 沙尘暴 沙漠迁移 锅炉炉膛内煤粉燃烧 鼓泡流化床反应器 循环流化床烟气脱硫装置 稀疏 稠密 湍流效应——流动非稳定性，影响 - 流动 - 燃烧 - 污染物形成 气固相互作用 流动极其复杂 概述 概述 预测 提供应用依据 寻求流动规律 合理描述 模型 数值模拟 深化认识 气固两相湍流模型 目的：正确描述两相湍流相互作用；合理预测燃烧过程和 污染物NOx排放规律 气固两相流动 力平衡关系——守恒方程组 (基于牛顿第二定律) t V 湍流： …… 瞬时 平均 脉动 未知关联项 湍流模型 在认识物理本质的基础上发展出的求解方法： 湍流模型系列 K-e -Kp模型——气相湍流动能-耗散-固相湍流动能 二阶矩模型——用更高阶的守恒方程描写和求解未知关联项。早于欧洲 统一二阶矩输运方程模型USM——各种不同的关联项用统一形式描述 USM＋气相湍流IPCM压力应变项封闭——改进气相湍流模拟的精度 USM拉氏改进模型——引入Lagrange概念分析颗粒惯性等作用 USM＋颗粒-壁面作用模型——计入壁面粗糙度的影响，提高颗粒雷诺应力预测精度 两相双尺度耗散USM模型——在湍流耗散项中引入2个湍流时间尺度，可根据实际流动变动 非线性NKP两相湍流代数模型——改进预测精度，存贮和计算量锐减，工程应用 ○ 实验 —— 二阶矩模型 (IPCM) GL … … IPCM＋壁面效应 图：旋转射流轴向速度分布（旋流数S＝0.53） ▲ 实验 —— NKP模型 USM模型 …… KP模型 颗粒相的轴向（上）和径向（下）脉动速度分布 0.0 0.1 0.0 0.1 0.2 0.0 0.1 0.2 0.0 0.1 0.2 1.0 0.8 0.6 0.4 0.2 0.0 r/R 1.0 0.8 0.6 0.4 0.2 0.0 0.0 0.1 0.2 0.0 0.1 0.2 0.0 0.1 0.2 0.0 0.1 0.2 r/R 2D矩形液柱模拟 平均垂直速度分布 (cm/s) （液相） （气泡相） 1100 液体 9个进气喷嘴 11.2 48.2 液相湍流动能(cm2/s2) 气泡体积分数 流化床曳力理论模型 曳力——流体对颗粒的拖拽力，气固两相相互作用 FD0 FD 根本原因——强烈的非均匀效应 浓相 稀相 气流 现状：国际上迄今为止未取得突破 （2002年） 曳力不准导致的计算结果 局部非均匀性——聚团流动 典型的非均匀流动 内循环 稀相区 外循环 密相区 整体非均匀性——上下和径向 局部非均匀性——颗粒团聚，曳力↓，颗粒可在高速气流中停留！ 计算域 气固速度滑移：Urel=Ug-Us，最大发生在环核交界处 现有曳力模型大多适于均匀流动。 团聚效应（实验） 模型4 ——唯一可计及团聚效应的经验 模型 团聚效应区：0es0.3， bmin 在 es=0.1，环核交界处 缺陷：受制于实验条件，无普 适性 β (Ns/m4) 颗粒体积分数es 模型1 模型2 模型3 模型4 欧拉方法暨“双流体”模型 曳力（Drag） 颗粒相 应力 脉动 团聚 壁面滑移 质量守恒 动量守恒 气相 质量守恒 湍流 动量守恒 动量守恒方程 时间项 对流项 静压项 扩散项 重力 相互作用项 用于单个孤立颗粒 - 用于多颗粒系统； - 经验常数，多来自于液固床，无普适性 守恒方程组没有完全封闭 曳力函数 综上所述 曳力是颗粒运动的根本驱动力； 曳力模型是准确描述稠密气固两相流动的关键； 曳力问题背后的机理：均匀与非均匀流动的差别； 非均匀流动的根本原因：颗粒团聚。 研究目标 发展新曳力模型 - 计及团聚效应； - 消除经验系数影响：模型普适性； - 真正封闭动量方程 研究思路 实验方法：气固困难 探索：计算流体力学(CFD) + 最小能量原理 ！ 分析非均匀流动特性 ? 采用理论方法合理描述曳力 最小能量原理 一种独立的能量分析理论 分析流体夹带颗粒运动和做功