低雷诺数下二维扩压叶型的设计的方法的研究.docVIP

低雷诺数下二维扩压叶型的设计的方法的研究 　　摘要：本文通过数值模拟方法，对低雷诺数叶型设计方法进行了初步探索。结果表明，对于低来流马赫数，在叶型弯角一定的前提下，适当增大中弧线前段的曲率，减小吸力面曲线前段的曲率，可以有效的抑制流动分离。 　　关键词：低雷诺数；高空低速；叶型设计 　　 　　前言 　　 　　当雷诺数小于某一临界值时，发动机性能急剧恶化。为了保持低雷诺数下发动机的性能参数，有必要对叶型设计方法进行研究，提出低雷诺数叶型的设计准则。 　　本文通过数值模拟方法，对低雷诺数下二维扩压叶型设计方法进行了初步探索。 　　 　　1计算方案及初始叶型 　　 　　本文采用多项式叶型造型方法。其中，中弧线采用了四次多项式函数来确定。吸力面与压力面曲线均以各自的厚度最大点分为前后两段，分别采用两段三次多项式样条函数来确定。初始叶型如图1所示，相关参数见表1。其中，e1和e2分别为吸力面和压力面的最大厚度在轴向弦长的位置。Cmax1和Cmax2分别为吸力面和压力面的最大厚度与轴向弦长之比。 　　 　　文中采用计算软件NUMECA，在17000m高空、0攻角、0.5来流马赫数条件下，通过改变叶型的各项几何参数，对不同叶栅的流场进行了数值模拟。探讨在低雷诺数、低来流马赫数下，如何选取各项叶型参数，以达到改善叶栅性能的目的。 　　计算选取定常粘性的雷诺平均Navier-Stokes方程为数学模型， S-A模型为湍流模型，AGS模型为转捩模型。边界条件给定入口总压、总温、0°攻角和出口反压。其中，入口总压、总温为17000m高空的大气压强和温度，调节出口反压保证来流马赫数为0.5。 　　 　　2中弧线曲率分布的设计 　　 　　2.1 中弧线及叶型的变化规律 　　中弧线的曲率分布是通过中弧线前端点的二阶导 来控制的。图2展示了叶型随中弧线的变化趋势。可以明显看出，当 =0时，整个叶形前半段比较直，曲率小；后半段比较弯，曲率大。随着 的减小，叶型前半段逐渐变弯，曲率变大；后半段逐渐变直，曲率变小。 　　2.2中弧线曲率分布对损失系数的影响 　　通过数值模拟发现，随着中弧线前端点二阶导的减小，总压损失系数逐渐减小，在 等于-40时达到最低值，后又逐渐增大。对于亚音速静叶栅来说，损失主要来自吸力面的附面层分离。可以推断：随着前段曲率的增大，附面层分离造成的损失先减小，后增大；适当安排曲率分布，可降低附面层损失。 　　2.3流场分析 　　为了进一步确认中弧线曲率的合理分布对附面层分离的抑制作用，以下选取若干典型流场进行分析。 　　图4(a)~(c)分别为 =0、-20、-40时的马赫数分布。可以明显看出，随着 的减小，在吸力面上，叶片表面曲线前半段对气流的加速能力增强；后半段逐渐变平缓后，附面层分离区范围在减小。 　　 　　3吸力面厚度分布的设计 　　 　　3.1叶型设计 　　取中弧线前端点二阶导 =-40，改变叶型吸力面最大厚度点的位置e1，探究该参数的最优值。图5对比了e1等于0.2和0.8两种情况下的叶型。 　　可以发现，随着最大厚度点的后移，吸力面曲线的前段曲率有变小的趋势，后段则逐渐向上抬。 　　3.2 对损失系数的影响 　　图6给出了随吸力面最大厚度位置的后移，总压损失系数的变化趋势。数值模拟结果显示，随着吸力面最大厚度位置由30%后移到80%，总压损失系数以接近线性的趋势，减小了33.78%。 　　3.3流场分析 　　图7(a)~(c)分别为e1=0.3、0.5、0.8时的马赫数、流线分布。可以明显看出，随着最大厚度点的后移，吸力面曲线前半段对气流的加速能力下降，因此，流经该位置气流的马赫数降低。同时，吸力面后半段在保持平缓的前提下略微上抬，对并不是轴向流动的气流起到了更好的导向作用，使分离区减小。 　　 　　4压力面厚度分布的设计 　　 　　4.1叶型设计 　　取中弧线前端点二阶导 =-40，吸力面最大厚度点的位置e1=0.8，改变叶型压力面最大厚度点的位置e2，探究该参数的最优值。图8对比了e2等于0.2和0.8两种情况下的叶型。可以发现，随着最大厚度点的后移，压力面曲线的前段曲率有变大的趋势，后段则逐渐向下沉。 　　4.2对损失系数的影响 　　图9给出了随压力面最大厚度位置的后移，总压损失系数的变化趋势。数值模拟结果显示，随着吸力面最大厚度位置由20%后移到60%，总压损失系数由0.0598逐渐减小到0.0580，减小了3.36%。由60%后移到80%的过程中，损失系数稳定在0.058左右。 　　 　　5优化结果 　　 　　为了使叶型损失最小，本文依次对初始叶型的中弧线曲率分布、吸力面最大厚度位置和压力面最大厚度位置进行了优化，最终选定 =