基于动网格的离心泵内部流场数值模拟_江帆

文章编号:　1005—0329(2007)07—0020—05

基于动网格的离心泵内部流场数值模拟

江　帆

1,2

,陈维平,王一军,区嘉洁

211

(1.广州大学,广东广州510006;2.华南理工大学,广东广州5100)

摘　要:　为了实时地模拟离心叶片泵的内部动态流场,将动网格技术引入到离心泵内部流场的数值模拟中,并分析了离心泵在叶轮旋转情况下的内部流场的动态变化及对水泵蜗壳结构改进计算。模拟结果显示,动网格技术能较好地模拟离心泵内部流场的动态变化,便于研究人员分析离心泵内部流场参数的瞬时变化情况。关键词:　离心泵;数值模拟;动网格

中图分类号:　TQ051.21　　　　文献标识码:　A

NumericalSmiulationofFlowFieldInsideofCentrifugalPumpBasedonDynamicsMesh

JIANGFan1,2,CHENWei-ping2,WANGYi-jun1,OUJia-jie1

(1.GuangzhouUniversity,Guangzhou510006,China;

2.SouthChinaUniversityofScienceandTechnology,GuangZhou5100,China)

Abstract:　Inordertoreallysimulatetheflowfieldinsideofcentrifugalpump,themovinggridtechnologyisintroducedtonu-mericalsimulationofflowfieldinsideofcentrifugalpump,andthevariationofinteriorflowfieldundertheconditionoftheimpel-lerrotatingisanalyzed.Theresultsshowthatdynamicsvariationoftheinteriorflowfieldiseasytobeobtainedbydynamicsmeshtechnology,anditprovidesawaytoanalyzethetransientstatusofparametersoftheinteriorflowfieldforresearchers.Keywords:　centrifugalpump;numericalsimulation;dynamicsmesh

1　前言

离心泵是常见的输水设备,在各个领域应用比较广泛。离心泵内部流场的数值模拟为泵叶片及泵腔的结构设计和优化提供详尽的参考数据,并代替部分实验,在离心泵设计中占有重要作用。

水泵内部流场模拟较多,对于叶片转动的动态模拟一般采用滑移网格进行瞬态模拟,这种模拟对叶片转动的处理是叶片单独划分一个区域,叶片的转动与周围的区域的数据交换是通过滑移面来传递的

[1]

(包括凸轮泵和齿轮泵)内流动。本文将动网格技术引入到离心泵内部流动的数值模拟中,动态分析离心泵内部流场随叶片转动时的流动变化,可以部分代替样机实验,帮助设计人员校核水泵的运行参数能否达到设计要求,以及为离心泵的结构优化及新型离心泵的设计提供参考。2　动网格技术

离心泵内叶片是整周转动,位移比较大,采用网格变形与局部重构相结合的方法实现动网格。基本做法为:当叶片转动较小时,叶片周围的网格的每个边看成一个个弹簧,随着叶片转动作微小变形;当变形较大时,叶片周围的网格变形超过一

[2]

。滑移网格处理泵内叶片的瞬态转

动具有一定好处,但不太精确,特别是滑移面附近

会出现数据不连续,同时,传统的模拟方法需要给定详细的水泵进出口参数。针对这个问题,国外

一些学者开始采用动网格技术进行泵内叶片转动定的限度,开始进行局部重划分,使网格达到要

[1,3]

的数值模拟,比较接近真实情况,如KrisRi-求。具体算法如下:emslagh等人采用动网格技术计算了旋转置换泵

收稿日期:　2006—08—10修稿日期:　2007—04—10

当小变形时,采用基于弹性变形的网格更新

2007年第35卷第7期　　　　　　　　　　流　体　机　械　　　　　　　　　模型,该模型将网格单元的边看作弹簧,平衡长度等于网格边的原长,故网格点的初始状态所受的合力为0。当边界网格点移动后,由虎克定律,网格节点i位移引起的力Fi为:

Fi=΢Kij(δj-δi)j=1

式中　Ni———与节点相邻节点的数量

Kij———连接节点i和j弹簧的弹性系数,

(N/m),Kij=1/lij

lij———节点i,j的边长

δ——节点i,j的位移距离,mi,δj—

根据力平衡原理,对于每一个节点,连接到其上的所有弹簧的合力必为0,则得到该节点的位移:

δi

k+1

Ni

Ni

21

du p u2 u v

ρ=ρfx++{μ[2-(+dt x x x3 x y

w u v

　　+)]}+[μ(+)]

z y y xw+ u)]　　+ [μ( z x z

dv p v2 u vρ=ρfy++{μ[2-(+dt y y y3 x y w v w　　+)]}+[μ(+)]

z z z y u v

　　+[μ(+)]

x y xdw p w2 u vρ=ρfz++{μ[2-(+dt z z z3 x y w w u　　+)]}+[μ(+)]

z x x z v w

　　+[μ(+)]

y z y

(5)

3.2　湍流模型

[5,6]

(1)

΢Kijδj

j=1=Ni

΢Kijj=1=xi

old

k

(2)

由式(2)可知,网格节点新的位置为:

xi

new

+δi(3)

当大变形时,采用局部网格重构的网格更新模型。对于给定的理想网格高度hideal,当h>(1+α时,网格将根据预定义的h)hideal(αh为高度系数)高度条件进行,这时,在层i中的网格面的高度将正好是理想高度hideal。相反,如果层j中的网格体积是被压缩的,当压缩到h<αhhideal时,这些被压缩的网格面将与邻近层的网格面合并成一个新的网格层。

3　控制模型及数值解法3.1　流动控制模型

(1)质量守恒模型

[4]

采用标准k-ε模型,形式如下:

μtdk kρ=[(μ+)]+Gk+Gb

dt xiσxik

-ρε-YM

μdε t εερ=[(μ+)]+C1ε(Gk

dt xiσxikε

ε

+C3εGb)-C2ερk

式中　Gk———平均速度梯度引起的湍动能,J

Gb———浮力影响引起的湍动能,J

YM——————可压缩湍流脉动膨胀对总的耗

散率响,J

kμ——湍流粘性系数,μρCμt—t=

ε

σ——湍动能k的湍流普朗特数,σ1.0k—k=σ2———耗散率ε湍流普朗特数,σ1.32=C1ε,C2ε,C3ε———常数,取C1ε=1.44,C2ε=

1.92,C3ε=0.093.3　基于动网格的流场控制模型的变化流场控制模型是基于流域固定形状的情况,

当采用动网格技术后,计算区域是变化的,所以要对上述流场控制模型进行改造,即考虑动边界移动的影响,表达式如下:

d　　Vsρ dV+ρ (u-ug)ndS Lsdt　　　=Γ ndS+q dV LsVs

(8)

[1,7]

2

2

(6)

(7)

ρ (ρu) (ρv) (ρw)

　+u+v+w

t x y z　　=0

式中　ρ———流体密度,kg/m

m/s

t———时间,s

x,y,z———三个坐标方向(2)动量守恒模型

[4]

3

(4)

u,v,w———流体在三个方向上的速度分量,

流体在三个方向上所受的体积力:fx,fy,fz

(N),压力:p(Pa),粘性系数μ(N s/m)。2

22　　　　　　　　　　　　　　FLUIDMACHINERY　　　　　　　　　　　Vol.35,No.7,2007式中　 ———通用变量

Vs———控制体积,m

Ls———控制体积的边界u———流体时均速度,m/sug———动网格边界移动速度,m/sn———表面Ls上方向朝外的法向单位向量Γ———扩散系数q ———源项3.4　边界条件及物性参数

进出水口设置为压力边界条件,叶片轮廓为

动边界,转速n为1450r/min,扬程H为50m,流量Q为400m/h。工作介质为水,密度ρ为998.2kg/m,粘性系数μ为0.001003N s/m。3.5　数值解法

在FLUENT平台进行数值计算,采用有限体积法求解,压力项用PRESTO格式离散,扩散项用中心差分格式离散,其余项用二阶迎风格式离散,压力速度耦合方程采用PISO算法求解。4　离心泵内部流场的动态模拟

模拟采用的离心泵结构是根据IS200-150-400水泵稍作改动后的结构,蜗壳基圆直径D为410mm,叶片厚度t为7mm,共6片,进口直径D1为150mm,出口直径DO为200mm。4.1　网格模型

离心泵内部流场计算区域采用三角形网格离散,其网格模型如图1所示。

图2　叶片转动0.0045s后的网格

3

2

3

3

动其转动。初始状况为3740个节点,个网格。在计算中,因为动边界的转动(即叶片转动),计算区域的网格会发生变化,如图2所示。

4.2　计算结果及讨论

离心泵内部不同瞬时的压力分布如图3和图4所示。不同的瞬时,压力分布有所变化,但是压力分布的大致规律相似,即压力从进口到出口递增。

从图3和图4中还可看到,叶片的压力面压力大于吸力面的压力,压力面压力分布比较均匀,而吸力面压力变化则比较大,并且在进口附近出现较大的负压,这是在叶轮进口处容易引起空化的主要原因。

图1　离心泵计算区域初始网格

图3　0.052s时泵内压力场分布

叶片部分定义为动边界,由UDF聚光函数驱

离心泵内部不同瞬时的速度矢量分布如图5、6所示,叶片的转动对水流速度影响较大,但内部2007年第35卷第7期　　　　　　　　　　流　体　机　械　　　　　　　　　流场分布规律接近一致,在出口会产生回流。离心泵内部不同瞬时的流线分布如图7、8所示,从图中看到在出口左侧有明显的回流。对比现有的泵内流场的计算结果,采用动网格技术计算与其它方法计算的区别就是叶片附近流速较大,主要是因为叶片动边界高速转动引起的,这个与实验

[8]

情况相符。

23

图6　0.062s时泵内速度场分布

图4　0.062s时泵内压力场分布

图7　0.052s时泵内流线

图5　0.052s时泵内速度场分布

基于动网格的内部流场模拟能较好地描述泵

内的瞬态现象,如图9所示,叶片进口处的压力随时间和叶片位置的变化,可以间接说明叶片进口空化的产生和扩散,为改进结构或运行参数以控制气蚀提供参考。

图8　0.062s时泵内流线

24　　　　　　　　　　　　　　FLUIDMACHINERY　　　　　　　　　　　Vol.35,No.7,2007

(a)　0.052s　　　　　　　　　　　　　(b)　0.053s　　　　　　　　　　　　　(c)　0.0s

图9　叶片进口处压力云图随时间和叶片位置的变化

rianfinite-volumemethodforthesimulationofrotary

5　结语

动网格技术对运动边界的描述与人们对运动边界认识是一致的,不象滑移网格技术通过交界面传递数据,给计算区域模型处理上带来不便,动网格直接建立计算区域模型,通过UDF或profile文件驱动动边界运动,给离心泵的流场模拟带来极大便利。

动网格技术针对计算区域的几何形状随时间变化的流体动力学进行瞬时动态模拟,可以捕捉离心泵内流场瞬时变化,特别是在分析泵内叶片间漩涡的生成、运动、合并、、脱落等水流特性,以及和多相流计算相结合,模拟泵内的空化产生、扩散等现象,为离心泵的设计和结构优化、及改进运行参数提供参考。

参考文献

displacementpumpflow[J].Mathematics,2000.32:419-433.

AppliedNumerical

[3]　MoureauV,LartigueG,SommererY,eta.lNumeri-calmethodsforunsteadycompressiblemulti-componentreactingflowsonfixedandmovinggrids[J].JournalofComputationalPhysics,2005,202:710-736.[4]　PotterMC,WiggertDC.Mechanicsoffluids[M].

ChinaMachinePress,2003.

[5]　江帆,陈维平,李元元,等.基于射流与两相流的射

流曝气器研究[J].流体机械,2005,33(6):18-21.[6]　JiangF,ChenWP,MaiMR,eta.lNumericalsimu-lationonwaterandsedimentstwo-phaseflowsinriver[C].

ProgressinSafetyScienceandTechnology,

2005,11,1634-1639.

[7]　杨建明,吴建华.动网格技术数值模拟挑流冲刷过

程[J].水动力学研究与进展,2001,16(2):156-161.

[8]　方毅,柯坚.离心泵流道和泵腔内流场的试验研究

[J].机械工程学报,2005,41(10):45-49.

[1]　江帆.悬浮载体材质-结构的生物特性及新型转笼

生物反应器[D].华南理工大学博士学位论文,2006.

[2]　KrisR,JanV,ErikD.AnarbitraryLagrangian-Eule-作者简介:江帆(1974-),男,博士,研究方向为机械CAE/CFD,汽车安全与环境,通讯地址:510006广东广州市广州大学信息与机电学院机电系。

(上接第9页)

[3]　王连泽,彦启森.旋风分离器减阻杆结构及减阻前

参考文献

后流场的测定与分析[J].实验力学,1998,13(4):469-476.

[1]　旋风分离器联合攻关组.PV型旋风分离的性能及

工业应用[J].石油炼制,1990,(1):37-42.

[2]　王连泽,彦启森,叶龙.旋风分离器压力损失及减阻

杆的研究[J].力学与实践,1998,20(4):32-36.

作者简介:张从智(1962-)男,副教授,博士研究生,通讯地址:257061山东东营市北二路271号中国石油大学(华东)多相流实验室。