基于fluen的天神离心叶轮湍流流动的数值研究.docxVIP

基于fluen的天神离心叶轮湍流流动的数值研究 1 离心叶轮流场模拟 离心叶片通道的流动是工程流量机械中最复杂的流动现象之一。该流量现象具有很强的三维特性。叶轮内湍流流动的研究为改进叶轮设计方法,提高叶轮效率提供了重要的依据,因此,近半个世纪以来,许多研究者对其进行了试验研究和数值模拟,1980年Johnson详细测量了叶轮从进口到出口5个截面的三维速度分布,并公开发表了该叶轮的几何型线,该叶轮被称为Ghost叶轮。在进行试验研究的同时,国内学者相继利用数值模拟方法对离心叶轮内部旋转流场进行了数值研究,分析了叶轮内的主流,二次流及出口“射流-尾迹”结构[1,2,3,4,5,6,7,8,9,10,11]。这些计算均采用k-ε模型,然而叶片通道是弯曲和旋转的,具有很强的各项异性,因此采用各项同性的湍流模型模拟的结果与试验相比差异较大。本文采用雷诺应力模型(RSM)对离心叶轮内的湍流流场进行了数值模拟,其中定子转子耦合模型采用混合平面模型(MP),得到了较好的计算结果,与试验结果相比基本一致,并分析了离心叶轮内部的粘性二次流和“射流-尾迹”现象以及尾迹区的位置。 2 数学的物理模式 2.1 叶轮入口长丝环境参数 文献提供的Ghost叶轮在圆柱坐标下的几何型线数据。Ghost叶轮的主要设计参数如下,叶轮宽径比b2/D2=0.08,出口安装角β2A=90°,叶片数Z=19,外径D2=909.3mm,轮盖处内径d=177.8mm,叶轮入口截面直径D1=571.5mm,设计转速n=500r/min。为了与试验叶轮接近,在叶轮入口前加装了导流锥;同时为了使计算易于收敛在叶轮的出口加了一个径向无叶扩压通道。生成的叶轮单叶片几何模型如图1所示。 2.2 结构网格划分 为了解决离心叶轮通道的计算域复杂而引起网格质量差的问题,本文采用了多块网格法,即将叶片流道求解区域划分成Ⅰ、Ⅱ、Ⅲ3个较为简单的块,分别划分结构网格如图2所示。 Ⅰ为叶片入口前的一个小三角区域,Ⅱ为贴近压力面的流道,Ⅲ为贴近吸力面的流道。沿叶片布置70个节点,周向布置44个节点,径向布置40个节点。叶片表面、轮盖和轮盘内表面采用边界层加密网格,出口扩压部分比较规则,可直接生成结构网格。 2.3 边界条件的确定 利用FLUENT软件,采用控制容积法离散控制方程,对流项采用QUICK格式,应用SIMPLEC算法对离散方程进行求解,其收敛标准为各个变量的相对残差均小于10-5。湍流模型采用RSM模型。转子和定子的耦合采用MP模型。 进口边界采用速度入口,速度大小按流量条件给定,流量m=3.3kg/s,该流量高于Ghost叶轮的设计流量2.5 kg/s。进口截面的湍动能设为进口截面平均动能的(0.5%～1.5%),耗散率ε按公式ηt=cμρk2/ε确定。出口边界条件采用局部单项化来处理。整个计算区域的两侧均设为周期性边界条件。固体壁面采用无滑移边界条件。 3 分析与讨论的结果 3.1 核心内流道的转子静压分布 文献给出了Ghost叶轮5个测量截面上的相对速度和主流速度分布,截面位置如图3所示。 图4～8给出了本文的计算结果,同文献的结果相比其结果是一致的。等值线表示垂直截面的主流速度,矢量表示截面内的二次流。 截面Ⅰ位于叶片流道的入口处,二次流速度存在由轮盘指向轮盖的分量,导致流道上半部分的主流速度较高。当流量大于设计值时,在压力面产生负冲角,在吸力面产生正冲角,尤其是在叶道上半部分冲角更大。由于冲击作用,轮盖/压力面角域形成了一个低转子静压区,在该区域内无量纲转子静压pr0.05,如图4(b)。而截面Ⅰ上转子总压分布比较均匀,因此在截面Ⅰ轮盖/压力面角域形成了一个主流高速区,如图4(a)。 截面Ⅱ位于诱导弯道的出口。旋转做功压力面的转子静压升高,如图5(b),由图4(b)和图5(b)可知,截面Ⅰ和截面Ⅱ之间流道的轮盖/压力面角域存在逆向的转子静压梯度,该角域处的流体减速。轮盖/吸力面角域则为正向转子静压梯度,该角域处流体被加速,因此主流高速区由截面Ⅰ的轮盖/压力面角域转向截面Ⅱ的轮盖/吸力面角域。受诱导弯道和轴向转向径向弯道的双重影响,截面Ⅱ轮盘附近产生一个与旋转方向相反的通道涡,如图5(a),该通道涡将轮盘边界层内的低能流体推向吸力面中部。 截面Ⅲ位于叶片流道的中间位置。轴向转向径向的弯道对该截面上的二次流的形态产生了显著的影响。弯道中出现了典型的二次流涡,如图6(a),将压力面与吸力面边界层内的低能流体推向轮盖面,同时受弯道旋转的影响,在轮盖面的边界层产生了由压力面指向吸力面的二次流,将轮盖边界层内的低能流体推向吸力面。由于低能流体的堆积,使轮盖面上的转子总压降低,并且轮盖/吸力面角域转子总压最低,如图6(b),因此在该角域处出现尾迹区,如图6(a)。 截面Ⅳ位于轴向转向径向的弯道的