黃書才，穆春玉，楊勤，沈飛，方兵

船用多級泵的水力仿真及汽蝕性能優(yōu)化分析

黃書才，穆春玉，楊勤，沈飛，方兵

為預(yù)測DCP200-4型船用多級泵的內(nèi)部流動規(guī)律及提高汽蝕性能，對其內(nèi)部流場進(jìn)行數(shù)值模擬，通過定性和定量分析對誘導(dǎo)輪進(jìn)行優(yōu)化設(shè)計。結(jié)果表明，數(shù)值模擬可直觀形象地分析多級泵內(nèi)部流動規(guī)律，優(yōu)化誘導(dǎo)輪可有效提高汽蝕性能。

多級泵；內(nèi)部流場；數(shù)值模擬；汽蝕性能

多級泵因揚(yáng)程高、占地小等優(yōu)點(diǎn)被廣泛應(yīng)用于石油、化工、電力、農(nóng)業(yè)、礦業(yè)等領(lǐng)域。近年來，隨著數(shù)值仿真技術(shù)的發(fā)展，對多級泵內(nèi)部流場進(jìn)行仿真分析的技術(shù)越來越成熟，眾多學(xué)者取得了一定的研究成果[1-5]。這些研究僅分析了內(nèi)部流場及外特性，有些甚至只針對某一級或某個部件(如葉輪、誘導(dǎo)輪、導(dǎo)葉等)進(jìn)行分析或優(yōu)化，而對全流道進(jìn)行仿真，以及對多級泵汽蝕性能的優(yōu)化研究則鮮見報道[6-10]。

本文研究的DCP200-4型船用多級泵為我公司自主研發(fā)，輸送介質(zhì)為液化石油氣(liquefied petroleum gas，LPG)。石油氣具有易燃、易爆，以及易汽化的特點(diǎn)，通常將其加壓為液態(tài)進(jìn)行運(yùn)輸。因此，該多級泵對汽蝕性能有很高的要求。

對多級泵的內(nèi)部流場進(jìn)行三維數(shù)值模擬，揭示內(nèi)部壓力場、速度場的變化規(guī)律，并通過定性和定量分析，對誘導(dǎo)輪進(jìn)行優(yōu)化設(shè)計，有效降低汽蝕余量，提高汽蝕性能。

1 基本參數(shù)

計算模型是一臺比轉(zhuǎn)速為140.8的立式、單吸、多級(四級)船用離心泵，其設(shè)計參數(shù)見表1。

表1 DCP200-4型船用多級泵的基本設(shè)計參數(shù)

2 數(shù)值計算

2.1 計算域及網(wǎng)格

計算域及網(wǎng)格劃分如圖1所示，包括吸水井、進(jìn)口段、誘導(dǎo)輪、進(jìn)口導(dǎo)葉、四級葉輪、四級正反導(dǎo)葉和出水管路7個部分。該計算域與實際工況完全一致，這樣可以使模擬更加逼近真實場景。對計算區(qū)域采用結(jié)構(gòu)與非結(jié)構(gòu)化網(wǎng)格相結(jié)合的方式進(jìn)行網(wǎng)格劃分，并分別在誘導(dǎo)輪、葉輪和導(dǎo)葉的葉片附近進(jìn)行局部加密。為排除網(wǎng)格因素對模擬結(jié)果的影響，在設(shè)計工況下對網(wǎng)格數(shù)進(jìn)行無關(guān)性檢查，發(fā)現(xiàn)當(dāng)網(wǎng)格總數(shù)為 260萬左右時，泵的揚(yáng)程、功率和效率等均趨于穩(wěn)定，沒有再增加網(wǎng)格的必要。

2.2 控制方程

泵內(nèi)部流動的介質(zhì)為LPG，可視為不可壓粘性流體，其流體運(yùn)動視為不可壓縮流體的湍流流動。整個過程可用連續(xù)性方程、雷諾平均N-S方程、湍動能k，以及湍動能耗散率ε的輸動方程等控制方程來描述[11]。

2.3 計算模型

運(yùn)用CFD技術(shù)對該型船用多級泵的內(nèi)部流場進(jìn)行數(shù)值模擬，計算中采用SIMPLEC算法求解控制方程，采用標(biāo)準(zhǔn)k-ε湍流模型來進(jìn)行數(shù)值模擬，方程的離散方式采用二階迎風(fēng)格式，收斂精度為10-5。

2.4 邊界條件

1)進(jìn)口邊界條件。采用質(zhì)量進(jìn)口邊界條件。

2)出口邊界條件。采用壓力出口邊界條件。

3)壁面邊界條件。絕熱且無滑移壁面。

2.5 性能參數(shù)的計算

表征泵的主要性能參數(shù)有流量、揚(yáng)程、功率、效率及汽蝕余量等。對該泵的揚(yáng)程、功率、效率及汽蝕余量進(jìn)行計算如下[12]。

揚(yáng)程的計算公式為

(1)

功率的計算公式為

(2)

式中：N為泵軸的轉(zhuǎn)矩，Nm；ω為角速度，rad/s。

效率的計算公式為

(3)

式中：Qout為流體的體積流量，m3/s；Hout為單位重量流體出口處的總能頭，m。

汽蝕余量的計算公式為

(4)

其中：u1為葉片進(jìn)口處的圓周速度；s1為葉片前緣點(diǎn)的厚度；s為葉片的最大厚度。

v1、w1和u1值由泵入口的速度數(shù)值計算結(jié)果決定，s1、s由幾何模型決定。那么在一定工況下，泵的汽蝕余量是一定的。式(4)表明，泵汽蝕余量NPSHr只與泵進(jìn)口部分的運(yùn)動參數(shù)(v1、w1、u1)有關(guān)，是由泵本身的條件決定的，而與裝置參數(shù)無關(guān)。NPSHr越小，表示壓力降小，因而泵的抗汽蝕性能越好。所以NPSHr的物理意義是表示液體在泵進(jìn)口部分的壓力下降程度，也就是為了保證泵不發(fā)生汽蝕，要求在泵進(jìn)口處單位重量液體具有超過汽化壓力水頭的富余能量。

判別泵的汽蝕性能的標(biāo)準(zhǔn)可以有2個，一個是上面提到的經(jīng)驗公式(4)，NPSHr的值越小泵的汽蝕性能越好，v1和w1的值可以在數(shù)值模擬的結(jié)果中提取，但系數(shù)μ和λ的取值不一樣，得到的汽蝕余量將有誤差。這個標(biāo)準(zhǔn)可以作近似的定量分析。另一個是根據(jù)相同計算邊界的數(shù)值結(jié)果得到的泵內(nèi)壓力最低處(一般都在位于葉片背面進(jìn)口稍后處)的壓力值，最低壓力值越大的，汽蝕性能越好，這個標(biāo)準(zhǔn)可以作定性分析。

3 結(jié)果分析

3.1 壓力場分析

圖2為額定工況時的壓力場。通過分析可知，由于每級葉輪的出口均與正反導(dǎo)葉的進(jìn)口相連，正反導(dǎo)葉包裹在葉輪外部，所以經(jīng)過葉輪增壓后，隨著級數(shù)的增加，正反導(dǎo)葉的靜壓值呈現(xiàn)均勻遞增的趨勢，并在末級正反導(dǎo)葉達(dá)到最大值。也就是說，隨著級數(shù)的增加，揚(yáng)程呈現(xiàn)規(guī)律性的增加，滿足了多級泵適用于高揚(yáng)程的設(shè)計需求。

3.2 流場分析

圖3為額定工況時的流線分布。通過分析可知，整體流動狀態(tài)是良好的，流線都比較光順，沒有明顯的折斷、沖擊、回流和漩渦等影響性能的現(xiàn)象，只是在進(jìn)口導(dǎo)葉和出口管路的入口處有較紊亂的局部流動。這是由其局部結(jié)構(gòu)變化較大造成的，并且不可避免。

3.3 誘導(dǎo)輪壓降分析

本文研究的多級泵對汽蝕性能有極高的要求，設(shè)計點(diǎn)的汽蝕余量只有1.3 m，而CFD數(shù)值模擬要比較精確地確定其汽蝕性能有一定難度，并且誘導(dǎo)輪設(shè)計的優(yōu)劣直接關(guān)乎整機(jī)的汽蝕性能。鑒于此，需對誘導(dǎo)輪進(jìn)行反復(fù)優(yōu)化，最終形成了下述3個誘導(dǎo)輪的設(shè)計方案，見圖4。

如前所述，誘導(dǎo)輪的最低壓力值可對汽蝕性能作定性分析，即最低壓力值越大的，汽蝕性能越好。由此可定性的判斷，方案002的汽蝕性能最優(yōu)，其次為方案001，方案003的汽蝕性能最差。

3.4 性能曲線分析

圖5為3種誘導(dǎo)輪的性能曲線。由圖5可見，方案001和方案003的H-Q曲線、η-Q曲線，以及P-Q曲線均具有較高的重合度，說明二者具有相似的水力性能，且曲線變化平滑，無突變；而方案002的H-Q曲線、η-Q曲線，以及P-Q曲線與其余2個方案相差較大，揚(yáng)程高、效率低、軸功率大。因此，單從水力性能考慮，應(yīng)優(yōu)選方案001和方案003。

通過對比3種方案的NPSHr數(shù)值(見表2)以及其對應(yīng)的NPSHr-Q曲線(見圖5)，可以發(fā)現(xiàn)，在設(shè)計工況點(diǎn)，3種方案的NPSHr值均達(dá)了設(shè)計要求。方案001的NPSHr最小，其次為方案002，方案003的NPSHr最大。由此可定量的判斷，方案001的汽蝕性能最優(yōu)。

表2 3種誘導(dǎo)輪方案在設(shè)計點(diǎn)的NPSHr

3.5 優(yōu)化結(jié)果

從定性分析的結(jié)果來看，方案002具有最佳的汽蝕性能，然而經(jīng)定量分析，方案002的汽蝕性能介于方案001和方案003之間，并且其揚(yáng)程、效率、軸功率均偏離設(shè)計點(diǎn)較遠(yuǎn)；從定量分析的結(jié)果來看，方案001的汽蝕性能最優(yōu)，從定性分析來看，其汽蝕性能也僅次于方案002。由此，綜合考慮定性分析和定量分析的結(jié)果，以及兼顧揚(yáng)程、效率和功率等外特性，最終選取方案001作為最優(yōu)方案。

4 結(jié)論

1)隨著多級泵級數(shù)的增加，正反導(dǎo)葉的靜壓值呈現(xiàn)均勻遞增的趨勢，該變化規(guī)律可為在設(shè)計階段較好地預(yù)測多級泵的流動性能提供理論依據(jù)。

2)進(jìn)口導(dǎo)葉和出口管路的入口處，有局部流動較紊亂，如何最大限度地減小該區(qū)域的紊亂現(xiàn)象，是未來的工作方向之一。

3)通過定性和定量分析，對誘導(dǎo)輪進(jìn)行了優(yōu)化設(shè)計，有效降低了汽蝕余量。該方法可為設(shè)計高汽蝕性能的多級泵提供理論依據(jù)

[1] 江見福,顧伯勤,邵春雷.多級泵內(nèi)部流動分析及其性能預(yù)測[J].南京工業(yè)大學(xué)學(xué)報:自然科學(xué)版，2012,34(5):94-98.

[2] 張學(xué)靜,楊軍虎.多級泵內(nèi)部流場的三維數(shù)值模擬及性能預(yù)測[J].流體機(jī)械，2011,39(8):24-28.

[3] 宋志光.多級離心泵導(dǎo)葉流動性能及優(yōu)化設(shè)計的研究[D].廣州：華南理工大學(xué),2014.

[4] 薛睿,李家文,唐飛.多級泵級間導(dǎo)葉的優(yōu)化與數(shù)值仿真[J].火箭推進(jìn)，2011,37(5):24-29.

[5] 黃思,桑迪科.多級離心泵三維流場數(shù)值模擬及性能預(yù)測[J].機(jī)械科學(xué)與技術(shù),2010,29(6):705-708.

[6] 石建偉,宋文武,金永鑫,等.多級離心泵葉輪水力模型優(yōu)化及性能分析[J].中國農(nóng)機(jī)化學(xué)報,2016,37(3):214-217.

[7] 王昌生，李志鵬，陳芳芳,等.多級離心泵首級葉輪內(nèi)部流動的數(shù)值模擬[J].水電能源科學(xué),2012,30(6):147-150.

[8] 張學(xué)靜,楊軍虎.多級離心泵首級水力模型的流場研究及優(yōu)化設(shè)計[J].中國農(nóng)村水利水電,2014,12:104-108.

[9] 周邵萍,胡良波,張浩.多級離心泵級間導(dǎo)葉性能優(yōu)化[J].農(nóng)業(yè)機(jī)械學(xué)報,2015,46(4):33-39.

[10] 叢小青,張嘯,袁丹青,等.多級中開式離心泵級間過水流道的流場分析[J].排灌機(jī)械工程學(xué)報,2015(10):840-845.

[11] 王福軍.計算流體力學(xué)分析：CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004.

[12] 關(guān)醒凡.現(xiàn)代泵理論與設(shè)計[M].北京：中國宇航出版社,2011.

(武漢船用機(jī)械有限責(zé)任公司，武漢 430084)

Hydraulic Simulation and Cavitation Performance Optimization of a Marine Multistage Pump

HUANG Shu-cai, MU Chun-yu,YANG Qin, SHEN Fei, FANG Bing

(Wuhan Marine Machinery Plant Co., Ltd, Wuhan 430084, China)

In order to predict the inner flow motion and raise cavitation performance of DCP200-4 marine multistage pump, numerical simulation was adopted to analyze its internal flow field, and optimize the design of the inducer. The results showed that numerical simulation could analyze the inner flow law of multistage pumps well, and the optimized inducer could sublimate its cavitation performance effectively.

multistage pump; inner flow; numerical simulation; cavitation performance

10.3963/j.issn.1671-7953.2017.01.014

2016-09-13

國家發(fā)改委項目(發(fā)改辦高技[2015]1409號)

黃書才(1986—)，男，碩士，助理工程師研究方向:流體機(jī)械設(shè)計研發(fā)

U664.5

A

1671-7953(2017)01-0058-04

修回日期：2016-10-24