三维槽道两相流颗粒运动的大涡模拟研究.pdfVIP

中国工程热物理学会 学科：多相流 学术会议论文 编号：056083 三维槽道两相流颗粒运动的大涡模拟 栗晶1柳朝晖2贺铸翁磊郑楚光 (华中科技大学煤燃烧国家重点实验室，湖北，武汉，430074)1 Tel：027—875424】7—8】07 Email：chongcaojin蛇000国163．COrn 摘要：利用基于动态粘性引稠格模型的大涡模拟方法，分别模拟了Re。=180的三维槽道湍流内的三种 颗粒，以及Re。分别为180和395时的塑料颗粒的运动状况。其中气相控制方程的求解采用分步投影方 法，泊松方程采用基十FFT的快速求解，颗粒的拉氏控制方程采用二阶龙格埚i塔积分。模拟统计得到 r流体和颗粒的平均速度．脉动速度等，不同的颗粒在流场中(特别是近壁区)表现出来了不同的行为， 并再现了颗粒在壁面附近局部富集的现象。艟Re提高，颗粒各方向上的速度脉动都有所增强。 关键宇：槽道，颗粒，两阳流，大涡模拟 1．引言 两相流动广泛存在于各种管道输运、化工反应器和燃烧装置中，如管道气力输运、 煤粉燃烧器、煤粉锅炉等，这类流动流场的准确预测对提高工程设计的水平极其重要。 在现实世界的工程中，流动通常是由壁面来限制和控制的。在近壁区域，流体一颗粒间 的耦合作用还会对流体的拟序结构及湍流发生机制产生深刻的影响。因此，在许多工程 系统中，详细了解壁面边界粒子的流动是很有必要的。 直接数值模拟(DNS)是对流场中的所有尺度进行求解，从而可以得到流场中的详 细信息，但是基于目前计算机技术的发展，DNS现在还仅限于低雷诺数的工况和规则几 何形状的流动问题。在大涡模拟方法(LES)中，对大尺度的物理量进行直接计算，仅对 小尺度量以及它们的相关量；}入亚网格模型，从而具有较好的准确性、通用性和较小的 计算量，因而大涡模拟正发展成为模拟湍流流动的有效方法。 本文将借助于大涡模拟方法，对槽道这一典型壁面约束流系统内发展中的气粒两相 和结果的正确性经行了证明，并详细了解和研究了壁面附近颗粒的运动特点，因此本文 的主要内容将是在其基础上进一步研究近壁区不同物性对颗粒分布的影响和不同雷诺 数对颗粒的运动的影啊。 2．数值模拟方法 2．1槽道中流体流动的大涡模拟 大涡模拟中，流体相的运动控制方程采用经过滤函数作用后的形式： 1第‘作者．博士研究牛 2联系作者，副教授，研究必趣为多相湍流反应流体力学及燃烧过程数值模拟 452 堕：o (1) 缸 鱼+百磐：一上皂+皇面(粤+§ (2) Ot j瓠。p ax|ax。 ax，axj 其中i’=ir+惫，i为分子动力粘性，_r为湍流粘性，由亚网格模型决定。 动态模式确定亚格子模式系数时，需要进行统计平均，在槽道流动中，当流动趋于 稳定以后，平{亍于壁酉的平面就可以看作是统计均匀方向，并且可以运用Meneveau等 提出的沿质点轨迹平均得动态确定模式系数的方法。根据动态模式确定的亚网格湍流粘 性为： i．：(e△)：√云忑其中(c／x)2= 踹 在壁面方向，由于流场速度变化较快，因而可以在壁面附近加密网格的数目，用以 捕捉到壁面流场的变化，同时在槽道中心使用比较稀疏的网格，这样就能控制总的网格 数，减少计算量。 2．2颗粒运动方程 流体拖曳力的影响(忽略重力的作用)，颗粒运动控制方程可以写为： 。