田俊良,鄭泉*,劉和明

(1.安徽農(nóng)業(yè)大學 工學院，安徽 合肥 230036；2.明光留香泵業(yè)有限公司，安徽 明光 239400)

基于雙向耦合的吸沙泵外殼與葉片摩損仿真分析

田俊良1,鄭泉1*,劉和明2

(1.安徽農(nóng)業(yè)大學 工學院，安徽 合肥 230036；2.明光留香泵業(yè)有限公司，安徽 明光 239400)

為研究吸沙泵磨損機理，建立吸沙泵的三維幾何模型和兩相流模型，通過仿真與試驗數(shù)據(jù)的對比驗證模型的準確性。分別采用單向和雙向流固耦合方法對建立的固液兩相流模型進行數(shù)值模擬，并對葉片及殼體相關區(qū)域的速度和壓力分布進行分析。結果表明：固相粒子普遍集中在泵體內(nèi)緣，出口處固相分布多于進口處；與單向耦合相比，雙向耦合更符合實際情況；摩擦磨損主要集中于泵體內(nèi)緣和葉片的工作面沿半徑方向外側(cè)；粒徑越大，固相分離現(xiàn)象越明顯，增加轉(zhuǎn)速使揚程降低，但是轉(zhuǎn)速過低會使效率急劇下降。

吸沙泵；多相流；雙向耦合；磨損

吸沙泵作為一種特殊用泵應用廣泛，但由于水、沙的磨損使其壽命縮短，因此有必要對復雜邊界條件下湍流結構的影響和高濃度固相粒子對摩擦磨損機理進行深入研究。目前國內(nèi)外對泵的研究主要集中在水利特性、內(nèi)部流場瞬態(tài)水力激振特性下的流固耦合以及動靜湍流場的非定常特性等方面，例如：文獻[1]采用發(fā)展的RN-S方程和液相/氣相界面跟蹤方法的單相空化模型數(shù)值求解技術，開展利用空化數(shù)對離心泵的水力性能影響特性的數(shù)值研究；文獻[2]基于代數(shù)滑移混合物模型(algebraic slip model,ASM)模型，對三維不可壓縮定常流動進行數(shù)值計算并采用凍結轉(zhuǎn)子法，確定最佳水力特性下轉(zhuǎn)子的轉(zhuǎn)動位置；文獻[3]分別考慮不耦合和耦合作用對泵流場內(nèi)部工作情況的影響；文獻[4]研究旋流泵輸送固液兩相流特性，證明旋流泵內(nèi)部兩相流動符合畸變速度原理；文獻[5-7]針對泵內(nèi)部流場瞬態(tài)水力激振特性對軸流泵葉頂區(qū)空化流氣液混合區(qū)域密度變化、軸流泵的流量-揚程曲線、空化特性及其誘導非定?？栈瘔毫γ}動以及水泵液壓系統(tǒng)瞬態(tài)特性分別進行研究；文獻[8]針對兩相流下葉片磨損情況的復雜性，考慮不同固相體積分數(shù)以及相對運動速度，得出葉輪磨損情況；文獻[9-10]使用二維渦方法研究離心泵葉輪和導葉(離心泵的轉(zhuǎn)能裝置)內(nèi)部的非定常壓力脈動。研究成果主要集中于水泵類產(chǎn)品，而對于吸沙泵這類特殊用泵以及其泵體內(nèi)部和葉片摩擦磨損分析等方面尚未深入研究。

本文采用ANSYS和CFX對吸沙泵進行耦合分析，研究不同工況下的速度場分布和壓強變化，并根據(jù)葉輪旋轉(zhuǎn)域以及殼體域的綜合分析得出相應葉片和殼體受力情況，對可能出現(xiàn)的摩擦磨損以及葉片表面剝離情況做出預測，為葉片葉形優(yōu)化和偏角優(yōu)化提供參考[11-12]。

1 流固兩相流理論

本文模型采用水沙固液兩相流且將水看成不可壓縮流體、而沙子默認為擬流體[13]，假設沙子外形成球形，粒徑固定不變且不發(fā)生相變。在 Eulerian 坐標系下，采用混合代數(shù)滑移模型(mixture algebraic slip model,MASM)[14]，其中兩相連續(xù)方程為：

(1)

(2)

(3)

式中：ρm為混合密度；ρk為k相密度；Umj為混合相的混合平均速度;Xi(i=1,2,3)為歐拉坐標;αk為k相體積系數(shù)；Ukj為k相的平均速度。

兩相動量方程為：

(4)

式中：Umi為混合相整體對流運動速度；pm為混合相動壓；τmi、τTmi、τDmi分別為平均粘性應力、紊動應力和擴散應力;gi為附加質(zhì)量力系數(shù);Fki為k相所受的作用力。

相對速度方程為：

(5)

式中：Ucpi為相間的相對速度；mp為混合相質(zhì)量。

2 吸沙泵模型與邊界條件

2.1幾何與網(wǎng)格模型

基于Catia軟件對吸沙泵進行實體建模，并將三維模型導入CFX軟件中建立流體模型[15-18]。吸沙泵的額定轉(zhuǎn)速為750 r/min、葉片數(shù)為4、進口直徑和出口直徑均為406 mm，揚程為14 m，流量為750 m3/h,配套電機功率為60 kW，設沙子密度為2 650 kg/m3，吸沙泵的模型如圖1a)所示，殼體域與旋轉(zhuǎn)域如圖1b)所示。基于ICEM軟件對流體域進行非結構網(wǎng)格劃分，外殼域選擇四面體網(wǎng)格，如圖2a)所示，旋轉(zhuǎn)區(qū)域采用了六面體網(wǎng)格劃分，如圖2b)所示，并在外殼域表層和細小尖角區(qū)域網(wǎng)格加密，共產(chǎn)生1 325 487個網(wǎng)格，對網(wǎng)格質(zhì)量檢查，殼體域與旋轉(zhuǎn)域的綜合網(wǎng)格平均值為0.789，標準差為0.188，符合網(wǎng)格質(zhì)量要求。

圖1 吸沙泵模型與內(nèi)部流模型示意圖 圖2 流體域網(wǎng)格

2.2邊界條件與算法

圖3 試驗與模擬的揚程-流量關系對比

外殼壁面與葉輪壁面采用無滑移的邊界條件，采用速度入口和自由壓力出口；在壓力與速度耦合上，采用RNGk-ε湍流模型[19-21]進行求解；對流區(qū)域的空間離散采用一階迎風格式；為了簡化計算，認為所有固相顆粒均是相同固定直徑的球形顆粒物；各個控制方程方根RMS殘差格式的收斂依據(jù)是10-5。在數(shù)值模擬計算中，控制方程采用MSAM模型；計算采用Double精度，基礎時間步長設置為2×10-4s；在收斂比較困難的部分，適當減小時間步長，加快收斂。

2.3外工作特性

為了驗證上述模型以及邊界條件的設置準確性，以單相清水為介質(zhì)，開展試驗測試，根據(jù)兩種耦合方式的數(shù)值模擬計算結果得出不同流量下?lián)P程外特性曲線，并與試驗揚程曲線進行對比[19-21]，如圖3所示。由圖3可知：數(shù)值模擬計算揚程與試驗結果誤差較小，均不超過6%，因此模型是可信的。

3 吸沙泵仿真分析

3.1不同粒徑下速度場分布

在設計轉(zhuǎn)速下，當工作介質(zhì)采用不同粒徑的沙子與水組成的固液兩相混合物時，泵腔內(nèi)的速度矢量分布如圖4所示(圖中數(shù)字單位為m/s)。

圖4 不同粒徑沙相速度矢量圖

由圖4可知：無論是單向還是雙向耦合，不同粒徑沙相的運行軌跡及運動狀態(tài)基本相似，且具有幾乎相等的最大速度；泵體內(nèi)輪緣處沙的速度比內(nèi)部輪轂速度快，且在高速運轉(zhuǎn)產(chǎn)生的強離心力作用下，泵體內(nèi)緣沙質(zhì)速度流線隨粒徑增大而集中，且轉(zhuǎn)速偏高區(qū)域的面積逐漸增大，這表明在該區(qū)域固相顆粒粒徑越大對于泵體的磨損越嚴重。出口處壓力大，水沙混合物隨葉輪運行到出口處受壓形成漩渦；隨著粒徑增大，沖擊能力增強，漩渦現(xiàn)象反而明顯減弱。通過對不同耦合方式的結果對比分析可以看出，渦旋位置分布以及出口段速度分布規(guī)律基本相同。隨著粒徑增加，葉輪工作背面速度大于工作正面，部分沙速度流線從葉輪工作背面向葉輪工作正面遷移，與工作背面相比，固相顆粒對于工作面磨損幾率明顯增大。

由此得知：粒徑的改變對于單向耦合和雙向耦合作用的整體差異影響不大，但在出口、隔舌、高低速運轉(zhuǎn)區(qū)和泵體輪緣的速度變化比較明顯；且與單向耦合作用相比，雙向耦合作用流線更為紊亂，這是因為考慮了葉輪振動以及流體對葉輪的作用反饋。

圖5 設計流量下葉輪壓強分布

3.2不同轉(zhuǎn)速下壓強分析

通過對速度場分析可知，在離心段外圈貼合泵體內(nèi)邊緣處速度比較大，則易存在較大的摩擦磨損；而葉片區(qū)域沙流速度變化不顯著。采用固相粒徑為0.01 mm沙子和水混合，在混合相中沙子的質(zhì)量分數(shù)為60%，其壓強分布結果如圖5所示(單位為MPa)。由圖5可知：工作背面的壓強普遍小于工作正面；沿著工作正面壓強分布往外擴展，而在入口處附近葉片的工作背面，則出現(xiàn)了小范圍的相對負氣壓，因而此處極易發(fā)生氣蝕現(xiàn)象，附近葉輪固體壁面在局部溫度和局部壓強的反復作用下被剝蝕。

圖6 兩種耦合方式下0.01 mm粒徑壓強分布圖

對葉片工作正面輪緣附近選一點(即圖5中的A點)進行壓強分析，如圖6所示。由圖6可知：雙向耦合方式下某點壓強隨轉(zhuǎn)速增加而增加(轉(zhuǎn)速和流量在效率不變、粘度不大的情況下成正相關)，而單向耦合得到的該點壓強在各個轉(zhuǎn)速下均大于雙向耦合。隨轉(zhuǎn)速增大，揚程下降，而沙泵的效率先增大后減小，且壓強逐漸增大，使葉片沿半徑方向工作表面磨損加劇，因而實際工作時均是選擇合適的轉(zhuǎn)速，這樣既保證了合適的揚程，又可以使沙泵的磨損情況比較理想。而對于粒徑為0.1 mm和1 mm下該點的壓強分布，也滿足圖6的分布規(guī)律。

4 結論

1)離心作用使固相粒子普遍集中在泵體內(nèi)緣，出口處固相分布多于進口處，出口易磨損；考慮雙向耦合的情形下，離心區(qū)域的流道分布更加紊亂，在出口、隔舌和離心區(qū)域產(chǎn)生的渦旋對于內(nèi)部流道產(chǎn)生較大影響。

2)離心泵轉(zhuǎn)速的提高導致流量增加，對于隔舌的最大壓強持續(xù)增大，且隨粒徑的增加最大壓強有所波動，不利于改善隔舌磨損；葉片工作面固相分布多于工作背面，對于進口處附近的葉片極易產(chǎn)生氣蝕，長期會導致葉片剝離，且壓強逐漸增大使得葉片沿半徑方向工作表面磨損加劇。

[1]李軍,劉立軍,李國君，等.空化數(shù)對離心泵水力性能影響的數(shù)值研究[J].工程熱物理學報，2010,31(5)：773-776. LI Jun, LIU Lijun,LI Guojun,et al. Numercial investigations on effects of cavitaion number on hydrodynamic performance of centrifugal pump[J]. Journal of Engineering Thermophysics, 2010, 31 (5) : 773-776.

[2]張玉良, 李昳, 崔寶玲,等. 兩相流離心泵水力輸送性能計算分析[J].機械工程學報, 2012, 48(14):169-176. ZHANG Yuliang, LI Yi, CUI Baoling, et al. Computational analysis of hydraulic transmission performance in two-phase flow centrifugal pump[J].Journal of Mechanical Engineering, 2012, 48 (14) : 169-176.

[3]袁壽其, 徐宇平, 張金鳳,等. 流固耦合作用對螺旋離心泵流場影響的數(shù)值分析[J].農(nóng)業(yè)機械學報, 2013, 44(1):38-42. YUAN Shouqi, XU Yuping, ZHANG Jinfeng, et al. Numercial analysis for effect of fluid-structure interaction on flow field in screw cventrifugal pump[J].Transactions of the Chinese Society for Agricultural Machinery, 2013, 44 (1) : 38-42.

[4]沙毅, 劉祥松. 旋流泵固液兩相流輸送特性試驗[J].農(nóng)業(yè)工程學報, 2013(22):76-82. SHA Yi, LIU Xiangsong. Performance test on solid-liquid two-phase flow hydrotransport of vortex pump[J].Transactions of the Chinese Society of Agricultural Engineering, 2013 (22) : 76-82.

[5]張德勝, 潘大志, 施衛(wèi)東,等. 軸流泵葉頂區(qū)的空化流場與葉片載荷分布特性[J].化工學報, 2014(2):501-507. ZHANG Desheng, PAN Dazhi, SHI Weidong, et al. Cavitation flow and blade loading characteristic in impeller tip region of axial flow pump[J].Journal of Chemical Industry and Engineering, 2014 (2) : 501-507.

[6]張德勝, 潘大志, 施衛(wèi)東,等. 軸流泵空化流及其誘導壓力脈動的數(shù)值模擬[J].華中科技大學學報(自然科學版), 2014,42(1):34-38. ZHANG Desheng, PAN Dazhi, SHI Weidong, et al. Numerical simulation of cavitation flow in axial flow pump and induced press[J].Journal of Huazhong University of Science and Technology(Natural Science Edition), 2014,42 (1) : 34-38.

[7]張德勝, 施衛(wèi)東, 潘大志,等. 基于數(shù)值模擬的特種混流泵水力性能優(yōu)化與試驗[J].機械工程學報, 2014, 50(6):177-184. ZHANG Desheng, SHI Weidong, PAN Dazhi,et al. Hydralic performance optimization and experiment of special mixed-flow pump based on numerical simulation[J].Journal of Mechanical Engineering, 2014, 50 (6) : 177-184.

[8]李仁年, 趙文杰, 李琪飛. 螺旋離心泵固液兩相流中的葉輪磨損分析[J].甘肅科學學報, 2014, 26(1):67-70. LI Rennian, ZHAO Wenjie, LI Qifei. Analysis of impeller abrasion by liquid-solid two-phase medium in screw centrifugal pumps[J].Journal of Gansu Science, 2014, 26 (1) : 67-70.

[9]WANG H，TSUKAMOTO H. The oretical study of pressure fluctuations downstream of diffuser pump impeller:part2:effects of volute,flow rate and radial gap[J].ASME Journal of Fluids Engineering,2001,123:737-747.

[10]WANG H. Experimental and numerical study of unsteady flow in a diffuser pump at off-design conditions[J].Journal of Fluids Engineering, 2003, 125(5):767-778.

[11]施衛(wèi)東, 李啟鋒, 陸偉剛,等. 基于CFD的離心泵軸向力計算與試驗[J].農(nóng)業(yè)機械學報, 2009, 40(1):60-63. SHI Weidong, LI Qifeng, LU Weigang, et al. Estimation and experiment of axial thrust in centrifugal pump based on CFD[J].Transactions of the Chinese Society for Agricultural Machinery, 2009, 40 (1) : 60-63.

[12]楊麗霞, 謝東, 田建濤,等. ANSYS 軟件在齒輪泵靜力分析中的應用[J].山東交通學院學報, 2014,22(3):75-78. YANG Lixia, Xie Dong ,TIAN Jiantao, et al. Application of ANSYS software tostatic analysis of gear pump[J].Journal of Shandong Jiaotong University, 2014 ,22(3) : 75-78.

[13]裴吉,袁壽其. 離心泵非定常流動特性及流固耦合機理[M].北京：機械工業(yè)出版社, 2014.

[14]趙琴,李嘉,安瑞冬.水庫渾水異重流的兩相流模型適用性研究[J].水動力學研究與進展,2010,25(1):76-84. ZHAO Qin, LI Jia, AN Ruidong. Two-phase flow models for turbid density current in a reservoir[J].Chinese Journal of Hydrodynamics, 2010,25(1):76-84.

[15]丁思云,蘇吉鑫,鄭國運,等. 雙葉片螺旋離心泵固液兩相流特性研究[J].流體機械, 2012, 40(5):26-29. DING Siyun, SU Jixin, ZHENG Guoyun, et al. Study on two-phase flow characteristics of double blade screw centrifugal pump solid-liquid[J].Fluid Machinery, 2012, 40 (5) : 26-29.

[16]王洋, 劉洋, 王維軍,等. 雙葉片排污泵固液兩相流的數(shù)值模擬[J].中國農(nóng)村水利水電, 2013(11):121-125. WANG Yang, LIU Yang, WANG Weijun, et al. A numerical simulation of solid-liquid two-phase flow in two-blade sewage pump[J].China Rural Water and Hydropower, 2013 (11) : 121-125.

[17]徐海良, 陳奇, 吳波. 深海采礦礦漿泵固液兩相流數(shù)值計算及分析[J].計算機仿真, 2014, 31(9):432-435. XU Hailiang, CHEN Qi,WU Bo. Numerical calculation and analysis of solid-liquid two-phase flow in slurry pump for deep-sea mining[J].Computer Simulation, 2014, 31 (9) : 432-435.

[18]劉勝柱, 郭鵬程, 羅興锜. 離心泵葉輪內(nèi)的網(wǎng)格生成與計算流體力學分析[J].農(nóng)業(yè)工程學報, 2004, 20(5):78-81. LIU Shengzhu, GUO Pengcheng, LUO Xingqi. Grid generation in the centrifugal pump impeller and computational fluid dynamics analysis[J].Transactions of the Chinese Society of Agricultural Engineering, 2004, 20 (5) : 78-81.

[19]耿少娟, 聶超群, 黃偉光,等. 不同葉輪形式下離心泵整機非定常流場的數(shù)值研究[J].機械工程學報, 2006, 42(5):27-31. GENG Shaojuan, NIE Chaoqun, HUANG Weiguang, et al. Numerical study of unsteady flow in centrifugal pump with different type of impellers[J].Journal of Mechanical Engineering, 2006, 42 (5) : 27-31.

[20]邵春雷, 顧伯勤, 陳曄. 離心泵內(nèi)部非定常壓力場的數(shù)值研究[J].農(nóng)業(yè)工程學報, 2009, 25(1):75-80. SHAO Chunlei, GU Boqin, CHEN Ye. Numerical simulation of unsteady pressure field in centrifugal pumps[J].Transactions of the Chinese Society of Agricultural Engineering, 2009, 25 (1) : 75-80.

[21]張德勝, 施衛(wèi)東, 陳斌,等. 低比轉(zhuǎn)速離心泵內(nèi)部流場分析及試驗[J].農(nóng)業(yè)工程學報, 2010, 26(11):108-113. ZHANG Desheng, SHI Weidong, CHEN Bin, et al. Turbulence analysis and experiments of low-specific-speed centrifugal pump[J].Transactions of the Chinese Society of Agricultural Engineering, 2010, 26 (11) : 108-113.

SimulationAnalysisofFrictionLossofPumpSuctionSandShellandBladeBasedonTwo-WayCoupling

TIANJunliang1,ZHENGQuan1*,LIUHeming2

(1.CollegeofTechnology,AnhuiAgriculturalUniversity,Hefei230036,China；2.MingguangLiuxiangPumpCo.,Ltd.,Mingguang239400,China)

In order to study the mechanism of the sand suction pump wear, the three-dimensional geometric model and two-phase flow model of the sand suction pump is established and the accuracy of the models through the comparison of simulation data and experimental data is verified. The unidirectional and bidirectional fluid-solid coupling methods are adopted respectively to make the numerical simulation of the established solid-liquid two-phase flow model, the velocity and pressure distribution of the relevant area of the blade and shell is analyzed. The results show as follows. The solid phase particles are generally focused on the inner edge of the pump body, more solid phase distribution at the exit than that at the inlet. Compared with the one-way coupling, the two-way coupling is more in line with the actual situation. The friction and wear are mainly concentrated on the inner edge of the pump body and the blade working face along the radius direction of the lateral. The bigger the particle is, the more obvious the solid phase separation phenomenon is. The rotating speed increase lowers the head, but the low speed drops sharply in efficiency.

sand suction pump; multiphase flow; two-way coupling; friction loss

TH117.1

：A

：1672-0032(2017)03-0088-06

(責任編輯：郎偉鋒)

2016-10-10

安徽省重點研究與開發(fā)計劃項目(1704a0902045)

田俊良(1989—)，男，合肥人，碩士研究生，主要研究方向為流體力學，E-mail:junliang106@163.com.

*通訊作者：鄭泉(1970—)，女，安徽蕪湖人，教授，碩士生導師，主要研究方向為機械設計，E-mail:wangww0618@163.com.

10.3969/j.issn.1672-0032.2017.03.014