李鵬程，蔡佑林，汲國(guó)瑞，李 寧

(1. 中國(guó)船舶及海洋工程設(shè)計(jì)研究院, 上海 200011；2. 噴水推進(jìn)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 上海 200011；3. 上海交通大學(xué) 海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 上海 200011)

基于三元可控速度矩法的水下航行體推進(jìn)泵參數(shù)化設(shè)計(jì)方法

李鵬程1,2，蔡佑林1,2，汲國(guó)瑞1,3，李 寧1,2

(1. 中國(guó)船舶及海洋工程設(shè)計(jì)研究院, 上海 200011；2. 噴水推進(jìn)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 上海 200011；3. 上海交通大學(xué) 海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 上海 200011)

隨著水下航行體對(duì)推進(jìn)器高性能和低噪聲要求與日俱增，迫切需要研究適用于水下航行體的推進(jìn)泵。本文選出影響水下航行體推進(jìn)泵水動(dòng)力性能的主要參數(shù)，基于三元可控速度矩方法，編寫Matlab程序?qū)ν七M(jìn)泵進(jìn)行參數(shù)化設(shè)計(jì)，并結(jié)合CFD數(shù)值模擬方法進(jìn)行設(shè)計(jì)優(yōu)化與試驗(yàn)驗(yàn)證。CFD計(jì)算結(jié)果表明模型泵水力效率能達(dá)到83%，達(dá)到設(shè)計(jì)目標(biāo)并準(zhǔn)確把握了設(shè)計(jì)工況點(diǎn)，驗(yàn)證了本方法的可靠性，為開發(fā)高性能水下航行體推進(jìn)泵奠定基礎(chǔ)。

可控速度矩；水下航行體推進(jìn)泵；前置導(dǎo)葉；CFD數(shù)值模擬

0 引 言

水下航行體推進(jìn)泵不同于傳統(tǒng)的軸流泵與混流泵，其輪轂為水下航行體尾部線型光順延伸，輪緣與輪轂構(gòu)成收縮型流道；為使結(jié)構(gòu)緊湊，降低葉輪導(dǎo)葉間脈動(dòng)壓力，常具有前置導(dǎo)葉。在國(guó)內(nèi)，水下航行體推進(jìn)泵的設(shè)計(jì)研究仍處于起步階段。由于軍事或商業(yè)保密原因，該泵的研究很少見(jiàn)諸報(bào)道。戴原星[1]嘗試過(guò)用升力法設(shè)計(jì)前置導(dǎo)葉軸流泵，結(jié)果表明升力法對(duì)前置導(dǎo)葉泵的設(shè)計(jì)效果不好。2005年，蔡佑林等[2]根據(jù)文獻(xiàn)中的混流泵可控渦設(shè)計(jì)理論，提出了混流泵設(shè)計(jì)的三元可控速度矩設(shè)計(jì)新方法。2009年，吳志旺[3]將三元可控速度矩法應(yīng)用于高比轉(zhuǎn)速混流泵設(shè)計(jì)，開展了幾種不同環(huán)量分布條件下的葉輪設(shè)計(jì)，試驗(yàn)結(jié)果表明該方法可有效用于混流泵設(shè)計(jì)。

本文繼承與發(fā)展了用于混流泵設(shè)計(jì)的三元可控速度矩方法，將其應(yīng)用到水下航行體推進(jìn)泵的設(shè)計(jì)中，編寫Matlab程序，形成了一套水下航行體推進(jìn)泵的參數(shù)化設(shè)計(jì)方法。該參數(shù)化設(shè)計(jì)方法可以根據(jù)設(shè)計(jì)要求靈活地調(diào)整各個(gè)參數(shù)，使得每次調(diào)整優(yōu)化都有很強(qiáng)針對(duì)性，大大縮減了設(shè)計(jì)與優(yōu)化時(shí)間。

1 參數(shù)化設(shè)計(jì)方法

1.1 數(shù)學(xué)模型

葉輪內(nèi)S2m流面上沿正交線n方向的速度梯度方程：

圖 1 兩類坐標(biāo)系Fig. 1 Two coordinate systems

對(duì)于推進(jìn)泵內(nèi)不可壓縮流體，流過(guò)泵內(nèi)各截面流量相等，得到連續(xù)性方程：

1.2 水下航行體推進(jìn)泵主要控制參數(shù)

1）S2m流面可控速度矩分布

為使速度矩沿徑向和流線變化均勻，在葉輪進(jìn)口邊滿足來(lái)流無(wú)沖擊條件，出口邊根據(jù)庫(kù)塔條件，Vur沿流線的梯度應(yīng)為 0。本文采用四次函數(shù)來(lái)擬合速度矩沿流線從進(jìn)口邊到出口邊分布，即可得到整個(gè)S2m流面的速度矩分布，如圖2所示。水平方向 2 個(gè)軸分別代表徑向和軸向，豎直方向軸為速度矩Vur?；B表示葉輪進(jìn)口邊速度矩分布，弧BC表示速度矩沿進(jìn)口邊到出口邊的變化規(guī)律。

對(duì)于前置導(dǎo)葉，其進(jìn)口處水流方向?yàn)檩S向，速度矩為 0，出口處存在為葉輪提供的負(fù)預(yù)旋，各剖面速度矩大小理論上應(yīng)為葉輪進(jìn)口速度矩。但考慮到導(dǎo)葉到葉輪間流體能量損失并且為便于調(diào)整導(dǎo)葉各剖面預(yù)旋大小，本文將通過(guò)導(dǎo)葉出口各剖面速度矩系數(shù)Ki來(lái)控制導(dǎo)葉出口速度矩分布，即

圖 2S2m流面速度矩分布圖Fig. 2 The velocity moment distribution ofS2mstream surface

式中，角標(biāo)i表示徑向第i條流線、葉片第i剖面。如此以來(lái)，能夠設(shè)計(jì)出與葉輪更匹配的前置導(dǎo)葉，使得葉輪出口周向速度得到更好的回收，從而提高泵效率。

2）葉輪、導(dǎo)葉軸面弦長(zhǎng)

葉片數(shù)一定的條件下，葉片弦長(zhǎng)大，摩擦損失增加，水力效率降低，葉片弦長(zhǎng)小，葉輪做功能力下降，為提供相同的揚(yáng)程，需增大葉片拱度或安放角，同時(shí)由于弦長(zhǎng)小，厚度比要增加，會(huì)導(dǎo)致空泡發(fā)生的危險(xiǎn)加大。葉根部做功能力差，主要考慮強(qiáng)度和低阻要求，適當(dāng)減少弦長(zhǎng)，葉梢附近剖面做功能力強(qiáng)可適當(dāng)增大弦長(zhǎng)，考慮到葉梢處流速大易發(fā)生空泡，故在葉梢處要求適當(dāng)減少弦長(zhǎng)。本文采用二次函數(shù)來(lái)控制葉輪和導(dǎo)葉沿葉根到葉梢的軸面弦長(zhǎng)。如圖3是Matlab運(yùn)行出的葉輪軸面弦長(zhǎng)分布結(jié)果。

3）導(dǎo)葉與葉輪匹配

圖 3 葉輪軸面弦長(zhǎng)分布圖Fig. 3 Axial chord length distribution of impeller

水下航行體推進(jìn)泵具有前置導(dǎo)葉，提供負(fù)預(yù)旋，導(dǎo)葉出口需要提供葉輪進(jìn)口所需速度矩。前置導(dǎo)葉與葉輪能否合理匹配大大影響到推進(jìn)泵水動(dòng)力性能，導(dǎo)葉與葉輪間距太大水力損失較大，但噪聲性能好；間距太小，則葉輪運(yùn)轉(zhuǎn)不穩(wěn)定。由于導(dǎo)葉安放得根據(jù)艇體尾部結(jié)構(gòu)，位置變動(dòng)范圍受限，本文通過(guò)調(diào)整葉輪位置來(lái)調(diào)整間距。觀察導(dǎo)葉出口處的液流。

1.3 數(shù)值方法

先根據(jù)設(shè)計(jì)參數(shù)用一元設(shè)計(jì)方法得到葉片初始形狀，得到葉片的初始安放角β。將葉輪子午面沿徑向均分為N條流線，軸向均分為M條近似的過(guò)水?dāng)嗝婢€，得到N×M個(gè)控制點(diǎn)，如圖4所示。

圖 4 葉輪軸面流道網(wǎng)格劃分示意圖Fig. 4 Axial flow passage numerical grid

將方程（2）離散，對(duì)于第j條過(guò)水?dāng)嗝婢€流量Qj為：

聯(lián)立方程（6）和方程（7），形成如下矩陣：

式中，ai=2πri,jdnj,i=1,2...N，

1.4 葉片構(gòu)型

解方程（8）得到葉輪與導(dǎo)葉處速度場(chǎng)后，即可通過(guò)葉片構(gòu)型微分方程（9）積分得到葉片展開面骨線型值。

由積分得出的葉片骨線形狀，用NACA66翼型厚度分布規(guī)律對(duì)骨線對(duì)稱加厚，形成葉片展開面的的吸力面與壓力面線。最后將加厚的葉片投影到錐面上就形成三維翼型，在ProE5.0中導(dǎo)入葉輪剖面線，如圖5所示。

圖 5 葉輪剖面線Fig. 5 The section lines of Impeller

1.5 迭代

葉輪和導(dǎo)葉的設(shè)計(jì)方法由基于S2m流面反問(wèn)題與CFD計(jì)算的正問(wèn)題構(gòu)成。得到葉片形狀再用CFD計(jì)算出流場(chǎng)信息，考察葉輪進(jìn)口處速度矩差異，若相差超過(guò)3%，那么將新速度矩分布代入方程（1），重復(fù)反問(wèn)題的計(jì)算。一般迭代 2 步就能滿足收斂條件，便認(rèn)為得到的葉片拱度線就能跟泵內(nèi)真實(shí)流線近似，由此正反問(wèn)題迭代計(jì)算得到水下航行體推進(jìn)泵模型。

2 計(jì)算實(shí)例與CFD模擬方法

現(xiàn)給定設(shè)計(jì)參數(shù)為：流量Q=0.46 m3/s；揚(yáng)程H=5 m；比轉(zhuǎn)速ns=1 074；泵效率η=82%。

根據(jù)本文所述參數(shù)化設(shè)計(jì)方法，在Matlab程序里輸入各主要參數(shù)，直接輸出用于ProE 5.0三維建模的流道，葉輪和導(dǎo)葉的ibl文件和用于TurboGrid的curve文件。形成一套推進(jìn)泵三維模型如圖6所示。紅色是9葉的葉輪，黃色為11葉前置導(dǎo)葉，灰色為輪轂，透明部分為輪緣殼體。

數(shù)值模擬方法選擇RANS方法，湍流模型選用κ-ωSST模型[6]，該模型集合了κ-ε和κ-ω兩種模型的優(yōu)點(diǎn)，在近壁面區(qū)域，采用κ-ω模型模擬；在湍流充分發(fā)展區(qū)域，采用κ-ε模型。為使進(jìn)口和出口為充分發(fā)展的湍流，推進(jìn)泵模型前后均加4D長(zhǎng)度的直管。借助商業(yè)CFD軟件TurboGrid劃分網(wǎng)格，并用CFX求解器進(jìn)行求解。

模擬采用單通道（只含一個(gè)葉片），全結(jié)構(gòu)化網(wǎng)格對(duì)進(jìn)流管、前置導(dǎo)葉、動(dòng)葉輪和出流管進(jìn)行網(wǎng)格劃分。在流場(chǎng)變化劇烈的區(qū)域如壁面、梢隙、葉片導(dǎo)邊隨邊附近對(duì)網(wǎng)格進(jìn)行加密。根據(jù)近壁區(qū)理論，在壁面附近添加邊界層網(wǎng)格，保證壁面網(wǎng)格y+值在30以內(nèi)。進(jìn)流管網(wǎng)格數(shù)量13萬(wàn)，導(dǎo)葉23萬(wàn)，葉輪28萬(wàn)，出流管18萬(wàn)。

圖 6 水下航行體推進(jìn)泵三維圖Fig. 6 The 3D model of the underwater vehicle water-jet pump

圖 7 葉輪和導(dǎo)葉網(wǎng)格Fig. 7 The grid of impeller and guide vanes

采用CFX-Pre13.0對(duì)劃分好的網(wǎng)格按照“進(jìn)流管-前置導(dǎo)葉-導(dǎo)葉-出流管”組裝，形成如圖8所示計(jì)算域。動(dòng)葉輪和導(dǎo)葉的動(dòng)靜耦合問(wèn)題采用多參考系模型，在導(dǎo)葉上建立靜止的慣性坐標(biāo)系，在葉輪上建立隨葉輪一起轉(zhuǎn)動(dòng)的旋轉(zhuǎn)坐標(biāo)系，將泵內(nèi)流場(chǎng)的非定常問(wèn)題轉(zhuǎn)化成定常問(wèn)題。從出口方向看，葉輪以1 450 r/min的轉(zhuǎn)速左旋。

圖 8 網(wǎng)格組裝圖Fig. 8 The mesh assembly

計(jì)算域進(jìn)口邊界條件采用總壓壓力進(jìn)口，給定進(jìn)口總壓為一個(gè)標(biāo)準(zhǔn)大氣壓，出口邊界條件為質(zhì)量流量出口。葉片表面、輪轂、輪緣處設(shè)為固壁邊界條件。在進(jìn)流管、導(dǎo)葉、葉輪、出流管與其他通道的臨界面上，設(shè)置周期性邊界條件。

本文選取設(shè)計(jì)工況點(diǎn)及其附近共 9 個(gè)流量點(diǎn)進(jìn)行數(shù)值計(jì)算，從0.38 m3/s到0.54 m3/s，增量0.02 m3/s，覆蓋0.8～1.2倍設(shè)計(jì)流量范圍，計(jì)算水下航行體推進(jìn)泵的揚(yáng)程、效率等力特性曲線。

3 CFD試驗(yàn)分析

在CFD-post中，在距葉輪進(jìn)出口2D處分別建立平面，分別計(jì)算靜壓和軸向速度。由伯努利方程得到揚(yáng)程計(jì)算公式為：

式中：Pout，V1m分別為葉輪出口2D處?kù)o壓和軸面速度；Pin，V2m為葉輪進(jìn)口2D處?kù)o壓和軸面速度；ρ為水密度，998.2 kg/m3；g為重力加速度，9.81 m/s2；H為泵揚(yáng)程，m。

將各流量點(diǎn)的效率和揚(yáng)程繪制成如圖9所示外特性曲線。

圖 9 0.8～1.2倍流量范圍的外特性曲線Fig. 9 The external characteristic curve between 0.8 to 1.2 times of design mass flow

分析可知，0.8~1.1倍設(shè)計(jì)流量范圍內(nèi)，效率曲線較平坦，高效范圍寬，以該模型為基礎(chǔ)的推進(jìn)裝置適應(yīng)變工況能力強(qiáng)且主機(jī)不易過(guò)載。重要的是，準(zhǔn)確把握了設(shè)計(jì)工況點(diǎn)且實(shí)現(xiàn)了設(shè)計(jì)指標(biāo)。

為考察葉輪后尾流場(chǎng)周向能量回收情況，在尾流后0.1D剖面上的做出周向投影速度圖，如圖10所示。

圖 10 尾流場(chǎng)速度矢量周向投影圖Fig. 10 The velocity vector of circumferential projection at tail flow field

由該投影圖可以看出，葉輪后的周向旋轉(zhuǎn)速度很小，約1 m/s，最大周向速度位于輪緣處，約為1.6 m/s，相比于近9 m/s的軸向流速為相對(duì)小量，周向速度回收效果很好。

4 結(jié) 語(yǔ)

本文形成了一套水下航行體推進(jìn)泵的基于三元可控速度矩法的參數(shù)化設(shè)計(jì)方法，經(jīng)CFD試驗(yàn)驗(yàn)證，設(shè)計(jì)結(jié)果準(zhǔn)確把握了設(shè)計(jì)工況點(diǎn)，效率達(dá)到83%，驗(yàn)證了該設(shè)計(jì)方法的可靠性，為高性能水下航行體推進(jìn)泵提供了優(yōu)秀的初步設(shè)計(jì)模型及初步設(shè)計(jì)方法，得到了如下結(jié)論：

1）混流泵設(shè)計(jì)中的三元可控速度矩方法可以用于水下航行體推進(jìn)泵的設(shè)計(jì)中，設(shè)計(jì)效果很理想。

2）本文方法設(shè)計(jì)的推進(jìn)泵，其葉片剖面不采用傳統(tǒng)設(shè)計(jì)方法用的NACA翼型拱度線，而是由泵內(nèi)流場(chǎng)直接計(jì)算得到的全新的翼型，更貼近真實(shí)流場(chǎng)。

3）只需輸入控制參數(shù)，運(yùn)行程序得到葉片型值后建模并進(jìn)行CFD計(jì)算，總耗時(shí)僅需一兩天，大大縮短了設(shè)計(jì)周期，并能通過(guò)CFD計(jì)算結(jié)果很針對(duì)性的對(duì)泵進(jìn)行控制、優(yōu)化。

[1]戴原星, 王立祥. 基于CFD的前置導(dǎo)葉軸流泵通用特性曲線預(yù)報(bào)[J]. 船舶, 2013, (5): 1–5.

[2]蔡佑林, 王立祥, 張新. 混流泵葉輪三元可控速度矩設(shè)計(jì)[J].流體機(jī)械, 2005, 33(11): 13–15.

[3]吳志旺. 基于環(huán)量分布的高比轉(zhuǎn)速混流泵設(shè)計(jì)[D]. 鎮(zhèn)江: 江蘇大學(xué), 2009.

[4]WU CH, A general theory of three-dimension flow in subsonic and supersonic turbo machines of axial, radial, and mixed-flow type, ASME Paper Number 50-A-79, ASME Transactions, 1952, or NACA TN 2604, 1952, 1–90.

[5]王海松. 軸流泵CAD-CFD綜合特性研究[D]. 北京: 中國(guó)農(nóng)業(yè)大學(xué), 2005.

[6]MENTER F R. Zonal two equation k-w turbulence models for aerodynamic flows[R]. AIAA-93–2906, 1993.

Parametric design method of underwater vehicle water-jet pump applying 3D controllable velocity moment theory

LI Peng-cheng1,2, CAI You-lin1,2, JI Guo-rui1,3, LI Ning1,2
(1. Marine Design &Research Institute of China, Shanghai 200011, China; 2. Laboratory of Science and Technology on Water-jet Propulsion, Shanghai 200011, China; 3. State Key Laboratory of Ocean Engineering, Shanghai Jiaotong University, Shanghai 200011, China)

With the increasing demand of underwater vehicles for the high performance and low noise of propulsion, it is urgent to study the water-jet pump which is suitable for underwater vehicles. According to 3D controllable velocity moment theory, this paper chooses the main parameters which greatly infect the hydrodynamic performance of the underwater vehicle water-jet pump and conduct parametric design through Matlab programming. Design and optimization are carried out with the help of CFD numerical simulation method. The CFD results show that the efficiency of the model pump can reach 83%, and all hydrodynamic performance achieve the design goal, moreover, the design operating point is grasped accurately thus testifying that this parametric design method is reliable and lay the foundation for high performance underwater vehicle water-jet pump.

controllable velocity moment；underwater vehicle water-jet pump；front guide vanes；CFD numerical simulation

U661.1

A

1672 – 7649(2017)07 – 0043 – 05

10.3404/j.issn.1672 – 7649.2017.07.009

2016 – 10 – 26；

2016 – 11 – 21

國(guó)家自然科學(xué)基金資助項(xiàng)目(51579145)

李鵬程(1991 – )，男，碩士研究生，專業(yè)方向?yàn)閲娝七M(jìn)。