利用STAR-CCM+对压气机叶型进行优化

       压气机作为发动机的重要部件，对发动机性能有重要影响；在高推重比涡轮风扇发动机中，高负荷跨声速压气机内部的流动情况复杂多变，叶型损失和二次流损失大大增加了压气机的设计难度；研究压气机三维叶片和流道对激波和复杂二次流动的影响机理及提升相应的控制技术，使得压气机优化设计成为高性能压气机设计过程中至关重要的环节。

   数值方法

研究模型

研究对象为跨声速压气机转子NASA Rotor 37转子， Rotor 37是NASA Glenn Research Center于 20 世纪 70 年代设计的四个高压压气机进口级之一，后来NASA采用激光测速仪和探针对孤立转子流场进行了详细的测量，如图。

rotor 37 参数

参数	值
叶片数	36
设计转速(rpm)	17188.7
设计流量(kg/s)	20.19
设计压比	2.106
叶尖速度(m/s)	454.14
设计转速下的叶尖间隙高度(cm)	0.0356
叶尖入口相对马赫数	1.48
叶根入口相对马赫数	1.13
叶尖稠度	1.29
展弦比(平均叶高/叶根轴向弦长)	1.56
转子进口轮毂比	0.7

Rotor 37测量位置

Rotor 37转子计算模型和网格如图，边界条件与文献相同，进口根据文献给定了径向分布的总温、总压，出口给定静压，叶片两侧给定周期性边界条件。对模型划分了网格，并进行了无关性验证，如图3，网格数量大于108万时满足无关性要求，通过划分网格和计算求解重复来调整Y Plus，满足无关性要求的Y Plus平均值约为0.74。

rotor 37计算域

网格无关性验证

数值方法及验证

       叶片流场的计算采用有限体积法求解粘性N-S方程的方法，采用二阶迎风格式对方程的对流项进行离散，采用中心差分格式对扩散项进行离散，离散方程的求解过程为连续方程和动量方程同时求解，湍流模型为SST模型，采用多重网格法加速收敛。

数值方法验证

   对Rotor 37转子进行了数值计算，计算的堵塞流量为 20.97kg/s，试验值为 20.93±0.14 kg/s，计算结果在试验误差范围内。图给出了转子特性曲线，可以看出数值计算得到的压比和效率曲线与试验结果的趋势基本一致。

Rotor 37压比特性

Rotor 37效率特性

   在98%堵塞流量下，对比了效率、压比、温比和气流角的出口节距平均值的计算结果与试验结果的差别，可以看出计算结果与试验结果基本一致。

Rotor 37出口压比分布

Rotor 37出口效率分布

Rotor 37出口温比分布

Rotor 37出口气流角分布

       通过对比Rotor 37 转子的计算结果与试验结果，计算结果和试验结果符合较好，该计算方法可以用来进行压气机优化计算。

优化

       在经过数值方法验证的98%堵塞流量工况下，对Rotor 37转子叶片进行了优化。转子叶片优化先进行初始流场求解，获得流场之后确定优化目标函数，为获得寻优方向，需求解优化目标函数对叶片的微分，即敏感性，计算方法采用伴随法，根据获得的敏感性利用自定义函数通过控制点控制网格变形改变叶片型面，最后求解新叶型的流动控制方程，计算出优化后的流场，根据目标函数是否满足预期要求再决定是否继续优化，叶片优化流程如图。

优化流程图

优化目标及设计变量

   叶片优化的目的是降低流动损失，例如激波损失、二次流损失等，因此将效率作为优化目标，为保证叶片维持设计的性能参数，需对压比和流量比进行约束，通过引入罚函数将压比和流量比约束引入到目标函数中，目标函数定义为：

 

   优化变量为控制网格变形的控制点。

敏感性分析

   在获得三维流场之后；求解伴随方程得到目标函数对型面的敏感性，如图，可以看出在叶片吸力面中部敏感性最大，提取10%、50%、90%三个截面敏感性矢量，如图，箭头方向代表使目标函数增大的叶片型面变形方向。  

a）吸力面              b）压力面

面的敏感性

敏感性矢量

网格变形

   压气机叶片是基于有限体积法对流场进行求解，网格的离散使梯度变化通常是不连续的，如果直接对叶片网格节点进行调整会导致型面变形量过大，不光滑，因此需要设置控制点，通过移动控制点的方式，使相应的网格节点的空间坐标发生变化，从而实现网格的变形；图为设置在叶片外的控制点，控制点均匀排布，用以控制叶片网格变形。

   根据敏感性计算控制点的位移量，控制点的移动量为控制点的累积位移量加上控制点上计算的网格敏感性得到的控制点所有方向上目标函数变化量的最大值乘以系数，系数的值根据经验取得，系数取得越小，所需的网格变形次数越多，系数取得越大，网格变形次数减少，但有可能变形过大，超过最优值。

控制点

约束条件

   叶片网格变形中需给定约束，一是叶片前缘和尾缘直径不能变，二是叶顶间隙不能变。为保证前后缘直径不变，在网格变形过程中，前后缘表面控制点的位移量通过自定义函数先进行周向平均然后赋给相应半径范围内的所有点，如图；为保证叶顶间隙不变，在变形过程中叶顶控制点的径向位移设置为0，切向和轴向随叶片变化而自由移动，如图。

a）前缘控制点         b）叶顶控制点

约束控制点

结果分析

性能分析

       优化前后设计点性能参数变化如表。优化后效率增加0.790%，压比提升0.475%，质量流量增加0.705%，出口熵增降低3.853%。

优化前后设计点性能参数

	流量/(kg•s-1)	总压比	效率
初始叶片	20.553	2.106	0.8730
优化叶片	20.698	2.116	0.8799

       图为优化前后出口性能参数对比，从图中可以看出出口效率在30%叶高以上提升，在30%叶高以下略有下降；总压比在整个径向高度全部提升。

出口效率对比

出口总压比对比

   图为效率特性和压比特性对比，从图中可以看出，叶片优化后效率在近堵点基本无提升，在近失速点提升较大，压比在近堵点提升较小，近失速点提升较大，说明叶片优化在逆压梯度越大时，优化效果越明显。总体来说叶片优化后提升了变工况性能。

效率特性对比

压比特性对比

叶片几何形状变化

       优化前后叶型对比如图，其中灰色为优化前叶型，红色为优化后叶型，叶型主要变化在前缘。

10%叶高叶片对比

50%叶高叶片对比

90%叶高叶片对比

流场分析

       叶片优化前后10%截面马赫数如图；从图中可以看出，叶片优化后叶片前缘马赫数降低，中部激波前马赫数降低，激波后分离区延后，尾缘分离区稍有增大。

a）初始叶片        b）优化叶片

10%叶高马赫数对比

   叶片优化前后50%截面马赫数如图；叶片优化后叶片前缘马赫数降低，中部激波变得更加倾斜，激波后分离区明显减小。

a）初始叶片        b）优化叶片

50%叶高马赫数对比

       叶片优化前后90%截面马赫数如图,叶片优化后前缘马赫数降低，中部激波更加倾斜，激波前马赫数降低。

a）初始叶片        b）优化叶片

90%叶高马赫数对比

    给出了叶片优化前后10%、50%、90%叶高叶片表面静压分布。从图中可以看出，优化叶片相比于初始叶片，叶片前部静压差增大，中部静压差减小，叶片中部负荷前移。10%叶高，优化后的叶片静压上升的位置延后，说明激波位置延后，在叶栅流动中激波与附面层相互干扰，产生流动分离，造成流动损失，激波位置的延后，能够减小分离区；在叶片后部，静压差稍有增加，叶片负荷增加。 50%叶高，静压上升位置延后，说明激波延后；叶片后部静压差减小，叶片负荷降低，在叶片后部由于激波和附面层的相互干扰，存在流动分离，流动情况较差，叶片后部负荷的降低，能够改善叶片后部的流动情况。在90%叶高，叶片中部激波位置延后，叶片后部静压差减小，负荷降低。

10%叶高静压对比

50%叶高静压对比

90%叶高静压对比

文章来源：starccm仿真学堂