基于非均相流模型的離心泵氣液兩相流動(dòng)數(shù)值研究

袁建平 張克玉 司喬瑞 周幫倫 唐苑峰 金中坤

(1.江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心， 鎮(zhèn)江 212013； 2.中國科學(xué)院聲學(xué)研究所， 北京 100190)

基于非均相流模型的離心泵氣液兩相流動(dòng)數(shù)值研究

袁建平1張克玉1司喬瑞1周幫倫1唐苑峰1金中坤2

(1.江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心， 鎮(zhèn)江 212013； 2.中國科學(xué)院聲學(xué)研究所， 北京 100190)

為研究離心泵氣液混輸狀態(tài)下的內(nèi)部流動(dòng)特性，基于Eulerian-Eulerian非均相流模型以空氣和水作為工作介質(zhì)，在不同進(jìn)口含氣率工況下對離心泵內(nèi)流場進(jìn)行定常和非定常數(shù)值計(jì)算，相間阻力作用采用Schiller Nauman模型，得到葉輪內(nèi)氣相分布情況以及氣液兩相的速度流線圖，探求氣液兩相在泵內(nèi)的流動(dòng)規(guī)律。計(jì)算結(jié)果表明：在葉片進(jìn)口邊壓力面氣相濃度較高處，會(huì)產(chǎn)生漩渦，說明葉輪流道內(nèi)漩渦的產(chǎn)生與氣體的聚集有很大的關(guān)系；增大進(jìn)口含氣率到10%時(shí)，葉輪流道內(nèi)靠近吸力面處已經(jīng)出現(xiàn)比較明顯的相態(tài)分離現(xiàn)象，氣相有沿著葉片吸力面向葉輪出口運(yùn)動(dòng)的趨勢；氣相在葉輪流道內(nèi)會(huì)沿流道中部向前后蓋板運(yùn)動(dòng)，隨著含氣率增大，靠近前后蓋板側(cè)的氣相濃度逐漸增大，靠近葉輪出口邊前蓋板側(cè)的氣相濃度增加較后蓋板更明顯，最終氣體可能會(huì)堵塞流道；在一個(gè)旋轉(zhuǎn)周期內(nèi)，葉輪出口壓力呈周期性變化，進(jìn)口含氣率從1%增至10%，葉輪出口壓力逐漸降低，監(jiān)測點(diǎn)在不同含氣率下壓力脈動(dòng)主頻在葉片通過頻率附近，進(jìn)口含氣率不超過10%時(shí)對監(jiān)測點(diǎn)壓力脈動(dòng)的主頻及次主頻影響不大；通過對比試驗(yàn)結(jié)果和數(shù)值計(jì)算結(jié)果，證明所采用的計(jì)算模型和方法基本可靠。

離心泵； 氣液混輸； Eulerian-Eulerian非均相流模型； 壓力脈動(dòng)； 數(shù)值模擬

引言

離心泵是重要的能量轉(zhuǎn)換裝置和流體輸送設(shè)備，在工程實(shí)際如化工流程、石油天然氣輸送、核泵失水事故等[1-2]情況下經(jīng)常遇到離心泵氣液混輸?shù)膯栴}。由于離心泵對工作介質(zhì)中的氣體含量十分敏感，氣液混輸條件下泵的運(yùn)行性能會(huì)發(fā)生較大改變，離心泵在進(jìn)口含氣率達(dá)到10%～15%時(shí)就容易產(chǎn)生氣液分離，導(dǎo)致?lián)P程陡降[3]。因此研究離心泵氣液混輸狀態(tài)下的內(nèi)部流動(dòng)特性對提高系統(tǒng)運(yùn)行性能及安全可靠性、完善泵設(shè)計(jì)理論具有重要意義[4]。

近年來由于計(jì)算流體力學(xué)的迅速發(fā)展，國內(nèi)外研究人員開始使用數(shù)值模擬手段對泵內(nèi)部氣液兩相流動(dòng)進(jìn)行研究[5-7]。CARIDAD等[8-10]采用雙流體模型分析了葉輪內(nèi)兩相流動(dòng)規(guī)律，卻忽略了轉(zhuǎn)子與定子之間相互作用的影響。BARRIOS等[11]采用Eulerian-Eulerian非均相流模型分析了電動(dòng)潛水泵葉輪內(nèi)氣液兩相流動(dòng)規(guī)律。潘兵輝等[12]采用Mixture模型研究了氣液兩相流對離心泵揚(yáng)程和效率的影響。余志毅等[13]基于細(xì)泡狀流動(dòng)假設(shè)，采用雙流體模型對混輸泵葉輪內(nèi)的氣液兩相流動(dòng)進(jìn)行分析，發(fā)現(xiàn)氣相旋渦是造成氣體局部聚集的主要因素。但在關(guān)于氣體對泵內(nèi)部壓力的瞬態(tài)影響規(guī)律方面研究卻還不夠。

本文基于Eulerian-Eulerian非均相流模型，在不同流量、不同進(jìn)口含氣率條件下，對模型泵進(jìn)行定常和非定常數(shù)值計(jì)算，得到葉輪內(nèi)氣相分布情況及氣液兩相的速度場，并得到含氣率對泵出口壓力的瞬態(tài)影響規(guī)律，探求氣液兩相在泵內(nèi)的流動(dòng)規(guī)律，最后通過試驗(yàn)驗(yàn)證計(jì)算結(jié)果的準(zhǔn)確性。

1 數(shù)值方法

1.1 Eulerian-Eulerian非均相流模型

Eulerian-Eulerian多相流模型有2種：均相流模型和非均相流模型。前者未考慮速度滑移，假設(shè)各相速度相同；后者不僅考慮速度滑移，也考慮相間質(zhì)量及動(dòng)量傳遞等[14]，非均相流模型中每相流體都有各自的流場并且通過相間傳遞單元進(jìn)行傳遞，即每相都有各自的速度場和溫度場，最后通過相間作用力和熱量傳遞使兩相速度和溫度得到平衡[15]，因此，更接近實(shí)際情況。本文采用非均相流模型，不考慮溫度場，液相為連續(xù)相，氣相為離散相，相間傳遞單元采用Particle模型，假定氣液兩相流流型為泡狀流，滿足質(zhì)量和動(dòng)量守恒。

1.2 控制方程

連續(xù)性方程和動(dòng)量方程[16]為

式中k——任意相(l代表液體，g代表氣體)ρk——k相密度pk——k相壓力αk——k相體積分?jǐn)?shù)μk——k相動(dòng)力粘度wk——k相流體相對速度Mk——k相所受相間作用力fk——與葉輪旋轉(zhuǎn)有關(guān)的質(zhì)量力

氣體和液體滿足

式中αg——進(jìn)口含氣率αl——進(jìn)口含液率Qg——?dú)怏w體積流量Ql——液體體積流量

液相采用SST湍流模型，該模型考慮了湍流剪切應(yīng)力的輸運(yùn)，可對逆壓梯度下的流動(dòng)分離現(xiàn)象給出更精確的預(yù)測[16]；氣相采用零方程理論模型。

1.3 相間作用力

兩相間的作用力只考慮了阻力作用，作用于液相的阻力計(jì)算公式為

其中

式中dB——?dú)馀葜睆絚D——阻力系數(shù)Re——?dú)馀堇字Z數(shù)

2 幾何模型

模型泵選用NKG 65-50-125/139型直聯(lián)式單級單吸離心泵，其設(shè)計(jì)工況參數(shù)及基本幾何尺寸如下：設(shè)計(jì)流量Qd=50.6 m3/h，設(shè)計(jì)揚(yáng)程H=20.2 m，額定轉(zhuǎn)速n=2 910 r/min，比轉(zhuǎn)數(shù)ns=132.2；葉輪進(jìn)口直徑D1=79 mm，葉輪出口直徑D2=140 mm，葉片出口寬度b2=15.5 mm，蝸殼基圓直徑D3=149 mm，泵出口直徑Dd=50 mm，葉片數(shù)Z=6。

采用商業(yè)軟件Pro/E對模型泵水體進(jìn)行三維造型，為了使流體能充分發(fā)展并避免進(jìn)出口回流的影響，對進(jìn)出口進(jìn)行了適當(dāng)?shù)难娱L，且為使蝸殼和泵出口交界處銜接更好，造型過程中將蝸殼水體與泵出口段水體造成一整體。計(jì)算域三維建模如圖1所示，分為泵進(jìn)口水體、葉輪水體、蝸殼水體。

圖1 模型泵三維幾何模型Fig.1 Three dimensional geometry model1.進(jìn)口水體 2.葉輪水體 3.蝸殼水體

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

3.1 網(wǎng)格劃分

采用ICEM對水體域進(jìn)行六面體結(jié)構(gòu)化網(wǎng)格劃分，結(jié)構(gòu)化網(wǎng)格的優(yōu)勢之一是可以靈活控制各節(jié)點(diǎn)的密度對邊界層進(jìn)行調(diào)整[17]。為了確保數(shù)值模擬的精度，對近壁區(qū)的網(wǎng)格進(jìn)行了加密。根據(jù)圖2網(wǎng)格無關(guān)性檢查結(jié)果可知，當(dāng)網(wǎng)格數(shù)大于1 433 589時(shí)，泵揚(yáng)程計(jì)算結(jié)果變化幅值小于1%。綜合網(wǎng)格數(shù)及計(jì)算時(shí)間，最終選擇網(wǎng)格數(shù)1 433 589、網(wǎng)格節(jié)點(diǎn)數(shù)1 496 732。葉輪和蝸殼網(wǎng)格如圖3所示，葉片表面Y+值(第1層網(wǎng)格質(zhì)心到壁面無量綱距離)如圖4所示。

圖2 網(wǎng)格無關(guān)性分析Fig.2 Analysis of grid independence

圖3 葉輪和蝸殼網(wǎng)格Fig.3 Meshes of impeller and volute

圖4 葉片表面Y+值Fig.4 Y+ value on blade surface

3.2 邊界條件設(shè)置

假設(shè)進(jìn)口氣液兩相分布均勻，且無空化現(xiàn)象發(fā)生。在定常模擬時(shí)，進(jìn)口邊界條件給定質(zhì)量流量進(jìn)口，并給定進(jìn)口含氣率，出口根據(jù)試驗(yàn)值設(shè)定出口平均壓力，葉輪水體壁面設(shè)置為旋轉(zhuǎn)，其余壁面設(shè)置為靜止，葉輪和蝸殼之間交界面設(shè)置為Frozen Rotor模式，壁面為無滑移光滑絕熱壁面。在非定常計(jì)算時(shí)，修改動(dòng)靜部件耦合模型為Transient Rotor Stator模式，總計(jì)算時(shí)間為0.206 s，葉輪每轉(zhuǎn)3°為一個(gè)時(shí)間步長，即葉輪每時(shí)間步內(nèi)旋轉(zhuǎn)3°，總共保存了葉輪旋轉(zhuǎn)10周的數(shù)據(jù)，并采用最后穩(wěn)定的3周數(shù)據(jù)進(jìn)行分析。葉輪轉(zhuǎn)頻f1=2 910/60=48.5 Hz，葉片通過頻率f2=6f1=291 Hz。

3.3 監(jiān)測點(diǎn)設(shè)置

為研究葉輪出口壓力變化情況，在葉輪出口設(shè)置監(jiān)測點(diǎn)p1～p5，如圖5所示。

圖5 監(jiān)測點(diǎn)位置Fig.5 Positions of monitoring points

4 計(jì)算結(jié)果與分析

4.1 兩相含氣率分布

圖6和圖7給出了設(shè)計(jì)流量下不同進(jìn)口含氣率時(shí)葉輪中截面上的含氣率分布及速度流線分布情況。由圖6可知，在靠近葉片進(jìn)口壓力面，氣相濃度很高，對比圖7的相對速度流線，發(fā)現(xiàn)在氣相濃度較高處，會(huì)對應(yīng)產(chǎn)生漩渦，從而說明葉輪流道內(nèi)漩渦的產(chǎn)生與氣體的聚集有很大的關(guān)系。在保證液相流量不變，增大進(jìn)口氣相濃度時(shí)，從圖中可以看出，葉輪流道中的液相分布越來越稀疏，氣相在流道內(nèi)所占比例越來越大，氣相會(huì)更多地聚集在葉片吸力面以及出口邊附近，說明氣相有沿著葉片吸力面向葉輪出口運(yùn)動(dòng)的趨勢，與MüLLER等[18]得到的結(jié)果相吻合，這是因?yàn)橐合嘣谶M(jìn)入葉輪后，受到較大的慣性力和離心力作用偏離正常流線軌跡而向葉片壓力面運(yùn)動(dòng)，而氣相受到的慣性力和離心力較小，再加上受液相排擠作用從而被迫偏向吸力面，因而吸力面含氣率較高[19]。在進(jìn)口含氣率達(dá)到10%時(shí)，葉輪流道內(nèi)靠近吸力面處已經(jīng)出現(xiàn)比較明顯的相態(tài)分離現(xiàn)象(MINEMURA等[20]認(rèn)為含氣率大于20%時(shí)會(huì)發(fā)生氣泡聚集，產(chǎn)生相態(tài)分離)。

圖6 設(shè)計(jì)流量下不同進(jìn)口含氣率時(shí)氣相體積分?jǐn)?shù)分布云圖Fig.6 Gas volume fraction distribution contours at different inlet gas volume fractions under design flow rate

圖7 設(shè)計(jì)流量下不同進(jìn)口含氣率時(shí)液相相對速度流線分布云圖Fig.7 Water relative velocity streamline distribution contours at different inlet gas volume fractions under design flow rate

圖8 不同進(jìn)口含氣率時(shí)葉片壓力面和吸力面氣相體積分?jǐn)?shù)分布云圖Fig.8 Gas volume fraction distribution and relative velocity streamline contours on pressure side and suction side at different inlet gas volume fractions1.前蓋板側(cè) 2.后蓋板側(cè) 3.壓力面 4.吸力面

圖9 進(jìn)口含氣率為5%時(shí)葉片壓力面和吸力面氣相相對速度分布云圖Fig.9 Gas relative velocity streamline contour on pressure side and suction side at inlet gas volume fraction of 5%

4.2 葉片表面氣相分布及速度流線分布

圖8給出了設(shè)計(jì)流量下不同進(jìn)口含氣率時(shí)葉片壓力面和吸力面氣相體積分?jǐn)?shù)分布云圖和氣相相對速度流線圖。由圖8可知，在壓力面上，進(jìn)口含氣率較小時(shí)，壓力面前后蓋板側(cè)的氣相體積分?jǐn)?shù)較其它地方高，而隨著進(jìn)口含氣率的增加，發(fā)現(xiàn)前后蓋板側(cè)的氣相體積分?jǐn)?shù)逐漸增大。葉輪出口邊前蓋板側(cè)的氣相體積分?jǐn)?shù)增加更明顯，最終可能會(huì)堵塞流道。吸力面靠近葉片進(jìn)口處的氣相體積分?jǐn)?shù)也會(huì)隨著進(jìn)口含氣率的增大而增大。對比圖9的氣相速度流線，可以看出兩相流體進(jìn)入葉輪流道后，氣相會(huì)沿著壓力面從流道中間向前后蓋板運(yùn)動(dòng)，流線的運(yùn)動(dòng)趨勢與氣相的聚集位置相對應(yīng)，而在葉片吸力面氣相分布與速度流線較工作面要均勻。

4.3 進(jìn)口含氣率對葉輪出口壓力的影響

圖10給出了泵在設(shè)計(jì)流量下進(jìn)口含氣率為5%時(shí)葉輪出口監(jiān)測點(diǎn)p1～p5的壓力脈動(dòng)時(shí)域圖。由圖可知，各監(jiān)測點(diǎn)壓力均呈周期性變化，在相同的時(shí)間范圍內(nèi)出現(xiàn)6次波峰和波谷，6次壓力脈動(dòng)的主要原因是葉輪和隔舌之間動(dòng)靜干涉作用所致。從p1到p4處靜壓逐漸減小，p5處靜壓介于p1和p2之間，p1和p5處幅值變化大于其余各點(diǎn)，這表明隔舌附近區(qū)域流動(dòng)不穩(wěn)定性較強(qiáng)。

圖10 監(jiān)測點(diǎn)p1～p5的壓力脈動(dòng)時(shí)域圖Fig.10 Pressure fluctuation in time domain of monitoring points p1～p5

圖11和圖12給出了泵在額定流量不同泵進(jìn)口含氣率工況下葉輪出口監(jiān)測點(diǎn)p3的壓力脈動(dòng)時(shí)域圖和頻域圖。由圖可知，隨著進(jìn)口含氣率的增大，葉輪出口壓力逐漸降低。在進(jìn)口含氣率不超過10%時(shí)，壓力波動(dòng)規(guī)律仍然呈周期性變化，并沒有受到破壞，說明在此含氣率條件下氣體對葉輪出口壓力的影響作用并不明顯。壓力脈動(dòng)峰值集中在葉頻及其倍頻段，隨著頻率的增大能量值逐漸降低。比較不同進(jìn)口含氣率下所對應(yīng)壓力脈動(dòng)幅值，發(fā)現(xiàn)隨著進(jìn)口含氣率增大，在主頻處壓力脈動(dòng)幅值有稍微減小，依次為10.3、9.95、9.79 kPa。不同含氣率下監(jiān)測點(diǎn)的壓力脈動(dòng)主頻均為289 Hz，與葉片通過頻率接近，說明泵進(jìn)口含氣率不超過10%時(shí)氣體對壓力脈動(dòng)主頻影響不大，與付強(qiáng)等[14]得到的結(jié)論一致。

圖11 監(jiān)測點(diǎn)p3的壓力脈動(dòng)時(shí)域圖Fig.11 Pressure fluctuation in time domain of monitoring point p3

圖12 監(jiān)測點(diǎn)p3的壓力脈動(dòng)頻域圖Fig.12 Pressure fluctuation in frequency domain of monitoring point p3

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

為了驗(yàn)證數(shù)值計(jì)算的可靠性，進(jìn)行了不同進(jìn)口含氣率下的外特性試驗(yàn)，采用的試驗(yàn)臺為開式試驗(yàn)臺(如圖13所示)，試驗(yàn)裝置和儀器設(shè)備滿足GB/T 3216—2005中Ⅱ級精度要求。通過保證液體流量在某一恒定值，通過氣體調(diào)節(jié)閥來控制氣體流量來得到不同含氣率的氣液混合物，從而進(jìn)行試驗(yàn)泵的氣液混輸試驗(yàn)。儲水箱4中的水通過電磁流量計(jì)后在氣液混合裝置中與壓縮機(jī)產(chǎn)生的空氣混合均勻后進(jìn)入試驗(yàn)泵，泵輸出的氣液混合物經(jīng)過電動(dòng)閥后最終流入儲水箱5，儲水箱4和5上方均開口并且是連通的，這樣能使氣體有效排出并保證進(jìn)入儲水箱4的水不再含有空氣。

模型泵在純水和氣液兩相流工況下的揚(yáng)程以及進(jìn)口含氣率滿足公式

式中Hl、Hm——純水、氣液兩相流工況下?lián)P程g——重力加速度pout——泵出口壓力pin——泵進(jìn)口壓力

圖13 試驗(yàn)臺Fig.13 Test rig1.壓縮機(jī) 2、6電磁流量計(jì) 3.氣體流量計(jì) 4、5儲水箱 7.氣液混合裝置 8.試驗(yàn)泵

模型泵純水工況下的外特性曲線如圖14所示。從圖中可以看出數(shù)值計(jì)算與試驗(yàn)結(jié)果在整個(gè)流量范圍內(nèi)保持較好的一致性，在小流量工況點(diǎn)，兩者非常吻合，在大流量工況點(diǎn)，模擬與試驗(yàn)結(jié)果存在一定偏差。在設(shè)計(jì)點(diǎn)工況，離心泵揚(yáng)程試驗(yàn)值為19.8 m，揚(yáng)程模擬值為20.66 m；效率試驗(yàn)值為73.3%，效率模擬值為73.8%，設(shè)計(jì)點(diǎn)揚(yáng)程和效率相對誤差都在5%以內(nèi)。造成試驗(yàn)和模擬結(jié)果的偏差可能是未考慮壁面粗糙度以及忽略了葉輪進(jìn)口口環(huán)泄漏損失導(dǎo)致的。圖15為在進(jìn)口含氣率為5%時(shí)試驗(yàn)泵揚(yáng)程數(shù)值計(jì)算與試驗(yàn)結(jié)果對比。由圖可知，除了在較小流量下(Q=20.24 m3/h)模擬值與試驗(yàn)值偏差較大，其余各點(diǎn)模擬值與試驗(yàn)值均比較接近，在允許的相對誤差范圍內(nèi)，但在相同進(jìn)口含氣率條件下，小流量工況下的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果存在偏差，說明現(xiàn)有模型在小流量工況點(diǎn)的計(jì)算還不夠準(zhǔn)確，有待進(jìn)一步研究。

圖14 數(shù)值模擬和試驗(yàn)結(jié)果對比Fig.14 Comparisons between simulated and experimental results

圖15 進(jìn)口含氣率為5%時(shí)揚(yáng)程數(shù)值模擬和試驗(yàn)結(jié)果對比Fig.15 Comparisons between simulated and experimental results at inlet gas volume fraction of 5%

5 結(jié)論

(1)在進(jìn)口氣相濃度較高處，會(huì)產(chǎn)生漩渦，說明葉輪流道內(nèi)漩渦的產(chǎn)生與氣體的聚集有很大的關(guān)系。進(jìn)口含氣率增大到10%時(shí)，葉輪流道內(nèi)靠近吸力面處已經(jīng)出現(xiàn)比較明顯的相態(tài)分離，氣相有沿著葉片吸力面向葉輪出口運(yùn)動(dòng)的趨勢。

(2)氣液混輸狀態(tài)下，氣相在葉輪流道內(nèi)會(huì)沿流道中部向前后蓋板運(yùn)動(dòng)，隨著含氣率增大，前后蓋板側(cè)的氣體濃度也會(huì)增大。葉輪出口邊前蓋板側(cè)的氣相濃度增加更明顯，氣體最終可能會(huì)堵塞流道，造成泵性能惡化。

(3)在一個(gè)旋轉(zhuǎn)周期內(nèi)，葉輪出口壓力呈周期性變化。進(jìn)口含氣率從1%增至10%，葉輪出口壓力逐漸降低，監(jiān)測點(diǎn)在不同含氣率下壓力脈動(dòng)主頻在葉片通過頻率附近，進(jìn)口含氣率不超過10%時(shí)氣相體積分?jǐn)?shù)對監(jiān)測點(diǎn)壓力脈動(dòng)的主頻影響不大。

(4)通過試驗(yàn)結(jié)果和數(shù)值計(jì)算所獲得的揚(yáng)程特性曲線，驗(yàn)證了所采用的計(jì)算模型和方法是基本可靠的，但鑒于在小流量工況下的誤差，說明現(xiàn)階段旋轉(zhuǎn)機(jī)械內(nèi)氣液兩相流模型和計(jì)算方法仍不夠成熟，需做進(jìn)一步研究。

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

2 MINEMURA K, UCHIYAMA T. Three-dimensional calculation of air-water two-phase flow in centrifugal pump impeller based on a bubbly flow model[J]. ASME Journal of Fluids Engineering, 1993, 115(4): 766-771.

3 盧金鈴, 席光, 祁大同. 離心泵葉輪內(nèi)氣液兩相三維流動(dòng)數(shù)值研究[J]. 工程熱物理學(xué)報(bào), 2003, 24(2): 237-240. LU Jinling, XI Guang, QI Datong. Numerical study on the gas-liquid two-phase 3-D flow in the impeller of a centrifugal pump[J]. Journal of Engineering Thermophysics, 2003, 24(2): 237-240. (in Chinese)

4 陳次昌, 劉正英, 劉天寶, 等. 兩相流泵的理論與設(shè)計(jì)[M]. 北京:兵器工業(yè)出版社, 1993.

5 TREMANTE A, MORENO N, REY R, et al. Numerical turbulent simulation of the two-phase flow (liquid/gas) through a cascade of an axial pump[J]. ASME Journal of Fluids Engineering, 2002, 124(2): 371-376.

6 CAMPO A, CHISELY E A. Experimental characterization of two-phase flow centrifugal pumps[J]. Journal of Bacteriology, 2010, 176(24): 7524-7531.

從第一模態(tài)的時(shí)間序列(圖2b)看出,感熱通量有明顯的年際和年代際變化,其中1983年、1998年和2012年有明顯的轉(zhuǎn)折,第一模態(tài)時(shí)間序列與氣候態(tài)時(shí)間序列變化較為相似,其相關(guān)系數(shù)為0.4(通過99%置信度的顯著性檢驗(yàn)),因此第一模態(tài)的時(shí)間序列可以很好地代表高原感熱通量隨時(shí)間變化。由于感熱通量具有明顯的年際和年代際變化,為此對第一模態(tài)的時(shí)間序列進(jìn)行了小波變換分析。從分析結(jié)果(圖2c)看,高原感熱通量具有4～5 a的主周期和8 a的副周期。

7 ZHU J, ZHANG H. CFD simulation of ESP performance and bubble size estimation under gassy conditions[C]∥SPE Annual Technical Conference and Exhibition, 2014:SPE-170727-PT.

8 CARIDAD J, KENYERY F. CFD analysis of electric submersible pumps (ESP) handling two-phase mixtures[J]. Journal of Energy Resources Technology, 2004, 126(2): 99-104.

9 CARIDAD J A, KENYERY F. Slip factor for centrifugal impellers under single and two-phase flow conditions[J]. Journal of Fluids Engineering, 2005, 127(2): 317-321.

10 CARIDAD J, ASUAJE M, KENYERY F, et al. Characterization of a centrifugal pump impeller under two-phase flow conditions[J]. Journal of Petroleum Science & Engineering, 2008, 63(1): 18-22.

11 BARRIOS L, PRADO M G. Modeling two phase flow inside an electrical submersible pump stage[J]. ASME Journal of Energy Resources Technology, 2009, 133(4): 227-231.

12 潘兵輝, 王萬榮, 江偉. 離心泵氣液兩相流數(shù)值分析[J]. 石油化工應(yīng)用, 2011, 30(12): 101-104. PAN Binghui, WANG Wanrong, JIANG Wei. Gas-liquid two-phase flow numerical simulation of centrifugal pump[J]. Petrochemical Industry Application, 2011, 30(12): 101-104. (in Chinese)

13 余志毅, 劉影. 葉片式混輸泵氣液兩相非定常流動(dòng)特性分析[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2013, 44(5): 66-69. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20130513&flag=1. DOI：10.6041/j.issn.1000-1298.2013.05.013. YU Zhiyi, LIU Ying. Characteristic analysis of unsteady gas-liquid two-phase flow in a multiphase rotodynamic pump[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(5): 66-69. (in Chinese)

14 付強(qiáng), 習(xí)毅, 朱榮生, 等. 含氣率對AP1000核主泵影響的非定常分析[J]. 振動(dòng)與沖擊, 2015, 34(6): 132-136. FU Qiang, XI Yi, ZHU Rongsheng, et al. AP1000 nuclear main pump internal unsteady analysis under gas-liquid two phase condition[J]. Journal of Vibration and Shock, 2015, 34(6): 132-136. (in Chinese)

15 BARRIOS L, PRADO M G, KENYERY F. CFD modeling inside an electrical submersible pump in two-phase flow condition[C]∥ASME 2009 Fluids Engineering Division Summer Meeting, 2009: 457-469.

16 張金亞, 蔡淑杰, 朱宏武, 等. 三級螺旋軸流式混輸泵可壓縮流場數(shù)值模擬[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2014, 45(9): 89-95. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=2014095&flag=1. DOI：10.6041/j.issn.1000-1298.2014.09.015. ZHANG Jinya, CAI Shujie, ZHU Hongwu, et al. Numerical investigation of compressible flow in a three-stage helico-axial multiphase pump[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(9): 89-95 (in Chinese)

17 周水清, 孔繁余, 王志強(qiáng), 等. 基于結(jié)構(gòu)化網(wǎng)格的低比轉(zhuǎn)數(shù)離心泵性能數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2011, 42(7): 66-69. ZHOU Shuiqing, KONG Fanyu, WANG Zhiqiang, et al. Numerical simulation for low specific-speed centrifugal pump with structured grid[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(7): 66-69. (in Chinese)

18 MüLLER T, LIMBACH P, SKODA R. Numerical 3D RANS simulation of gas-liquid flow in a centrifugal pump with an Euler-Euler two-phase model and a dispersed phase distribution[C]∥Proceedings of 11th European Conference on Turbomachinery Fluid Dynamics & Thermodynamics, 2015: 760-772.

19 劉建瑞, 蘇起欽. 自吸泵氣液兩相流數(shù)值模擬分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2009, 40(9): 73-76. LIU Jianrui, SU Qiqin. Numerical simulation on gas-liquid two-phase flow in self-priming pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(9):73-76. (in Chinese)

20 MINEMURA K, UCHIYAMA T. Three-dimensional calculation of air-water two-phase flow in a centrifugal pump based on a bubbly flow model with fixed cavity[J]. JSME International Journal, Series B, Fluids & Thermal Engineering, 1994, 37(4): 726-735.

Numerical Investigation of Gas-Liquid Two-phase Flow in Centrifugal Pumps Based on Inhomogeneous Model

YUAN Jianping1ZHANG Keyu1SI Qiaorui1ZHOU Banglun1TANG Yuanfeng1JIN Zhongkun2

(1.NationalResearchCenterofPumps,JiangsuUniversity,Zhenjiang212013,China2.InstituteofAcoustics,ChineseAcademyofSciences,Beijing100190,China)

In order to study the flow characteristics of centrifugal pumps when transporting the gas-liquid mixture, water and air were chosen as the working medium, numerical simulation was conducted on a centrifugal pump under different conditions of inlet gas volume fraction (IGVF) based on the Eulerian-Eulerian inhomogeneous model. The drag force was approximated by the Schiller Nauman model. No other interfacial forces were considered. The gas distribution and velocity streamline in the impeller were obtained to discuss the gas and water flow characteristics of the pump. The results showed that gas concentration was high at the inlet pressure side of the blade, where the vortex would exist, indicating that the gas concentration had a great relationship with the vortex aggregation in the impeller passages. When the IGVF was increased to 10%, phase separation had obviously appeared at the suction side of the blades in the impeller passages, and gas had a movement trend along the suction side to the outlet of the impeller. Gas was moved from the center of the passages to the front and back shroud, and with IGVF increased, gas fraction in the back shroud near the leading edge and back shroud would be increase. Gas fraction in the front shroud near the trading edge was increased significantly than that in the back shroud, and the gas may eventually choke the passages. Within a rotation period, the outlet pressure of impeller was cyclically changed, when IGVF was increased from 1% to 10%, the outlet pressure of impeller was gradually decreased, and pressure pulsation frequency of the monitors was near the blade passing frequency at different IGVF values. When the IGVF was not higher than 10%, IGVF had no obvious effects on the pressure pulsation frequency and time frequency of the monitoring points. By comparing the experimental results with the numerical results, the reliability of the mathematical model and calculation methods was confirmed.

centrifugal pumps; gas-liquid transportation; Eulerian-Eulerian inhomogeneous model; pressure pulsation; numerical simulation

10.6041/j.issn.1000-1298.2017.01.012

2016-04-26

2016-08-12

“十二五”國家科技支撐計(jì)劃項(xiàng)目(2015BAD20B01)、國家自然科學(xué)基金項(xiàng)目(51509108)、江蘇大學(xué)高級人才基金項(xiàng)目 (15JDG048)和江蘇省自然科學(xué)基金項(xiàng)目(BK20150516)

袁建平(1970—)，男，研究員，博士生導(dǎo)師，主要從事流體機(jī)械內(nèi)部流動(dòng)及優(yōu)化設(shè)計(jì)研究，E-mail： yh@ujs.edu.cn

TH311； O359+.1

A

1000-1298(2017)01-0089-07