04-9-粒径对离心泵叶片摩擦和碰撞磨损强度影响分析修改稿

粒径对离心泵叶片摩擦和碰撞磨损强度影响分析*程效锐,董富弟,杨从新,赵伟国,张楠(兰州理工大学能源与动力工程学院,甘肃兰州730050)摘要：采用雷诺时均N-S方程、RNGk-ε模型和SIMPLE算法，以含沙水为介质，基于代数滑移混合物模型(Algebraicslipmixturemodel,ASME)对一台单级双吸式离心泵内固液两相流动进行了全三维不可压缩定常流动数值模拟，其中转子与定子之间耦合方式采用“冻结转子法”实现。通过对比清水及含沙水介质时泵外特性试验数据与数值模拟结果，验证了数值计算方法的可靠性。对固相颗粒直径分别为0.019、0.036、0.076mm时叶片工作面和背面摩擦磨损强度和碰撞磨损强度进行了预测。结果表明：在同一粒径时，从叶片进口至出口碰撞磨损强度逐渐增大，且工作面大于背面，摩擦磨损强度呈现先增大后减小，又逐渐增大的趋势；随着粒径的增大，叶片表面碰撞磨损强度和摩擦磨损强度逐渐增大，摩擦磨损强度沿着整个叶片均大于碰撞磨损强度。关键词：离心泵；固液两相流；数值模拟；磨损预测；磨损模型中图分类号：TH31Analyzethefrictionandcollisionwearstrengthtocentrifugalpumpbladeswithdifferentparticlesize*1CHENGXiao-rui,DONGFu-di,YANGCong-xin,ZHAOWei-guo,ZHANGNan(CollegeofEnergyandPowerEngineering，LanzhouUniv.ofTech.，Lanzhou730050,China)Abstract：BasedonAlgebraicslipmixturemodel,usingSIMPLEalgorithmandRNGk-εturbulencemodelandfrozenrotormethod,thesolid-liquidtwo-phaseflowindouble-suctionscentrifugalpumpwassimulated.TheUnsteadyReynoldsAveragedNavier-Stokes(RANS)approachhasbeenappliedtosolvetheunsteady,incompressible,three-dimensionalturbulent.Thewearratesoffrictionandcollisionatthepressureandsuctionsurfacesofbladeswithdifferentsanddiameters(0.019mm、0.036mm、0.076mm)wereforecasted.Bycomparisonsoftheexternalcharacteristicaboutthetestandnumericalsimulationdatesinadouble-suctioncentrifugalpumpwhenpumpedclearwaterandsolid-liquidtwo-phaseflow,thereliabilityofnumericalcalculationmethodhasbeenverified.Theresultsshowedthatwiththesamesanddiameter,fromthebladesinlettooutlet,thecollision收稿日期：2014-04-09基金项目：国家自然科学基金（51269011），甘肃省自然科学基金（2011GS04252）作者简介：程效锐(1972-),男,甘肃兰州人,博士,副教授.E-mail：[e-mail]cesisincreasesandthepressuresurfacesisseriousthansuctionsurfaces,andthefrictionwearrateisfirstincreasesthendecreases,andagainincreases.Withsanddiameterincreasing,thewearratesoffrictionandcollisionatthepressureandsuctionsurfacesofbladesareincreasesandthewearrateoffractionisseriousthancollision.Keywords：centrifugalpump；solid-liquidtwo-phaseflow；numericalsimulation；predictionoferosion；wearmodel双吸式离心泵由于叶轮的对称结构几乎不产生轴向力，且便于维修，被广泛应用于引黄灌区的提水灌溉中[1]。黄河中的悬移质泥沙造成泵过流部件材料磨损、使泵扬程和效率降低。因此研究双吸式离心泵在含沙水流下的流动规律及沙粒属性变化对材料磨损特性具有重要意义。目前针对离心泵内固液两相流动的研究主要是借助计算流体动力学（computationalfluiddynamics，CFD）理论，结合部分试验方法对泵在含沙水流下的流动及磨损特性进行综合分析，并已取得了一定成果[2-8]。但鉴于泵内固液两相流动的复杂性和现有试验手段的制约，导致对离心泵叶轮在固液两相流作用下的磨损机理和定量磨损预测方面的研究还不够深入。本研究运用CFD数值模拟方法，结合颗粒摩擦和碰撞模型[9]，对双吸式离心泵内的固液两相流动进行了全三维数值计算。对叶轮叶片摩擦磨损和碰撞磨损进行定量预测，并且对影响这两种磨损的相关参数在固液两相流动中的变化规律做进一步研究分析，得到了固相颗粒直径(以下简称粒径)对双吸式离心泵叶片磨损特性的影响规律，为双吸式离心泵磨损机理和定量磨损预测方面的研究及优化设计提供依据。1数学模型的建立1.1基本假设以含沙水流为研究介质，由于双吸式离心泵内部流动复杂，为了深入研究固液两相流内部流动及磨损机理，现做出以下假设：1）整个流道内介质为定常流动，转子和定子之间耦合方式采用”冻结转子法”实现；2）固相为直径均匀的球形颗粒，不考虑相变；3）固相为连续介质，液相为不可压缩流体，且每相的物理特性均为常数。1.2数值计算方法和边界条件采用“冻结转子法”，建立相对坐标系下的雷诺时均连续方程和全三维不可压缩时均N-S方程，通过RNG-k模型来简化和封闭方程组。两相流采用代数滑移混合物模型（algebraicslipmixturemodel，ASMM）[4]来描述。速度和压力之间耦合方式采用SIMPLE算法。对流项和扩散项的空间离散均采用一阶迎风格式，源项的空间离散采用线性化标准格式。各个控制方程收敛精度均为10-5。假设壁面绝热且满足无滑移条件，近壁面区域采用标准壁面函数处理，壁面材料取为Ni-Cr白口铸铁（ρ=7600kg/m3）。采用速度进口边界条件，假定进口处浓度分布均匀且等于固相输送浓度。出口采用自由出流条件。1.3磨损模型图1为固相颗粒沿叶片运动的磨损示意图。将颗粒对叶片的磨损分为两种形式：一种是颗粒与叶片摩擦运动造成的摩擦磨损；另一种是颗粒与叶片碰撞运动造成的碰撞磨损。为便于分析比较本研究将摩擦磨损率WSL、碰撞磨损率WI(μm/h)[9]分别定义为摩擦磨损强度和碰撞磨损强度。图1固相沿叶片的碰撞磨损和摩擦磨损Fig.1Thesandcollisionandfrictionwearalongtheblade1）摩擦磨损：nkkSL1SpNkCuWEd(1)式中：sl8Spslslsl10nEdAdCB为单位摩擦磨损能，J/m3，nC是固相浓度；kdudl为剪切应力，Pa；其中μ为流体动力粘度，而dudl表示任一点的速度梯度；ku为固相切向速度，m/s；d为粒径，μm；其他参数是与材料相关的试验常数，对于Ni-Cr白口铸铁，取值如下[9]：14sl4.23610A；sl180B；sl490C；sl3.861n。2）碰撞磨损：3knkI1Ikk,NkCVWEd(2)式中：IkkadjIk,EdCE；22kkkVuv；IkO41EE；imadjimimim1.65nCAdCB；ku和k分别为固相速度kV在x和y方向上的速度分量（叶轮绕z轴旋转）；为固体颗粒碰撞角，1kkktanu，rad；k为固体颗粒密度，等于2500kg/m3；IkE为单位碰撞磨损能，J/m3;OE为当2时磨损表面系数，和粒径及壁面材料有关，J/m3；Cadj为自调整因子，其他参数是与材料相关的试验常数，对于Ni-Cr白口铸铁，取值如下[9]：15O2.7710E；10im=8.57810A；2im=3.64510B；im=490C；im3.668n。1.4计算模型图2为本研究双吸式离心泵的三维模型。研究的双吸式离心泵设计工况下的性能参数为:流量qv=400m3/h，扬程H=33m，转速n=1480r/min。叶轮主要参数见表1。图2双吸式离心泵三维模型Fig.23Dmodelofthedouble-suctioncentrifugalpump表1叶轮主要参数Tab.1Themainparametersoftheimpeller参数数值叶轮出口直径D2/mm355轮毂直径Dh/mm100叶轮进口直径D1/mm200叶片进口安放角β2(°)23叶片出口安放角β1(°)38叶片出口宽度B2/mm22.5叶片数6叶片包角φ(°)1641.5网格划分吸水室、压水室和叶轮三部分计算区域均采用非结构化网格及不同网格尺度进行离散，叶片表面进行局部加密处理。通过对比分析数值模拟结果，进行了网格无关性检查。最终确定计算域网格总数为15875640，吸水室、蜗壳及叶轮区域的网格数分别为6147840、6709830、3017970。2结果与分析2.1数值计算可靠性验证为验证数值模拟结果的可靠性，采用开式实验台对双吸式离心泵进行外特性试验。经多年实测分组粒径资料统计，黄河上游流域多年平均粒径为0.036mm[10]。图3为清水及含沙水（固相颗粒属性：Cv=10%，d=0.036mm，ρk=2500kg/m3）介质下泵的扬程和效率在八个工况点处（0.4qv～1.4qv）的数值模拟值和试验值的比较。0.40.60.81.01.21.4151821242730333639扬程试验值（清水）效率试验值（清水）扬程计算值（固液两相）扬程试验值（固液两相）效率计算值（固液两相）效率试验值（固液两相）455055606570758085q/qvη/%H/m图3.数值模拟与试验外特性曲线Fig.3Theexternalcharacteristiccurveofnumericalsimulationandexperimental数据显示模型泵扬程曲线没有出现“驼峰”现象，效率具有较宽的“高效区”。含沙水介质时泵扬程和效率的试验及计算值明显小于清水介质时的外特性试验数据，可见固相介质的存在造成了严重的能量损失，影响了泵的外特性。固液两相流介质时计算结果与试验值吻合较好，最大扬程误差不超过2%，最大效率误差不超过3.8%，在允许误差范围内。证明了代数滑移混合物模型在描述固液两相流流动时具有良好的适用性。2.2颗粒直径对叶片磨损特性的影响在设计流量及固相浓度为Cv=10%，密度ρk=2500kg/m3情况下，分别对粒径为0.019mm、0.036mm和0.076mm时叶片的磨损特性进行计算分析。为便于分析比较，将叶片平均中心流线从进口至出口距离定义为相对长度L=1。2.2.1粒径对摩擦磨损特性的影响从式(1)可以看出，固相切向速度、剪切应力和固相浓度分布是影响摩擦磨损强度的主要因素。图4为不同直径固相切向速度沿叶片变化规律。从图中可以看出，固相切向速度不随粒径的变化而