单转子轴流压气机涡动力学失稳机理研究

为了研究单转子轴流压气机的涡动力学失稳机理,采用基于Shear Stress Transport(SST)湍流模型的尺度自适应雷诺平均/大涡(RANS/LES)混合模拟的方法对低速单转子轴流压气机进行了非定常数值模拟。研究结果表明:在设计转速3kr/min条件下,叶顶泄漏涡、二次泄漏涡以及诱导涡破碎引起的叶顶区域的堵塞是触发单转子轴流压气机内部流动失稳的主要因素。压气机由近堵塞工况点向小流量工况点逼近的过程中,叶顶泄漏涡轨迹与轴向的夹角由70°增加到76°,二次泄漏涡起始点位置前移加速叶顶泄漏涡向转子前缘移动。近失速工况点叶顶泄漏涡的轴向动量与主流的轴向动量之间存在一种平衡,叶顶泄漏涡稳定在转子前缘。压气机进一步节流主流的轴向动量减小,对叶顶泄漏涡轴向动量的抑制能力减弱,叶顶泄漏涡的位置不再稳定,诱发尖脉冲型失速先兆。     

三维超声速不等密度平面混合层的小激波形成与演变

基于大涡模拟方法与高精度混合格式,对两种密度比(1.0与0.138)条件下的高对流马赫数(Mac=1.1)三维空间发展超声速平面混合层进行了数值模拟。数值结果清晰描述了扰动激励下超声速平面混合层的演变过程,讨论了三维混合层中小激波的结构及其形成过程。分析了混合层两侧密度变化对混合层失稳及小激波形成的影响,发现当上、下层密度不等时,混合层失稳加速,小激波强度减弱,而密度较低侧的小激波会更弱。     

端区流动对跨声速压气机失速影响的数值研究

为了揭示近失速工况下动叶间隙和静叶角区的复杂流动结构,探索诱发压气机失速的主要因素,对高负荷跨声速压气机开展数值研究,获得了整级压气机和单转子的特性曲线,进而对压气机的流场进行了详细分析。结果表明:对于研究对象,失速最先起始于静叶上角区。上角区完成由开式分离到闭式分离的转变过程,可以认为是静叶发生失速的标志。在近失速工况下,动叶叶尖泄漏涡发生"泡式破裂"堵塞流道,泄漏涡破裂引起的堵塞作用一直从叶顶延伸至70%叶高左右。     

冲击距与气膜孔方位角对旋流气膜冷却性能影响

为了获得冲击距Hi与气膜孔方位角α对旋流气膜冷却性能的影响规律,以六边形供气腔圆形气膜孔平板气膜冷却结构为研究对象,对五种冲击距参数(0.74D,1.14D,1.54D,1.94D,2.34D)(D为气膜孔直径)和五种气膜孔方位角参数(0°,10°,15°,20°,25°)进行了三维数值计算研究,得到了绝热壁面气膜冷却效率、展向平均气膜冷却效率、流场空间无量纲浓度分布等随冲击距与气膜孔方位角的变化规律,分析了肾形涡对旋流气膜冷却性能的影响机理。结果表明:冲击距对绝热壁面气膜冷却效率展向分布规律影响不大,而方位角增加能够明显提高绝热壁面气膜冷却效率及展向气膜覆盖面积,方位角0°模型展向气膜冷却效率最大值为0.42,方位角25°模型的最大值为0.48,相比前者增加14.3%;绝热壁面同一流向位置的展向平均气膜冷却效率随冲击距的增加而增大,随方位角的增加而增大,Hi=2.34D时的展向平均气膜冷却效率最佳,α=20°时的展向平均气膜冷却效率曲线最佳。方位角的增加能够明显破坏流场中存在的肾形涡结构。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

基于 LES 方法的增升装置气动噪声特性分析

在气动噪声数值计算中,流场的求解精度对涡流扰动的细节计算以及声学的求解结果有着重要的影响。本文应用 LES 方法对增升装置的流场进行数值模拟,采用可穿透积分面的 Ffcows Wil1iams-Hawkings(FW-H)积分方法进行远场噪声计算。采用圆柱绕流算例对本文的数值计算方法进行了验证,验证结果表明：本文所使用的LES 方法能准确地捕捉到涡脱落、流动分离等非定常流动现象,可为远场气动噪声的计算提供精确的近场流动的数值解;基于 FW-H 的声类比方法能够精确高效求解远场气动噪声。在此基础上,对增升装置噪声产生的流动特性、远场特性、风速影响等进行了数值模拟研究。结果表明：缝翼产生气动噪声的主要原因是,流动在缝翼和主翼之间的凹槽形成的不稳定波以及缝翼钝后缘的小脱落涡;襟翼产生气动噪声的主要原因是,襟翼附近由于流动分离产生的高频的小尺度不稳定涡和低频的大尺度涡。     

基于Broyden算法的航空发动机气路故障诊断

针对基于Kalman的故障诊断算法响应速度慢、多故障诊断及非设计点诊断精度低的问题,提出一种基于改进Broyden算法求解方程组的航空发动机气路故障诊断方法。针对涡轴发动机,以模型输出跟踪发动机输出为准则确定3个方程,结合发动机模型中的2个平衡方程,构建气路故障诊断方程组,通过改进Broyden算法求解方程组以获得部件性能退化因子及模型猜值。数字仿真结果表明,所提出的基于Broyden算法求解方程组的航空发动机气路故障诊断方法,在包线内的单故障和多故障诊断稳态误差均小于0.35%,且诊断过程算法单步运行最大耗时小于2ms,具有良好的实时性,远优于Kalman滤波方法,验证了算法的先进性。     

高度欠膨胀声速射流的自相似特性

构建了高压气体射流模型和高分辨率的计算网格,对喷压比（NPR）为5.60、流动雷诺数Re为105量级的高度欠膨胀射流进行了三维大涡模拟（LES）计算.讨论了时均的射流近场结构,发现大涡模拟成功捕捉到了高度欠膨胀射流近场的典型波系结构,并与文献结果吻合较好.研究着重定量考察流向速度和氮气质量分数的分布规律,以揭示高度欠膨胀射流的自相似特性.结果表明:在射流核心区之后,不同流向位置上流向速度和氮气质量分数沿径向的分布呈现出与亚声速射流类似的自相似特性.但射流流场开始呈现自相似的位置仍然为超声速,这是高度欠膨胀射流自相似特性的独有特征.提出了高斯拟合公式,在靠近和远离射流中心线的区域内均可较好地表征射流的自相似特性.此外,还考察了射流剪切层的发展特征,比较了计算得到的射流剪切层增长率与先前实验测量结果的差异,并分析了其中可能的原因.      

扰动形状对钝头体非对称流动的影响

通过在钝头体头部施加人工扰动块可以得到确定的大攻角下的非对称背涡结构。为了研究扰动块形状对非对称背涡结构的影响,本文在攻角50°、雷诺数Re_D=1.54×105的条件下,利用数值模拟对周向角90°、子午角10°的扰动位置的半球形、D型及方形3种扰动块形状分别进行了研究。研究发现在同一扰动位置,半球形扰动主控下的背涡结构为右涡型,而D型扰动和方形扰动主控下的背涡结构呈现左涡型,且方形扰动主控下的背涡结构的非对称性弱于其他2种扰动主控的非对称背涡。通过分析发现扰动块所引起的微流动直接影响钝头体非对称背涡结构。因此为了更精准地通过施加人工扰动得到确定的非对称背涡结构,应尽量选择形状简单、表面平滑过渡的扰动块形状。      

基于大涡模拟的平流层浮空器气动特性分析

为研究不同球体间距对“球形囊体型”平流层浮空器(SLV)气动特性的影响,采用有限体格式、结构网格和大涡模拟(LES)计算方法求解不可压缩的Navier-Stokes(N-S)方程,对超临界雷诺数“球形囊体型”平流层浮空器绕流进行数值模拟,并对不同球体间距下的数值计算结果进行详细的分析比较.通过对比试验数据,单球体数值模拟的阻力系数时均值与Achenbach的试验数据一致,验证了计算方法分析超临界雷诺数球体绕流问题的准确性.研究不同间距的双球体阻力变化规律以及振动频谱特性,随上下游球体间距G的增加,合阻力先增大后减小,上游球体的阻力占优振动频率逐渐减小;G=1.5D(D为球体直径)和G=2D时,上下游球体有相同的占优振动频率.随间距G的增加,两球体相互作用与上游球体对下游球体的尾涡结构影响逐渐减弱.