两个典型离心式压缩机级静子元件中流动的数值模拟

1、两个典型离心式压缩机级静子元件中流动的数值模拟 周晓东 徐 忠/西安交通大学流体机械研究所 孙玉山/沈阳鼓风机（集团）有限公司 摘要：使用NUMACA 软件对两个不同流量系数的离心式压缩机级静子元件中的流动进行了数值模拟。分析了两个级无叶扩压器中的流场及其中一个级弯道和回流器中的流场。对离心式压缩机级设计提供了改进信息。关键词：离心式压缩机 静子元件 数值模拟中图分类号：TH452 文献标识码：B 文章编号：1006-8155(2006)02-0001-05Numerical Simulation of Flow between the Static Components of 2 Stage

2、s of Typical Centrifugal CompressorAbstract: Numerical simulation of flow is carried out on static components of 2 stages of centrifugal compressor whose flow coefficient are different using NUMACA software. It analyzes the Flow field of vaneless diffuser of 2 stages and the bend and return channel

3、for 1 stage. The information how to improve the stage design of centrifugal compressor stage is put forward.Key words: Centrifugal compressor Static component Numerical simulation 1 引 言 离心式压缩机的静子元件通常有吸气室、扩压器、弯道、回流器及蜗壳（或排气室）等。其作用主要是把气体从前一元件引到后一元件中去，并使之具有一定的速度和方向。静子元件对整个压缩机的工作效率有相当大的影响，随着人们对叶轮研究的不断深入，

4、叶轮的效率已经达到了很高的水平。相对而言，人们对静子元件的研究还很不深入，因此近年来人们对静子元件的研究越来越重视。本文主要涉及的是扩压器、弯道和回流器这3个静子元件。一般情况下，离开离心式压缩机级叶轮的气流动能为总输入功的3040。这些动能需要在扩压器中进行回收。在扩压器中气流流道面积增大使平均速度减小，并因此使动能减小，静压增大，同时伴随着能量损失，气流分离和失速。在弯道和回流器中亦同样会有能量损失和气流分离。无叶扩压器具有结构简单，性能曲线平坦，稳定工况范围宽等特点。因此，在多级离心式压缩机中得到广泛应用。本文就两个带有无叶扩压器（其中一个为末级，另一个为中间级带有弯道和回流器）的不同流

5、量系数的离心式压缩机级进行数值模拟。分析了两个级无叶扩压器中的流场及一个级弯道和回流器中的流场。以其对离心式压缩机级提供改进信息，从而提高级的设计水平。2 数值模拟方法 采用NUMECA软件的FINE/TURBO模块进行数值模拟。该软件采用时间相关法求解雷诺平均NS方程，中心节点的有限体积离散，显式龙格库塔时间积分，全多重网格初场处理，以及多重网格迭代加速。计算中采用SpalartAllmaras一方程湍流模型，对流项离散格式采用中心差分并结合四阶耗散项，时间推进选用四阶显示Runge-Kutta方法，CFL数取值2。分析了流量系数为10.08和0.16的两个离心式压缩机无叶扩压器级：（1）1

6、0.08级（以下称为级A），D2450mm，2A45°，D4720mm，回流器叶片进口角-61°，由叶轮、扩压器、弯道和回流器组成，计算的进口位置在叶轮的进口处，出口位于回流器的出口；（2）10.16的级（以下称为级B），D2240mm，2A65°，长短叶片，D4384mm，由叶轮、扩压器和蜗壳组成，计算的进口位置在叶轮的进口处，出口位于扩压器的出口。为计算方便，两个级均在进、出口位置增加了延伸段。计算网格：级A的网格分为两大部分，第一部分为叶轮和扩压器，I型网格，分为两块，第一块为叶片流道（41×41×193），第二块为叶片尾缘钝头下游区域（

7、17×41×41）；第二部分为弯道和回流器 (41×41×193)，H型网格，总共约62万个网格。级B的网格分为8块，第一、第二块为主流道，I型网格，网格数均为25×49×193，长叶片顶部间隙采用蝶型网格（分为两块，网格数分别为9×9×193和17×9×89），短叶片顶部间隙采用蝶型网格（分为两块，网格数分别为9×9×193和17×9×89），长、短叶片尾缘钝头下游区域分为两块（网格数均为17×49×33），总共约59万个网格。图1和

8、图2 分别为两个级的网格图。边界条件：进口给定总温、总压和气流角（轴向进气），出口给定质量流量，叶片及内外环壁为无滑移、绝热壁面。收敛准则：进出口流量误差小于0.5，残差水平下降到104以下。 图1 级A网格图 图2 级B网格图3 结果分析3.1 级特性参数分析 图3给出了级A特性参数曲线，从图中看出，该级的工况范围为92.5125；图4给出了级B特性参数曲线，从图中可以看出，该级的工况范围为80%130%。级B有比较理想的工况范围，而级A的设计点偏向小流量区。图3 级A多变效率及压力曲线图4 级B多变效率及压力曲线3.2 设计工况下级内部三维粘性流场分析 图5和图6给出了两个级的子午平均速度

9、分布，从图中看出级A无明显低速区，而级B在扩压器内靠近盖侧发生分离并产生倒流，形成较大的漩涡，气流大部分被压缩到盘侧，使盘侧气流形成较大的速度。 图5 级A子午平均速度 图6 级B子午平均速度（右图为扩压器分离区放大图）图7图12给出了两个级的横向平面速度分布，从图中看出，级A在回流器叶片吸力面有一个分离区，该区域在轴盘侧较小，位于叶片的出口，向盖盘侧逐渐增大，在叶高约50的位置可以明显看到较大的分离，该分离区在近盖侧已从叶片出口处向进口延伸了约2/3的叶片长度；级B在扩压器盖侧有一个低速区。A图7 级A 25%叶高位置横向平面速度（右图为A处局部放大图）B图8 级A 50%叶高位置横向平面速

10、度（右图为B处局部放大图B）C图9 级A 75%叶高位置横向平面速度（右图为C处局部放大图）图10 级B 25%叶高位置横向平面速度 图11 级B 50%叶高位置横向平面速度D 图12 级B 75%叶高位置横向平面速度（右图为D处局部放大图）图13和图14给出了两个级的子午平均静压分布，图中显示出静压梯度在两个级的扩压器中分布都比较均匀，级A在弯道和回流器中压力变化较小。从图15和图16的子午平均绝对总压分布图中看出级A扩压器出口处的轴盘侧绝对总压较低，级B在扩压器进口处的盖盘侧绝对总压较低。从图17和图18熵的子午平均分布图中可以看出相应绝对总压较低的位置熵值较大，表明该处损失较大。 图13

11、 级A子午平均静压分布图 图14 级B子午平均静压分布图 图15 级A子午平均绝对总压分布图 图16 级B子午平均绝对总压分布图 图17 级A子午平均熵分布图 图18 级B子午平均熵分布图图19和图20给出了两个级在50叶高处的静压分布。两个级在叶轮部分压力梯度分布比较均匀。级A在无叶扩压器里的压力梯度仍较均匀，在弯道、回流器中压力变化较小；而级B扩压器中的静压在前半部分梯度较小，这是由于此处存在较大的分离区而使通道有效面积减小、扩压效果减弱的缘故，后半部分压力梯度较大且比较均匀。图19 级A 50叶高位置横向平面静压分布 图20 级B 50叶高位置横向平面静压分布图21和图22给出了两个级的

12、扩压器子午平均绝对气流角分布，从图中可以明显看到两个扩压器中气流角分布有明显的差别。级A扩压器盖侧气流角在进口处沿扩压器宽度的分布比较平缓，轴盘侧的气流角略大于盖盘侧，流道中部气流角随直径的增大而逐渐增大；级B的气流角沿流道宽度变化较大。扩压器气流角在靠近盖盘侧20宽度附近气流角最小，向轴盘方向气流角逐渐增大，流道中部气流角变化不大。1.2D21.4D21.05D21.6D21.6D21.4D21.2D21.05D2 图21 级A扩压器子午平均绝对气流角 图22 级B扩压器子午平均绝对气流角4 结论本文对两个流量系数的离心式压缩机级进行了数值模拟，从中看出，两个级都有较高的效率，内部流场总体上都比较好，但仍存在以下问题需要加以改进。（1）级A叶轮气流转弯较均匀，且b2/D2较小，因而扩压器中流场的分布较好；而级B由于叶轮盖侧气流转弯较快，同时由于b2/D2较大，因而在扩压器盖盘侧产生分离引起较大的损失。可考虑改进扩压器设计，使b3<b2，加速此处的气流速度，同时将叶轮出口到无叶扩压器进口之间做成圆滑过渡的收缩型，以使气流转弯均匀减少扩压器盖盘侧的分离。（2）在扩压器中，级A的气流角沿扩压器宽度较为平缓，而级B的气流角在盖盘侧和轴盘侧差别很大；在扩压器盖侧级A的气流角较大，而级B的气流角较小，这与级B在此处有大的分离有关；叶道中部，级A的气流