旋涡风机内部流场特性及涡分布特征研究.pdf

书书书设 计 研 究 年 期（总第 期）旋涡风机内部流场特性及涡分布特征研究陈少华，张弘，谢军龙（华中科技大学能源与动力工程学院，湖北 武汉 ）摘要：旋涡风机是一种特殊回转式气体压缩机械，有效补充了容积式和离心式压缩机械的小流量和高压比应用间各自不足。风机内部涡流动机理对于其整体性能和可靠性具有至关重要影响。建立了旋涡风机数值模型，通过实验验证了模拟可靠性。通过旋涡风机速度流线分析，发现叶片边缘处流动特性较为复杂。通过分析叶片边缘的湍动能和静压分布曲线，揭示了旋涡风机内部气体流动规律。利用 准则对旋涡风机叶轮流域不同叶高处的涡分布情况进行定性分析。研究结果表明：叶片边缘处涡的分布直接关系到风机性能。关键词：旋涡风机；流动机理；叶片边缘；涡识别中图分类号：文献标志码：文章编号：（）收稿日期：，（，）：，：；引言旋涡风机主要用于气体增压、输送和循环，其独特之处在于能够在小流量下实现高压比。旋涡风机兼具容积式机器的可控性和稳定性，和离心式机器低振动、高可靠性、无需润滑和机械部件较少等典型优点 。由于旋涡风机工作原理和几何尺寸复杂，其内部流场呈现高度湍流 状态，导致流体传设 计 研 究 年 期（总第 期）输中损失能量较多，其多变效率无法达到 以上。尽管与其他压缩机相比，旋涡风机效率较低，但在旋转速度相同的情况下，它具有在较大流量下（相对容积式）的高压比（相对离心式）特性，且具有制造成本相对低、简单易用、自吸能力强等优点，在许多应用中可以提供很好的替代方案 。世纪 年代至 年代，部分学者基于湍流理论 和动量交换理论 解释了旋涡风机工作原理。通过实验和理论研究发现，动量交换理论在描述风机内部流场方面更为合适 。根据这个理论，叶轮中的气流受离心力的作用，在周向形成一系列螺旋状的运动 。随着叶轮旋转，壳体内部流体受叶片剪切力作用形成周向流，叶轮内部流体沿叶片及径向方向做旋转运动而获得动能。流体随后进入壳体流域并将动能传递给其中的周向流。气流受压差作用速度得到增加，重新进入到叶片根部。经过多次循环，气体在风机内以周向螺旋运动不断实现压力的增加。然而，动量交换理论并没有对整个场进行综合描述，意味着我们在研究整个风机时需要考虑更多的因素。旋涡风机的内部流场受到众多因素的影响。通过响应面分析研究得出，叶片出口角对压力系数的影响要大于进口角。唐照付 研究得出叶片侧边的形状影响风机的性能 。研究了轮毂形状对由叶轮旋转产生的径向流的引导作用。研究指出叶轮区叶片边缘处的压力分布增加明显，其中沿叶轮旋转方向中间位置的压力增加最为明显。叶片边缘处流动特性直接影响着气流动能转化为静压能的效率，进而影响到真空风机的性能和使用寿命。综合上述，叶片边缘位置的流动特性对风机的性能有重要影响，但其起作用的内流机理尚未被充分解释，需要对风机叶片边缘的气流流动特性开展深入研究。本文采用 建模分析方法对一款旋涡风机原型机模型验证并通过外特性试验验证对比，在此基础上对风机进行优化设计。在确认模型的可靠性后，对整体进行了流线分析，以确定风机整体流动状态。其中，叶片边缘处气流的运动状态和涡旋结构具有重要的影响，因此需要进行详细的流动特性分析。最后，在典型位置的涡识别分析中，主要关注涡旋的形态、强度和位置等参数，并借助准则等流场分析方法进行定性判断。通过对涡旋的准确识别和分析，能够更好地揭示气流的运动规律和机制，为更深入的流场分析和优化设计提供理论依据。数值方法 几何模型和边界条件旋涡风机主要由壳体、叶轮、进气口、出气口等部件组成。原型机的额定流量为 ，转速为 ，额定流量下的全压升可达 。在容积式压缩机中，间隙尺寸是影响泄漏量的重要因素 。为探究轴向间隙对本文旋涡风机流动性能的影响，在 计算中设置了 个轴向间隙（、），模拟发现轴向间隙为 时，泄漏量对压力曲线影响在可考虑范围内，这也与 的研究结果一致。旋涡风机的计算域划分如图 所示。为了方便网格划分和对各不同区域设置计算参数，将整个流体区域分成叶轮旋转区、壳体静止区和间隙区三部分。求解方法本文的数值模拟采用了有限体积的离散方法，内部湍流模型采用标准方程。通过对质量、能量和动量的交换来模拟流体在系统中的运动状态。设置了速度入口和静压出口边界条件，并使用普通连接模型和冻结转子算法进行坐标变换。为保证计算结果的准确性，在考虑数值模拟的收敛情况时，要求残差收敛曲线均小于 ，并限制监测的物理量波动范围在 以内。这样可以有效地提高数值模拟的计算精度和可靠性。网格无关性验证将旋涡风机计算区域分为旋转域与静止域，两者通过 连接，网格连接方式为 。旋涡图 风机计算域设 计 研 究 年 期（总第 期）风机网格划分如图 所示。为保证计算精度，对叶轮域和间隙域进行网格加密处理。原型机网格无关性验证如图 所示。采用了四组网格数据，网格数量分别为 万、万、万、万。设计流量下，不同网格数量的风机全压监测结果如图 所示。当网格数量从 万增加到 万时，风机全压从 变化到 ，压力波动小于 。综合考虑数值计算结果准确性与经济性，采用 万网格进行后续计算。试验验证旋涡风机原型机试验装置示意图如图 所示。风机进出口接风管到室外，进风口管道采用蝶阀控制流量，在进出口管上设置智能压力风量风速测试仪（系列，测量范围为 ，精度为 ）用于测量压力和流量。试验测试时，利用进口风阀调整流量至额定工况，采用测试仪监测风机进出口压力。为了验证数值计算的准确性，将试验数据与数值计 算 数 据 进 行 对 比，如 图 所 示。流 量 从 变化到 的过程中，全压模拟值及试验值的变化趋势相近，最大偏差不超过，验图 网格划分图图 网格无关性验证图 旋涡风机试验装置示意图图 试验验证证了所采用的数值模型和计算方法是可信的。结果与讨论 流线分析旋涡风机在额定工况下的流线图如图 所示。为更好地描述旋涡风机内部流动状态，在主视图分别取截面 、截面 来分析进口管、出口管流场，在俯视图取截面 、截面 、截面 、截面 个径向截面分析流动过程中壳体区和叶轮区的流场。观察图（）截面流线图，进口管截面 的流线分布相对均匀，没有明显的涡旋出现；而出口管截面 的流线却形成了两个反向旋转的涡流，涡区域内流体的速度较慢，压力较高，容易引起阻力和能量损失。个径向截面均呈现纵向旋涡特性，上下 个主涡附近均出现了一个副涡，位于叶片（出口）边缘附近，其对称流动性较差，导致流动损失增加。因此需要对叶轮区叶片（出口）边缘位置处的流动特性进一步分设 计 研 究 年 期（总第 期）图 旋涡风机总体和截面流线图析，探究其流动损失的主要原因。叶片边缘处流动特性叶轮区叶片边缘位置在整个风机中呈现为三维圆柱面，无法直观分析其流场特征，因此将其展开成二维平面，以便后续分析。沿圆柱面展开的叶片边缘湍动能耗散云图如图所示，可以看出湍动能耗散主要出现在 区域、叶片附近和出口管附近。（）区域：当流体通过 区域时，流体将发生强烈的扭曲、旋转和变形，并产生较大的速度梯度和剪切力，导致局部的湍流强度较高。（）出口管附近：图 沿圆柱面展开的叶片边缘湍动能耗散云图气流流经出口管区域时发生急转，高速气流急速冲击出口管的内侧（图中管的右侧），进而形成一个相当大的旋涡，阻碍了大部分流道，使主流区速度显著增大，产生了较大的湍流能量损失。（）叶片附近：叶片附近的湍动能耗散主要受到空气流动状态和叶片形状的影响。叶片吸力面附近的湍动能耗散相对较小，而叶片压力面（顶部和底部）附近的湍动能耗散则相对较大。因此，在设计风机时需要考虑叶片形状和气流流动状态对湍动能耗散的影响，以降低湍动能损失，提高风机性能。不同叶高位置处沿圆柱面展开的叶片边缘湍动能分布曲线如图 所示。由图可知：（）当旋涡风机叶片开始旋转时，气流开始在叶片表面形成作用力。由于叶片的阻力，气流逐渐在边缘处集聚，进而产生旋涡。旋转的涡流会对叶片边缘产生强烈的扰动，导致湍动能增加。（）叶高位置的湍动能大小明显高于 叶高和 叶高。原因是壳体底部流道截面面积较大，叶轮流道和壳体流道中的流体所受离心力的差值也较大，导致纵向旋涡增强，湍动能增大。沿圆柱面展开的叶片边缘静压云图如图 所示。受进口负压影响，气流沿进气口被吸入壳体。随叶轮旋转，气流静压逐渐升高。气流沿叶片表面流动时，叶片的吸力面的气流流速加快，静压降低。而叶片的压力面的气流流速减慢，静压升高。吸力面与压力面产生的静压差将推动气流从叶片底部向图 不同叶高处沿圆柱面展开叶片边缘湍动能分布曲线图 沿圆柱面展开的叶片边缘静压云图设 计 研 究 年 期（总第 期）顶部流动。不同叶高位置处沿圆柱面展开的叶片边缘静压分布曲线如图 所示。由图可知，叶高处的静压分布显著高于 与 叶高处的静压。原因是在叶轮轮毂形状（中部凸起）的影响下，气流在叶片间的旋涡大多分布在叶片底端和叶片顶端附近，与 叶高处的湍动能高，而叶片中部流动湍动能低。对比叶片边缘处湍动能耗散分布和静压分布可知，湍动能对旋涡风机流场分布影响较大，而湍动能又与涡分布息息相关，为更好描述风机内部流场，需要进一步分析其涡分布特征。通过对叶轮区涡分布特征的详细分析，进一步研究风机内部流动特性。叶轮区涡分布特征了解涡识别准则需要从速度梯度张量入手，而速度梯度张量可以分解成两部分速度张量可以分解成两部分珗（珗 珗）（珗 珗）（）其中，对称矩阵记作 ，指变形引起的应变速率张量；反对称矩阵记作 ，指与整个流场涡度相关的旋转速率或涡量张量。引入矩阵 和矩阵 的 范数平方，如式（）（）所示图 沿圆柱面展开的叶片边缘静压分布曲线（不同叶高）（）（）（）（）（）（）值定义为速度梯度张量的第二不变量，如式（）所示。表示涡存在 。（）（）准则下的旋涡风机涡量分布如 所示。为进一步分析旋涡风机涡分布，采用 准则对叶轮流域不同叶高处的涡分布云图定性分析。通过计算流场的旋度与速度梯度的比值，结合 准则来判断涡旋的存在和强度。当 时，表示该处存在着明显的涡旋结构。使用无量纲形式来计算 的值 ，如式（）（）所示 （）（）（）（）（）准则下旋涡风机叶轮流域不同叶高处的涡分布云图如图 所示。旋涡风机的叶轮内部涡主要分布在叶片边缘位置。叶高位置处，气流运动空间较小，气流的运动状态较为平稳，缺乏充分的能量聚集和扰动，导致涡旋生成的几率较小；图 准则下旋涡风机涡识别（速度分布）图 准则下旋涡风机叶轮流域不同叶高处的涡分布设 计 研 究 年 期（总第 期）叶高和 叶高位置，流域横截面积较大，气流速度和压力变化较大，容易形成不稳定的气流扰动，从而促进了涡旋的形成和发展。结语本文分析了旋涡风机内部流动特性和涡分布特征，主要结论如下：（）旋涡风机流场的径向截面呈现纵向旋涡特性，在叶片边缘位置会出现不规则涡分布，导致流动损失增加。（）叶 片 边 缘 的 湍 动 能 耗 散 主 要 出 现 在 区域、叶片附近和出口管附近。另外，受流道截面面积影响，叶高位置的湍动能较高处纵向旋涡增强，湍动能较大。（）利用准则定性分析旋涡风机叶轮流域涡分布情况，在 叶高位置，漩涡发展空间受限，气流的运动状态平稳；而在 叶高和 叶高位置，流域截面积增大，涡分布明显增加。参考文献：，：，：，（）唐照付 叶片侧边型线对旋涡风机性能影响的研究 华东理工大学，唐照付，聂波，张俊林，等 旋涡风机叶片侧边型线的研究 流体机械，（）：，（）：，（）：李敏，王亚洲，刁安娜，等 螺杆压缩机内密封齿间隙泄漏量的数值仿真 压缩机技术，（）：，：，（），（）作者简介：陈少华（），女，硕士研究生，现就读于华中科技大学能源与动力工程学院。：檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪殏殏殏殏 欢迎订阅 年 压缩机技术杂志