离心式通风机内稀疏气固两相流整机三维数值模拟.docVIP

离心式通风机内稀疏气固两相流整机三维数值模拟 田贵昌 李宝宽/东北大学 本文档由【中文word文档库】提供，转载分发敬请保留本信息； 中文word文档库免费提供海量范文、教育、学习、政策、报告和经济类word文档。 陈中才 李泰勋/ 沈阳鼓风机通风设备有限公司 摘要：利用颗粒轨道模型对通风机内部稀疏气固两相流动进行了数值模拟，重点分析了最优工况下不同粒径粒子轨迹分布情况。模拟得到了不同直径的固体颗粒在风机中的运动轨迹，积尘分布。这些结果为研究风机磨损，延长使用寿命提供理论基础。 关键词：离心式通风机；三维数值模拟；气－固两相流 中图分类号：TH432 文献标识码：B 文章编号：1006-8155（2007）05-0009-07 3-D Numerical Simulation of the Sparseness Gas-solid Two-phase Flow Field in Centrifugal Fan Abstract: In this paper, particle trajectory model is applied to numerically simulate the sparseness gas-solid two-phase flow field in fan. Especially, the tracks’ distribution of particles with different diameters is analyzed under the best operation. The result provides theory basis for studying the abrasion and prolonging service life of centrifugal fan. Key words: centrifugal fan; 3-D numerical simulation; gas-solid two-phase flow 0　引言 叶轮机械中的通风机、鼓风机等用途比较广泛，受工作环境的影响，这些叶轮机械在运行中经常有固体颗粒渗入。据我国风机行业协会统计，输送气固两相混合物的风机约占我国风机年产量的40%，所以对通风机气固两相流的研究有重要的意义。 据调查，“磨损比较严重的风机在半个月以内就要更换叶轮，一般运行3个月或者半年以后，就要进行补焊，才能维持大修周期，方能使风机在恶劣的环境下工作时能保证较高的效率和安全系数” [1]。国外对此研究较早，A.Hamed[2]在1982年采用三维气体流场计算了汽轮机定子中的粒子运动轨迹。M. Fathy和Tabakoff E[3]在1983年提出了反弹模型理论，但是并没有见到在现实中的应用。国内的费瑞乾[4]考察了边界层对粒子轨迹的影响。姜晓敏[9-10]。 1.2 颗粒轨道模型的基本方程 忽略颗粒相自身各变量的脉动，并且忽略流体相的密度脉动及变质量源脉动，则可得到确定轨道模型的湍流两相流基本方程。 流体相连续方程为 颗粒相连续方程为 流体相动量方程为 颗粒相动量方程为 流体相能量方程为 颗粒相能量方程为 流体组分方程为 其中；；；；； 式中为颗粒流速垂直于流向断面的分量；为有效粘性系数；为常数；为组分质量分数；为相变中组分的贡献分数；为流体相中组分反应率；为密度；为时间；为速度；为温度；下标为颗粒相的变量；带上标的变量为湍流脉动值。在动量方程中为考虑浮力影响的重力项。在能量方程中为单位体积中流体相由于变质量所造成的能源项；为流体辐射热；为颗粒辐射热；为流体相单位体积中反应放热；为颗粒表面热效应；为颗粒的数密度；为每个颗粒的质量；为Magnus力；为单位体积中平均的物质源；为驰豫时间。 2　计算对象及边界条件 2.1 风机模型参数 分析对象为4-73№10D离心通风机，如图1所示，由进气室、集流器、叶轮和蜗壳组成。在Pro/E中建立模型，为解决问题的方便，在整机的装配中让绝对坐标和相对坐标处于同一位置，原点位于叶轮后盘中心外壁上，X轴负方向为蜗壳出口方向；Y轴负方向为蜗壳的进气方向；Z轴正方向为进气室进口方向。叶片后倾，进、出口角分别为32°、45°，叶轮内径720mm,叶轮外经1000mm,叶片进口宽350mm,叶片出口宽250mm,进气室吸风口为1300mm×600mm,蜗壳宽650mm,出风口为900mm×650mm，叶片12个，转速1200r/min。 2.2　网格划分 在GAMBIT中对流道区域划分网格，如图2所示。由于风机结构较复杂，采取四面体和六面体网格相结合的方式划分，网格共计676045。叶轮流动区域采用旋