嚴(yán)文澤，姚寶運(yùn)，姚會(huì)茹

（大連深藍(lán)泵業(yè)有限公司，大連 116031）

隨著我國(guó)工業(yè)裝置大規(guī)模發(fā)展，在深海油氣輸送、石化、核電、航空航天、冶金等領(lǐng)域內(nèi)，離心泵可輸送介質(zhì)呈現(xiàn)多樣化，特別是在輸送易汽化介質(zhì)和氣體低溫液化輸送（液氧、液氮等）等應(yīng)用領(lǐng)域的不斷拓展，對(duì)離心泵穩(wěn)定運(yùn)行可靠性、輸送經(jīng)濟(jì)性、汽蝕性能都提出了更高要求。

眾多學(xué)者對(duì)離心泵抗汽蝕性能的深入研究，已提出一些行之有效的技術(shù)解決方案，其中，離心泵入口前置誘導(dǎo)輪已成為常用措施，國(guó)外技術(shù)先進(jìn)的泵企業(yè)通過(guò)配置高抗汽蝕誘導(dǎo)輪，使得泵汽蝕比轉(zhuǎn)速達(dá)到了6 000 ～ 7 000。匹配幾何尺寸理想的誘導(dǎo)輪，不但可以顯著改善泵的汽蝕性能[1]，同時(shí)也可以減弱泵汽蝕而導(dǎo)致的振動(dòng)、噪聲及材料磨損及破壞[2]。崔寶玲等[3]通過(guò)CFD 計(jì)算手段對(duì)等螺距誘導(dǎo)輪的內(nèi)部流動(dòng)進(jìn)行數(shù)值模擬，揭示了等螺距誘導(dǎo)輪內(nèi)部壓力場(chǎng)及流場(chǎng)的分布；黎義斌等[4]對(duì)口環(huán)間隙變化對(duì)誘導(dǎo)輪離心泵汽蝕性能影響進(jìn)行研究；孫健等[5]研究了變螺距誘導(dǎo)輪的設(shè)計(jì)參數(shù)選擇，總結(jié)了一些關(guān)鍵參數(shù)的選擇。潘中永[6]對(duì)誘導(dǎo)輪的設(shè)計(jì)進(jìn)行系統(tǒng)論述并給出計(jì)算方 法。

鑒于上述學(xué)者研究以及實(shí)際泵產(chǎn)品性能規(guī)劃及開(kāi)發(fā)，大連深藍(lán)泵業(yè)采用CFD 水力優(yōu)化設(shè)計(jì)手段，理論與試驗(yàn)驗(yàn)證相結(jié)合的方法，使得配置誘導(dǎo)輪的泵產(chǎn)品汽蝕性能達(dá)到國(guó)際領(lǐng)先水平，且其他性能指標(biāo)（效率，振動(dòng)噪聲，曲線(xiàn)形狀等）均滿(mǎn)足要求，同時(shí)形成了誘導(dǎo)輪的水力設(shè)計(jì)規(guī)范及水力設(shè)計(jì)計(jì)算程序。本文以某LNG 國(guó)產(chǎn)化首臺(tái)套泵研發(fā)項(xiàng)目為依托開(kāi)展高抗汽蝕誘導(dǎo)輪設(shè)計(jì)論述。

1 高抗汽蝕誘導(dǎo)輪的設(shè)計(jì)策略及原則

1.1 高抗汽蝕誘導(dǎo)輪的設(shè)計(jì)思路及流程

誘導(dǎo)輪是軸流式葉輪，其具備軸流葉輪的汽蝕特性，在輪緣處產(chǎn)生的氣泡，在沿軸向向前運(yùn)動(dòng)的過(guò)程中，由于輪轂側(cè)液體受離心力作用向外壓的原因，使氣泡被壓縮在輪緣局部區(qū)域并在此凝結(jié)，不會(huì)造成整個(gè)流道堵塞。另外，選擇較大的葉柵稠密度、進(jìn)口過(guò)流面積以及較薄的葉片進(jìn)口厚度和較小的進(jìn)口角，使得誘導(dǎo)輪抗汽蝕性能又整體優(yōu)于軸流葉輪。高抗汽蝕誘導(dǎo)輪的設(shè)計(jì)思路可以概括為：“通過(guò)控制葉片進(jìn)口角及過(guò)流面積保證高抗汽蝕性能、出口角及輪緣直徑保證必須的揚(yáng)程；通過(guò)協(xié)調(diào)設(shè)計(jì)，使得葉柵稠密度處在合理范圍內(nèi)”。

通常情況下，誘導(dǎo)輪的計(jì)算如式（1）、圖1 所示。

圖1 誘導(dǎo)輪與主葉輪的速度三角形示意Fig.1 Velocity triangle between the inducer and the main impeller

在一般情況下，誘導(dǎo)輪的設(shè)計(jì)流程如圖2 所示。

1.2 高抗汽蝕誘導(dǎo)輪結(jié)構(gòu)形式的確定原則

根據(jù)誘導(dǎo)輪汽蝕比轉(zhuǎn)速的限定要求，在設(shè)計(jì)過(guò)程中推薦使用的誘導(dǎo)輪結(jié)構(gòu)形式如表1 所示。

在滿(mǎn)足性能要求，節(jié)約成本的前提下，在設(shè)計(jì)過(guò)程中推薦選用的誘導(dǎo)輪結(jié)構(gòu)形式的優(yōu)先級(jí)為：

圖2 誘導(dǎo)輪設(shè)計(jì)流程Fig.2 Design flow of induction wheel

表1 優(yōu)先推薦使用的誘導(dǎo)輪結(jié)構(gòu)形式Table 1 Preferentially recommends the use of the inducer wheel structure

等螺距誘導(dǎo)輪＞圓柱型變螺距誘導(dǎo)輪＞圓錐型變螺距誘導(dǎo)輪

1.3 誘導(dǎo)輪的優(yōu)選運(yùn)行范圍及進(jìn)口流量系數(shù)的匹配

基于額定流量設(shè)計(jì)的誘導(dǎo)輪，其優(yōu)選的運(yùn)行工作范圍為80% Qbep～ 120% Qbep，當(dāng)流量Q＞1 000 m3/ h時(shí)，不建議誘導(dǎo)輪偏大流量工況使用；在一般情況下，不建議超工況范圍使用，這是由誘導(dǎo)輪的自身空化特性所決定的。圖3 為一般情況下誘導(dǎo)輪的汽蝕特性曲線(xiàn)：

圖3 誘導(dǎo)輪汽蝕特性曲線(xiàn)Fig.3 Cavitation characteristic curve of inducer

存在最優(yōu)的進(jìn)口流量系數(shù)φt，獲得最低的NPSHr，即：

布倫菲爾德（Brumfield）準(zhǔn)則：

式中 λ—— 誘導(dǎo)輪進(jìn)口壓降系數(shù)，通常取0.01 ～ 0.02。

由于誘導(dǎo)輪進(jìn)口沖角的影響，進(jìn)口流量系數(shù)φt匹配關(guān)系間如表2 所示。

表2 進(jìn)口流量系數(shù)的匹配關(guān)系Table 2 Matching relation of inlet flow coefficient

汽蝕比轉(zhuǎn)速C 與流量系數(shù)的關(guān)系式如式（4）所示。

已經(jīng)有大量的文獻(xiàn)資料對(duì)其他關(guān)鍵參數(shù)（葉柵稠密度、輪緣間隙等）進(jìn)行分析研究，在此不再贅述。

2 高抗汽蝕誘導(dǎo)輪的設(shè)計(jì)實(shí)現(xiàn)

2.1 誘導(dǎo)輪水力設(shè)計(jì)及三維模型

按照上述設(shè)計(jì)理論及原則，針對(duì)某LNG 國(guó)產(chǎn)化首臺(tái)套高壓泵，開(kāi)展高抗汽蝕誘導(dǎo)輪的水力開(kāi)發(fā)，其設(shè)計(jì)參數(shù)如表3 和表4 所示。

表3 某國(guó)產(chǎn)化首臺(tái)套LNG 高壓泵水介質(zhì)設(shè)計(jì)參數(shù)Table 3 Design parameters of water medium of the first LNG high-pressure pump

表4 某國(guó)產(chǎn)化首臺(tái)套LNG 高壓泵高抗汽蝕誘導(dǎo)輪部分設(shè)計(jì)幾何參數(shù)Table 4 Design parameters of the first LNG high pressure pump with high cavitation resistance induction wheel

根據(jù)已知設(shè)計(jì)參數(shù)，選取水力模型的葉輪NPSHr= 8 m，并不能滿(mǎn)足泵汽蝕余量的設(shè)計(jì)要求，且相差甚遠(yuǎn)；通過(guò)葉輪汽蝕性能優(yōu)化，亦無(wú)法滿(mǎn)足，因此，采用前置誘導(dǎo)輪的設(shè)計(jì)方法，期望滿(mǎn)足設(shè)計(jì)要 求。

按照上文敘述的高抗汽蝕誘導(dǎo)輪設(shè)計(jì)思路及設(shè)計(jì)流程，通過(guò)高抗汽蝕誘導(dǎo)輪設(shè)計(jì)計(jì)算程序，確定誘導(dǎo)輪設(shè)計(jì)參數(shù)；應(yīng)用SolidWorks 對(duì)設(shè)計(jì)出的不同方案水力模型進(jìn)行三維建模。三維模型如圖4～5所 示。

2.2 網(wǎng)格劃分

采用ANSYS ICEM 網(wǎng)格劃分軟件，采用適應(yīng)性較強(qiáng)的非結(jié)構(gòu)化四面體網(wǎng)格類(lèi)型，對(duì)于葉片進(jìn)口、出口、圓角等重要位置進(jìn)行了局部加密，葉片表面采用棱柱層網(wǎng)格處理，誘導(dǎo)輪與葉輪計(jì)算域總網(wǎng)格數(shù)量約250 萬(wàn)。

當(dāng)模型網(wǎng)格數(shù)大于250 萬(wàn)時(shí)，網(wǎng)格數(shù)對(duì)性能影響不大，而且計(jì)算耗時(shí)比較合理，因此選擇該網(wǎng)格數(shù)進(jìn)行數(shù)值計(jì)算。

2.3 汽蝕模型和邊界條件

圖4 誘導(dǎo)輪三維模型Fig.4 3D model of induction wheel

圖5 誘導(dǎo)輪和首級(jí)葉輪三維模型Fig.5 3D models of the induction wheel and head impeller

空化模型選用ANSYS CFX 軟件提供的均相流（Homogeneous）模型，此模型認(rèn)為蒸汽速度與液體速度相同，采用Rayleigh-Plesset 方程描述空泡形成和潰滅時(shí)液相與氣相之間質(zhì)量傳遞的過(guò)程。誘導(dǎo)輪離心泵內(nèi)部流動(dòng)為三維不可壓黏性湍流流動(dòng)，控制方程為連續(xù)性方程和雷諾時(shí)均N-S 方程，湍流模型選用能夠較好地處理高應(yīng)變率及流線(xiàn)彎曲程度較大流動(dòng)的RNG k-ε 模型。采用高階求解差分格式，收斂精度為10-5。壁面滿(mǎn)足無(wú)滑移邊界條件，近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)法；以單相流動(dòng)數(shù)值計(jì)算的收斂結(jié)果作為空化數(shù)值計(jì)算的初始流場(chǎng)，通過(guò)不斷減小進(jìn)口總壓力值進(jìn)行定?？栈瘮?shù)值計(jì)算。

圖6 計(jì)算模型網(wǎng)格數(shù)-揚(yáng)程變化曲線(xiàn)Fig.6 Grid-head variation curve of the calculated model

圖7 計(jì)算模型網(wǎng)格示意Fig.7 Grid diagram of calculation model

2.4 數(shù)值模擬結(jié)果與評(píng)價(jià)

根據(jù)表5 汽蝕性能計(jì)算結(jié)果，認(rèn)為NPSHr= 20 m時(shí)，誘導(dǎo)輪與葉輪內(nèi)沒(méi)有發(fā)生汽蝕，此時(shí)的揚(yáng)程134.2 m 作為未發(fā)生汽蝕的衡量指標(biāo)，按API 610 標(biāo)準(zhǔn)要求，以揚(yáng)程下降3%時(shí)的NPSHr值為考核點(diǎn)。當(dāng)NPSHr= 0.9 m 時(shí)，若只按揚(yáng)程指標(biāo)來(lái)看，此時(shí)仍未發(fā)生汽蝕，需結(jié)合流場(chǎng)做進(jìn)一步分析，其曲線(xiàn)見(jiàn)圖8。

表5 額定工況汽蝕性能預(yù)測(cè)結(jié)果Table 5 Cavitation performance prediction results at rated conditions

圖8 汽蝕性能預(yù)測(cè)曲線(xiàn)Fig. 8 Cavitation performance prediction curve

圖9、圖10 流線(xiàn)圖中顯示的誘導(dǎo)輪內(nèi)部紅色高亮部分為氣體體積分?jǐn)?shù)ξ ≥0.5 所占用的體積空間。當(dāng)NPSHr= 1.2 m 時(shí)，氣泡被制約在誘導(dǎo)輪葉片1/3 位置上。當(dāng)NPSHr= 1 m 時(shí)，氣泡消失于誘導(dǎo)輪葉片2/3 位置上，并未對(duì)主流形成過(guò)多干擾。而當(dāng)NPSHr= 0.9 m 時(shí)，氣泡布滿(mǎn)整個(gè)誘導(dǎo)輪葉片占據(jù)大量空間，對(duì)主流形成明顯的堵塞作用，誘導(dǎo)輪為液體提供能量的作用減弱，大量氣泡聚集在葉輪葉片入口背面位置上，影響葉輪做功能力。

圖9 不同NPSHr值下氣相體積分?jǐn)?shù)分布圖Fig.9 Distribution of gas phase volume fraction under different NPSHrvalues

圖10 不同NPSHr值下的氣相分布圖Fig.10 Gas phase distribution under different NPSHrvalues

圖11～圖13 為不同NPSHr 值下，誘導(dǎo)輪與首級(jí)葉輪不同流面（Span=0.1、0.5、0.9）下的氣相分布?？梢钥闯觯S著入口總壓不斷降低，NPSHr不斷降低，在3 個(gè)流面上，氣相體積占比不斷增加的；而在這3 個(gè)流面中，NPSHr相同時(shí)，氣相分布規(guī)律明顯：從進(jìn)輪轂流面（Span=0.1）至近輪緣流面（Span=0.9），氣相占比逐漸增加，這也與上文分析相互對(duì)應(yīng)。

圖11 不同NPSHr下，誘導(dǎo)輪+葉輪近輪轂流面（Span=0.1）的氣相分布Fig.11 Gas phase distribution of the induction wheel and impeller near the hub flow surface (Span=0.1) under different NPSHr

圖12 不同NPSHr下，誘導(dǎo)輪+葉輪中間流面（Span=0.5）的氣相分布Fig.12 Gas phase distribution of inducer and impeller middle flow surface (Span=0.5) under different NPSHr

圖13 不同NPSHr下，誘導(dǎo)輪+葉輪近輪緣流面（Span=0.9）的氣相分布Fig.13 Gas phase distribution on the flow surface (Span=0.9) of the induction wheel and impeller at different NPSHr

3 試驗(yàn)驗(yàn)證

該LNG 高壓泵在大連深藍(lán)泵業(yè)有限公司的I 級(jí)試驗(yàn)臺(tái)以及LNG 介質(zhì)試驗(yàn)臺(tái)分別進(jìn)行性能試驗(yàn)。在試驗(yàn)過(guò)程中，泵頭插入深度L=1.8 m，以首級(jí)揚(yáng)程下降3%時(shí)的汽蝕值作為臨界汽蝕值。測(cè)試結(jié)果如下表6 和圖14 所示。需要說(shuō)明的是，入口取壓點(diǎn)位于液面以下，測(cè)試的NPSHr值是以試驗(yàn)泵安裝平面為基準(zhǔn)的，需要將NPSHr值換算值誘導(dǎo)輪入口中心處，因此：

表6 額定工況汽蝕性能測(cè)試結(jié)果Table 6 Test results of cavitation performance at rated conditions

圖14 額定工況水介質(zhì)下，NPSHr性能測(cè)試曲線(xiàn)Fig.14 NPSHrperformance test curve under rated working condition

NPSHr實(shí)際= NPSHr測(cè)試+ L = -0.95 + 1.8 = 0.85 m（5）式中 L——泵頭插入深度，m。

在隨后的LNG 介質(zhì)下，NPSHr實(shí)際= 0.5 m，很好地實(shí)現(xiàn)了設(shè)計(jì)預(yù)期。

對(duì)比NPSHr預(yù)測(cè)結(jié)果與最終結(jié)果不難發(fā)現(xiàn)，按 照高抗汽蝕誘導(dǎo)輪設(shè)計(jì)思路及結(jié)構(gòu)選取原則，借助CFD 分析手段，使得該高抗汽蝕誘導(dǎo)輪設(shè)計(jì)開(kāi)發(fā)非常成功，設(shè)計(jì)吻合性?xún)?yōu)異。同時(shí)，也可以看出，在離心泵的設(shè)計(jì)領(lǐng)域內(nèi)，無(wú)論是重大項(xiàng)目研發(fā)，還是日常的水力性能優(yōu)化中，CFD 正扮演著越來(lái)越重要的角 色。

4 結(jié)論

本文提煉出若干離心泵高抗汽蝕誘導(dǎo)輪的設(shè)計(jì)準(zhǔn)則，并以某LNG 高壓泵國(guó)產(chǎn)化首臺(tái)泵項(xiàng)目為依托，采用CFD 和試驗(yàn)相結(jié)合的方法，很好地實(shí)現(xiàn)了高抗汽蝕誘導(dǎo)輪的水力開(kāi)發(fā)及驗(yàn)證，并從中得到以下結(jié) 論：

（1）借助CFD 分析手段，在高抗汽蝕誘導(dǎo)輪的水力設(shè)計(jì)中，可以很好地與理論設(shè)計(jì)相結(jié)合，提高了水力設(shè)計(jì)準(zhǔn)確率，節(jié)約經(jīng)濟(jì)成本和縮短設(shè)計(jì)周期。

（2）從實(shí)踐角度出發(fā)，本文在一定范圍內(nèi)總結(jié)了離心泵高抗汽蝕誘導(dǎo)輪的水力設(shè)計(jì)思想及形式選取原則，可對(duì)相關(guān)從業(yè)人員在設(shè)計(jì)時(shí)進(jìn)行參考。

（3）當(dāng)NPSHr值相同時(shí)，氣體相從近輪轂流面到近輪緣流面逐漸增加，主要?dú)馀葑罱K聚集在輪緣處。

（4）在提高離心泵汽蝕性能的設(shè)計(jì)措施中，前置誘導(dǎo)輪是一種行之有效的設(shè)計(jì)方法，良好的設(shè)計(jì)匹配，能夠很好地提升離心泵的汽蝕性能。