颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响.pdf

颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响王凯，柳涵宇，王李科，卢金玲（西安理工大学省部共建西北旱区生态水利国家重点实验室，西安710048）摘要：固液两相流条件下半开式叶轮离心泵中颗粒冲击、泄漏涡发展和颗粒轨迹之间存在紧密交互作用，导致过流部件的磨损行为复杂多变。该研究结合双向耦合欧拉-拉格朗日方法和颗粒磨损 Finnie 模型，对不同颗粒体积浓度下半开式叶轮离心泵固液两相流场进行求解，分析了颗粒体积浓度对泄漏涡结构特征、颗粒运移轨迹和磨损特性的影响，揭示了颗粒体积浓度、叶顶间隙泄漏涡和过流部件表面磨损规律的关联机制。结果表明：随着颗粒体积浓度的增加，颗粒的频繁撞击加剧了叶片压力面进水边和后盖板磨损程度，叶片吸力面出水边的磨损范围向进水边方向延伸；颗粒体积浓度小于 1%时，颗粒的轴向运动和叶顶间隙泄漏涡的阻碍作用导致颗粒易与叶片前缘靠近叶根处和吸力面出水边靠近叶顶的区域发生撞击，诱发严重磨损，且呈现点状磨损；当颗粒体积浓度大于 3%时，叶轮后盖板的整体磨损强度大于叶片，颗粒体积浓度的增加造成流入叶顶间隙层的颗粒数增加，颗粒对叶顶间隙泄漏涡的冲击导致涡流的破碎、分离、再融合，加剧不稳定流动，泵的扬程和效率均明显下降。该研究可为固液两相半开式叶轮离心泵优化设计和安全稳定运行提供理论参考。关键词：离心泵；模型；磨损；固液两相流；泄漏涡doi：10.11975/j.issn.1002-6819.202305243中图分类号：TK730文献标志码：A文章编号：1002-6819(2023)-16-0044-10王凯，柳涵宇，王李科，等.颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响J.农业工程学报，2023，39(16)：44-53.doi:10.11975/j.issn.1002-6819.202305243http:/www.tcsae.orgWANG Kai,LIU Hanyu,WANG Like,et al.Effects of particle volume concentration on the leakage vortex and erosioncharacteristicsofsemi-opencentrifugalpumpsJ.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2023,39(16):44-53.(inChinesewithEnglishabstract)doi:10.11975/j.issn.1002-6819.202305243http:/www.tcsae.org0引言与闭式叶轮离心泵相比，半开式叶轮以其结构简单紧凑、维护方便、便于输送黏性较大流体等优点，在石油化工、水利工程等领域越来越受到青睐1。但是，叶片与泵体形成的叶顶间隙会导致在叶片压力面和吸力面压差的作用下产生流动分离、二次流和脱流等不稳定流动，使泵内流动损耗较高2-4。中国河流含沙量较大，离心泵长期运行于此类泥沙浓度较大的流体环境中时，水中的固体颗粒会对过流部件产生撞击和切削，使过流表面产生变形，发生磨损破坏，导致机组运行稳定性变差，效率降低5-7。目前国内外研究学者对颗粒浓度对水力机械的磨损速率的影响进行了大量的研究。LEI 等8基于 SST 湍流模型，结合欧拉-拉格朗日方法模拟泵内泥沙颗粒的运动轨迹和磨损机理，从而建立颗粒参数与泵性能及转轮叶片磨损情况的关联。SHEN 等9基于欧拉-拉格朗日的离散相模型，模拟了不同固相质量浓度工况下的螺杆离心泵内部流场，分析了固相颗粒的能量变化、轨迹特征及磨蚀规律。LAI 等10采用双向耦合欧拉-拉格朗日方法对离心泵进行磨蚀计算并预测磨蚀机制，得到了过流部件磨损严重的区域分布。张自超等11-12考虑相间阻力和泥沙扩散系数两方面的影响，基于改进的 Euler-Euler 算法对双吸离心泵流场进行数值模拟，得出湍流强度和颗粒动态尺度对磨损严重的块状磨损区分布。张玉良等13基于代数滑移混合物模型对离心泵内部流场进行定常流动数值计算，发现在靠近隔舌的叶轮出口处的三种速度组合的双剪切层射流-尾流结构和固相体积分数对叶片表面的磨损程度影响。汪家琼等14基于 Particle 模型和非均相模型，分析了离心泵过流部件壁面处固体颗粒体积浓度与磨损程度间的关系。周世杰等15基于相关试验改进了固相粒子动量方程，结合 Alert 磨损经验公式对抽黄用泵进行数值模拟计算和分析，固相体积浓度递增加剧了叶轮进口回流现象和部件表面磨损程度。半开式叶轮叶顶间隙泄漏流研究方面，张青扬等16对半开式叶轮离心泵进行全流道数值模拟，发现叶顶间隙影响叶轮流道中的压力分布，产生的回流、涡流、扰流等不稳定流动严重影响泵内部流场性能参数。AYAD等17通过数值模拟研究发现不同流速条件下，发现不同流量工况下离心泵扬程和效率的变化与叶顶间隙变化量近似为回归效应，增加叶顶间隙宽度，二次流产生的涡旋会阻碍核心流动。MOUSMOULIS 等18采用了一种齐次方法分别对半开式径向叶轮离心泵在两相条件下进行数值模拟，分析了叶顶间隙区域与叶轮回流和靠近叶片前缘的二次流间的关联性。贾晓奇等19对不同工况下半开式叶轮离心泵内部三维湍流流场进行数值模拟，揭示收稿日期：2023-05-30修订日期：2023-08-08基金项目：国家自然科学基金面上项目（52376148）；国家自然科学基金重 点 项 目（51839010）；陕 西 省 教 育 厅 重 点 实 验 室 科 研 计 划 项 目（19JS045）作者简介：王凯，博士，讲师，研究方向为水力机械磨损-腐蚀预测及防护。E-mail：通信作者：卢金玲，博士，教授，研究方向为水力机械流动安全评估与防护。Email：jinling_第39卷第16期农 业 工 程 学 报Vol.39No.16442023年8月TransactionsoftheChineseSocietyofAgriculturalEngineeringAug.2023了不同流量工况下叶顶间隙泄漏涡的流动特性、叶顶间隙总压与相对速度分布规律。现有针对半开式叶轮离心泵的研究多集中在单相流条件下的不稳定流动特性，关于固液两相流磨损规律的研究多为颗粒参数与叶片磨损速率的外在关联。但是，从动力学机理的内在角度来讲，泄漏涡对流道的阻碍作用会改变颗粒的运动趋势，包括颗粒运移方向和速度，同时颗粒穿过涡会破坏涡结构，因此颗粒撞击、泄漏涡变化和颗粒运动轨迹之间的相互作用加剧了磨损行为的复杂性和多变性。而对于固液两相流条件下半开式叶轮离心泵在不同颗粒体积浓度时泄漏流特性、颗粒运动特性和磨损特性之间内在关系的研究较为匮乏，尤其缺乏颗粒运移对流体扰动的影响。因此，本文采用双向耦合的欧拉-拉格朗日方法，考虑流体-颗粒间的相互作用，同时结合 Finnie 磨损模型，通过数值模拟探究不同颗粒相体积浓度对半开式叶轮离心泵内固液两相流场及过流部件磨损规律的影响，深入分析颗粒运移、泄漏涡结构和过流表面磨损之间的交互作用机制，拟为固液两相流半开式叶轮离心泵优化设计和安全稳定运行提供理论参考。1数学模型1.1液相本文中的颗粒浓度较小，将半开式叶轮离心泵内液相均视为三维不可压缩流体，连续性方程和动量方程分别如式（1）和式（2）所示，动量交换源项为式（3）20：uixi=0（1）(uit+ujuixj)=xj(uixj+ujxi)pxi+gi+Fi（2）Fi=(FD+FB+FM+FP)Vcell（3）式中 u 为流体速度，m/s；p 为静压，Pa；为流体密度，kg/m3；为流体的动力黏度，Ns/m2；考虑固液两相的交互作用，通过增加一个固液两相动量交换的源项Fi（i=1,2,3）表示颗粒对流体的影响，FD为曳力，N；FB为浮力，N；FM为虚拟质量力，N；FP为压力梯度力，N；Vcell是一个计算单元体积。1.2固相固体离散相模型粒子轨迹的计算是独立的，单个粒子在特定时间间隔内被追踪运动轨迹21。半开式叶轮离心泵中颗粒运动轨迹的改变受重力、曳力、虚拟质量力、Basset 力、Saffman 升力和 Magnus 升力等主导，由于连续相为水相时 Basset 力、Saffman 升力和 Magnus 升力相对于曳力非常小，且 Saffman 升力和 Magnus 升力仅对直径处于 110m 的颗粒有明显影响22-23。本文的颗粒粒径为 0.5mm，且连续相为水相，因此仅对曳力、重力、浮力、虚拟质量力和压力梯度力进行分析。控制方程为广义牛顿第二定律，可表示为24-25mpdupdt=FD+FG+FB+FM+FP（4）式中 mp为粒子质量，kg；up为粒子速度，m/s。FD=mp18CDRep24pdp2(uup)（5）Rep=dp?uup?（6）p式中 CD为曳力系数；Rep为颗粒雷诺数；为颗粒密度，kg/m3；dp为颗粒直径，mm。重力 FG和浮力 FB可表示为FG+FB=mp(1p)g（7）虚拟质量力 FM可表示为FM=mppCM(dudtdupdt)（8）其中，CM为虚拟质量系数 0.5。压力梯度力 FP可表示为FP=mppdudt（9）1.3磨损模型颗粒对部件过流表面进行撞击和切割，导致部件变形，材料发生疲劳破坏。产生的磨损是与颗粒特性、冲击条件和壁材特性相关的函数。国内外学者常采用的磨损模型有：Finnie 模型，DNV 模型，Oka 模型以及Ahlert 磨损模型。DNV 模型、Oka 模型的模型修正数据均来自气固两相流动，与本文的液固两相流环境不同，其中 DNV 模型未考虑不同角度下磨损机制的变化26，且 Oka 模型参考的基准冲蚀速率为冲击角度 90时的工况，与离心泵常见的倾斜角度冲击机制不同27，而Ahlert 模型根据经验和大量试验数据进行拟合得出，并未从磨损过程中材料去除机制出发26。Finnie 模型从材料去除的颗粒位移和切削机制出发，并进行不同角度下磨损机制的区分，特别适用于离心泵等水力机械倾斜冲击角度下的磨损机制28，更符合本文的研究条件。因此，本文采用 Finnie 磨损模型研究固液两相条件下半开式叶轮离心泵过流部件的磨损规律。Finnie 提出的经典磨损模型如下29：E=kVnpf()（10）f()式中 E 为磨损量无量纲参数；k 为获得无量纲冲蚀因子的常数；Vp为固相颗粒撞击速度，m/s；n 为速度指数，与壁面材料属性相关，常用金属材料取值范围为 2.32.5；是关于颗粒冲击角的无量纲函数。f()=13cos2tan13sin(2)3sin2tan13（11）其中，=18.42是区分滑动磨损和冲击磨损的临界值，18.42，则冲击磨损占主导地位30。2计算模型及网格划分2.1模型描述研究对象为一台半开式叶轮离心泵，主要设计参数第16期王凯等：颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响45为：流量 Qd=47m3/h，扬程 Hd=7m，转速 n=980r/min。叶轮主要几何尺寸为：叶片进口安放角 1=16.6，叶片出口安放角 2=18.6，包角=99.4，叶轮进口直径 D1=100mm，叶轮出口直径 D2=232mm。计算域为进口管、叶轮和无叶扩压器（图 1a）。a.计算域出口Outlet进口Inleta.Calculation domain b.网格12345b.Mesh1.进口管2.叶片后缘3.叶轮4.叶片前缘5.无叶扩压器1.Inletpipe2.Bladetrailingedge3.Impeller4.Bladeleadingedge5.Vanelessdiffuser图 1半开式叶轮离心泵计算域和网格Fig.1Calculationdomainandmeshofsemi-opencentrifugalpump2.2网格划分及无关性验证采用 ANSYSICEM 软件对半开式叶轮离心泵进口管、叶轮和无叶扩压器模型进行六面体结构化网格划分。为消除网格数量对计算结果的影响，采用美国机械工程协会（AmericanSocietyofMechanicalEngineers，ASME）推荐的网格收敛指数 GCI 进行网格离散误差的估计31-33，对网格进行无关性验证。该方法需要设置三套不同数目的网格，分别为粗糙网格、中等网格和细密网格。本文选择叶片扭矩和叶轮出口总压作为关键变量，消除变量之间带来的运算误差。计算设计工况下 GCI 网格离散误差的估计方法中的各参数，得到叶片扭矩的近似相对误差（中等-粗糙和细密-中等）分别为 2.442%和0.436%；总压的近似相对误差分别为 2.175%和 0.199%；叶片扭矩的外推相对误差分别为 0.512%和 0.075%；总压的外推相对误差分别为 0.213%和 0.015%；叶片扭矩的不确定度分别为 0.643%和 0.093%；总压的不确定度分别为 0.267%和 0.019%，3 套网格以渐进形式收敛。综合考虑计算精度和效率，最终确定选取中等网格，各部分网格数为叶轮5561412，进口管427276，无叶扩压器976080，网格数共 6964768。半开式叶轮离心泵网格如图 1b 所示。2.3数学计算方法及边界条件设置本文基于 SSTk-两方程湍流模型，采用双向耦合欧拉-拉格朗日方法求解半开式叶轮离心泵内的固液两相流动，将液相视为连续相，固体颗粒视为离散相。叶轮设置为旋转域，进口管和无叶扩压器均设置为固定域，壁面处设置为无滑移壁面；流体介质为水，温度为 25，密度为 997.05kg/m3；进口处给定总压，假定颗粒为球形，颗粒密度为 2650kg/m3，进口处颗粒和流体速度相等，且颗粒在进口面上均匀分布；出口给定质量流量。2.4数值模拟的验证2.4.1试验设备及方法试验测量工作在离心泵闭式试验台上进行，试验系统的详细描述可参考我们的前期工作34。转速通过变频器进行调节，流量通过在出口管路上安装的电磁流量计测量，压力通过在进出口管路上布置的压力传感器测量，扭矩通过安装在泵与电机之间的扭矩转速仪测量。试验测量时通过数据采集设备将测量参数（转速、流量、进口压力、出口压力、扭矩）转换为电流信号，经信号线传送至数据采集系统，数据采集系统将其转换为可识别的数字信号，进行数据采集，计算得到离心泵的扬程和效率，对比清水工况下半开式叶轮离心泵的扬程和效率的计算值和试验值，验证数值模拟的可靠性。2.4.2外特性验证为了验证对半开式叶轮离心泵性能模拟的准确性，对清水工况下离心泵的扬程和效率的模拟和试验结果进行比较。如图 2 所示，数值模拟结果与试验值分布一致，扬程最大相对误差 7.52%，最小相对误差 0.96%，平均相对误差 3.77%；效率最大相对误差 7.78%，最小相对误差 5.61%，平均相对误差 6.41%。数值模拟与试验结果具有相同的规律，相对误差在可接受范围内。扬程 Head H/m效率 Efficiency/%无量纲流量 Dimensionless flow rate7.06.56.05.55.04.54.03.53.00.20.40.60.81.01.21.41.690807060504030模拟扬程Simulated head试验扬程Experimental head试验效率Experimental efficiency模拟效率Simulated efficiency图 2试验与数值模拟外特性结果对比图Fig.2Comparisonoftheheadandefficiencybetweenexperimentalandnumericalresults2.5试验设计及方法根据中国泥沙问题严重的河流沿岸水力机械运行规律，平均固相体积浓度超过 2%，最大固相体积浓度超过 4%，同时参考了国内外众多学者9,14,30在研究不同固相浓度工况下离心泵的内部流场时所选择的固相浓度范围（多为 1%10%）。本文选择在设计转速 980r/min、设计流量 1.0Qd、dp=0.5mm 颗粒粒径的条件下，对比分析 4种颗粒相体积浓度（Cv=1%、3%、5%、7%）对半开式叶轮离心泵泄漏涡特征和磨损特性的影响。为了便于分析半开式叶轮离心泵内的固液两相特性、泄漏涡特性及磨损特性，采用 ANSYS 软件对计算结果进行后处理，得到叶轮内部的的压力分布、泄漏涡结构、颗粒轨迹、颗粒平均体积分数分布及磨损分布。采用 Omega 准则捕捉阈值 R 为 0.52 时各颗粒体积浓度下的泄漏涡结构，众多国内外学者证明阈值 R 为0.52 时具有阈值不敏感特性，并能同时有效识别强、中、弱多种涡结构35。为了定量识别涡结构，采用标准螺旋度 Hn 着色，Hn 表示涡矢量与相对速度矢量夹角的余弦值，二者方向相同时 Hn 为 1，正交时为 0，相反时为46农业工程学报（http:/www.tcsae.org）2023年134。为了能深入理解泄漏涡的运动机理，采用涡量输运方程对叶顶间隙区域的泄漏涡结构进行更详细的分析，涡量输运方程如下36-38：ddt=()W(W)2(W)+p2+v2（12）W式中为涡矢量；为相对速度矢量；为旋转角速度，rad/s；为流体运动黏度，Pas。方程左边表示涡量的变化速率，右边第一项为相对涡量拉伸项（Relativevortexstretching，RVS），由相对速度梯度引起涡的拉伸和弯曲；第二项为涡膨胀项，反映流体可压缩性引起的涡量变化情况；第三项为科氏力项（Coriolisforce，CORF），与旋转运动有关；第四项为斜压扭矩项，表示压力梯度与密度梯度对涡量的影响；最后一项黏性扩散项（Viscosdiffusion，VISD）主要是由于流体的黏性引起的。对于不可压缩流体，忽略第二项和第四项。由于叶顶间隙的存在，叶轮内部尤其是靠近叶顶区域的不稳定流动复杂，已有工作对叶轮进出口环量进行了分析2，发现叶顶间隙对叶轮的影响主要体现在靠近叶顶侧，因此选取 0.93 叶高（Span=0.93）（图 3）对清水和固液两相工况的流场和涡系特征（RVS、CORF 和 VISD）的分布规律进行分析。0.93叶高周向截面Circumferentialsection of 0.93 span0 Span 1叶轮进口Impeller inlet图 3叶高（Span）示意图Fig.3Diagramofthebladeheight(Span)3结果与讨论3.1半开式叶轮离心泵固液两相流动特性3.1.1颗粒体积浓度对泵性能的影响由图 4 可知，由于固体颗粒密度远大于液体的密度，加入较小体积浓度的颗粒时，固液混合介质的密度增加，泵的有效功率增加，同时颗粒的存在抑制了边界的流动分离，减少了涡流脱落造成的损失，因此泵的扬程和效率相比于清水工况均增大。随着颗粒体积浓度的增大，颗粒之间相互干扰，发生更多的碰撞和摩擦，流体输送颗粒所需的能耗增加，因此泵的扬程和效率明显下降。6.05.95.85.75.65.55.45.35.25.15.001234567884828078767472706866颗粒体积浓度 Particle volume concentration Cv/%扬程 Head H/m效率 Efficiency/%扬程 Head效率 Efficiency图 4不同颗粒体积浓度外特性变化曲线Fig.4Headandefficiencycurvewithdifferentparticlevolumeconcentration3.1.2叶轮内部流场特性图 5 为设计流量下 0.93 叶高截面的压力分布图。由于叶片对吸入介质做功，流场压力沿着叶轮进口向出口方向逐渐增大。清水工况下，由于叶顶间隙的存在，相对液流角的增大造成叶轮进口形成较大的负冲角，结合叶片压力面与吸力面的压差，二者共同驱动形成泄漏流，导致叶片吸力面进水边出现面积较大的低压区。相比于清水工况，颗粒的堆积导致叶轮流道内的压力值普遍升高，最大压力上升 4%。部分颗粒在叶轮进口处因流道变化和叶轮旋转的离心力，运动方向由轴向变为径向，与叶片前缘发生碰撞，抑制了边界层的不稳定流动，低压区面积减小，但颗粒碰撞后速度减小，流场内部分压能又转化为驱动颗粒继续运动的动能，产生低压，因此低压区面积仅轻微减小。3.2颗粒相体积浓度对泄漏涡特性的影响3.2.1泄漏涡结构特征各颗粒体积浓度下的泄漏涡结构如图 6 所示。颗粒体积浓度为 1%时，少量颗粒的加入对泵内涡流结构的影响较弱，流道中泄漏涡、分离涡和通道涡与清水工况相似。随着颗粒体积浓度的增加，在曳力、离心力和惯性力的耦合作用下，颗粒对流体的撞击频率和跟随流体流入叶顶间隙的颗粒数量增加，流动的不均匀性增强，涡流的破碎程度逐渐加剧，破碎分离的涡流间通过周向运动进行融合。同时，流体旋转消耗的能量增大，抑制了泄漏涡、叶顶分离涡、通道涡的发展，其中二次泄漏涡与主流卷吸形成的二次泄漏流强度明显减弱。3.2.2泄漏涡系特征图 7 为叶顶间隙层的相对涡量拉伸项（RVS）分布，随着颗粒相体积浓度的增加，颗粒对叶片前缘的碰撞更加频繁，驱动叶片前缘的流动分离，因此叶片前缘的高RVS 区逐渐增大。颗粒相体积浓度的改变对高 RVS 区域在叶片吸力面侧的带状分布趋势影响较小，但颗粒的冲击导致吸力面侧破碎分离出更多的漩涡。科氏力项（CORF）分布如图 8 所示，叶轮的旋转产生科氏力并导致在叶片压力面和吸力面形成压力梯度，泄漏流受到压力梯度的驱动而形成。当颗粒相体积浓度为 1%时，固液两相工况下的 CORF 分布与清水工况相似，说明较小颗粒浓度对叶轮旋转产生的压力梯度影响作用较小。随着颗粒相体积浓度的增大，叶片前缘受颗粒撞击脱落的涡旋强度增加，叶片吸力面侧的高 CORF分布区波动均增大，面积减小，反映出流场不稳定性增第16期王凯等：颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响47强。当颗粒相体积浓度增至 7%时，区域 A 中吸力面中部的高 CORF 区明显波动、分裂，这是由于颗粒的不断累积，在这一区域的流体驱动其偏转跟随运动消耗的能量增多，造成涡量出现一定程度的减少。86 00092 50088 70087 500100 80093 90097 000107 700114 600121 500128 400135 300142 200149 100156 00094 400101 300108 200115 100122 000128 900135 800142 700149 600156 50085 630102 560109 490116 420123 350130 280137 210144 140151 070158 00099 000105 500112 000118 500125 000131 500138 000144 500151 000压力 Pressure/Pa压力 Pressure/Pa压力 Pressure/Pa88 30095 170102 040108 910115 780122 650129 520136 390143 260150 130157 000a.Cv=0 b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%压力 Pressure/Pa压力 Pressure/Pa图 5不同颗粒体积浓度（Cv）0.93 叶高截面压力分布图Fig.5Pressuredistributionof0.93spanwithdifferentparticlevolumeconcentration(Cv)1.000.780.560.330.110.110.330.560.781.00标准螺旋度 Standard helicity12345涡流破碎Brcakdown流动分离Flow separntion泄漏涡1TLVI+二次泄漏涡Second TLV流动融合Flow fusiona.Cv=0 b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%1.叶顶分离涡2.尾迹涡3.叶顶泄漏涡4.通道涡5.二次叶顶泄漏涡1.Tipseparationvortex(TSV)2.Wakevortex(WV）3.Tipleakagevortex(TLV)4.Passagevortex(PV)5.Secondtipleakagevortex(TLV2)图 6不同颗粒体积浓度泄漏涡结构Fig.6Structureofleakagevortexwithdifferentparticlevolumeconcentration0e+001e+052e+053e+054e+055e+056e+057e+058e+059e+0510e+06相对涡量拉伸项 Relative vortex stretching/S2流动方向Flow directiona.Cv=0 b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%图 7不同颗粒体积浓度相对涡量拉伸项分布Fig.7Distributionofrelativevortexstretchingitemwithdifferentparticlevolumeconcentration黏性扩散项（VISD）分布如图 9 所示，在叶顶间隙层，叶片压力面与吸力面的压差形成微射流，微射流的运动方向与叶片旋转方向相反，因此与主流和固定的前盖板形成相对运动，导致在壁面处产生涡量。随着颗粒体积浓度的增加，固液两相介质粘性作用的变化及颗粒的运动影响涡流向流体内部的扩散，涡旋难以稳定发展，因此高 VISD 区的波动更加明显，同时伴随着剧烈的涡旋分离、融合。0e+001e+052e+053e+054e+055e+056e+057e+058e+059e+0510e+06科氏力项 Coriolis force/S2流动方向Flow direction区域 ARegion Aa.Cv=0 b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%图 8不同颗粒体积浓度科氏力项分布Fig.8DistributionofCoriolisforceitemwithdifferentparticlevolumeconcentration3.3颗粒相体积浓度对磨损特性的影响不同颗粒体积浓度下颗粒速度及轨迹如图 10 所示。当 Cv3%时，随着颗粒体积浓度的增加，流入叶顶间隙的颗粒数量增多，运动至叶片吸力面出水边靠近叶顶处堆积的颗粒增多。0e+001e+052e+053e+054e+055e+056e+057e+058e+059e+0510e+06黏性扩散项Viscos diffusion/S2流动方向Flow directiona.Cv=0 b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%图 9不同颗粒体积浓度黏性扩散项分布Fig.9Distributionofviscosdiffusionitemwithdifferentparticlevolumeconcentration同时，颗粒体积浓度越大，颗粒间的排挤导致颗粒逐渐占据整个流道，在叶片吸力面出水边靠近叶根处堆积的颗粒数量增多，且明显多于在叶顶附近聚集的颗粒，说明大部分颗粒在流道中运动时靠近后盖板，叶顶间隙泄漏涡扩散到叶根处时对颗粒的运动已无明显影响。不同颗粒体积浓度下叶片表面及后盖板的平均颗粒体积分数分布情况如图 11 所示。当颗粒体积浓度增大，同等粒径下颗粒数量增多，由于颗粒本身粒径较大，则质量和惯性力较大。颗粒速度 Particle velocity/(ms1)012345678910b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%图 10不同颗粒体积浓度下颗粒速度和轨迹Fig.10Particlevelocityandtrajectorieswithdifferentparticlevolumeconcentrationa.后盖板a.Hubb.叶片表面b.Blade颗粒平均体积分数 Particle averaged volume fraction/%012345678910颗粒平均体积分数 Particle averaged volume fraction/%012345吸力面Suction surface压力面Pressure surface吸力面Suction surface压力面Pressure surface吸力面Suction surface压力面Pressure surface吸力面Suction surface压力面Pressure surfaceb.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%b.Cv=1%c.Cv=3%e.Cv=7%d.Cv=5%图 11不同颗粒体积浓度下平均颗粒体积分数分布Fig.11Particleaveragedvolumefractionwithdifferentparticlevolumeconcentration第16期王凯等：颗粒体积浓度对半开式叶轮离心泵泄漏涡和磨损的影响49在进入叶轮后颗粒保持原有的轴向运动趋势的能力较强，容易较早的偏转对后盖板进行冲击。同时，压力梯度力和曳力的耦合作用限制了颗粒跟随流线移动，驱使其产生明显的向叶片压力面侧运动的趋势，因此随着颗粒体积浓度增大，聚集在后盖板和叶片压力面侧进水边颗粒明显增多，且分布在后盖板的颗粒多于叶片。图 12 所示，随着颗粒相体积浓度的增大，固相颗粒与过流部件的接触和碰撞几率增加，过流部件的磨损情况逐渐加重。颗粒在流道中运动时，径向的惯性力不能抵消叶轮旋转的离心力，导致后盖板磨损严重区域发生在曲率最大处及叶片曲率最大处的交界处。粒径 0.5mm的颗粒属于粒径较大颗粒，自身质量大，惯性大，流体对颗粒的裹挟能力差，因此叶片压力面侧磨损区域主要集中于进水边至中部。a.后盖板a.Hubb.叶片表面b.Blade surfaceCv=1%Cv=3%Cv=5%Cv=7%Cv=1%吸力面Suction surface压力面Pressure surfaceCv=5%吸力面Suction surface压力面Pressure surfaceCv=7%吸力面Suction surface压力面Pressure surfaceCv=3%吸力面Suction surface压力面Pressure surface0.0e+001.0e062.0e063.0e064.0e065.0e066.0e067.0e068.0e069.0e061.0e05磨损率 Erosion rate/(kgm2s1)0.0e+001.0e072.0e073.0e074.0e055.0e076.0e077.0e078.0e079.0e071.0e06磨损率 Erosion rate/(kgm2s1)图 12不同颗粒体积浓度下磨损分布Fig.12Erosiondistributionwithdifferentparticlevolumeconcentration当 Cv3%时，颗粒相体积浓度增大，流入叶顶间隙的颗粒增多，吸力面出水边靠近叶顶的点状磨损加剧；流道内颗粒之间的排挤严重，部分颗粒偏离原有的运行轨迹，跟随流体偏向吸力面侧运动，对吸力面侧进行冲击，撞击角较小，因此吸力面侧出水边靠近后盖板的位置也出现严重磨损，磨损区域呈现线状磨损且向叶片进水边方向延伸。由图 13 可知，颗粒相体积浓度的增大，叶片表面及后盖板磨损率均逐渐增大，但后盖板磨损率增速更快，当颗粒相体积浓度大于 3%时，后盖板平均磨损率超过叶片平均磨损率。1234567颗粒体积浓度 Particle volume concentration Cv/%1010686420平均磨损率Average erosion rate/(kgm2s1)叶片 Blade后盖板 Hub图 13不同颗粒体积浓度下叶片表面及后盖板的平均磨损率Fig.13Averagederosionrateofbladeandhubwithdifferentparticlevolumeconcentration4结论本文对半开式叶轮离心泵输送含沙水进行数值计算，分析了不同颗粒体积浓度条件下模型泵内泄漏涡特性和过流表面磨损特性的关联性，主要结论如下：1）设计流量工况下，半开式叶轮离心泵内压力由叶50农业工程学报（http:/www.tcsae.org）2023年轮进口向出口方向逐渐增大，相比于清水工况，固体颗粒的堆积造成流道内压力普遍升高，同时抑制了叶片边界层部分不稳定流动，叶片前半部分低压区面积减小。2）少量的颗粒可以抑制边界的流动分离，有利于泵的性能提升；当颗粒相体积浓度大于 1%，颗粒的堆积和频繁碰撞增加了泵的水力损失，泵的扬程和效率明显下降。3）颗粒相体积浓度为 1%时，固液两相工况下叶顶间隙泄漏涡、分离涡、通道涡结构与清水工况条件下相似；颗粒相体积浓度由 3%增至 7%时会造成泄漏流的破碎、分离、再融合，加剧流动的不稳定性。4）受到泄漏涡的阻碍，流入叶顶间隙的颗粒易绕过流道低速区偏向叶片吸力面运动，因此磨损严重区域主要分布于叶片前缘、叶片吸力面出水边、叶片压力面进水边及叶轮后盖板；体积浓度的增大会加剧叶轮整体磨损，颗粒具有明显向后盖板和叶片压力面的运动趋势，后盖板的平均磨损率上升速度更快。参考文献WANGL,LUJ,LIAOW,etal.Numericalinvestigationonvibration and pressure fluctuation characteristics in acentrifugalpumpunderlowflowrateJ.ProceedingsoftheInstitution of Mechanical Engineers,Part C:Journal ofMechanicalEngineeringScience,2021,235(19):3888-3899.1王李科.半开叶轮离心泵叶顶泄漏涡的流动特性及其抑制方法研究D.西安：西安理工大学，2020.WANG Like.Investigation of the Characteristics andSuppressionMethodofTipLeakageVortexinASemi-openCentrifugalPumpD.Xian:XianUniversityofTechnology,2020.(inChinesewithEnglishabstract)2司乔瑞，郭勇胜，田鼎，等.半开式叶轮离心泵气液两相条件下内部流动特性分析J.农业工程学报，2021，37(24)：30-37.SIQiaorui,GUOYongsheng,TIANDing,etal.Studyontheinner flow characteristics of an inside unshroud impellercentrifugal pump under gas-liquid two phase conditionJ.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2021,37(24):30-37.(inChinesewithEnglishabstract)3张华，陈斌，王炳祺，等.叶顶间隙大小对螺旋离心泵内部压力脉动的影响J.农业工程学报，2017，33(1)：84-89.ZHANGHua,CHENBin,WANGBingqi,etal.Influenceoftip clearance on internal pressure fluctuation of screwcentrifugalpumpJ.TransactionsoftheChineseSocietyofAgricultural Engineering(Transactions of the CSAE),2017,33(1):84-89.(inChinesewithEnglishabstract)4PARSI M,NAJMI K,NAJAFIFARD F,et al.Acomprehensivereviewofsolidparticleerosionmodelingfor5oil and gas wells and pipelines applicationsJ.Journal ofNatural