梁月娟
陽江市新力工業有限公司,廣東 陽江 529500
利用三維建模軟件Pro-E建立模擬所用的葉輪模型和泵外殼模型,如圖1所示。葉輪的葉片數為7,設計為后彎葉片式,導葉的葉片數為10,設計為徑向式結構。利用FLUENT的前處理軟件Gambit對模型進行網格的劃分,劃分后的模型網格如圖2所示。

圖1 葉輪和泵體物理模型

圖2 模型網格示意圖
選取FLUENT軟件中的3D模型求解。在MODELSOLVER中選取Pressure-Based求解,采用標準的k-ε湍流模型,各值保持默認設置。為使計算加快收斂,在SOLVECONTROLS中,采用SIMPLEC算法來計算壓力-速度耦合(Pressure-Velocity Coupling)。在差分離散方法中,壓力方程選用標準格式(Standard)。為使單位表面計算值有二階精度,采用二階迎風格式(Second Order Upwind)計算動量方程、湍動能方程與耗散率輸運方程。在迭代計算的過程中,打開殘差監視器,通過殘差監視器判斷計算是否收斂。各項的收斂標準設為0.001。
計算中采用了多重參考系,葉輪為移動參考系,其余為絕對參考系。
對于本計算模型,共需設置三個邊界條件,分別是進口邊界,出口邊界和壁面邊界。
進口邊界設置:根據泵使用的一般情況,把進口邊界設置為質量入口邊界(Mass-Flow-inlet)。由實際的工況決定表中的質量流速和質量流量,總溫,壓力和流動方向等具體數值。湍流計算選取水力直徑和湍流強度,其數值大小根據理論計算給出。
出口邊界設置:因不知道出口的速度和壓力,因此出口邊界設置為自由出流邊界條件(Outflow),保證流動是完全發展的。
壁面邊界條件:壁面邊界設置為無滑移的固壁條件(Wall)。
泵的工況及流體性質參數如表1所示,流體為清水。

表1 泵工況及流體性質參數表

圖3 葉輪壓力分布圖
圖3為模擬計算后,葉輪處的壓力分布圖。由圖中可看出,在流動方向上,葉輪內的壓力是逐漸增加的。壓力最大值出現在泵體內的周圍附近,壓力在導葉后的區域數值略有下降。由此分析,在流動方向上,葉輪內的壓力的增加,這是由于葉輪以旋轉的方式做功。在泵體內,流體的能量變化由動能變為勢能,出現最大壓力值。在導葉后,流體因水力損失,導致壓力值有所下降。
圖4為模擬計算后,葉輪處的流場分布圖。從圖中可以看出在葉輪處流體的三維空間速度大小分布和流向趨勢。由圖中可看出,葉輪的絕對速度最大為16.7m/s。在葉輪出口附近,靠近葉片吸力面和葉輪前盤附近,特別是在這兩個面的交匯處附近,流體的相對速度比較低,形成尾跡區。在壓力面附近的流體的相對速度比較高,形成了射流區,出現射流-尾跡結構。同時,由于導葉葉頭處的干擾,葉輪的出流在該處出現了與葉輪運動方向相反的逆流。在導葉內,流體的速度在凸面略高于凹面。

圖4 葉輪的流場分布圖
對于離心泵,其揚程,泵級軸功率,泵效率可用下列公式計算。
泵揚程計算公式[1]:

式中:Pout為流體在出口處單位面積上的平均總壓力,Pa;
Pin為流體在入口處單位面積上的平均總壓力,Pa;
ρ為流體密度,kg/m3;
g為重力加速度,m/s2。
泵級軸功率計算公式[2]:

式中:M為作用在葉輪上的轉軸的力矩,N·m;
n為泵的轉速,rpm。
泵效率計算公式[1]:

式中,Q為泵的體積流量,m3/s。
根據有關文獻[3],在計算揚程和效率時,容積效率ηV取0.96,圓盤摩擦效率ηm取0.96。
根據上述公式,結合模擬數值,得出泵的流量-揚程,流量-軸功率和流量-效率預測曲線。同時,為把模擬計算得到的預測結果與實際的實驗結果作對比,在陽江市新力工業有限公司的水泵測試站對離心泵的性能進行實際測試,得到實際測試的性能曲線。
圖5,圖6,圖7為模擬數值計算結果與實際測試結果的對比。

圖5 流量-揚程曲線圖

圖6 流量-軸功率曲線圖

圖7 流量-效率曲線圖
從以上圖可見,模擬預測值和實際測量值吻合十分好。模擬預測值能比較準確地反應泵的實際工作狀況。
1)使用FLUENT模擬軟件能有效地描述離心泵葉輪的工作狀況;
2)通過模擬計算預測曲線和實際測試數值十分吻合,說明在泵設計的參數優化過程中,模擬計算軟件能起到很好的輔助作用。
[1]王秀勇,王燦星.基于數值模擬的離心泵性能預測[J].流體機械,2007,35(10):9-13.
[2]張淑佳,朱保林,林鋒,等.基于仿真分析的離心泵特性曲線計算[J].中國機械工程,2006(8):46-48.
[3]關醒凡.泵的理論與設計[M].北京:機械工業出版社,1987.