趙以奎,許建強(qiáng)
(1.合肥通用機(jī)械研究院,安徽 合肥 230031;2.江蘇大學(xué)流體機(jī)械工程技術(shù)研究中心,江蘇 鎮(zhèn)江 212013)
固液兩相流離心泵廣泛地應(yīng)用于礦山、冶金、電力和化工等領(lǐng)域,是進(jìn)行精礦、尾礦、灰渣、泥沙等固體物料水力輸送的關(guān)鍵核心設(shè)備。由于固液兩相流離心泵內(nèi)部漿體對(duì)過流部件的沖蝕、磨損作用,使葉輪和蝸殼的使用壽命成為重要的考核應(yīng)用指標(biāo)[1]。隨著計(jì)算機(jī)技術(shù)及計(jì)算流體動(dòng)力學(xué)(CFD)技術(shù)的飛速發(fā)展,固液兩相流泵的內(nèi)部流動(dòng)特性和磨損機(jī)理已經(jīng)取得了較好的研究成果[2-5]。Roco M C等[6-8]分析顆粒的運(yùn)動(dòng)規(guī)律后得出過流部件的磨損原因?yàn)轭w粒的直接沖擊、紊亂沖擊和滑動(dòng)摩擦。劉娟[9]對(duì)固液兩相流泵的磨損進(jìn)行了數(shù)值計(jì)算,研究發(fā)現(xiàn)大顆粒沿著葉片工作面運(yùn)動(dòng),易與葉片頭部發(fā)生撞擊,在多次撞擊中葉片受到嚴(yán)重磨損;而小顆粒運(yùn)動(dòng)趨勢(shì)則與之相反。李昳[10]通過對(duì)固液兩相流泵的研究,揭示了漿體流動(dòng)對(duì)泵過流部件磨損特性的影響,研究發(fā)現(xiàn)葉片和隔舌處分別發(fā)生滑動(dòng)磨損和沖擊磨損。
該基于CFD技術(shù),從固相顆粒屬性及其他外部因素出發(fā)研究了顆粒的運(yùn)動(dòng)規(guī)律,通過數(shù)值計(jì)算獲得離心泵葉輪和蝸殼內(nèi)固相體積濃度分布和固液兩相流動(dòng)特性,定量分析固體顆粒對(duì)過流部件磨損的影響,為高性能的固液兩相流離心泵的研制提供指導(dǎo)。
該文研究的模型泵設(shè)計(jì)參數(shù)流量(Q)為2440 m3/h,揚(yáng)程(H)為40 m,轉(zhuǎn)速(n)為 990 r/min,比轉(zhuǎn)速(ns)為187,葉片數(shù)(z)為6,葉輪結(jié)構(gòu)形式為半開式。工作介質(zhì)為含有石灰石顆粒的溶液,顆粒的平均直徑為60~100 μm。在對(duì)模型泵的內(nèi)部流場(chǎng)數(shù)值模擬之前,采用三維造型軟件PRO/E分別對(duì)葉輪、壓水室進(jìn)行建模,得到脫硫泵整體的三維實(shí)體造型,葉輪實(shí)體造型見圖1,模型泵的流道實(shí)體造型見圖2。
圖1 葉輪實(shí)體造型Fig.1 Solid model of impeller
圖2 模型泵全流道實(shí)體造型Fig.2 Full flow field model of pump
考慮計(jì)算模型較為復(fù)雜,網(wǎng)格劃分采用含有六面體、四面體、錐體和楔形體的混合網(wǎng)格。葉輪網(wǎng)格節(jié)點(diǎn)數(shù)為184434,蝸殼網(wǎng)格節(jié)點(diǎn)數(shù)499157。
數(shù)值計(jì)算采用三維不可壓縮N-S方程描述泵內(nèi)部流動(dòng),選用標(biāo)準(zhǔn)k-ε湍流模型封閉方程組;采用SIMPLEC算法求解壓力速度耦合方程組,對(duì)基本方程的離散差分主要用如下格式:對(duì)流項(xiàng)的離散采用二階迎風(fēng)格式,擴(kuò)散項(xiàng)、源項(xiàng)采用二階中心差分格式。多相流模型采用Mixture模型,考慮相與相之間的滑移速度,液固交換系數(shù)采用Wen-Yu模型,顆粒碰撞歸還系數(shù)設(shè)為0.9。葉輪與壓水室的耦合采用多重參考坐標(biāo)系法。
定義水為第一相,石灰石為第二相且全部為單一同種顆粒,各相為連續(xù)不可壓縮相,固體粒徑均勻且不發(fā)生相變。假定進(jìn)口速度軸向均勻分布,兩相之間均勻分布且速度等于固體輸送速度;采用自由出口邊界條件;在固壁處采用無滑移邊界條件,近壁處采用標(biāo)準(zhǔn)壁函數(shù)。
設(shè)計(jì)流量下顆粒體積分?jǐn)?shù)(Cv)為11%,顆粒直徑d分別為0.06,0.08和0.10 mm時(shí),葉輪中截面固相體積分?jǐn)?shù)分布見圖3a;設(shè)計(jì)工況下,顆粒直徑(d)為0.06 mm,顆粒體積分?jǐn)?shù)分別為7%,9%和11%時(shí),葉輪中截面固相體積分?jǐn)?shù)分布見圖3b。
圖3 葉輪中截面上固相體積分?jǐn)?shù)分布Fig.3 Solid phase volume concentration distribution of middle section
從圖3a中可以看出,當(dāng)顆粒直徑從0.06 mm增大到0.10 mm時(shí),固體顆粒主要分布在葉輪工作面所在區(qū)域,說明顆粒直徑的變化對(duì)固相的離析作用有較大影響。當(dāng)顆粒直徑較小時(shí),固體顆粒有向葉輪背面逐漸靠近的運(yùn)動(dòng)趨勢(shì);而顆粒直徑較大時(shí),固體顆粒不斷向葉輪的工作面運(yùn)動(dòng),因而造成葉輪工作面上的固相體積分?jǐn)?shù)不斷增大,固體顆粒在葉輪工作面的積聚和增加,必將導(dǎo)致該區(qū)域發(fā)生切削磨損。
從圖3b中可以看出,在相同顆粒直徑下,葉輪內(nèi)的固相體積分?jǐn)?shù)呈現(xiàn)這樣的變化趨勢(shì):在葉輪進(jìn)口處,葉輪工作面的固相體積分?jǐn)?shù)較高;從葉輪進(jìn)口到出口,固相體積分?jǐn)?shù)逐漸降低;在葉輪出口處達(dá)到最低。同時(shí),葉輪工作面、背面的固相體積分?jǐn)?shù)都隨著初始固相體積分?jǐn)?shù)的增大而增大。顆粒直徑對(duì)這種變化趨勢(shì)沒有太大影響,但葉輪工作面、背面的固相體積分?jǐn)?shù)在大粒徑工況下的升高速度明顯大于小粒徑工況下。
結(jié)合葉輪固相相對(duì)速度矢量圖可以發(fā)現(xiàn),葉輪區(qū)域磨損主要出現(xiàn)在葉輪進(jìn)口處、出口段。葉輪進(jìn)口處破壞主要是由于顆粒的直接撞擊造成切削磨損,切削磨損導(dǎo)致葉輪的進(jìn)口半徑增大,從而加大液流沖角,由此造成液流能量轉(zhuǎn)換能力相對(duì)下降。與實(shí)驗(yàn)結(jié)果對(duì)比,實(shí)驗(yàn)中葉輪出口邊緣被磨成鋸齒狀,計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果相符。
設(shè)計(jì)流量下顆粒體積濃度為11%時(shí),不同顆粒直徑下壓水室中截面固相體積分?jǐn)?shù)分布見圖4。從圖4中可以看出,壓水室內(nèi)的固相濃度是從壓水室的進(jìn)口開始,隨著半徑的增大而逐漸增大,在壓水室最大半徑處達(dá)到最大值。同時(shí),壓水室內(nèi)的固相濃度分布分別隨固體顆粒直徑和初始固相濃度的增大而增大。在小粒徑、低濃度工況下,由于離心力較小,壓水室內(nèi)的固相濃度分布相對(duì)均勻,影響相對(duì)較小;在大粒徑、高濃度工況下,由于漿液從葉片流道中流出進(jìn)入壓水室后,固相速度高于液相速度,固體顆粒在離心力的作用下與液相在該區(qū)域發(fā)生嚴(yán)重分離現(xiàn)象,顆粒直徑越大,其離心力也越大,致使壓水室中固相濃度分布不均勻,大部分固體顆粒呈現(xiàn)向壓水室壁面偏移的趨勢(shì),在壓水室近壁區(qū)形成高濃度區(qū)域;該區(qū)域同時(shí)還存在較大的速度梯度且離壁面越近速度越小,由此導(dǎo)致一部分固體顆粒沉積在壓水室壁面上。
圖4 壓水室中截面上固相體積分?jǐn)?shù)分布Fig.4 Solid phase volume concentration distribution in middle section of pumping chamber
(1)顆粒直徑對(duì)固相的離析作用有較大影響,隨著顆粒直徑的不斷增大,固體顆粒呈現(xiàn)向葉輪工作面偏移的趨勢(shì),由此造成葉輪工作面的切削磨損;固相體積濃度的增加對(duì)顆粒在葉輪流道內(nèi)的分布略有影響;葉輪的磨損主要發(fā)生在葉輪進(jìn)口處和葉輪的出口段。
(2)壓水室內(nèi)的固相體積濃度分布隨著固體顆粒直徑和初始固相濃度的增大而增大,固體顆粒在離心力的作用下與液相在壓水室發(fā)生分離現(xiàn)象,造成壓水室內(nèi)的固相濃度分布不均勻,壓水室壁面附近的速度梯度比較大,造成壁面磨損。
[1]崔巧玲,李昳,金守泉.非定常固液兩相流對(duì)雙流道泵磨損性能的影響[J].浙江理工大學(xué)學(xué)報(bào),2012,29(4):543-548.
[2]劉建瑞,徐永剛,王董梅,等.離心泵葉輪固液兩相流動(dòng)及泵外特性數(shù)值分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,41(3):86-90.
[3]張玉良,李昳,崔寶玲,等.兩相流離心泵水力輸送性能計(jì)算分析[J].機(jī)械工程學(xué)報(bào),2012,48(14):169-176.
[4]許海明,莊衛(wèi)將,邰俊,等.旋流泵固液兩相流數(shù)值模擬[J].中國農(nóng)村水利水電,2011(6):86-88.
[5]劉棟,楊敏官,李輝,等.化工泵葉輪內(nèi)部固液兩相流場(chǎng)的研究[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)固體顆粒運(yùn)動(dòng)規(guī)律與磨損的數(shù)值模擬[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,39(6):54-59.
[10]李昳,何偉強(qiáng),朱祖超,等.脫硫泵固液兩相流動(dòng)的數(shù)值模擬與磨損特性[J]. 排灌機(jī)械,2009,27(2):124-128.