，，，

(海軍工程大學(xué) 動力工程學(xué)院，武漢 430033)

噴水推進(jìn)在高速高性能船舶領(lǐng)域得到廣泛應(yīng)用[1-3]。噴水推進(jìn)器主要由噴水推進(jìn)泵(以下稱為噴泵)和進(jìn)水流道兩部分構(gòu)成。噴泵是噴水推進(jìn)核心部件，噴泵對效率、揚程、流量、抗空化等方面有較高要求，設(shè)計難度較大[4-5]。近年來，計算機技術(shù)和計算流體力學(xué)(computational fluid dynamics，CFD)發(fā)展很快，與傳統(tǒng)水力模型試驗相比，其周期短，耗資較少，CFD在船舶噴水推進(jìn)技術(shù)研發(fā)、預(yù)報、設(shè)計中的應(yīng)用日趨廣泛。在以前葉片進(jìn)口邊位置和形狀主要憑經(jīng)驗來設(shè)計，對性能曲線的影響只有試驗后才知曉[6]。在理論上，通過泵氣蝕余量公式，認(rèn)為修圓葉片進(jìn)口外圓部分可以提高葉片的抗汽蝕性能[7]。楊軍虎等[8]利用CFD軟件FLUENT對低轉(zhuǎn)速軸流泵進(jìn)行過葉片外緣修圓對性能的影響進(jìn)行過研究，重點對葉頂間隙的流場進(jìn)行了分析。國外Angelo Cervone、Giovanni Pace等利用水力試驗，研究發(fā)現(xiàn)導(dǎo)邊幾何形狀對空化特性影響較大。對于高轉(zhuǎn)速軸流式噴水推進(jìn)泵，空化是非常重要的考慮因素，葉片導(dǎo)邊形狀對該空化性能和水力性能之間的影響研究報道仍比較少。為了研究導(dǎo)邊形狀對水力效率和空化性能的影響情況，文中運用CFD方法，通過改變某高速摩托艇噴水推進(jìn)泵導(dǎo)邊形狀，研究導(dǎo)邊形狀對該類噴水推進(jìn)泵敞水效率和空化性能的影響,為高轉(zhuǎn)速泵的設(shè)計、性能優(yōu)化提供基本參考。

1 控制方程與湍流模型

1.1 控制方程

流體流動須遵守質(zhì)量守恒規(guī)律，動量守恒定律和能量守恒定律?？刂品匠碳礊檫@些守恒定律的數(shù)學(xué)描述，包括連續(xù)性方程和動量方程(也稱Navier-Stokes方程，簡稱N-S方程)。在定常條件下，選用不可壓縮的三維連續(xù)方程和RANS方程模擬噴水推進(jìn)泵內(nèi)流場。

1)質(zhì)量守恒方程。

▽·(ρV)=0

(1)

2)動量守恒方程(N-S方程)。

(2)

3)能量守恒方程。

(3)

1.2 湍流模型

利用有限體積法對RANS偏微分方程進(jìn)行離散，進(jìn)而對其進(jìn)行數(shù)值求解。所有算例都采用SST(shear stress transport)模型來封閉方程組。SST模型結(jié)合了κ-ε模型和κ-ω模型的優(yōu)點, 在近壁面區(qū)域，采用κ-ω模型模擬，湍流耗散較小，收斂性好。在湍流充分發(fā)展區(qū)域，采用κ-ε模型，計算效率高，對復(fù)雜流場的適應(yīng)性更好。兩種模型的結(jié)合能夠較好地模擬粘性底層的流動[9]。SST模型的最大優(yōu)點就在于考慮了湍流剪切應(yīng)力，從而不會對渦流黏度造成過度預(yù)測[10]。研究采用高階SST模型。

1.3 CFD方法驗證

在對該高速摩托艇軸流式噴泵進(jìn)行CFD計算前，先以KaMeWa71SⅡ型噴水推進(jìn)器的混流式噴水推進(jìn)泵為例對其水動力特性進(jìn)行數(shù)值模擬，并與廠商提供的噴水推進(jìn)器數(shù)據(jù)進(jìn)行比較，以驗證所用方法和模型的可信性。

KaMeWa71SII型噴水推進(jìn)器的混流式噴水推進(jìn)泵幾何模型見圖1。

圖1 KaMeWa71SII幾何模型

該噴泵葉輪進(jìn)口直管處直徑為710 mm，葉輪有6片，導(dǎo)葉11片。該噴泵流體性能數(shù)值計算的幾何分為噴口、導(dǎo)葉、葉輪以及本文為計算敞水性能所添加的進(jìn)口直管。噴口和進(jìn)口直管都用ICEM軟件，劃分六面體網(wǎng)格。導(dǎo)葉部和葉輪部用TurboGrid軟件劃分六面體網(wǎng)格，對葉輪進(jìn)行網(wǎng)格無關(guān)性計算，計算結(jié)果見圖2。

圖2 網(wǎng)格無關(guān)性計算

經(jīng)計算葉輪網(wǎng)格至少需要70萬，本研究該葉輪網(wǎng)格實取98.3萬。其他部件的網(wǎng)格數(shù)見表1。

邊界條件為總壓進(jìn)口，靜壓出口。時間步按CFX推薦設(shè)置，設(shè)為物理時間1/ω(ω為轉(zhuǎn)速，(°)/s)。

噴水推進(jìn)泵軸功率計算結(jié)果見表2。由表2可見，計算值與廠商提供數(shù)據(jù)相對誤差均在1.5%以內(nèi)，滿足工程要求。故從工程應(yīng)用的角度來看，該方法計算水動力特性是可信的。

表1 各部分網(wǎng)格數(shù)

表2 CFD計算值與試驗值對比

2 CFD模擬計算

研究對象為某高速摩托艇噴水推進(jìn)器的噴水推進(jìn)泵，其轉(zhuǎn)速較高，額定轉(zhuǎn)速為7 500 r/min。僅以該軸流式泵體的敞水性能為研究對象。該噴水推進(jìn)泵為單級軸流泵，通過改變導(dǎo)邊的修圓情況，研究導(dǎo)邊對該泵性能的影響。

2.1 幾何建模

根據(jù)摩托艇噴水推進(jìn)器幾何，利用CAD三維建模軟件UG進(jìn)行幾何建模。噴泵分為葉輪體、導(dǎo)葉體和噴口段3部分。噴泵葉輪葉片為3片，導(dǎo)葉為6片，噴口段也對應(yīng)設(shè)置了6片二級導(dǎo)葉以加強整流效果。該軸流式噴泵尺寸較小，葉輪直徑為155.5 mm，葉頂間隙為0.5 mm。導(dǎo)葉體和噴口段幾何模型見圖3，網(wǎng)格劃分見圖4。

圖3 噴水推進(jìn)泵幾何模型

圖4 噴水推進(jìn)泵網(wǎng)格

對于該噴泵葉輪的建模分為修圓導(dǎo)邊幾何與未修圓導(dǎo)邊的幾何。原噴泵幾何為葉輪a。葉輪a較未修圓幾何葉片(以下稱為葉輪c)在導(dǎo)邊形狀上有較大改變，導(dǎo)邊形狀為“c”字形。葉輪b修圓程度介于葉輪a、c之間，見圖5。

圖5 葉輪a、葉輪b、葉輪c單通道幾何

2.2 網(wǎng)格劃分

葉輪和導(dǎo)葉用ANSYS 12.0軟件中專為旋轉(zhuǎn)機械劃分網(wǎng)格的CFX-TurboGrid軟件。由于葉輪扭曲程度較大，故采用J型拓?fù)?，?dǎo)葉采用H型拓?fù)浣Y(jié)構(gòu)。噴口由于包含了二級導(dǎo)葉，幾何較為復(fù)雜，用ICEM軟件，四面體網(wǎng)格進(jìn)行劃分，在壁面邊界用菱柱形網(wǎng)格進(jìn)行加密，對細(xì)小部件用密度盒加密，確保對幾何描述的準(zhǔn)確性。在對葉輪葉片網(wǎng)格無關(guān)性分析中，通過對葉輪a單通道網(wǎng)格節(jié)點數(shù)設(shè)定變化，對計算的軸力矩進(jìn)行監(jiān)測，力矩變化波動見圖6。

圖6 網(wǎng)格無關(guān)性計算

經(jīng)計算葉輪總網(wǎng)格至少需70.8萬節(jié)點，本研究實取80.5萬以上。其他部分取較大的網(wǎng)格數(shù)以保證網(wǎng)格滿足要求。具體網(wǎng)格數(shù)和節(jié)點數(shù)見表3。

表3 具體網(wǎng)格數(shù)和節(jié)點數(shù)

2.3 邊界條件設(shè)置

對噴泵敞水性能進(jìn)行計算時，為保證均勻流動，在葉輪進(jìn)口前加直管進(jìn)流道，管徑與葉輪進(jìn)口直徑相等，長度在2倍直徑以上。采用邊界條件為流量進(jìn)口，靜壓出口邊界條件，葉輪葉片與輪轂為相對靜止壁面條件，泵殼設(shè)為絕對靜止壁面條件。在計算空化性能時，在原有基礎(chǔ)上加入空化模型進(jìn)行計算。

3 計算結(jié)果分析

噴泵的敞水性能指孤立泵在均勻進(jìn)流條件下的性能，即不加進(jìn)水流道的水力性能。模擬計算噴泵在7 500 r/min的情況下，噴泵葉輪分別采用葉輪a、b和c計算流量190～230 kg/s的軸功率、揚程和效率，見圖7～9。

圖7 軸功率對比

圖8 揚程對比

圖9 效率對比

分析發(fā)現(xiàn)，葉輪a、b和c的功率和揚程都有微小變化。葉輪a的功率最大，葉輪c的功率最小。因為該軸流式噴泵的葉片包角較大，修圓導(dǎo)邊后造成葉片做功面積減小，但進(jìn)流角相對于葉輪c會變大，在這兩種作用的綜合作用下，造成該泵的功率在葉輪a的情況下所需功率最大。

分析還發(fā)現(xiàn)修圓導(dǎo)邊對該裸泵的工作點效率影響微小，但在小于工況點的流量時，效率明顯低于未修圓導(dǎo)邊。分析認(rèn)為葉輪外緣修圓后，泄漏流動加強，降低了做功效率。

為進(jìn)一步計算導(dǎo)邊形狀對空化性能的影響，在原計算模型的基礎(chǔ)上加入空化模型，分析葉片在最高效率工況點，流量為210 kg/s時的空化情況見圖10。

圖10 空化對比

從空化的對比圖中可以發(fā)現(xiàn)，葉片空化都發(fā)生在葉稍和導(dǎo)邊部分，葉輪c的空化面積明顯最多，葉輪a空化面積最小，這與用理論分析的情況一致。一方面，流進(jìn)葉輪的液體先與導(dǎo)邊中部接觸，獲得能量后一部分液體沿著導(dǎo)邊向葉稍方向流動，旋轉(zhuǎn)能量增加了，壓力也增加了，延遲了空化的發(fā)生；另一方面通過后處理葉輪的流場矢量，可以發(fā)現(xiàn)在葉稍處流動非常劇烈，葉輪a的弧形導(dǎo)邊相對于葉輪c的直線導(dǎo)邊可以明顯地減少來流的沖擊，削弱了二次流的影響，見圖11。

圖11 葉輪流場矢量

葉頂部分是做功和流動最強烈的地方，葉片外緣導(dǎo)邊部分修圓以后，造成葉片葉稍卸載，減小了空化發(fā)生面積。

4 結(jié)論

1)對于高轉(zhuǎn)速的軸流泵，修圓導(dǎo)邊的葉片對軸流泵的功率和揚程有微小的影響，會造成小于工作點流量時的效率降低；

2)修圓導(dǎo)邊的葉片在空化性能上優(yōu)于直導(dǎo)邊的葉片，空化面積明顯減少，有利于葉片葉稍卸載。

[1] BULTEN N.Numerical analysis of a waterjet propulsion system[D]. Eindhoven, The Netherlans:The Eindhoven University of Technology,2006.

[2] 劉承江,王永生,丁江明等.現(xiàn)代噴水推進(jìn)裝置的演變[J].艦船科學(xué)技術(shù)，2006,28(4):8-12.

[3] ALLISON J L.Marine waterjet propulsion[J].Transaction of SNAME,1993,101:275-335.

[4] ALLISON J L,JIANG C B.Modern tools for waterjet pump design and recent advances in the field[C]∥Internationa Conference on Waterjet Propulsion Ⅱ.RINA,Amsterdam,1998:128-142.

[5] 湯方平.噴水推進(jìn)軸流泵設(shè)計及紊流數(shù)值分析[D].上海：上海交通大學(xué)，2006.

[6] 李世煌.葉片泵的非設(shè)計工況及其優(yōu)化設(shè)計[M].北京：機械工業(yè)出版社，2006.

[7] 關(guān)醒凡.軸流泵和斜流泵[M].北京：中國宇航出版社，2009.

[8] 楊軍虎，馬靜先，張人會,等.軸流泵葉片外緣修圓對泵性能的影響[J].蘭州理工大學(xué)學(xué)報，2007,33(5):60-63.

[9] 靳栓寶.基于CFD噴水推進(jìn)泵性能分析及優(yōu)化設(shè)計[J].船海工程，2010，39(2)：125-129.

[10] 謝龍漢，趙新宇，張炯明.ANSYS CFX流體分析及仿真[M].北京：電子工業(yè)出版社，2012.