用二阶矩颗粒壁面相互作用新模型对突扩两相湍流地研究.pdfVIP

中国工程熟物理学会 多丰H流 第卜届年会 9 编l，：01601 用二阶矩颗粒壁面相互作用新模型 对突扩两相湍流的研究1 张夏周力行 ‘清毕大学丁=程力学系．僦酌清洁燃烧国家重点实验室．100084) 摘要车史吸收了硬流体模型用PDF积分和颗粒拉氏模型全面考虑伺f撺过程的优点．攫…一阶币顺粒 壁面相互作用新模型，考虑壁面粗糙度对两相湍流的影响。对轴对称突靠流动进行r数值模拟，着荤 探讨r不同颗粒相边界条件的影响。结果表明，壁画处颗粒相边界条件的选取足卉恰兰j，小便蟛啊壁 面Ⅲ近，而且会影响到全场的两相湍流．颗粒越^，受到的这种影响也越人。由十考虑，再方向雷诺 应力之间的相互转化，雷诺应力从平均运动中得到能量，以及壁面对颗粒运动的裒减作用等…翥．也： 文提出的包括摩擦系数、恢复系数、壁面粗糙度等影响的新模型得到的结果≮实验符台得史好．111J以 前使用的梯度为零的颗粒壁耐扛界条件则给出完全失真的模拟结果。 关键词f颗粒与壁面作用y突扩流动，两相湍流 1．引言 现有的气闽两相流动的各种二阶矩两相湍流模型，都能昭合理地预报山两相时平 均速度场，但不能较好地预报两相湍流脉动，两相雷诺应力再分擐的模拟结果与实验 值之间尚存在较明显的筹别『11。一个重要原因是忽略了颗粒肇面碰撞机理，目前往拄 麻瑚十分简化的颗粒相边界条件，例如，设罐面处颗粒法向速度为零，浓度、切向速 度和雷诺应力的法向梯度为零等。实验中已经发现，颗粒碰鼙后，特别是在弯道中， 颗粒湍流脉动显著增大。本文作者在吸收颗粒拉氏碰壁模型考虑鼙面粗糙度和烈流体 模型j}j概率密度积分法处理颗粒碰壁模型优点的基础上，提出了二阶矩两相湍流模刑 中颗粒壁面相互作用新模璎．引入了壤面粗糙度对碰擘颗粒湍流脉动影响的机理[2]。 本文以轴对称突扩两相湍流流动为例，进行了数值模拟，探讨包括新模刑往内的各种 颗粒相边界条件对两相湍流特性模拟结果的影响。 2．模拟对象、基本方程和数值方法 本文模拟的突扩通道几何形状及尺寸如图1所示。空气进口速度为27．5rigs，颗粒 1．1 um和95m。轴对称柱坐标系中二阶矩两相湍流模刑 质量载荷比为O．005，粒径为50 气粒两相流动方程组可表示成如下形式： (1) 导(pv中)+一1i8(pr矿中)：昙(r。罢)+一Ii0(r。r罢)+氏 (2) 7一x(nYpO，)+7l磊J(～r_。，)=嘉(～套+÷暑(r卸r争+％ 式(1)、(2)分别代表气相和颗粒相，由、L、矗、由，、r却_手lI％的具体内容见文献 [2l。采俐非均匀交错网格，廿点数取24X46。用SIMPLEC算法，收敛标准为最人相 PC机上运行一个I：况需10分钟。 对余源和小T-104。在PENTIUM．11-266 本文为国家重点基础研究专项经费资助项目(G1999-0222-08)研妃成果 139 3．二阶矩颗粒墅面相互作用新模型 考虑拉氏模型颗粒壁面碰撞受摩擦系数、恢复系数和壁面粗糙度等的影响，运用双 流体模型中PDF积分方法，得到如下颗粒浓度、速度和雷诺应力的壁面边界条件： ％=撕+驰击等) ④ ％=帆，一巧．力fl一扣2I (4) 匕=÷∥(咋。+o，，)