离心压气机近失速工况叶尖流场特征及泄漏涡面涡模型

为研究航空微型涡轮发动机离心压气机的失速机理,采用数值模拟的方法,对压气机在设计点和近失速点的流场特征进行了研究,并对叶尖泄漏涡进行了建模。结果表明:在近失速工况下,更多的来自轮毂处的低能流和二次流会汇入叶轮通道近壁面的低速区,使得叶尖的堵塞加重,泄漏涡轨迹的偏转及其卷吸作用可能是造成这种现象的原因。通过在压气机机匣外侧引入虚拟镜像涡的方法,应用机翼涡对的不稳定性理论和面涡模型对泄漏涡建立了模型。应用该模型对本文研究的离心压气机以及三个国内外具有典型代表性的压气机进行了预测,表明模型得到的频率同计算和实验结果具有较好的吻合性,对德国宇航院压气机的预测频率与其实际频率的误差仅为2.8%。总体而言,当间隙高度较大时,模型的预测更为准确。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段;在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波;第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波;稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s;斜激波的存在使燃烧室出口平面流场产生了较大波动。      

预旋喷嘴径向角度对预旋特性影响的数值研究

为了降低低位进气预旋流路的气动损失,针对带有不同径向角度（0°~30°）预旋喷嘴的预旋系统进行了数值仿真,并对流动特性、温降特性和比熵增特性进行了分析。结果表明:随着预旋喷嘴径向角度的增大,预旋系统无量纲温降先增大后减小,流动阻力减小,预旋系统的流量随之增大。旋转雷诺数为2.3×107时,预旋喷嘴带径向角度的预旋系统无量纲温降比传统喷嘴最大可提高18.3%,存在某一角度使预旋温降特性达到最好。预旋系统内的耗散主要发生在预旋腔和共转盘腔内,径向角度为10°时其比熵增变化量分别占整个预旋系统总体比熵增的42.4%和30.2%;合理设计预旋喷嘴的径向角度,能改善预旋腔内气流的流动效果,并且可以减少整个预旋系统的不可逆损失。      

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

CHN-T1标模2.4米风洞气动特性试验研究

大飞机标模可用来校验风洞流场品质,检验和提高大型飞机风洞试验数据质量,标模的外形特点及其气动特性能否反应现代大飞机设计特点尤为重要。气动中心前期已完成了一套大展弦比飞机标模CHN-T1的设计研制,为了验证设计结果,在2.4米跨声速风洞中进行了一期验证试验,试验Ma数范围0.40～0.90,模型名义迎角-6°～15°,侧滑角-3°～12°,雷诺数Re=(3.3～7.5)×10~6。试验内容包括纵横航向基本特性试验、重复性精度试验、变雷诺数试验、转捩对比试验、流谱观察试验和变形测量试验。结果表明,该飞机外形具有良好的升阻特性,符合现代大展弦比飞机的典型气动特征,可用于2.4m跨声速风洞大型飞机标模试验数据体系建设。     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

基于充气前缘技术的旋翼翼型动态失速抑制

动态失速的发生会在直升机旋翼桨叶和桨毂上产生高的交变扭转振动载荷,并限制直升机高速重载状态下的使用包线。本文利用计算流体力学(CFD)方法对基于充气前缘(ILE)技术的SC1095旋翼翼型动态失速抑制进行研究,分析了ILE抑制动态失速的控制机理,获得了ILE结构布置和充放气方式对动态失速的影响规律。研究表明:ILE可以有效抑制动态失速的发生;ILE最大膨胀程度越大,其抑制动态失速的效果越好,但膨胀程度过大后抑制效果开始减弱;ILE在翼型上仰至最大迎角时恰好达到最大膨胀状态,其对动态失速的抑制效果最好;ILE保持最大膨胀状态的时间长短对抑制效果影响不大;在翼型上仰至不同迎角时开始对ILE充气会对动态失速抑制有较大影响;ILE整流段与翼型连接位置对动态失速抑制有很大影响,整流段越长,抑制效果越好。     

小功率氮电弧推力器的性能

采用自然辐射冷却结构的小功率电弧推力器,实现了以氮气为推进剂的长时间稳定运行。采用间接测力方法得到推力,利用铜镜反射法拍摄喉道处的放电状态,结合测量的弧电压、弧电流和气流量数据以及导出的比冲和推力效率,对推力器运行性能和放电特性进行了研究。结果显示:在气流量为100~700mL/min,输入功率为35~55W的条件下,最大推力约为24mN,最大比冲接近175s,当弧电流为80~120mA范围内变化时,弧电压变化范围为420~520V,并且弧电压随气流量的增加呈现出先下降后上升的趋势。      

基于载荷当量安装误差的面齿轮振动特性

以面齿轮传动系统为研究对象,考虑载荷作用下面齿轮传动系统中支撑结构变形和轮齿弹性变形,通过将传动系统变形等效到面齿轮和直齿轮安装误差方向,建立包含当量安装误差的面齿轮多自由度耦合振动分析模型;采用4阶Runge-Kutta法求解了面齿轮传动系统运动微分方程,得到了当量安装误差对面齿轮振动加速度和法向动态啮合力的影响;开展了面齿轮动态特性测试试验,试验结果表明:偏置距误差和轴交角误差引起的面齿轮沿x方向振动加速度大于沿z方向振动加速度;偏置距误差对面齿轮x方向振动加速度的影响大于轴交角误差。