颗粒在大尺度涡结构中弥散.pdfVIP

中国工程熟物理年会 多相流 2003年学术会议论文 编号：036056 颗粒在大尺度涡结构中的弥散 王兵，张会强，王希麟 清华大学工程力学系，北京100084 摘要：本文采用大涡模拟方法和颗粒轨道模型，分别研究了二维和三维后台阶气粒两相湍流流动的 大尺度涡结构，进而对不同入流速度的颗粒在大涡结构中的弥散进行了研究，重点研究了高速释放的 大颗粒在流动中的弥散特性。三维流动中大尺度涡结构具有明显的脱离、发展、台并和破碎过程，而 二维流动中大尺度涡没有表现出破碎过程。小颗粒在二维流动与三维流动中的弥敖稍有不同，二维呈 现出颗粒与流体处于平衡态，颗粒完全受大涡结构控制，三维流动中由于大涡的破碎，颗粒与大涡之 问不完全平衡，颗粒在大尺度涡的边缘呈现稍高的富集。高速释放到流场中的大颗粒受惯性影响大． 保持在其原有动量方向上运动。 关键词：颗粒弥散；大涡模拟：颗粒轨道模型 湍流两相或多相流动广泛存在于自然界以及工程应用中，诸如含尘空气流、含沙水 流、生物体的血液流动，粉料或粒料的气力或液力管道输运、粉尘的分离和收集、喷雾 干燥与喷雾冷却、流化床等等。颗粒在湍流流场中弥散机制复杂，因此，研究湍流两相 流动中颗粒的弥散机制，对于：弄清其运动规律以及强化混合和提高燃烧效率具有重要的 学术价值和工程价值。 关于不同入流速度下颗粒在大涡结构中的弥散行为的信息比较少。然而对于广泛存 在于自然界和生产、工程中的机器设备中的常见的扩散器，具有分离流动的螺旋桨，突 扩燃烧器，下水道处理场中的沉积池等等，需要深入了解颗粒在流场中的弥散行为．尤 其需要弄清当颗粒进入流场的入流速度不同时，颗粒的弥散行为。对于大部分两相流动， 流场中的颗粒大小不同，可以从几个微米到几百个微米，那么对于不同大小的颗粒在不 同入流速度时在流场中的弥散行为的研究也尤为需要展开。 本文采用大涡模拟研究气相场，用轨道法模拟颗粒相运动，在研究气相人涡结构的 基础上，重点揭示颗粒在大涡结构中的运动规律，并展开了不同大小，不同入流速度的 颗粒在三维流场大涡结构中的弥散行为的研究。 1基本方程 1．1连续相 进行直接模拟，对小于滤波宽度的小尺度涡运动则引入亚格子模型进行模拟。以后台阶 高度作为长度尺度，入口速度作为速度尺度，进一步对滤波方程进行无量纲化，得至Ⅱ连 续相大涡模拟的控制方程如下： 连续性方程： 国家自然科学基金资助项目(No 堕：0 Ox。 动量方程： 詈+毒c甄巧，=一古善+毒M考+詈卅鲁 。’ 式中，上标一’’表示大尺度分量，o是亚格子应力a对于亚格子应力本文采用了涡粘 性假设确定，如下式 乃刊要+争 (3) 其中，h是亚格子涡粘性系数，本文采用Smagorinsky模式，定义 蜥；p(CsA)2阼p(GA)2√雨 (4) ‘ 其中 墨=÷c要+争 (5) 模型中的常数C：取为O．17。 1．2离散相 本文研究中采用Lagrangian方法模拟颗粒的行为。几条基本假设为：颗粒为稀相： 颗粒对流体的作用和颗粒间的相互作用忽略不计；流场中的颗粒均为球形，且密度Pp和 直径瓯相同；忽略压差力、Saffman力；因为重点研究颗粒和大涡之间的作用，因此本 文目前忽略重力的影响。基于以上假设，颗粒的控制力主要为拖曳力，那么沿颗粒运动 轨道的无量纲方程为： 下dxp(t)=V， (6) 盟： !