王永生，彭云龍，劉承江

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

0 引 言

仿生學(xué)研究不僅是模仿生物具體的自然結(jié)構(gòu)，而是強調(diào)要充分利用在自然界中發(fā)現(xiàn)的特獨現(xiàn)象和功能［1］。所以最廣泛地運用類比或模擬方法是仿生學(xué)研究方法的精髓［2］。本研究利用類比或模擬方法來指導(dǎo)高速巡邏艇噴水推進(jìn)器進(jìn)水流道的結(jié)構(gòu)設(shè)計，獲得了適合研究對象的流體性能優(yōu)良的進(jìn)水流道。

噴水推進(jìn)主要部件是進(jìn)水流道和噴水推進(jìn)泵［3］(簡稱噴泵)(見圖1)。進(jìn)水流道為噴泵提供了從船底吸水、將作功介質(zhì)水流輸送至噴泵進(jìn)口處的傳輸通道，噴泵把驅(qū)動軸傳遞給它的機械功通過葉輪做功和導(dǎo)葉整流變成作功介質(zhì)水流的壓能和動能［4］。高壓高速大流量的水流從噴泵噴口射出時的反作用力即為船舶推力。

圖1 噴水推進(jìn)器的基本組成Fig.1 Components of waterjet

要保持噴水推進(jìn)高效率，不僅進(jìn)水流道內(nèi)的流動損失要小，噴泵效率要高，而且“船體-進(jìn)水流道-噴泵”流體匹配性能要最佳。本文介紹用仿生學(xué)原理指導(dǎo)流體匹配性能最佳的適合于高速巡邏艇的進(jìn)水流道的空間形狀(流體結(jié)構(gòu))設(shè)計。

1 研究對象及設(shè)計任務(wù)

設(shè)計任務(wù)是為某高速噴水推進(jìn)巡邏快艇提供航速不低于80 km/h 的雙軸噴水推進(jìn)系統(tǒng)。在該航速下噴水推進(jìn)巡邏艇已處于滑行狀態(tài)，噴水推進(jìn)器軸線已遠(yuǎn)高于靜浮狀態(tài)時的水線。

依據(jù)船模試驗所得到的船體阻力、主機功率和轉(zhuǎn)速，運用艦船動力裝置和噴水推進(jìn)原理，對噴水推進(jìn)泵進(jìn)行選型，得到設(shè)計航速時噴泵的選型參數(shù):流量、揚程、比轉(zhuǎn)速、汽蝕比轉(zhuǎn)速、噴速比、推進(jìn)效率和預(yù)報航速。根據(jù)上述這些噴泵選型參數(shù)，運用三元理論進(jìn)行噴泵的水動力設(shè)計(見圖2);依據(jù)船體設(shè)計提供型線和限制條件，如流道高度、流道長度，運用參數(shù)化設(shè)計方法，設(shè)計進(jìn)水流道幾何結(jié)構(gòu)(見圖3)。

圖2 用三元方法設(shè)計得到泵噴三維幾何結(jié)構(gòu)Fig.2 Geometry of the pump designed by 3D method

圖3 用參數(shù)化方法設(shè)計的進(jìn)水流道幾何Fig.3 Geometry of duct designed by parametric method

運用已被驗證的計算流體力學(xué)方法(CFD)數(shù)值計算了噴泵和進(jìn)水流道的敞水性能，計算得到噴泵和進(jìn)水流道的敞水效率分別為91.5%和96.1%。將噴泵和進(jìn)水流道安裝于該巡邏艇船體上并對“船體-進(jìn)水流道-噴泵”系統(tǒng)進(jìn)行數(shù)字自航并使泵噴產(chǎn)生推力與船體阻力相等。對數(shù)值計算結(jié)果進(jìn)行后處理以深入分析在設(shè)計航速、設(shè)計轉(zhuǎn)速時“船體-進(jìn)水流道-噴泵”流場特點、水力性能和推進(jìn)性能。圖4 為設(shè)計航速、設(shè)計轉(zhuǎn)速時進(jìn)水流道和噴泵內(nèi)的流動情況。

2 問題分析及解決辦法

從圖4 可以看到，從圖右的船頭方向沿船底流至進(jìn)水流道內(nèi)的水流，沒有充滿流道背部。從進(jìn)水口進(jìn)入進(jìn)水流道的水流只沖到軸線高度后就進(jìn)入了噴泵。由于流體沒有充滿整個進(jìn)水流道，致使流道背部產(chǎn)生低壓區(qū)，流體流動不均勻并存在大量渦流，故流道流動損失顯著增大;又由于流量不足(遠(yuǎn)低于設(shè)計流量)，故噴泵做功能力顯著下降，揚程、功率和推力遠(yuǎn)小于設(shè)計值，快速性預(yù)報結(jié)果發(fā)現(xiàn)該船未達(dá)設(shè)計航速。

圖4 流道主流區(qū)曲線化截取Fig.4 Main stream area curve extraction

分析上述現(xiàn)象即水流未能充滿進(jìn)水流道的原因，可以概括為3 點:1)進(jìn)水流道軸線偏高;2)進(jìn)水流道縱向長度偏短;3)流道背部形狀不適合航速為80 km/h 滑行艇的高速進(jìn)流。該艇航速80 km/h時已處在滑行狀態(tài)，此時噴泵軸線高度已遠(yuǎn)高于船靜浮狀態(tài)時的水線，且船體有3°左右的縱傾角，所以從船底前駐點到進(jìn)水口之間的有效濕表面積的邊界層縱向長度很短，造成進(jìn)水流道進(jìn)水口處的進(jìn)流幾乎無邊界層。換句話說，水流幾乎是直接沖入進(jìn)水口的。由于進(jìn)流速度很高，基于高速流動慣性的原因，又由于滑行狀態(tài)時進(jìn)水流道軸線高出水面而進(jìn)水流道縱向長度又偏短，故沖入進(jìn)水口的水流無法或來不及在垂直方向充滿流道背部時就從流道出口流出了。所以，要增大流量提高流體做功能力的措施之一是降低流道高度，增大流道縱向長度。

第3 個原因即“流道背部形狀不適合高速進(jìn)流”的問題如何解決呢?仔細(xì)觀察進(jìn)水流道內(nèi)主體水流最上方流線形狀(圖4 中標(biāo)注的a-b-c- d線)。從船頭駐點沿船底向后流動的水流，流入進(jìn)水流道后，沿著流線a-b-c-d 的流動應(yīng)該是該航速下流動距離最短、流動損失最小的流動軌跡。如果流道背部形狀和該流線一致，流道內(nèi)的流動有可能既使流體能充滿整個流道且流動損失可最小。

依據(jù)上述設(shè)想作者重新設(shè)計了進(jìn)水流道并按以下步驟完成了進(jìn)水流道背部形狀的類比設(shè)計或形狀模擬:1)在確定的坐標(biāo)系中用曲線取點數(shù)字軟件完成了進(jìn)水流道內(nèi)主體水流最上部的流線a-b-c-d 取值;2)確定步驟1 中背部形狀曲線的放大倍數(shù)等于進(jìn)水流道出口截面頂部高度H/主體水流流出進(jìn)水流道出口時的頂部流體高度h(見圖4)。該放大倍數(shù)的意義是在流道高度保持不變時增大流道長度且流道背部曲線按該船高速流體進(jìn)入流道后的流線形狀來設(shè)計。

3 仿生流道水動力性能數(shù)值預(yù)報

本文用計算流體力學(xué)方法來預(yù)報利用仿生原理所得到的進(jìn)水流道(簡稱仿生流道)的流體動力性能。為了說明所用CFD 數(shù)值方法的可信性，先用一個有驗證數(shù)據(jù)的“國外某噴泵+國外某進(jìn)水流道”系統(tǒng)的流體性能數(shù)值計算結(jié)果的校驗來證明。然后用該CFD 方法分析國外某進(jìn)水流道的敞水性能，仿生流道的敞水性能，以及“仿生流道+某巡邏艇噴泵”的推進(jìn)性能。

3.1 “國外噴泵+國外進(jìn)水流道”流體性能數(shù)值模擬

本校核對象為標(biāo)稱直徑0.71 m 的國外某混流式噴水推進(jìn)泵，包括葉輪、導(dǎo)葉和噴口三部分，葉輪葉片數(shù)6 片，導(dǎo)葉葉片數(shù)11 片，葉輪輪轂前與泵軸銜接處加有導(dǎo)流帽以改善進(jìn)流平滑過渡?；炝魇絿娝七M(jìn)泵幾何如圖5所示。

圖5 某國外混流式噴水推進(jìn)泵幾何形狀Fig.5 Geometry model of an overseas mixed waterjet pump

由于六面體網(wǎng)格具有較高的求解精度和較精細(xì)刻劃邊界層的優(yōu)勢，故采用全六面體網(wǎng)格對噴泵計算域進(jìn)行網(wǎng)格劃分。葉輪和導(dǎo)葉的幾何結(jié)構(gòu)是周期性的，所以劃分網(wǎng)格只需進(jìn)行單通道［5］。葉輪和導(dǎo)葉幾何采用Ansys 中的專業(yè)旋轉(zhuǎn)機械網(wǎng)格劃分模塊TurboGrid 劃分，葉輪扭曲度較大采用J 型拓?fù)浣Y(jié)構(gòu)，導(dǎo)葉采用H 型拓?fù)浣Y(jié)構(gòu)，葉片附近用O-Grid 捕捉葉片表面邊界層(圖6 ～圖7)。該混流式噴水推進(jìn)泵的總網(wǎng)格數(shù)為250 萬。

圖6 葉輪網(wǎng)格Fig.6 Mesh of rotor

圖7 導(dǎo)葉網(wǎng)格Fig.7 Mesh of stator

根據(jù)廠家提供幾何圖紙，對與該混流式噴水推進(jìn)泵配套使用的進(jìn)水流道進(jìn)行幾何建模，并與噴泵裝配成整體的“噴水推進(jìn)泵-進(jìn)水流道-船體”流體系統(tǒng)以進(jìn)行數(shù)字自航計算，本研究對艇體作了如下簡化:用平板代替船底。

對噴水推進(jìn)器進(jìn)行流場數(shù)值模擬，需要考慮來流速度和船底邊界層對進(jìn)流的影響。通常用的方法是在流道進(jìn)水口下方建立一個適當(dāng)大小的區(qū)域模擬船底流場，根據(jù)經(jīng)驗一般長度為噴泵進(jìn)口直徑20倍，寬度為10 倍，深度為8 倍［6］(見圖8)。

圖8 噴水推進(jìn)器數(shù)值計算區(qū)域Fig.8 Numerical computational domain of waterjet

流道網(wǎng)格在ICEM 中劃分網(wǎng)格，對流動變化劇烈區(qū)域進(jìn)行網(wǎng)格加密。近壁面處采用O-grid，一方面提高網(wǎng)格質(zhì)量，另一方面方便調(diào)整近壁面網(wǎng)格厚度以精確模擬近壁面流動。第一層網(wǎng)格厚度設(shè)置為0.2 mm，壁面Y+值保持在250 以內(nèi)以滿足湍流模型要求。流道網(wǎng)格如圖9所示。湍流模型與噴水推進(jìn)泵裸泵數(shù)值計算一樣選取SST 模型，整個計算域網(wǎng)格總數(shù)353 萬。計算所得推力采用壁面積分法獲取，與廠商提供噴水推進(jìn)器推力和功率比較結(jié)果見表1(進(jìn)行歸一化處理)。誤差在工程允許范圍內(nèi)，證明本文所用數(shù)值模擬方法的可信性。

圖9 進(jìn)水流道網(wǎng)格Fig.9 Mesh of waterjet duct

表1 推力計算值與噴泵制造廠商提供的數(shù)據(jù)的比較Tab.1 Comparison of calculation result and manufacture data

圖10 是噴水推進(jìn)器在設(shè)計航速時噴泵和進(jìn)水流道內(nèi)流線圖，進(jìn)水流道內(nèi)流通順暢，無流動分離或漩渦。

圖10 設(shè)計航速時噴水推進(jìn)器流線Fig.10 The streamline of waterjet at the design speed

3.2 國外某進(jìn)水流道的敞水性能的數(shù)值計算

實際情況中，水流流過噴水推進(jìn)器時由于水流粘性、變速變向、壁面摩擦等因素存在流動損失的(稱為流道損失)。而流動損失大小客觀反映了該流道性能的好壞，因此引入流道效率的概念。流道效率是衡量進(jìn)水流道對來流能量利用程度的物理量，定義為流道出口面總能量與進(jìn)流面總能量之比［7］:

式中:E1、E3分別為流道進(jìn)流面和出口面的總能量。故流動損失的定義為ξ=(E1－E3)/E1。

流道流動特性另一個重要指標(biāo)是流道出口不均勻度系數(shù)ξ，它的作用是評判進(jìn)水流道出流(即噴泵進(jìn)流)的不均勻度，其定義式為［8］:

式中:Q 為流道出口流量;u 為CFD 各流體單元出口處的軸向速度;為出口處的平均軸向速度。不均勻度ξ 越小，表示流道出流越均勻，流動性越好。

利用上述網(wǎng)格和數(shù)值計算模型對國外某流道進(jìn)行敞水性能數(shù)值模擬，知該進(jìn)水流道敞水效率93.2%，流道出口出流不均勻度0.26。圖11 為流道內(nèi)部流線圖，可以看到流道內(nèi)基本充滿且流動順暢無漩渦。

圖11 敞水流道流線圖Fig.11 The streamline of duct

3.3 仿生流道敞水性能數(shù)值計算及與國外流道對比

根據(jù)類比或模擬原理來設(shè)計的仿生流道，其網(wǎng)格如圖12所示。同樣用CFD 方法計算其流道敞水性能，與國外某流道敞水性能對比結(jié)果如表2，可見仿生流道敞水性能優(yōu)于從國外進(jìn)口的進(jìn)水流道，不論在流道效率或是出口不均勻度上都有顯著提高。

圖12 仿生流道網(wǎng)格Fig.12 The mesh of a bionicduct

表2 國外某與仿生流道性能對比Tab.2 Comparison of the overseas duct and the bionic duct

3.4 “艇+仿生流道+噴泵”流體數(shù)值計算及對比

將為某巡邏艇設(shè)計的噴泵與仿生流道安裝到某高速巡邏艇的實尺度數(shù)值模型上，進(jìn)行數(shù)字自航以預(yù)報推進(jìn)性能，并與原流道性能進(jìn)行對比。建立控制域如圖13所示，整個計算域長度為7 倍船長，船頭前部1倍船長，船尾后5 倍船長，計算域?qū)挾葹? 倍船長。采用VOF 的方法來追蹤自由液面，湍流模型采用SST模型，對流項離散選用二階離散格式［9］。采用四面體和六面體網(wǎng)格對不同的區(qū)域進(jìn)行網(wǎng)格劃分，整體計算區(qū)域網(wǎng)格量總數(shù)為2.2 千萬。由于2 臺浸沒式噴水推進(jìn)器布置相對船體中心面對稱，為減小計算時間，計算時只對一半結(jié)構(gòu)進(jìn)行流場計算。

表3 為同一噴泵兩型流道(原流道直徑D=0.5 m，流道高度H=0.577 m，流道長度L=2.5 m;現(xiàn)仿生流道直徑D=0.5 m，流道高度H=0.46 m，流道長度L=3.3 m)推進(jìn)性能對比結(jié)果。由表中可知，仿生流道流量、功率、推力較原流道均有顯著增加，有利于實現(xiàn)噴泵設(shè)計的推進(jìn)性能。

圖13 “船+噴泵+流道”計算控制域Fig.13 Whole computational domain of“hull-pump-duct”

表3 推進(jìn)系統(tǒng)性能對比Tab.3 Propulsive properties of the original and the bionic duct

4 結(jié) 語

用常規(guī)設(shè)計方法為航速80 km/h 的某高速艇設(shè)計了敞水效率分別為91.5%和96.1%的噴泵和進(jìn)水流道。將它們和艇體進(jìn)行設(shè)計航速、設(shè)計轉(zhuǎn)速時的數(shù)字自航，發(fā)現(xiàn)進(jìn)水流道和噴泵內(nèi)未能充滿水流，造成噴泵流量、揚程、功率和所產(chǎn)生的推力遠(yuǎn)低于設(shè)計值。為什么部件性能均優(yōu)(艇體阻力小、噴泵效率高、流道流動損失小)但系統(tǒng)集成后推進(jìn)性能反而不優(yōu)呢?究其原因，是它們的流場匹配性能不優(yōu)。具體是:流道高度過高、流道長度過短、流道背部形狀這三者不適用于高速進(jìn)流。在降低流道高度、增大流道長度的同時，應(yīng)用仿生學(xué)原理——類比或模擬方法——設(shè)計流道背部形狀以及它與船底的銜接，最終結(jié)果是噴水推進(jìn)性能顯著改善:流量增加35%，功率增大55%，推力增大61%?？偨Y(jié)本研究工作有2 點:1)在進(jìn)行系統(tǒng)設(shè)計時，不僅追求部件性能優(yōu)良，更重要追求系統(tǒng)性能優(yōu)異。性能優(yōu)良的部件構(gòu)成系統(tǒng)后并不能保證系統(tǒng)性能一定滿足要求，原因是系統(tǒng)部件之間的匹配性能至關(guān)重要。2)觀察和利用自然界中某些特殊流動現(xiàn)象，運用仿生學(xué)原理即類比或模擬方法來輔助性能優(yōu)良的進(jìn)水流道形狀或結(jié)構(gòu)設(shè)計，是一項有意義和實用價值的研究和工程應(yīng)用。

［1］http://en.wikipedia.org/wiki/Bionics.

［2］http://www.baike.com/wiki.

［3］金平仲.船舶噴水推進(jìn)［M］.北京:國防工業(yè)出版社，1988.

JIN Ping-zhong.Marine waterjet propulsion［M］.Beijing:National Defence Industry Press，1988.

［4］靳栓寶.高性能噴水推進(jìn)器和泵噴推進(jìn)器的設(shè)計與性能分析［D］.武漢:海軍工程大學(xué)，2013:23-24.

JIN Shuan-bao.Dseign and analysis of high performance waterjet and pumpjet［D］.Wuhan:Naval University of Engineering，2013:23-24.

［5］彭云龍.基于三元理論和CFD 方法的浸沒式噴水推進(jìn)泵水力設(shè)計研究［D］.武漢:海軍工程大學(xué)，2013:28-29.

［6］PENG Yun-long.Hydraulic design of submerged waterjet pump based on three-dimensional thory and CFD method［D］.Wuhan:Naval University of Engineering，2013:28-29.

［7］劉承江，王永生，張志宏.噴水推進(jìn)器數(shù)值模擬所需流場控制體的研究［J］.水動力學(xué)研究與進(jìn)展，2008，23(5):592-595.

LIU Cheng-jiang，WANG Yong-sheng，ZHANG Zhi-hong.Study on flow control volume in numerical simulation of waterjet propulsor［J］.Chinese Journal of Hydrodynamics，2008，23(5):592-595.

［8］KEEGAN D，MARTIN D，MICHAEL E，et al.Use of RANS for waterjet analysis of a high-speed sealift concept vessel［C］.First International Symposium on Marine Propulsors.Trondheim，Norway:SMP′09，2009.

［9］Norbert Willem Herman Bulten.Numerical analysis of a waterjet propulsion system［D］.Eindhoven:Eindhoven University，2006.

［10］ARASH E.Investigation of Waterjet/Hull Interaction Effects［D］.Gothenburg Sweden:Chalmers University of Technology，2012.