气固两相圆柱绕流直接数值模拟.pdfVIP

中国工程热物理学会 多相流 2002年学术会议 026040 气固两相圆柱绕流的直接数值模拟 樊建人胡桂林，姚军，岑可法 浙江大学热能工程研究所，杭州，310027 摘要：本文运用谱元方法对气固两相圆柱绕流进行了直接数值模拟．在获得高精度计算流场信息基础 上，进行颗粒扩散运动的研究．通过研究颗粒大小，流动雷诺数两个关键因素对颗粒扩散运动的影响 来刻画气固两相圆柱绕流的物理特征。研究时发现，由于圆柱尾流独特的涡结构特征，颖粒在流场中 心轴线附近区域的扩散机理主要是决定于吸力作用，该作用会驱使小颗粒卷吸进入圆柱尾部近壁区， 甚至与圆柱后壁发生碰撞，这种扩散机理与混合层中颗粒扩散机理明显不同． 关键词；潜元方法，直接数值模拟，气固两相，圆柱绕流，扩散． 1)( 1．引言 颗粒在钝体尾迹中的扩散分布是自然界中～种普遍现象，特别是颗粒在圆柱尾迹中 的运动分布更为普遍。它们有着广泛的工业应用和技术背景，譬如，煤粉气流喷入带圆 形钝体稳焰装置的燃烧室中，此时煤粉在圆形钝体尾流中的运动就属于典型的固体颗粒 在圆柱尾迹中的运动。又如大型循环流化床锅炉的撞击圆管式惯性分离器；含灰的气流 冲刷锅炉对流受热面的管束等。因此，对颗粒在圆柱尾流中运动扩散的物理机理进行研 究具有重要的工程价值和实际意义。 曾经有研究者”…对颗粒圆柱绕流运动问题用不同的数值方法进行计算。但因为计算 所用方法精度不高，计算结果与实验数据相差较大．此外，他们并没有对(不同大小) 颗粒在圆柱绕流拟序结构中扩散分布物理特征进行研究。近些年来，应用高精度计算方 法求解圆柱绕流流场的报道很多，但研究气固两相圆柱绕流流动直接数值模拟的论文未 曾发现。本文研究目的是应用谱元方法对圆柱绕流问题进行直接数值模拟．获得高精度 的收敛结果。在此基础上对颗粒在圆柱绕流拟序结构中扩散分布物理特征进行研究。 鲫 图1计算区域网格单元 326 2． 数学方法及求解 本文采用如图1所示的计算网格。计算区域为40x28．圆柱半径R为1，整个计算区 域分为50个单元，每个单元采用8阶插值函数。计算采用时间步长为0．005。 2．1气相控制方程 原始变量形式的无量纲化不可压N．S方程可写为： 詈一即+去￡(叭砌) 在。内 (1) v·V=0 在Q内 (2) 其初始条件和边界条件视具体问题而定。其中v为速度．P为压力(密度P一1)。雷 诺数Re=Ud／v， U为流动的特征速度，d为流动的特征长度， v为运动学粘度。线 性扩散算子和非线性对流算子分别为： L(V)=V‘V N(V)=一÷Iy·Vy+V·(yy)l 对流项表示成斜对称形式以减小走样误差。 2．2求解方法 N．S方程采用混合显式隐式算子分裂格式进行离散I”，方程的求解分三步完成 尹一笙qP。一 q=O = ：=f △f 岛 _Ⅳ -y ¨∑唧 !二 竺兰；一W一 (4) ‘ △f ： 业：：二受Re--v：眇 (5) At ．： o为适当选择的加权系 其中矿，旷为中间速度变