离心泵全三维流场的大涡数值模拟

离心泵全三维流场的大涡数值模拟第34卷第4期2019年4月华南理工大学学报 (自然科学版)JournalofSouthChinaUniversityofTechnology(NaturalScienceEdition)Vol.34No.4April2019文章编号:10002565X(2019)0420191204离心泵全三维流场的大涡数值模拟黄 思 12(1. 华南理工大学工业装备与控制工程学院 , 广东.)摘 要:, 计算单级蜗壳式离. 发现泵叶轮内各通道的流量、流速及 , 流动呈现明显的非对称性 , 泵内流动旋涡一般出现在叶轮叶片工.文中还将泵性能的预测值与实测值作了对比,验证了计算结果的有效性.关键词:大涡模拟;离心泵;三维计算;FLUENT软件中图分类号:TH311.022文献标识码:A离心泵在农业及民用等许多部门中有着广泛的应用. 实现离心泵的流动数值模拟 , 对于泵的优化设计以达到增效节能的目的具有重要的现实意义.离心泵系统一般导入管、叶轮及泵壳等元件.随着近年流体力学计算技术的迅速发展,离心泵的三维数值ε等模拟已成为现实 . 但离心泵湍流模拟多见于 κ2[128]双方程湍流模型的使用对流动中的较大尺度旋涡进行

.

大涡模拟技术(LES) 是近年来才发展起来的湍流模型Navier2Stokes(N2S) 方程直接计算(DNS),

, 它对较小尺度旋涡进行平均近似处理

(RANS).

研究表明大涡模拟要优于通常的

κ2ε湍ε模型那样计算出流动的流模型

, 它不仅能够像

κ2时均规律, 而且能够较准确地描述流动中的不对称特性 . 受计算机硬件的限制 , 以往大涡模拟常见于简单几何体内的三维湍流计算中,近年来一些研究者尝试使用大涡模拟以较稀疏的计算网格应[9210]用于工程计算中. 为此本文选取一低比转速 (ns=45) 离心泵作为研究对象 , 以流动软件FLUENT为工具,使用大涡模拟湍流模型对包括导入管、叶轮及泵壳在内的全流场进行设计工况下的三维计算.行网格过滤, 从而得到较大尺度旋涡的基本方程组 .连续性方程(ρ(1)+u=0i)9t9xi大涡的动量方程组 (ρ)(ρ)u+u=iiuj9t9xj-++,i,j=1,2,3(2)9xi9x9xj9

目前使用最广泛的模型方程为-

τμtSij+

τ(3)ij=-2kk

δij3其中Siji(4)Sij=+

定义为29xj9

对μt

使用

Samagorin2Lilly

模型,

此模型方程为2μt=

ρ(5)LsS

式中S=SijSij1/3--(6))Ls 使用的计算公式为Ls=min( κd,CsV理论计算方法1.1 基本方程组大涡模拟的思想是在流动区域内对 N2S方程进收稿日期:2019204201作者简介:黄思(19622),男,博士,副教授,主要从事流体机械流动理论与设计研究.E2mail:huangsi@(7)(FLUENT 软件中取值 Cs=0.1, κ=0.4187).1.2 计算区域、网格的生成计算区域由离心泵导入管、叶轮及泵壳组成 . 使用FLUENT软件中的前处理程序 Gambit生成计算区域几何体,再进行网格划分,得到如图1所示的四面体网格单元.为降低计算成本,实现大涡模拟计算的实用化、工程化,本文具体计算时采用了与文献华南理工大学学报(自然科学版)第34卷[9210] 同等级网格尺寸的计算网格 . 四面体网格单元的数目如下:入水管29098单元; 叶轮76853单元; 泵壳95017单元.为此采取减小与速度压力相关的松弛因子 , 当迭代计算到一定阶段时将松弛因子设为0.1 或更低. 图2给出残差随松弛因子减小而收敛的过程 . 使用上述-4方法,可使所有参数的残差控制在10以内.Fig.ononinthepump1.3 针对含有运动固体边界的流动问题 ,FLUENT软件中提供了多重参考坐标系.具体计算时将叶轮区域设在运动坐标系,导入管和泵壳区域设在固定坐标系.计算使用了如下的边界条件:进口边界条件.按入口质量流量值设定,具体数值由泵设计工况给出.入口湍流取值按水力直径大小及湍流强度 IUa图2计算残差随松弛因子变化而收敛的过程Fig.2Convergencyprocessofcomputationalresidualwiththevariationofrelaxationfactor计算结果及其分析图3至图7给出了流动趋近稳态后 (t=3600s) 的部分计算结果. 图3给出泵全三维的绝对流速矢量图 . 由图可见, 绝对流速大小在 40m/s 以内, 较高的流速出现在泵壳区域内. 为了了解泵内出现旋涡、回流的情况 , 除了将绝对流速矢量图作局部放大观察外 ,也可以观察如图 4所示的流体质点迹线图 . 由图4可见, 较大尺度旋涡一般出现在叶轮叶片工作面上.图5给出的是叶轮中心回转面上的相对流速矢量图 . 值得注意的是叶轮 6个通道内部流场明显地表现出不一致 (即非对称性), 靠近泵出口通道的流速明显高于其它通道的流速 , 其中最靠近泵出口的通道流速最高.k 给定(I 取5%);3出口边界条件.取流动充分发展条件,即假设所有变量的扩散通量为0;3)壁面边界条件.FLUENT软件在应用于大涡模拟的具体计算中,为使近壁面网格间距不受严格的限制,使用了如下的壁面近似方法,即假定与壁面相邻的网格单元的质心处于边界层的对流区域,其方程可表达为ρu(8)=ln τκuμτρ式中,E 为常数9.793,u. τ=w/1.4 流体物性及泵操作工况参数设定计算中所使用的泵工况、流体物性等参数见表 1.表1 泵工况参数及流体物性参数Table1 Parametersofpumpoperatingpointandfluidproperties泵流量-1泵转速-1泵比转速s泵入口操作压力p/Pa流体密度流体粘度 μρ/(kg ?m-3)/(kg ?(m?s)-1)998.20.001003/(kg ?s)/(r ?min)551450451013251.5 计算收敛性及控制方法求解旋转坐标系流动问题遇到的主要困难是动[9]量方程的高度耦合,另外区域内一般还存在回流,这种耦合及回流现象会导致求解过程的不稳定.图3 离心泵内流场的绝对流速矢量图Fig.3 Vectorsofabsoluteflowvelocityinsidethecentrifugalpump第4期黄 思等:离心泵全三维流场的大涡数值模拟 113增加. 在泵壳的一段范围内 , 液流的动能转换成势能使得泵级内的静压值达到最大 . 后因沿程出现的水力损失使静压值有所降低.Fig.4ofthepu图7 离心泵中心面上的静压分布Fig.7 Staticpressuredistributiononcentralsurfaceofthepump图5 叶轮中心面上的相对流速矢量图Fig.5 Relativeflowvelocityoncentralsurfaceoftheimpeller鉴于此发现, 将叶轮各通道所通过的流量值计算出并示于图 6中. 由图6可见, 最高流量值(通道1) 与最低流量值(通道3) 相差几倍. 上述现象的原因是因为通道 2～5与泵壳壁面相距较近且离泵出口较远 , 叶轮出流被泵壳壁面堵塞使过流量减少 . 相反, 与泵壳壁面相距较远或离泵出口较近的通道 , 流量及流速就比较高.由模拟流场的结果可算出泵出口位置 b到入口位置 a的总压差即为该水泵的扬程H.22H=dA-dA++AbAbpg2AaAa

ρg2(9)为验证本文计算结果的有效性 , 需要将计算结果与实测结果进行对比尚缺乏同型号水泵的内流场实测数据 , 因此本文只进行离心泵的外特性对比分析算出该泵的扬程 H=35.8m, 接近产品样本提供的数据 (H=38m). 由此可见大涡模拟湍流模型模拟离心泵的全三维湍流场是现实可行的 .

.,

由于目前. 由式(9)本文使用3

结论通过使用FLUENT软件的大涡模拟湍流模型模拟计算离心泵的全三维流场 , 计算结果表明泵叶轮各通道的流量、流速及压力分布等表现出明显的非对称性 , 其中最高值与最低值相差数倍. 泵内流动大尺度旋涡一般出现在叶轮叶片工作面上 . 本文还将泵性能的预测值与实测值作了对比 , 验证计算结果是正确的 . 符号说明:图6 叶轮各个通道内的流量分布Fig.6 Flowratedistributionineachchanneloftheimpeller—过流断面的面积,m;Aa—截面a的面积,m;Ab—截面b的面积,m;g—重力加速度,m/s;I—湍流强度,%;Ls—网格的混合长度,m;2222Cs—Samagorin常数;H—水泵的扬程,m;k—湍流动能,m/s;ns—离心泵比转速;22图7给出叶轮与泵壳中心面上的静压分布.由图可见,叶轮各通道内静压场也表现出非对称性;由于叶轮旋转做功,叶轮内的静压值随流动方向逐渐—网络到壁面的最近距离,m;E—计算常数;114p—静压,Pa;t—时间,s;Ua —截面a的平均流速,m/su 粘性摩擦速度,m/s; τ——计算单元的体积,m;y—网格点到壁面的距离,m;3华南理工大学学报 (自然科学版)Sij(i,j=1,2,3) —旋转张量第34卷马福喜,李志伟.大涡模拟水环境中污染物团的运动分量,m/s2;ui,uj(i,j=1,2,3)

—流速分

量,m/s;xi,xj(i,j=1,2,3)

—空间几规律[J]. 水利学报,2002(9):55260.MaFu2xi,LiZhi2wei.Largeeddysimulationofpollutantmovementinwater[J].JournalofHydraulicEngineer2ing,2002(9):55260.MajidK.Numericalstudyofunsteadyflowinacentrifugalpump[J].JournalofTurbomachinery,2019,127(2):3632.[7GtoA.forpumpsbasedon3CFD,[J].JournalonsoftheASME,2002,.郭鹏程,罗兴奇,刘胜柱.离心泵内叶轮与蜗壳间耦合δ单位张量ij(i,j=1,2,3) —何坐标,m; ε—湍能耗散率,m2/s3;κ—VonKarman 常数; 分量;μ—流体动力粘度,kg/(m ?s); μt —亚网格湍流粘度 , ρ—流体密度,kg/m3;τ应力张量ij(i,j=1,2,3) — 分量, kg/(m ?s);τ,w—唐辉,.[J].水泵技流动的三维紊流数值模拟[J].农业工程学报,2019,21(8):125.GuoPeng2cheng,LuoXing2qi,LiuSheng2zhu.Numericalsimulationof3Dturbulentflowfieldsthroughacentrifugalpumpincludingimpellerandvolutecasing[J].Transac2tionsoftheChineseSocietyofAgricultureEngineering,2019,21(8):125.杨建明,刘文俊,吴玉林.用大涡模拟方法计算尾水管术,2002(3):328.TangHui,HeFeng.Numericalsimulationofflowsinacentrifugalpump[J].PumpTechnology,2002(3):328.[2] 王彦辉, 康志忠, 徐鸿, 等.32Sh219离心泵流场的全三维数值模拟[C]∥2002年FLUENT中国用户年会论文集 . 上海:[s.n.],54260.YanJ,SmithDG.CFDsimulationof32Dimensionalflowinturbomachineryapplications[R] ∥[s.l.]:Turboma2chineryFlowPrediction Ⅷ.ERCOFTACWorkshop,2000:342 38.杨建明,刘文俊,吴玉林.用大涡模拟方法计算尾水管内非定周期性湍流 [J]. 水利学报,2001(8):79284.YangJian2ming,Liuwen2jun,WuYu2lin.Thelargeeddysimulationmothodforcomputationof32Dunsteadyperi2odicturbulentflowindrafttube[J].JournalofHydraulicEngineering,2001(8):79284.马福喜,李志伟.大涡模拟水环境中污染物团的运动内非定常周期性湍流[J].水利学报,2001(8):79284.YangJian2ming,LiuWen2jun,WuYu2lin.Thelargeeddysimulationmethodforcomputationof32Dunsteadyperi2odicturbulentflowindrafttube[J].JournalofHydraulicEngineering,2001(8):79284.规律[J]. 水利学报,2002(9):55260.MaFu2xi,LiZhi2wei.Largeeddysimulationofpollutantmovementinwater[J].JournalofHydraulicEngineer2ing,2002(9):55260.Large2EddyNumericalSimulationofThree2DimensionFlowinCentrifugalPumpHuangSi WuYu2lin12(1.CollegeofIndustrialEquipmentandControlEngineering,SouthChinaUniv.ofTech.,Guangzhou510640,Guangdong,China;Dept.ofThermalEngineering,TsinghuaUniv.,Beijing100084,China)Abstract:Athree2dimensionturbulent