高浓度颗粒二阶矩方法及提升管流化特性的模拟研究.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：096171 高浓度颗粒二阶矩方法及提升管流化 特性的模拟 1 2 1 1 2 1 孙丹 ，王建志 ，N. L. MBOUANA , 王帅 ，陈明 ，陆慧林 （1.哈尔滨工业大学，能源科学与工程学院，哈尔滨市 150001 ；2. 中船重工第703 研究所，哈尔滨市 150036） (Tel:0451E-mail: sundan@hit.edu.cn ) 摘要：采用颗粒动理学方法，考虑颗粒速度脉动各向异性，建立颗粒相二阶矩模型。应用初等输运理 论，对三阶关联项进行模化和封闭。考虑颗粒与壁面之间的能量传递和交换，建立颗粒相边界条件模 型。数值模拟提升管内气固两相流动特性，预测颗粒浓度分布与 Benyahia 等实测结果相吻合，模拟得 到提升管内气固两相流动典型的环－核结构，模拟结果表明提升管内颗粒相湍流脉动具有明显的各向 异性。 关键词：二阶矩，各向异性，颗粒动理学，提升管，数值模拟 0 前言 流态化技术被广泛用于化工、能源、环境以及颗粒干燥等领域。循环流化床锅炉， 石油催化裂化装置等快速流态化设备，都是密相提升管工业应用的代表。提升管内的气 固两相流动中，颗粒在管中心处浓度较底，颗粒向上流动，而在壁面附近颗粒浓度较大， 颗粒向下流动，在提升管内形成颗粒的内循环流动。这种环－核结构造成提升管内气固 两相流动的不均匀性，增大了提升管内气固两相流动控制难度，是工业应用的难点，也 是提升管内气固两相流研究的核心。 Ibsen 等[1]在 12MW 循环流化床锅炉的 1/9 比例的试验台上，使用激光多普勒测速 仪测量了炉内颗粒的运动规律，实验统计所得颗粒雷诺应力在循环流化床内，尤其是锅 炉底部，具有明显的各向异性。Tartan 等[2]应用离子图像测速法技术测量了提升管内的 颗粒速度分布，测量结果表明，对于颗粒脉动速度二阶矩和雷诺应力型二阶矩，轴向分 量远远大于其在径向及切向上的二阶矩分量。Bhusarapu[3]采用放射性粒子自动跟踪技术 测量了提升管内的颗粒运动，测量结果表明，快速流态化下颗粒的脉动的各向异性明显 高于稀疏气粒输送状态。快速流态化状态下，颗粒的雷诺应力二阶矩轴向和径向的比值 为 1.5-2.5，而切向雷诺应力二阶矩较小，仅占径向雷诺应力二阶矩的 1/10 左右。 稠密气固两相流动的提升管内，颗粒运动呈现明显的各向异性。目前所应用的颗粒 动理学模型是在颗粒脉动为各向同性的假设基础上建立的[4]，无法预测颗粒脉动的各向 异性。本文基于稠密气体分子动力学，考虑颗粒间的碰撞作用导致的脉动各向异性，建 立高颗粒浓度下颗粒动理学模型；提出颗粒相的动量和速度二阶矩的输运方程；采用输 *基金项目： 自然科学基金项目的资助 运现象的初等原理对高阶的关联项进行封闭。应用各向异性颗粒动理学模型数值模拟提 升管内气固两相流动，研究提升管内颗粒的流动规律及其脉动的各向异性行为。 1 数学模型和边界条件 1.1 各向异性颗粒动理学 仿照分子动力学理论，定义颗粒的速度分布函数f (c, x, t)表示 t 时刻在点 x 附近dx 邻域内出现速度范围为(c, c + dc) 的几率。定义C = c - u 为颗粒脉动速度，则颗粒脉动速 度的 N 阶矩为： N 1 (1) M x t CCL C CCL C f c x t dc