CFDforcentrifugalpumps:areviewofthestate-of-the-artSRShah,SVJain,RNPatel,VJLakheraAbstract:离心泵内部流动的分析非常复杂,由于3维流动结构涉及到湍流,二次流,汽蚀和其他不稳定的因素。计算资源和进步的数值方法的准确性日益普及带来了从纯理论研究工作的涡轮机械计算流体动力学(CFD)进军竞争激烈的工业市场。本文是就利用CFD分析离心泵在未来发展的趋势所做的综述。CFD技术已被研究者应用于离心泵的研究中,这些研究包括在设计工况与非设计工况泵的性能的预测,泵的参数研究,汽蚀分析,扩散泵分析以及运转在涡轮泵工况下性能等。非定常的雷诺平均纳维斯托克斯方程和k两方程的湍流模型进行离心泵的CFD分析是适合的。蜗壳流动的研究以及蜗壳与叶轮之间的相互作用是一个有趣的研究领域,也是未来泵性能研究的方向,最活跃的研究和发展领域是两相流,泵在非牛顿流体下的运行和流固耦合。Keywords:离心泵;CFD;流动分析;数值模拟;综述1.Introduction离心泵是工业,农业和国内最常用的泵。它的叶片的设计需要详细的了解在部分负荷和额定负荷下叶轮内部的流动。设计经济高效泵,在制造之前预测泵的性能是至关重要的,这一预测是需要了解泵体内部不同部分之间的流动规律。模型试验是泵性能预测的一条途径,但它太单一,且耗时耗财。相反的,理论途径只是给出了一个值不能够决定引起泵特性差的根源。近年来,由于软件设计的改进,CFD在成功预测泵与涡轮机械方面发挥着重要的作用。CFD在水泵和涡轮机械设计中的应用始于30多年以前,这一起始恰逢将将有限元法应用于CFD,并以把全3维的流动求解简化为准3维的欧拉求解为特征。多年来经过了多个发展阶段:从3维的欧拉求解,到利用雷诺时均法对单个叶片流道利用有限容积法的模拟,再到整个机械的稳态模拟,直到今天利用先进的湍流模型求解非稳态雷诺时均方程。现在最为活跃且研究最多的是两相流和流固耦合。计算资源和进步的数值方法的准确性日益普及带来了从纯理论研究工作的涡轮机械计算流体动力学(CFD)进军竞争激烈的工业市场。市场上有很多的软件可以用于旋转机械的数值分析。Fluent(UKandUS),CFX(UKandCanada),Fidap(US),Polyflow(Belgium),Phoenix(UK),StarCD(UK),Flow3d(US),ESI/CFDRC(US),SCRYU(Japan)。最近,Kaupertetal.,PottsandNewtonandSunandTsukamoto分别利用CFX-TASCflow,FLUENT和STARCD来研究泵在非设计工况下了性能。通常情况下,CFD提供了提供了三种叶轮机械内部流动的计算方法:多参考系模型,混合平面和滑移网格。前两种方法是基本的定常流方法。在多参考系方法中,旋转轴被保持在一个固定的位置以及用于在旋转参考系旋转轴的控制方程求解,包括科里奥利力和离心力,静止部分是在绝对参考系中求解的。Takamara、Goto、Goto等人已广泛利用CFD分析低比转速和中比转速泵;Sedlar和Mensik利用CFD分析径向流泵。本文中,进行的是针对于CFD在离心泵的应用的综述。阐述CFD的优点和局限性,讨论CFD作为数值模拟的工具分析离心泵内部流动的当下趋势和长远发展。2.Literaturereview泵的设计者不断受到挑战,使得机械运行的更高效,平静,更经济。许多研究者将CFD作为数值模拟的工具来研究离心泵。本节介绍了研究人员通过CFD方法针对离心泵所进行的研究工作。2.1.Performancepredictionatdifferentoperatingconditions离心泵被广泛应用,因此泵系统需要在不同的应用场合有很宽的流动变化范围。最早的数值模拟聚焦于设计或接近设计工况的泵。很少会去研究非设计工况泵的性能,此时泵的性能会下降。在CFD方法使用之下,可以研究在不同的运行条件下通过泵体内部不同部分的复杂流动,这有益于改进泵在非设计工况下的性能。Mentzoset进行了关于后弯式叶片离心泵的数值模拟。使用MRF方法考虑叶轮和蜗壳的相互作用是不可行的,由于其固定的耦合方程。瞬态分析被看做是一个真正能够分析叶轮和蜗壳之间相互作用的工具。3维离心泵计算模型如下图(a)所示。(a)离心泵的3维计算模型(b)泵的静压云图Mentzos利用有限体积法,采用结构画网格来求解离散控制方程,进行了离心泵叶轮内部流动的模拟。CFD技术被用来预测流动参数,压力分布和水头-流量曲线。有文章称,尽管网格的尺寸对于局部边界层变量的研究不适合,但整体的信息却能很好地扑捉到。Shah利用雷诺时均方程进行了离心泵在200m3/hr流量下的稳态数值模拟。在非设计工况下,泵的各部分流场是不均匀的,这导致了泵的效率的下降。与RNGk湍流模型相比,kSST湍流模型可以得到更好的结果。通过数值模拟预测的工作特性曲线与模型试验的结果进行了对比,发现吻合良好。额定流量下泵的静压云图如图1(b)所示。2.2.ParametricstudyCFD可以在水利机械制造之前预测各个部分的流动性能,避免修改现有系统,可直接对计算模型进行修改,在制造之前可预测结果。CFD分析有助于研究各种参数以及通过形成独立的无量纲组对泵的性能的影响。Bacharoudis分析了保持出口叶片直径相同情况下,改变叶片出口的角度对泵性能的影响。进行了三维的数值模拟不可压缩的纳维斯托克斯方程是利用一个商业的CFD有限体积法程序求解的。在额定流量下,当叶片出口角从20度增加到50度,水头增加超过6%,但是效率了4.5%。但是,当流量很大的时候,叶片出口角度的增加会使得水力效率显著上升。为求解雷诺方程,Anagnostopoulos发展了模拟离心泵叶轮3维湍流的数值模型。叶轮的几何形状是由多个可控的设计变量表示,方便修改叶轮形状和测试不同的配置。参数研究结果表明,优化叶轮几何形状可以显著增加叶轮的水力效率。修改叶轮形状和其对于水头及效率的影响如图2(a)所示。(a)(b)PatelandRamakrishnan研究改变毂轮廓曲线和定子角在混流泵在工作点和在部分负荷对泵性能的影响。研究的结果如下:(1)水头与功率曲线和标准工况下混流泵的相似;(2)泵的效率在额定工况下效率的+5%范围内变化(3)通过对比非设计工况下的变量,在匹配了改变定子角度和轮毂曲线后,效率增加了1%。实测和CFD分析的泵的效率和功率曲线如图2(b)所示。2.3.Cavitationanalysis汽蚀发生在泵内压力低于对应温度下汽化压力的地方。汽蚀破坏的机理已经研究了100多年,但至今仍然没有统一的汽蚀破坏理论来描述泵叶轮汽蚀速率和预测泵在设计阶段的汽蚀性能。Medvitz采用多相流的CFD模型分析离心泵汽蚀时的性能。混合物的动量方程和连续性方程中所用到的两相雷诺时均方程和气泡体积分数方程耦合求解。叶片的汽蚀,甚至是汽蚀断裂,都是与实验观测进行定性分析。Nohmiet利用可压缩空气-蒸汽-液体两相介质(TE模型)和恒焓汽化(CEV)模型研究了低比转速离心泵的汽蚀。研究表明,在大流量下,气泡出现在压力面的进口边,水头下降很缓慢。TE模型能够预测水头的缓慢下降,但是计算发现是不稳定的,CEV模型不能预测水头的缓慢变化。在这两种模型中,需作进一步的修改才能获得稳定而又准确的结果。Caridad进行了离心泵叶轮输送汽水混合物的数值模拟,这与有汽蚀的流动很相似。进行了气泡间隙和气泡直径之间的相关分析。叶轮出口水头和相对出口流动角的变化被描述为液相体积流速和气相体积分数的函数。研究发现,更大直径的气泡更容易导致叶轮出口水头的波动。数值模拟的结果与试验数值吻合良好。2.4.Investigationsoninteractingcomponenets叶轮和蜗壳之间的相互移动,不仅会影响泵的整体性能,而且还会引起压力波动。与蜗壳相互作用产生的压力波动会对机械部件产生动态作用(主要是非定常力),这是产生震动和噪音的一个很重要的原因。Gonzalez利用数值模拟得出了离心泵内由于叶轮和蜗壳相互作用而产生的动力和非恒定流的影响。运用粘性的纳维斯托克斯方程和滑移网格技术求解叶轮与蜗壳之间的相互作用。蜗壳壁面的压力波动如图3(a)所示。WangandTsukamoto发展的数值方法能够很好的预测考虑扩散泵在改变运行工况点时由于定子和转子之间的相互作用产生的压力波动。图4显示了在失速情况下叶轮叶片流道1和2附近的速度矢量。运用涡分析法求解了2维非定常不可压缩流的压力波动,此时壁面产生的涡是基于动量方程。2.5.Axialthrustincentrifugalpumps卧式泵很有挑战性的研究是估算受泵的特性,工作条件和内部压力场影响的轴向推力。求解单级泵这一问题很简单,但对于多级泵,却是很复杂的。Salvadori进行了多级卧式泵的数值模拟并通过研究内部流动和压力场给出了各个部分对于轴向的负荷。该方法简易为了确保多级泵的工作靠靠性,可以增加合适尺寸的推力轴承。更重要的是做了漏流罩室和平衡鼓泵轴向力平衡分析。所提出的技术是考虑泵了泵的工作和磨损条件对轴向推力的影响。Gatta提出了一种技术可以预测多级泵的参与轴向力,讨论了不同变量对于轴向推力的影响。混合盘面技术可以用来研究静止部分和转动部分之间的相互作用,这一技术考虑了不同的边界条件和进口速度分布。认为侧室对于轴向力的分布有很大的影响,压力分布主要受空腔内部泄漏质量流量和相对转速影响。2.6.Non-newtonianfluidhandlingpumps泵广泛应用于采矿业,化工,冶金和煤炭行业进行泥浆的输送。典型的浆液,固体颗粒具有不同的直径和浓度,这取决于浆料的类型和具有特定应用。固体颗粒对于泵的水利和汽蚀性能有很大影响。Pagalthivarthi利用FLUENT6.1中的欧拉多相流模型模拟了浓浆流过离心泵的情况。采用一阶迎风格式离散动量方程。叶片曲率,套管宽度,颗粒的入口浓度对于壁面速度和壁面的压力分布有影响。分析的结论是,固体浓度和固体壁剪切应力从叶片的上游到下游的腹部区域单调增加。Yu从实验和数值计算两个方面研究了叶轮几何形状对于离心式血泵性能的影响。在设计工况点测试了4个叶轮,利用血液作为流体介质。结果发现,叶片通道内的应力水平一般在发生广泛的红细胞损伤阈值150N/m2水平以下。对靠近叶片的前缘部分的局部区域进行观察,其应力水平为60%以上的阈值水平。2.7.Pumprunninginturbinemode一个小泵可在微水电厂运行在涡轮机械模式下,具有很多优势。许多研究者将CFD分析泵称为涡轮泵。NatanasabapathiandKshirsagar利用数值模拟的方法研究涡轮泵的性能。利用PRO/E创建3维的计算模型,然后导入到CFX中计算如图3(b)所示。最初,模拟是用非结构化网格做,水头下降和实验值相差较大。(a)(b)Fernandez利用滑移网格进行了涡轮泵的3维数值模拟。数值计算的结果与实验数据如图5(a)所示。获得的总水头相对误差为9%,叶轮边缘平均静压的相对误差为5%。Barrio利用商业软件进行了离心泵工作于泵和离心泵状态下的内部流动,使用非恒定雷诺时均方程和有限容积法求解的纳维斯托克斯方程进行数值模拟,数值模拟的结果和实验结果误差在3-5%内。2.8.Mini/micropumpanalysis由于旋转泵是液体输送的关键机械,其进一步发展是很有必要的。叶轮直径在5到50毫米之间的泵被定义为微型泵。Liu利用实验和数值模拟对微型泵叶轮几何形状进行了研究。观察在雷诺数大于1.0×105的范围内的泵的相似特性的规律。由叶轮几何尺寸的缩减带来的叶尖间隙的减小对于泵的影响和较大的叶片出口角对泵性能的影响一样。结果基于雷诺时均方程和k湍流模型的3维数值模拟来研究微型泵是可以的。TsuiandLu利用集总系统和CFD的方法分析了微泵中阀的非恒定流场。移动膜