01 研究背景

什么是压气机的旋转失速问题？

旋转失速是沿压气机周向的非均匀流动状态。失速现象一般首先发生在叶轮处，当离心式或轴流式压缩机的操作工况发生变动时，气流会在叶片的凹面附近形成气流漩涡，气流漩涡的聚集会阻碍通道内的气流流通，减少通道内的有效流通面积，形成气流堵塞团，不但会使发动机性能（推力、经济性）大为恶化，限制发动机的工作范围，更严重的可能会引起发动机突然熄火，或引起压气机叶片剧烈振动以致叶片断裂而造成整台发动机的损坏。

失速现象可能对压缩机造成严重损害，因此预测在何种工况下会出现旋转失速的现象就显得十分有必要。在本次介绍的算例中，就使用法国电力（EDF）开发的通用计算流体力学（CFD）求解器 code_saturne进行数值模拟，并将仿真结果同实验测量数据以及FLUENT商业软件的仿真结果进行比较，证明了code_saturne进行可压缩流场仿真计算的可信性。

code_saturne中的可压缩算法和湍流模型的改进

压气机失速现象涉及到对旋转可压缩气流流场不稳定性的高精度的捕捉，针对本案例，对code_saturne内的可压缩流求解部分的源代码做出了以下的优化和改进。

Ø 对流项的离散格式增加至3阶

Ø 使用2阶Crank-Nicholson时间差分格式

Ø 植入并试验了由Čada和Torrilhon提出的3阶限制器

同时针对code_saturne中的k-ω SST湍流模型进行了旋转和曲率修正，使之更加适用于旋转气流仿真。

为了对改进后的k-ω SST湍流模型，在code_saturne中测试了旋转管流的经典验证算例，模拟绕轴向旋转的管道中的气流流场。

算例中流场对应的雷诺数（Re）为5800，测试不同罗斯贝数(Ro)下的模拟结果并同DNS仿真结果以及FLUENT商业流体力学软件的仿真结果进行比较，验证了改进后的k-ω SST湍流模型对于旋转流场湍流特征的捕捉精度。

02 应用案例

NASA Rotor 67多级压气机风扇转子

NASA Rotor 67 是可压缩流体的经典验证案例之一，整个结构由22个风扇叶片组成，其设计工况为转速16043 rpm，流量为33.25kg/s。本计算案例采用y+平均值为35的网格量为230万的网格进行计算。计算使用 S-A湍流模和k-ω SST模型分别进行计算，并使用1阶时间差分格式和2阶MUSCL格式。

使用code_saturne仿真计算所得出的结果同公开的实验测量数据以及使用Fluent商业计算软件的仿真结果进行了对比，code_saturne和FLUENT的计算结果十分接近，且都与实验所得数据基本吻合。

NASA CC3离心式压缩机的仿真

案例中模拟的NASA CC3离心式压缩机的设计运转工况为：

压缩机结构由15个叶片和15个分流叶片组成，整个压缩机分别使用平均y+值为30的2180万网格量的网格，和平均y+值为5的7000万网格量的两套网格进行了仿真模拟。模拟使用改进的 k-ω SST 湍流模型，1阶时间差分格式 和2阶 MUSCL格式。

使用code_saturne仿真计算所得出的数据与使用Fluent进行模拟的结果相比十分接近，同时也更为接近实验数据。

通过观察在轴面上的速度云图，也可以在叶轮与扩压器之间捕捉到旋转失速的现象，在叶轮和扩压器中发现了四个脱流点，叶轮中脱流点处的速度是叶轮中旋转速度的0.0587倍。还可以观察到压力扰动向上移动至叶轮的入口，同时高速气流团向下游的的扩散器移动。

03 结论

经过改进后的code_saturne算法采用了更高阶的离散格式，可以准确地预测到离心式压缩机中所产生的失速现象，仿真结果同Fluent模拟结果相比十分接近，且更符合实验数据。如果想要得到更加精确的仿真结果，可以采用更加精细的网格或者LES湍流模型。但是由于压缩机中气流流场的速度较大，需要的计算资源也会更多。

更多资讯可登录格物CAE官方网站

https://cae.yuansuan.cn/

或关注微信公众号【远算云学院】

远算在bilibili、头条、知乎、仿真秀定期发布课程视频等内容

敬请关注