嚴(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)噪聲,曲線形狀等)均滿足要求,同時(shí)形成了誘導(dǎo)輪的水力設(shè)計(jì)規(guī)范及水力設(shè)計(jì)計(jì)算程序。本文以某LNG 國(guó)產(chǎn)化首臺(tái)套泵研發(fā)項(xiàng)目為依托開(kāi)展高抗汽蝕誘導(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 所示。
根據(jù)誘導(dǎo)輪汽蝕比轉(zhuǎn)速的限定要求,在設(shè)計(jì)過(guò)程中推薦使用的誘導(dǎo)輪結(jié)構(gòu)形式如表1 所示。
在滿足性能要求,節(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)輪
基于額定流量設(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)輪的汽蝕特性曲線:
圖3 誘導(dǎo)輪汽蝕特性曲線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)行分析研究,在此不再贅述。
按照上述設(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,并不能滿足泵汽蝕余量的設(shè)計(jì)要求,且相差甚遠(yuǎn);通過(guò)葉輪汽蝕性能優(yōu)化,亦無(wú)法滿足,因此,采用前置誘導(dǎo)輪的設(shè)計(jì)方法,期望滿足設(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所 示。
采用ANSYS ICEM 網(wǎng)格劃分軟件,采用適應(yīng)性較強(qiáng)的非結(jié)構(gòu)化四面體網(wǎng)格類型,對(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ì)算。
圖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)變率及流線彎曲程度較大流動(dòng)的RNG k-ε 模型。采用高階求解差分格式,收斂精度為10-5。壁面滿足無(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)程變化曲線Fig.6 Grid-head variation curve of the calculated model
圖7 計(jì)算模型網(wǎng)格示意Fig.7 Grid diagram of calculation model
根據(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)一步分析,其曲線見(jiàn)圖8。
表5 額定工況汽蝕性能預(yù)測(cè)結(jié)果Table 5 Cavitation performance prediction results at rated conditions
圖8 汽蝕性能預(yù)測(cè)曲線Fig. 8 Cavitation performance prediction curve
圖9、圖10 流線圖中顯示的誘導(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í),氣泡布滿整個(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)下的氣相分布。可以看出,隨著入口總壓不斷降低,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
該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è)試曲線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ì)吻合性優(yōu)異。同時(shí),也可以看出,在離心泵的設(shè)計(jì)領(lǐng)域內(nèi),無(wú)論是重大項(xiàng)目研發(fā),還是日常的水力性能優(yōu)化中,CFD 正扮演著越來(lái)越重要的角 色。
本文提煉出若干離心泵高抗汽蝕誘導(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ì)匹配,能夠很好地提升離心泵的汽蝕性能。