石保虎，劉金龍，李晨昊，張高正，王李科

（1.國家石油天然氣管網(wǎng)集團(tuán)有限公司華南分公司，廣東省廣州市 510623；2.西安航天泵業(yè)有限公司，陜西省西安市 710077；3.西安理工大學(xué)，陜西省西安市 710048）

0 引言

雙吸泵以其流量大、軸向受力小和檢修便捷等優(yōu)點(diǎn)，被廣泛用于工農(nóng)業(yè)生產(chǎn)中。但由于實(shí)際生產(chǎn)中泵往往會偏離設(shè)計(jì)工況點(diǎn)在小流量下運(yùn)行，不利的入流條件會導(dǎo)致泵的性能下降、并改變其運(yùn)行穩(wěn)定性。因此探究小流量工況下泵的性能變化，以及工況對運(yùn)行穩(wěn)定性的影響對指導(dǎo)生產(chǎn)活動以及雙吸泵設(shè)計(jì)具有重要研究意義。

數(shù)值模擬作為目前重要的研究手段，已經(jīng)被廣泛應(yīng)用于離心泵的性能預(yù)測及內(nèi)部流動分析。劉健瑞等[1]基于商用軟件CFX 中提供RNGk-ε模型研究了交錯(cuò)葉片對雙吸泵性能的影響，并指出交錯(cuò)葉片角度為30°時(shí)的效果最好。黃凱杰等[2]也通過數(shù)值模擬的方法對一雙吸泵的性能以及內(nèi)部流動進(jìn)行了分析研究，并通過分析結(jié)果對結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)，使得泵效率得到了有效提升。朱榮生等[3]則同樣采用RNGk-ε分析了重心偏移對雙吸泵葉輪受力特性的影響，并通過與試驗(yàn)數(shù)據(jù)的對比驗(yàn)證了數(shù)值模擬的準(zhǔn)確性。馮憲等[4]則采用數(shù)值模擬結(jié)合試驗(yàn)的方法對一雙吸泵的空化特性進(jìn)行了分析，不僅驗(yàn)證了數(shù)值模擬的可靠性，并且基于模擬結(jié)果提出了相關(guān)優(yōu)化建議。熵產(chǎn)理論作為熱力學(xué)原理的一個(gè)延伸，目前已經(jīng)在水力機(jī)械能量損失分析方面得到了廣泛應(yīng)用[5，6]。苗森春等[7]采用熵產(chǎn)理論探究了雙吸泵做透平過程的能量變化機(jī)理，并指出湍流熵產(chǎn)和壁面熵產(chǎn)基本占據(jù)了總損失的95%以上，而葉輪和尾水室則是損失產(chǎn)生的主要部件。李德友[8]也采用熵產(chǎn)理論對水泵水輪機(jī)泵工況駝峰區(qū)的能量變化進(jìn)行了研究。任蕓等[9]對離心泵內(nèi)的流動損失進(jìn)行了分析，并將引起葉輪損失的類型劃分為由于動靜干涉引起的損失和水流沖擊引起的損失。李東闊等[10]也采用熵產(chǎn)理論對水泵水輪機(jī)泵工況的能量變化進(jìn)行了討論，指出葉輪和導(dǎo)葉內(nèi)的水力損失是泵損失的主要來源。

綜合來看，許多研究者已經(jīng)采用數(shù)值模擬結(jié)合熵產(chǎn)的方法在泵的能量特性分析和受力預(yù)測等方面進(jìn)行了應(yīng)用，并取得了不錯(cuò)的結(jié)果。因此，本文同樣采用數(shù)值模擬的方法，基于熵產(chǎn)理論探究小流量工況下雙吸泵的特性能量特性及其運(yùn)行穩(wěn)定性。研究結(jié)果可為雙吸式離心泵優(yōu)化設(shè)計(jì)提供理論依據(jù)。

1 計(jì)算模型與數(shù)值方法

1.1 控制方程

式中：ρ——流體密度，kg/m3；

ui——i方向流體的流速；

p——流體靜壓，Pa；

μeff——有效黏滯性系數(shù)，μeff=μ+μt，μt為湍流黏度；

Si——源項(xiàng)。

1.2 幾何模型

本文以一單級雙吸離心式輸油泵為研究對象，該泵額定轉(zhuǎn)速n=1480r/min，額定流量Qd=650m3/h，額定揚(yáng)程Hd=82m。主要過流部件由進(jìn)水管、葉輪和蝸殼三部分組成，具體結(jié)構(gòu)如圖1所示。

圖1 雙吸離心泵幾何模型Figure 1 Geometric model of double suction centrifugal pump

1.3 網(wǎng)格劃分和邊界條件設(shè)置

由于雙吸泵幾何結(jié)構(gòu)復(fù)雜，因此參考文獻(xiàn)[7]的研究處理方法，對結(jié)構(gòu)簡單的進(jìn)口管采用商用軟件ICEM 進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分，而對于結(jié)構(gòu)復(fù)雜的蝸殼和葉輪則采用非結(jié)構(gòu)化網(wǎng)格。具體網(wǎng)格如圖2所示。同時(shí)為了節(jié)約計(jì)算資源并保證計(jì)算精度，采用4 組不同數(shù)量的網(wǎng)格模型對泵的設(shè)計(jì)工況進(jìn)行了試算結(jié)果如圖3所示，發(fā)現(xiàn)當(dāng)網(wǎng)格數(shù)量大于330 萬后，計(jì)算獲得的泵揚(yáng)程值變化幅度小于2%，因此最終選擇網(wǎng)格數(shù)量為334 萬進(jìn)行后續(xù)研究。采用CFX 進(jìn)行模擬計(jì)算，忽略水的熱交換和可壓縮性，進(jìn)口邊界條件設(shè)為總壓進(jìn)口，出口邊界條件設(shè)為質(zhì)量流量出口。湍流模型選用RNGk-ε模型。動靜交界面間質(zhì)量傳遞采用凍結(jié)轉(zhuǎn)子技術(shù)，壁面設(shè)為絕熱無滑移壁面，并用標(biāo)準(zhǔn)壁面函數(shù)法對近壁區(qū)域進(jìn)行處理。采用二階離散格式對計(jì)算域進(jìn)行離散，收斂精度為平均殘差小于10-5。

圖2 網(wǎng)格結(jié)構(gòu)圖Figure 2 Mesh grid

圖3 網(wǎng)格無關(guān)性驗(yàn)證Figure 3 Mesh independence

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

試驗(yàn)在大型工業(yè)泵水力試驗(yàn)系統(tǒng)上完成，系統(tǒng)主要由蓄水池、前置泵、試驗(yàn)臺、閥門和各種儀器儀表及控制系統(tǒng)組成。圖4 為試驗(yàn)和數(shù)值模擬獲得的外特性曲線對比圖，從圖中可以看出在額定工況附近數(shù)值模擬獲得的外特性與試驗(yàn)值基本一致，雖然偏離額定工況后誤差略有增大，但整個(gè)流量范圍內(nèi)其揚(yáng)程相對誤差不超過3%，效率誤差不超8%，證明本文采用的仿真方法是可靠的。

圖4 試驗(yàn)與數(shù)值模擬外特性對比Figure 4 Comparison of external characteristics between experiment and numerical simulation

2 模擬結(jié)果分析

2.1 能量損失分析

由熱力學(xué)第二定律可以知道，不可逆熱力過程中熵的微增量總是大于零的，因此在旋轉(zhuǎn)機(jī)械做功的過程中也伴隨有熵增現(xiàn)象。文獻(xiàn)[8]將旋轉(zhuǎn)機(jī)械中的熵增現(xiàn)象分為兩種：一種是由于流體脈動引起的黏性耗散，更為具體的又可以分為由于時(shí)均運(yùn)動引起的耗散和由于速度脈動引起的湍流耗散，另一種則是由于壁面摩擦引起的內(nèi)能耗散。熵產(chǎn)計(jì)算公式見文獻(xiàn)[8]。圖5 為不同流量下雙吸泵各部件的熵產(chǎn)損失分布，由于平均速度引起的耗散熵產(chǎn)很小，因此本文將其歸于脈動熵產(chǎn)。由圖可知葉輪和蝸殼內(nèi)的損失是泵損失的主要來源，整個(gè)流量范圍內(nèi)葉輪內(nèi)損失占比在20%～43%之間，而蝸殼的損失占比在54%～73%之間，進(jìn)水管的損失占比最大不超過7%。葉輪內(nèi)脈動損失隨流量增大逐漸減小，這表明越靠近設(shè)計(jì)工況，葉輪內(nèi)流動狀態(tài)越好，而壁面熵產(chǎn)也隨著流量的增大略有減小，但變化不劇烈。蝸殼內(nèi)整體損失變化與葉輪基本相同，但不同的是葉輪內(nèi)脈動損失占比更高其在0.4Qd～0.8Qd之間，脈動損失基本保持在64%左右，這表明葉輪內(nèi)損失的主要來源是不穩(wěn)定的流動形態(tài)。而蝸殼內(nèi)壁面損失的占比則更高，整個(gè)流量范圍內(nèi)其壁面損失占比保持在53%～57%之間。

圖5 各部件熵產(chǎn)損失分布Figure 5 Entropy loss distribution of each component under different flow rate

2.2 葉輪內(nèi)流場分析

從前文可以看出，葉輪內(nèi)的流動損失主要來源于流動的脈動，因此對0.4Qd、0.6Qd和0.8Qd流量下葉輪內(nèi)的熵產(chǎn)分布和流場分布進(jìn)一步分析。由于雙吸泵為對稱結(jié)構(gòu)，因此選取一邊葉輪作為研究對象進(jìn)行分析，圖6 為葉輪0.5 倍葉高處截面上的熵產(chǎn)率分布（定義后蓋板到前蓋板無量綱距離為0～1）。從圖中可以看出高損失區(qū)域主要集中在轉(zhuǎn)輪出口和葉片壓力面附近，尤其是出口靠近隔舌和隔板尾端的區(qū)域。對比三種流量下的熵產(chǎn)率分布發(fā)現(xiàn)，隨著流量增大，整體流道內(nèi)的高熵產(chǎn)區(qū)域逐漸減小。但由于靠近隔舌區(qū)域，流道4和流道5 的高熵產(chǎn)區(qū)域減小程度要明顯小于其他流道。這表明隨著流量逐漸靠近設(shè)計(jì)點(diǎn)葉輪內(nèi)的流動狀態(tài)雖然得到了改善，但受隔舌的影響，流道5 出口的流動狀態(tài)仍較為復(fù)雜，且這種復(fù)雜流動會逆著旋轉(zhuǎn)方向傳播至流道4。

圖6 葉輪中間截面脈動熵產(chǎn)損失分布Figure 6 Distribution of entropy generation loss in the middle section of impeller

圖7 為相應(yīng)截面上的流場分布，對比熵產(chǎn)分布圖可以發(fā)現(xiàn)，雖然葉輪流道內(nèi)出現(xiàn)了大面積渦流，但渦流結(jié)構(gòu)相對穩(wěn)定、且強(qiáng)度較低，因此并未誘發(fā)大的速度脈動產(chǎn)生高損失。而高損失區(qū)域的流動主要表現(xiàn)為由于動靜干涉引起的不穩(wěn)定流動和葉片表面流動分離流動。具體分析各流道內(nèi)的流動變化可以看出流道1 由于受隔板頭部的影響，在0.4Qd時(shí)出口壓力面處形成了顯著渦流，但隨著流量的增大該渦流逐漸與葉片壓力面的流動分離區(qū)域融合，在葉片背面形成較為穩(wěn)定的分離流動。而且從整個(gè)流量范圍來看，流道1 出口的不穩(wěn)定流動會影響流道2 內(nèi)的流動，從圖中可以看出，當(dāng)流道1 出口區(qū)域的渦流消失后，流道2 內(nèi)的大面積渦流也消失了，且出口處的高速區(qū)域也得到了改善。而流道4 與流道5 由于受蝸殼隔舌的影響，可以發(fā)現(xiàn)其流道內(nèi)的流動整體比較復(fù)雜，但變化趨勢與1、2 一致，只是在0.8Qd下其出口處的流動仍不穩(wěn)定。從1、2、4 和5 四個(gè)流道的流動變化可以看出，流量對動靜干涉的影響極為顯著，且蝸殼隔舌引起的動靜干涉明顯強(qiáng)于隔板頭部的影響，這種差異可能要是由于隔舌與隔板頭部與葉輪出口的距離不同引起的。

圖7 葉輪中間截面流場Figure 7 Middle section flow field of impeller

2.3 蝸殼內(nèi)流場分析

蝸殼中間橫截面的熵產(chǎn)分布如圖8所示，可以看出蝸殼內(nèi)高脈動熵產(chǎn)區(qū)域主要集中在靠近葉輪出口的區(qū)域、隔板頭部以及蝸殼隔舌附近。隨著流量的增大，高熵產(chǎn)區(qū)域逐漸縮小，且數(shù)值也逐漸降低，但對于0.8Qd可以發(fā)現(xiàn)在蝸殼外壁面和靠近葉輪出口的隔板壁面區(qū)域出現(xiàn)了顯著的高熵產(chǎn)區(qū)，這可能是由于流量增大后，蝸殼過流面積不足，慣性力和離心力迫使流體向壁面區(qū)域集中而形成的擾動所導(dǎo)致的。

圖8 蝸殼內(nèi)脈動熵產(chǎn)損失分布Figure 8 Distribution of fluctuating entropy production loss in volute

圖9 為相應(yīng)的流場分布，從圖中可以看出隨著流量的增大，由于隔板右側(cè)的低速區(qū)域逐漸向隔板壁面靠近，流場得到顯著改善，這使得蝸殼內(nèi)的出流更加順暢，可以發(fā)現(xiàn)此時(shí)隔板左側(cè)的流動相較小流量工況更加平穩(wěn)，而對比熵產(chǎn)圖也可以發(fā)現(xiàn)此時(shí)相應(yīng)位置的熵產(chǎn)明顯減弱。而對比三種流量下隔舌位置的流動可以發(fā)現(xiàn)，隨著流量增大隔舌位置的渦流面積逐漸減小，且逐漸靠近隔舌。但即便是減小的渦流也仍會壓縮過流通道影響蝸殼的出流，使得隔舌上游區(qū)域出現(xiàn)高速流動，產(chǎn)生較為明顯的損失。

圖9 蝸殼中間截面流場Figure 9 Middle section flow field of volute

為了進(jìn)一步探究蝸殼內(nèi)高壁面熵產(chǎn)的原因，圖10 給出了蝸殼內(nèi)壁面熵產(chǎn)率的分布，可以看出不同流量下蝸殼壁面熵產(chǎn)的分布差異不大，高熵產(chǎn)率區(qū)域主要是位于蝸殼漸擴(kuò)段和隔板右側(cè)壁面附近。由于隔板的存在壓縮了蝸殼的流道，使得隔板右側(cè)區(qū)域蝸殼流道內(nèi)部出現(xiàn)了明顯的高速流動，進(jìn)而導(dǎo)致在蝸殼壁面以及隔板右側(cè)壁面上均出現(xiàn)了顯著的高熵產(chǎn)。因此綜合來看，隔板在小流量工況下不僅引起了頭部與葉輪出口的不穩(wěn)定流動，加劇了脈動損失，同時(shí)在整個(gè)流量范圍內(nèi)也壓縮了過流通道，使得在蝸殼內(nèi)產(chǎn)生了較高的壁面損失，因此應(yīng)對其結(jié)構(gòu)進(jìn)行優(yōu)化，以增加其在整個(gè)工作范圍內(nèi)的適用性。

圖10 蝸殼內(nèi)壁面熵產(chǎn)分布Figure 10 Distribution of entropy generation rate caused by wall shear stress of volute

2.4 葉輪受力特性分析

雙吸泵由于軸向結(jié)構(gòu)對稱，所以軸向受力很小，因此影響雙吸泵穩(wěn)定性的主要是其徑向力。而且從前文也可以看出，小流量工況下葉輪內(nèi)的流動分布并不對稱，因此極有可能在葉輪上產(chǎn)生較大的徑向力，威脅泵的運(yùn)行穩(wěn)定性。數(shù)值模擬中，葉輪所受徑向力由其x軸方向和y軸方向的合力計(jì)算得到，其計(jì)算如式（3）所示。由于通過公式計(jì)算的徑向力在小流量工況下遠(yuǎn)大于數(shù)值模擬值，但兩種方法獲得的徑向力變化規(guī)律較為一致。本文主要討論流量對葉輪徑向力的影響，為了方便探究流量對徑向力變化的影響，本文定義無量綱徑向力系數(shù)Δf［見式（4）］[11]來探究流量對葉輪徑向力的影響，結(jié)果如圖11所示。從圖11 可以看出隨著流量的增大，葉輪所受徑向力基本呈先減小后增大的一個(gè)趨勢，其最小值在設(shè)計(jì)點(diǎn)附近。產(chǎn)生這種現(xiàn)象的原因主要是由于偏工況下葉輪內(nèi)流動分布不對稱，使得作用葉輪在葉輪上的壓力在某一方向集中，故而產(chǎn)生了較大的徑向力。而越靠近額定工況，其各流道內(nèi)分布越一致，所以作用相互抵消，故徑向力較小。

圖11 葉輪徑向力隨流量變化曲線Figure 11 Variation of Δf with flow rate

式中：F——徑向力大?。?/p>

Fx——徑向力在x軸方向的分量；

Fy——徑向力在y軸方向的分量。

式中：Fd——設(shè)計(jì)工況點(diǎn)葉輪所受徑向力。

3 結(jié)論

本文以一雙吸式離心泵為研究對象，采用數(shù)值模擬的方法，對不同流量下泵內(nèi)的損失分布、流動特性和葉輪的受力特性進(jìn)行了研究，主要研究結(jié)論如下：

（1）小流量工況下泵的損失主要來源于葉輪和蝸殼的脈動損失和壁面損失，兩部件內(nèi)的損失占據(jù)了泵總損失的95%以上；葉輪內(nèi)脈動損失占比更高基本保持在64%左右，而蝸殼內(nèi)壁面損失的占比更高，基本保持在55%左右。

（2）葉輪內(nèi)脈動損失占比更高，且高損失區(qū)域主要位于葉輪出口靠近隔板頭部和蝸殼隔舌的區(qū)域，表明由于動靜干涉引起的不良流動是損失產(chǎn)生的主要來源；蝸殼內(nèi)壁面損失占比更高，且在隔板附近出現(xiàn)明顯的高熵產(chǎn)區(qū)，說明小流量工況下隔板的存在會加劇損失，帶來不利影響。

（3）不同流量下葉輪所受徑向力呈先減小后增加的趨勢，額定工況下葉輪所受軸向力最小。