新一代战斗机非定常流动数值研究综述

新一代战斗机强调超机动能力和强隐身性，其中大攻角下的静态失速、动态失速及内埋弹仓绕流是与高机动和强隐身密切相关的、极具挑战的几类典型的非定常流动，它们对数值方法提出了极高的要求。为了高精度地仿真流场、清楚地揭示流动机理，有效地控制非定常流动，非常有必要发展高精度且高效率的RANS-LES混合方法体系，包含RANS-LES混合方法本身、与RANS-LES混合方法匹配的高精度自适应耗散格式、基准湍流模式、高质量计算网格、高精度时间推进方法、非定常量的统计方法等，具有极强的紧迫性。提出、发展、验证并应用该类方法数值仿真新一代战斗机（包括单独部件、组合体、甚至全机）的非定常流动，数值预测结果与风洞实验数据吻合良好；此类方法可为新型战斗机设计提供理论依据和分析手段。     

跨声速串列转子失速机制的数值研究

为了解决跨声速串列转子的低裕度问题,就必须了解跨声速串列转子的流场结构与失速机制.设计了叶尖切线速度为450m/s,负荷系数为0.56的高负荷跨声速串列转子.基于数值模拟的结果,分析了该串列转子在0.5mm叶尖间隙下的叶尖流场结构与失速机制,并在此基础上分别探讨了叶尖间隙和前、后排叶片周向位置对串列转子特性的影响和失速机制的变化.结果表明:前排叶片的叶尖区域是 影响串列转子稳定性的关键;随着叶尖间隙的增加,串列转子的失速机制也发生变化,从前排叶片叶尖区域的尾迹与径向潜流堵塞后排叶片通道转变为前排叶片叶尖泄漏流堵塞;在较大周向相对位置(后排叶片压力面周向远离前排吸力面)的情况下,串列转子获得最好的效果,随着周向相对位置(PP)的增加,失速部位从后排叶片转移至前排叶片.      

基于非接触式测量的旋转叶片动应变重构方法

基于叶端定时非接触式测量和振动响应传递比的概念,开展高速旋转叶片动应变重构方法的研究。在频域内推导了叶片任意测点位移与任意测点动应变的传递比,给出了单模态共振下响应传递比关于位移和应变模态振型的解析表达式;建立旋转叶片的三维(3D)有限元模型,开展考虑旋转预应力效应的叶片模态分析,提取位移和应变模态振型,获得任意转速下叶端位移与叶根关键点动应变的传递比。开展高速旋转叶片叶端定时非接触式测量实验,采用周向傅里叶算法对叶端定时信号进行处理,获得叶片在不同转速单模态共振下的叶端位移,结合响应传递比,重构5个旋转叶片的关键点动应变。结果表明:旋转叶片在9000r/min和13000r/min转速下发生1阶共振时,与应变片实测结果相比,叶根处应力最大点、次大点和边缘点3个关键点的动应变平均重构误差均小于15%,验证了旋转叶片动应变重构方法的有效性。     

航空发动机适航审定喘振与失速影响因素

将气动热力基础理论与适航领域工程问题有机结合,对典型的喘振和失速影响因素开展了系统的研究,发展了一套多影响因素作用下航空发动机稳定裕度快速预估方法;在此基础上以JT9D发动机为研究对象,开展了喘振和失速主要影响因素的影响规律研究,结果显示对风扇部件稳定裕度影响权重最大的是进气畸变,对增压级和高压压气机影响权重最大的是寿命期内结构衰变。研究结果进一步明确了喘振和失速适航条款符合性验证活动的关键要素,完善了条款适航符合性验证流程,为进一步提升我国民用发动机适航审定能力提供了技术支撑。     

POD和DMD方法分析不同间隙压气机旋转不稳定性特性

对一台单级低速轴流压气机旋转不稳定性(RI)进行分析。在前期的旋转不稳定性特性、流场成因和不同间隙影响的研究基础上,进一步基于全通道数值模拟,采用两种流场降阶方法分析不同叶顶间隙条件下旋转不稳定性各阶模态对应的流场成分、模态能量占比和稳定性等特性。计算结果表明,前几阶模态能量占比大,对流场起主要贡献,并且稳定性相对较强。随着间隙的减小,旋转不稳定性的能量占比和稳定性下降,各阶的叶片转动成为流场的主要成分。      

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

旋转对曲率表面气膜冷却效率影响的数值研究

通过对旋转状态下曲率通道的气膜冷却现象的流动和换热进行数值模拟,得到了不同旋转数下凸表面和凹表面的冷却效率分布。计算选用剪切应力输运(SST)湍流模型,主流雷诺数ReD=4 797,吹风比M=0.4,旋转数Rt=0~0.023 9。研究结果表明,旋转数对冷却效率的影响明显：旋转数的增大使得气膜的轨迹偏转越明显,并且凸表面上气膜轨迹的偏转程度高于凹表面的。凸表面的冷却效率随着旋转数的增加而逐渐减小,而凹表面的冷却效率随着旋转数的增加而递增,并且旋转数的增加会弱化凸表面和凹表面上冷却效率的差别。     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之一:理论模型及算法分析

本文提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier－Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架。对压气机不稳定边界的判断通过计算和检验特征值的方式来实现。本文给出了有限差分法和基于Gauss－Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性。     

固定小翼颤振激励系统试飞验证

针对一种新研制的带开缝旋转圆筒固定小翼激励系统,对其工作特性和激励力特性进行了试飞验证。试验结果表明,该系统工作特性基本达到了最初的设计指标要求,激励力特性理论计算与实测结果基本吻合,为该系统进一步的应用奠定了基础。     

轴流压气机动态失速过程三维计算

基于三维非定常欧拉方程和三维激盘模型发展了一种用于研究轴流压气机动态失速过程的三维计算方法.利用该计算方法研究了某型高压压气机第一级的动态失速过程,并就转子总静压升特性的三种不同径向分布对失速过程的影响进行了分析.计算结果表明该方法能够反映出压气机的三维动态失速过程、失速团的三维空间结构,并能有效地表现出转子沿径向的特性变化对压气机气动稳定性的影响,而且还能够用于判断压气机的危险截面.