趙以奎，許建強

(1.合肥通用機械研究院，安徽 合肥 230031;2.江蘇大學(xué)流體機械工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

固液兩相流離心泵廣泛地應(yīng)用于礦山、冶金、電力和化工等領(lǐng)域，是進行精礦、尾礦、灰渣、泥沙等固體物料水力輸送的關(guān)鍵核心設(shè)備。由于固液兩相流離心泵內(nèi)部漿體對過流部件的沖蝕、磨損作用，使葉輪和蝸殼的使用壽命成為重要的考核應(yīng)用指標(biāo)[1]。隨著計算機技術(shù)及計算流體動力學(xué)(CFD)技術(shù)的飛速發(fā)展，固液兩相流泵的內(nèi)部流動特性和磨損機理已經(jīng)取得了較好的研究成果[2-5]。Roco M C等[6-8]分析顆粒的運動規(guī)律后得出過流部件的磨損原因為顆粒的直接沖擊、紊亂沖擊和滑動摩擦。劉娟[9]對固液兩相流泵的磨損進行了數(shù)值計算，研究發(fā)現(xiàn)大顆粒沿著葉片工作面運動，易與葉片頭部發(fā)生撞擊，在多次撞擊中葉片受到嚴重磨損;而小顆粒運動趨勢則與之相反。李昳[10]通過對固液兩相流泵的研究，揭示了漿體流動對泵過流部件磨損特性的影響，研究發(fā)現(xiàn)葉片和隔舌處分別發(fā)生滑動磨損和沖擊磨損。

該基于CFD技術(shù)，從固相顆粒屬性及其他外部因素出發(fā)研究了顆粒的運動規(guī)律，通過數(shù)值計算獲得離心泵葉輪和蝸殼內(nèi)固相體積濃度分布和固液兩相流動特性，定量分析固體顆粒對過流部件磨損的影響，為高性能的固液兩相流離心泵的研制提供指導(dǎo)。

1 計算模型

該文研究的模型泵設(shè)計參數(shù)流量(Q)為2440 m3/h，揚程(H)為40 m，轉(zhuǎn)速(n)為 990 r/min，比轉(zhuǎn)速(ns)為187，葉片數(shù)(z)為6，葉輪結(jié)構(gòu)形式為半開式。工作介質(zhì)為含有石灰石顆粒的溶液，顆粒的平均直徑為60～100 μm。在對模型泵的內(nèi)部流場數(shù)值模擬之前，采用三維造型軟件PRO/E分別對葉輪、壓水室進行建模，得到脫硫泵整體的三維實體造型，葉輪實體造型見圖1，模型泵的流道實體造型見圖2。

圖1 葉輪實體造型Fig.1 Solid model of impeller

圖2 模型泵全流道實體造型Fig.2 Full flow field model of pump

考慮計算模型較為復(fù)雜，網(wǎng)格劃分采用含有六面體、四面體、錐體和楔形體的混合網(wǎng)格。葉輪網(wǎng)格節(jié)點數(shù)為184434，蝸殼網(wǎng)格節(jié)點數(shù)499157。

2 計算方法

數(shù)值計算采用三維不可壓縮N-S方程描述泵內(nèi)部流動，選用標(biāo)準(zhǔn)k-ε湍流模型封閉方程組;采用SIMPLEC算法求解壓力速度耦合方程組，對基本方程的離散差分主要用如下格式:對流項的離散采用二階迎風(fēng)格式，擴散項、源項采用二階中心差分格式。多相流模型采用Mixture模型，考慮相與相之間的滑移速度，液固交換系數(shù)采用Wen-Yu模型，顆粒碰撞歸還系數(shù)設(shè)為0.9。葉輪與壓水室的耦合采用多重參考坐標(biāo)系法。

定義水為第一相，石灰石為第二相且全部為單一同種顆粒，各相為連續(xù)不可壓縮相，固體粒徑均勻且不發(fā)生相變。假定進口速度軸向均勻分布，兩相之間均勻分布且速度等于固體輸送速度;采用自由出口邊界條件;在固壁處采用無滑移邊界條件，近壁處采用標(biāo)準(zhǔn)壁函數(shù)。

3 數(shù)值計算結(jié)果與分析

3.1 固體顆粒對葉輪磨損的影響

設(shè)計流量下顆粒體積分數(shù)(Cv)為11%，顆粒直徑d分別為0.06，0.08和0.10 mm時，葉輪中截面固相體積分數(shù)分布見圖3a;設(shè)計工況下，顆粒直徑(d)為0.06 mm，顆粒體積分數(shù)分別為7%，9%和11%時，葉輪中截面固相體積分數(shù)分布見圖3b。

圖3 葉輪中截面上固相體積分數(shù)分布Fig.3 Solid phase volume concentration distribution of middle section

從圖3a中可以看出，當(dāng)顆粒直徑從0.06 mm增大到0.10 mm時，固體顆粒主要分布在葉輪工作面所在區(qū)域，說明顆粒直徑的變化對固相的離析作用有較大影響。當(dāng)顆粒直徑較小時，固體顆粒有向葉輪背面逐漸靠近的運動趨勢;而顆粒直徑較大時，固體顆粒不斷向葉輪的工作面運動，因而造成葉輪工作面上的固相體積分數(shù)不斷增大，固體顆粒在葉輪工作面的積聚和增加，必將導(dǎo)致該區(qū)域發(fā)生切削磨損。

從圖3b中可以看出，在相同顆粒直徑下，葉輪內(nèi)的固相體積分數(shù)呈現(xiàn)這樣的變化趨勢:在葉輪進口處，葉輪工作面的固相體積分數(shù)較高;從葉輪進口到出口，固相體積分數(shù)逐漸降低;在葉輪出口處達到最低。同時，葉輪工作面、背面的固相體積分數(shù)都隨著初始固相體積分數(shù)的增大而增大。顆粒直徑對這種變化趨勢沒有太大影響，但葉輪工作面、背面的固相體積分數(shù)在大粒徑工況下的升高速度明顯大于小粒徑工況下。

結(jié)合葉輪固相相對速度矢量圖可以發(fā)現(xiàn)，葉輪區(qū)域磨損主要出現(xiàn)在葉輪進口處、出口段。葉輪進口處破壞主要是由于顆粒的直接撞擊造成切削磨損，切削磨損導(dǎo)致葉輪的進口半徑增大，從而加大液流沖角，由此造成液流能量轉(zhuǎn)換能力相對下降。與實驗結(jié)果對比，實驗中葉輪出口邊緣被磨成鋸齒狀，計算結(jié)果與實驗結(jié)果相符。

3.2 固體顆粒對壓水室磨損的影響

設(shè)計流量下顆粒體積濃度為11%時，不同顆粒直徑下壓水室中截面固相體積分數(shù)分布見圖4。從圖4中可以看出，壓水室內(nèi)的固相濃度是從壓水室的進口開始，隨著半徑的增大而逐漸增大，在壓水室最大半徑處達到最大值。同時，壓水室內(nèi)的固相濃度分布分別隨固體顆粒直徑和初始固相濃度的增大而增大。在小粒徑、低濃度工況下，由于離心力較小，壓水室內(nèi)的固相濃度分布相對均勻，影響相對較小;在大粒徑、高濃度工況下，由于漿液從葉片流道中流出進入壓水室后，固相速度高于液相速度，固體顆粒在離心力的作用下與液相在該區(qū)域發(fā)生嚴重分離現(xiàn)象，顆粒直徑越大，其離心力也越大，致使壓水室中固相濃度分布不均勻，大部分固體顆粒呈現(xiàn)向壓水室壁面偏移的趨勢，在壓水室近壁區(qū)形成高濃度區(qū)域;該區(qū)域同時還存在較大的速度梯度且離壁面越近速度越小，由此導(dǎo)致一部分固體顆粒沉積在壓水室壁面上。

圖4 壓水室中截面上固相體積分數(shù)分布Fig.4 Solid phase volume concentration distribution in middle section of pumping chamber

4 結(jié)論

(1)顆粒直徑對固相的離析作用有較大影響，隨著顆粒直徑的不斷增大，固體顆粒呈現(xiàn)向葉輪工作面偏移的趨勢，由此造成葉輪工作面的切削磨損;固相體積濃度的增加對顆粒在葉輪流道內(nèi)的分布略有影響;葉輪的磨損主要發(fā)生在葉輪進口處和葉輪的出口段。

(2)壓水室內(nèi)的固相體積濃度分布隨著固體顆粒直徑和初始固相濃度的增大而增大，固體顆粒在離心力的作用下與液相在壓水室發(fā)生分離現(xiàn)象，造成壓水室內(nèi)的固相濃度分布不均勻，壓水室壁面附近的速度梯度比較大，造成壁面磨損。

[1]崔巧玲，李昳，金守泉.非定常固液兩相流對雙流道泵磨損性能的影響[J].浙江理工大學(xué)學(xué)報，2012，29(4):543-548.

[2]劉建瑞，徐永剛，王董梅，等.離心泵葉輪固液兩相流動及泵外特性數(shù)值分析[J].農(nóng)業(yè)機械學(xué)報，2010，41(3):86-90.

[3]張玉良，李昳，崔寶玲，等.兩相流離心泵水力輸送性能計算分析[J].機械工程學(xué)報，2012，48(14):169-176.

[4]許海明，莊衛(wèi)將，邰俊，等.旋流泵固液兩相流數(shù)值模擬[J].中國農(nóng)村水利水電，2011(6):86-88.

[5]劉棟，楊敏官，李輝，等.化工泵葉輪內(nèi)部固液兩相流場的研究[J].水泵技術(shù)，2007(3):16-19.

[6]Roco M C.Modeling erosion wear in centrifugal slurry pumps[C].Symposium of 9th Int1.Conf.on Hydraulic Transport of Solid in Pipes，1984:291-316.

[7]Schilling Rudo l.Numerical simulation of two-phase flow in centrifugal pump impellers[C]Procedings of the 2002 ASME JointU. S- European Fluids Engineering Conference，Simposia，2002.

[8]Mehta M，Kadambi J R，Sastry S，et al.Study of particulate flow in the impeller of a slurry pump using PIV[C].ASME 2004 Heat Transfer/Fluids Engineering Summer Conference(HT-FED2004)，2004:489-499.

[9]劉娟，許洪元，唐澍，等.離心泵內(nèi)固體顆粒運動規(guī)律與磨損的數(shù)值模擬[J].農(nóng)業(yè)機械學(xué)報，2008，39(6):54-59.

[10]李昳，何偉強，朱祖超，等.脫硫泵固液兩相流動的數(shù)值模擬與磨損特性[J]. 排灌機械，2009，27(2):124-128.