張 宇， 管仁偉， 安 偉*， 徐峰

(1．江南大學(xué)機(jī)械工程學(xué)院，江蘇無錫 214122;2．錫泵制造有限公司，江蘇無錫 214000)

目前我國廣大地區(qū)的水利抽排水工程、城市防汛工程，以及鄉(xiāng)村防洪抗旱工程等，大多采用固定泵站。其特點(diǎn)是容量大，投資大，建成后固定在確切位置來完成送水功能。但有許多防汛抗旱或需臨時排水場合沒有固定泵站，這時機(jī)動靈活的移動泵站是對固定泵站最好的補(bǔ)充，在很多情況下具有固定式泵站無法比擬的優(yōu)勢。但應(yīng)急搶險移動泵站相對于固定式泵站受到了很多的限制，不能發(fā)揮其所有性能。因此，對應(yīng)急搶險移動泵站用輕載型軸流泵的設(shè)計(jì)與優(yōu)化，研究軸流泵內(nèi)葉輪導(dǎo)葉的流場機(jī)理，對其效率、揚(yáng)程等性能進(jìn)行深入研究，以提高其性能指標(biāo)，不僅具有理論意義，更具有實(shí)際的應(yīng)用價值。但是目前國內(nèi)關(guān)于軸流泵葉輪葉型優(yōu)化的文獻(xiàn)相對較少，僅是對泵內(nèi)流場進(jìn)行模擬。因此，提出泵葉輪葉型優(yōu)化方法，完善軸流泵的設(shè)計(jì)方法，成為當(dāng)前需要解決的難題。

1 數(shù)值模擬及可靠性驗(yàn)證

1．1 建模及網(wǎng)格劃分

利用三維建模商業(yè)軟件Pro/E建立輕載型軸流泵的模型，該模型的結(jié)構(gòu)參數(shù)符合應(yīng)急搶險移動泵站用輕載型軸流泵的要求。圖1為根據(jù)其結(jié)構(gòu)的幾何參數(shù)和葉片木模圖生成的簡化軸流泵實(shí)體。

圖1 模型Pro/E造型Fig．1 Pro/E modeling

計(jì)算網(wǎng)格采用分區(qū)分塊網(wǎng)格生成的方法，將復(fù)雜的幾何區(qū)域分割成若干簡單的小塊，然后將每一小塊區(qū)域單獨(dú)生成合適的網(wǎng)格，最后將所有的網(wǎng)格耦合起來。葉輪導(dǎo)葉區(qū)域利用Numeca中專門針對葉輪旋轉(zhuǎn)機(jī)械部分網(wǎng)格生成器Autogrid5進(jìn)行劃分，設(shè)置葉輪網(wǎng)格屬性為泵用模板的拓?fù)浣Y(jié)構(gòu)、軸流，設(shè)置轉(zhuǎn)速、葉片數(shù)、第一層的網(wǎng)格厚度等參數(shù)。彎管、喇叭管區(qū)域在Numeca子模塊igg中劃分，采用蝶形網(wǎng)格。最后在igg中將三部分進(jìn)行耦合。見圖2。

圖2 Numeca中生成的模型Fig．2 Numeca modeling

圖3為葉輪導(dǎo)葉三維固壁網(wǎng)格，圖4為全流道的三維結(jié)構(gòu)網(wǎng)格。

圖3 葉輪導(dǎo)葉三維固壁網(wǎng)格Fig．3 3D mesh of impeller and guide vane

圖4 全流道三維結(jié)構(gòu)網(wǎng)格Fig．4 3D structured mesh modelling

其中葉輪導(dǎo)葉劃分單元的網(wǎng)格總數(shù)892 221，彎管和喇叭管劃分單元的網(wǎng)格總數(shù)412 548，網(wǎng)格總數(shù)1 304 769，后面的數(shù)值模擬收斂狀況良好，網(wǎng)格質(zhì)量也較好。

1．2 數(shù)值計(jì)算方法和邊界條件

利用CFD軟件——Numeca中的Fine/turbo建立了基于真實(shí)幾何模型的輕載型軸流泵內(nèi)部的三維湍流數(shù)值模擬的計(jì)算平臺，使用時間相關(guān)法求解Reynolds時均方程。為了快速計(jì)算各種復(fù)雜流場，該軟件采用多塊/多重網(wǎng)格的計(jì)算技術(shù)，并有多種可供選擇的運(yùn)動方程和湍流模型。計(jì)算工質(zhì)為清水，運(yùn)動方程為湍流N-S方程，對流項(xiàng)采用基于TVB量差分分裂方法的高精致迎風(fēng)格式，其中心格式采用了Jameson人工粘性顯格式，上風(fēng)格式采用基于TVD與通量差分分裂方法的高精致格式，方程求解采用多階Runge-Kutta法，并應(yīng)用多重網(wǎng)格與隱式殘差光順法等加速收斂技術(shù)，湍流模型采用Spalart-Allmaras模型。計(jì)算時給定進(jìn)口流體速度，工質(zhì)流入方向，靜溫;出口給定靜壓;固壁hub，葉片blade設(shè)定轉(zhuǎn)速，其它固壁設(shè)置靜止。

1．3 數(shù)值計(jì)算結(jié)果及可靠性分析

圖5為不同工況下軸流泵葉輪的總體性能模擬數(shù)據(jù)和泵性能試驗(yàn)數(shù)據(jù)的對比結(jié)果。

圖5 不同工況下水泵揚(yáng)程及效率曲線Fig．5 Pump lift and efficiency curve under the different conditions

從對比結(jié)果可以看出，揚(yáng)程和總效率都相差很小，其設(shè)計(jì)工況在工程允許的誤差范圍5%范圍內(nèi)。但由于計(jì)算模型中忽略了葉頂間隙的影響，使泵的效率和揚(yáng)程在設(shè)計(jì)工況均有不同程度的增加。在流量較小的工況，不容易收斂，誤差相對較大，但已經(jīng)可以充分說明數(shù)值模擬的正確性和合理性［1］。

1．4 葉輪導(dǎo)葉內(nèi)部的靜壓分析

圖6給出了設(shè)計(jì)工況下輪轂、輪蓋和葉輪導(dǎo)葉的靜壓分布云圖。

圖6 葉輪導(dǎo)葉區(qū)域靜壓云圖Fig．6 Static pressure diagram

由圖6可以看出，靜壓分布在圓周方向，具有較好的對稱性;工質(zhì)從進(jìn)口到出口，壓力呈現(xiàn)增大趨勢，在出口處壓力明顯得到提升。流體進(jìn)入葉輪以后，軸向速度增加，導(dǎo)致葉片吸力面一側(cè)的靜壓降低，尤其是葉輪葉頂一側(cè)的速度最快，流體繞過葉片的進(jìn)口端時，由于急速轉(zhuǎn)彎，使流速進(jìn)一步加快，在葉片吸力面進(jìn)口邊葉頂處形成全流場壓力的最低值，該位置與在實(shí)際運(yùn)行時發(fā)生汽蝕的位置是相吻合的，由此也可以證明計(jì)算的正確性和合理性［2-4］。如圖7 所示。

圖7 葉輪葉片的汽蝕Fig．7 Cavitation performance of pump

1．5 葉輪內(nèi)部的流動分析

為了研究速度矢量沿著Z軸的變化趨勢，在設(shè)計(jì)工況下從葉輪和導(dǎo)葉體出口端之間做了4個垂直于Z軸的截面，從圖8中可以大體了解到，進(jìn)口處的速度以旋轉(zhuǎn)流動為主，旋轉(zhuǎn)流動基本上順著流道沿徑向減弱，到出口處，軸向速度大于徑向速度。在靠近葉輪葉片出口邊的截面2處可以看到流動分離和漩渦的跡象。

圖8 不同截面處的相對速度矢量分布Fig．8 Relative velocity vector distribution in different section

文中研究的軸流泵葉輪葉片屬于扭曲比較嚴(yán)重的葉片，葉根附近的流動狀況混亂、復(fù)雜，葉片的壓力面附近的根部存在流動分離，葉根的出口處產(chǎn)生了明顯的漩渦，近壁面處的流動分離已經(jīng)嚴(yán)重影響到葉片流道主區(qū)域的流動，導(dǎo)致?lián)p失增大，效率降低，見圖9。

圖9 速度極限流線圖Fig．9 Speed limit flow diagram

可見，葉根出口處明顯的沖擊和馬蹄渦的存在，嚴(yán)重影響了流道的近壁面流動，葉輪葉根尾部結(jié)構(gòu)設(shè)計(jì)可能存在不合理性［5］。分析原因，可能是由于葉片的出口角和葉型壓力面型線的彎曲不符合流動規(guī)律，因而導(dǎo)致馬蹄渦的形成，致使馬蹄渦后的分支在葉片兩側(cè)產(chǎn)生了漩渦和分離流動［6-7］。針對上述的流動問題，文中對旋轉(zhuǎn)葉輪提出葉型修正的方案——優(yōu)化葉輪葉片的吸壓力面的型線，以期能夠消除或者減少不必要的流動損失，提高葉輪的水力性能。

2 葉型優(yōu)化

2．1 葉片幾何參數(shù)化擬合

根據(jù)提出的修正方案，只針對葉輪葉片的吸壓力面的型線進(jìn)行優(yōu)化，其它的結(jié)構(gòu)均不變。具體控制條件設(shè)定:輪轂線、輪蓋線設(shè)置20點(diǎn)控制的B-Spline樣條曲線;流面設(shè)置從左到右的默認(rèn)迭代;積疊線位置設(shè)置在進(jìn)口邊，葉片子午面定位設(shè)置為10點(diǎn)控制的Bezier曲線;葉片的進(jìn)口邊和出口邊均設(shè)置成圓角;周向定位設(shè)置為10點(diǎn)控制的Bezier曲線;中弧線設(shè)置為3點(diǎn)控制的Bezier曲線;葉片的壓力面和吸力面型線均采用3點(diǎn)控制的Bezier曲線。圖10為葉輪葉根處葉型最后的擬合結(jié)果。

圖10 葉輪葉根處葉型參數(shù)化擬合結(jié)果Fig．10 Leaf type parametric fitting results of blade root

2．2 數(shù)據(jù)庫生成

在生成數(shù)據(jù)庫之前，首先要清楚選擇什么幾何參數(shù)作為優(yōu)化的幾何約束(自由參數(shù))，盲目選擇自由參數(shù)只會為數(shù)據(jù)生成工作和后期的優(yōu)化工作帶來極大的難度。本課題中的葉片共7個截面，設(shè)置截面1和截面7的共6個控制點(diǎn)作為可變參數(shù)。對于中間截面，將其和截面1和截面7的變化關(guān)聯(lián)起來，用關(guān)系式表達(dá)，使從根部到頂部逐漸變化。這樣既可以減少自由參數(shù)的數(shù)量，又不至于由于各個截面自由變化而導(dǎo)致葉片變成不切實(shí)際的形狀。最后設(shè)置隨機(jī)數(shù)取樣模式，樣本數(shù)選取30。數(shù)值模型利用igg多級腳本文件，采用和原型泵相同的模型。

2．3 優(yōu)化方法和過程

OPtimization優(yōu)化過程實(shí)際上是對數(shù)據(jù)庫中的一系列的葉型所對應(yīng)的結(jié)果進(jìn)行近似函數(shù)擬合，并利用ANN(人工神經(jīng)網(wǎng)絡(luò))的理論尋找擬合曲線的最優(yōu)值，并對最優(yōu)值對應(yīng)的幾何形狀進(jìn)行流場分析，實(shí)現(xiàn)循環(huán)尋優(yōu)的過程。在此過程中，它不斷地向已經(jīng)存在的數(shù)據(jù)庫文件中添加新的樣本。對每一個幾何形狀都要進(jìn)行幾何建模、網(wǎng)格生成、流場分析以及后處理計(jì)算等步驟。

對數(shù)據(jù)庫樣本進(jìn)行優(yōu)化時，選擇葉輪的水力效率和進(jìn)出口壓差作為目標(biāo)函數(shù)，將水力效率值設(shè)置為1，即效率的目標(biāo)是1，進(jìn)出口的壓差設(shè)置為不小于29 066．00 Pa，數(shù)據(jù)庫尋優(yōu)步數(shù)設(shè)置為10。

2．4 葉輪葉片優(yōu)化的結(jié)果

下文提取了具有代表性的4個樣本的結(jié)果進(jìn)行參考，最后一個為最優(yōu)葉輪樣本，見圖11。

圖11 優(yōu)化樣本Fig．11 Optimization sample

圖12和圖13為葉根、葉頂截面上的葉片優(yōu)化前后的葉型圖。

從以上兩圖可以看出，壓力面和吸力面的型線發(fā)生了變化。

表1為優(yōu)化前后原型葉輪和最優(yōu)葉輪總體性能參數(shù)的對比。

表1 原型葉輪、最優(yōu)葉輪對比Tab．1 Numerical simulation results comparison between original impeller and best impeller

可以看出，在設(shè)計(jì)工況下，優(yōu)化設(shè)計(jì)后葉輪的水力效率相對提高了2．73%，揚(yáng)程提高了8．83%，表明最優(yōu)葉輪的水力性能要大大好于原型葉輪。

改進(jìn)葉型后，原葉輪葉片在吸力面根部形成的流動分離和漩渦基本消失，流道中的通道渦和出口處沖擊受到抑制，葉根進(jìn)口處的相對速度矢量分布趨于合理，近壁區(qū)域流動趨于良好，葉輪的水力效率得到很大的提升，說明改型措施正確有效［8-10］。文中的優(yōu)化設(shè)計(jì)證明，控制葉輪壓力面和吸力面的型線可以控制葉片出口邊的馬蹄渦，消除近壁面流動分離和漩渦，減小流動中的通道渦強(qiáng)度和影響范圍，提高軸流泵整體的水力效率。

3 結(jié)語

1)運(yùn)用 Numeca軟件中的前處理模塊 igg、autogrid5，計(jì)算模塊Fine/trubo，對輕載型軸流泵進(jìn)行了全流道的三維數(shù)值模擬，數(shù)值模擬的結(jié)果與試驗(yàn)數(shù)據(jù)較好吻合，證明了文中針對該型軸流泵所選用的數(shù)值模擬方法和計(jì)算模型的準(zhǔn)確性和可靠性。

2)根據(jù)模擬的結(jié)果，分析了流動中存在的問題，探討了葉輪葉型結(jié)構(gòu)的不合理處，提出了改進(jìn)的方法，利用design3D對葉型進(jìn)行了優(yōu)化設(shè)計(jì)，優(yōu)化效果明顯。在水力效率相對提高2．73% 的同時，水泵的揚(yáng)程也大大增加，原型葉輪中存在的漩渦和壁面脫離現(xiàn)象基本消失，葉根進(jìn)口沖擊和流道渦團(tuán)影響范圍減弱，葉輪的水力性能顯著提高。

3)研究表明，漩渦和流動分離是高比轉(zhuǎn)速、高扭曲軸流式葉輪損失的主要原因，驗(yàn)證了通過改變壓力面和吸力面的葉型形狀，可以控制流道內(nèi)各個渦團(tuán)的形成和結(jié)構(gòu)，減小或者消除流道內(nèi)不必要的流動損失，提高軸流泵葉輪的水力性能。

4)針對當(dāng)前國內(nèi)外對輕載型軸流泵葉輪葉型優(yōu)化的相關(guān)研究較少的情況，提出泵葉輪葉型優(yōu)化的方法，完善了軸流泵的設(shè)計(jì)和優(yōu)化方法。

［1］楊軍虎，張煒，王春龍．潛水軸流泵全流道三維湍流數(shù)值模擬及性能預(yù)估［J］．排灌機(jī)械2006(4):105-107．YANG Jun-hu，ZHANG Wei，WANG Chun-long．Numerical simulation of turbulent flow through full passage of submersible axial pump and performance predictions［J］．Drainage and Irrigation Machinery，2006(4):105-107．(in Chinese)

［2］郎繼興．軸流泵輪軸面速度與環(huán)量最佳分布規(guī)律的初探［J］．流體工程，1990(12):31-36．LANG Ji-xing．Axial flow pump shaft surface speed and circulation optimal distribution of the trial［J］．Journal of Fluids Engineering，1990(12):31-36．(in Chinese)

［3］楊昌明，陳次昌，千金諾．軸流泵的三維造型與數(shù)值模擬［J］．西南交通大學(xué)學(xué)報(bào)，2004，19(2):139-141．YANG Chang-ming，CHEN Ci-chang，QIAN Jin-ro．The 3 d modeling and numerical simulation of axial flow pump［J］．Journal of Southwest Jiaotong University，2004，19(2):139-141．(in Chinese)

［4］Goto A，Nohmi M，Sakumi T，et al．Hydrodynamic design system for pumps based on 3-D CAD，CFD，and inverse design method［J］．Journal of Fluids Engineering，2002，124(2):329-335．

［5］陸林廣．立式軸流泵裝置的三維湍流流動數(shù)值模擬［J］．排灌機(jī)械，2007，25(1):29-32．LU Lin-guang．Numerical simulation of 3D turbulent flow vertical axia-flow pump system［J］．Drainage and Irrigation Machinery，2007，25(1):29-32．(in Chinese)

［6］彭玉成．基于CFD數(shù)值解析的軸流泵優(yōu)化設(shè)計(jì)［J］．水泵技術(shù)，2004(5):24-26．PENG Yu-cheng．Based on the CFD numerical analytical axial flow pump optimization design sail［J］．Pump Technology，2004(5):24-26．(in Chinese)

［7］趙斌娟，王澤．離心泵葉輪內(nèi)流數(shù)值模擬的現(xiàn)狀和展望［J］．農(nóng)機(jī)化研究，2002(3):34-36．ZHAO Bin-juan，WANG Ze．Centrifugal pump impeller internal flow numerical simulation of the present situation and prospect［J］．Journal of Agricultural Mechanization Research，2002(3):34-36．(in Chinese)

［8］談明高，劉厚林，王勇，等．葉輪外徑對離心泵內(nèi)流影響的 CFD分析［J］．排灌機(jī)械，2009，27(5):314-318．TAN Ming-gao，LIU Hou-lin，WANG Yong，et al．CFD analysis on effects of impeller outlet diameter on flow field in centrifugal pump［J］．Drainage and Irrigation Machinery，2009，27(5):314-318．(in Chinese)

［9］Anagnostopoulos J S．A fast numerical method for flow analysis and blade design in centrifugal pump impellers［J］．Computers and Fluids，2009，38(2):284-289．

［10］戴哲敏，李玉仁．基于CFD的軸流泵內(nèi)部流場模擬［J］．中國陶瓷工業(yè)，2008，15(5):35-37．DAI Ze-min，LI Yu-ren．CFD-Base internal flow field simulation of axial flow pump［J］．China Ceramic Industry，2008，15(5):35-37．(in Chinese)