航空发动机中冰晶结冰的研究进展

一般认为经过风扇压缩后空气温度高于冰点，压气机内不会发生结冰。然而近年来，国外研究者对若干起发动机推力损失故障的研究表明冰晶能够引起压气机内部结冰，并在试验室中再现了温度高于冰点时压气机内结冰这一现象。我国在航空发动机冰晶结冰研究方面刚刚起步，为准确掌握国内外研究现状，从数值和试验两方面对国内外已开展的压气机内结冰研究进行了分析与总结，讨论了当前冰晶结冰的主要研究成果和存在的局限，提出了开展冰晶结冰研究需重点关注的方向，为我国航空发动机冰晶摄入结冰研究和适航审定研究提供一定参考。     

S形进气道流动控制数值模拟研究

采用CFD技术,结合试飞数据,对某S形进气道进行了加涡流发生器的流动控制数值模拟研究.着重分析了三个不同位置加涡流发生器后,进气道内部二次流的发展;之后比较了不加涡流发生器及不同位置加涡流发生器时进气道出口总压恢复、畸变等情况.结果表明涡流发生器明显地影响着进气道内部二次流的发展变化,涡流发生器对进气道出口周向稳态总压畸变有较大程度改善,但是对于提高总压恢复效果不明显.     

旋转径向孔流动特性计算模型

为满足空气系统一维仿真分析的需要,深入分析了影响径向孔流动特性的主要因素,建立了针对旋转径向孔的计算模型。模型考虑进气预旋和轴向错流的影响,并利用Simulink建立其仿真模型,计算径向孔的流量系数和空气流量。利用文献中的试验数据对径向孔流动特性计算模型进行验证。结果表明:在相同条件下,模型的结果与试验结果的最大误差不超过6.2%,最小误差约为0.3%。径向孔模型准确可靠,可以作为空气系统一维计算中的通用模型使用,能有效地提高建模和计算效率。     

S弯扩压器中四种湍流模型的比较

采用Realizableκ-ε,S-A(Spalart-Allmaras),SST(Shear stress transport剪切应力输运),RSM(Reynolds stress equation model雷诺应力模型)四种湍流模型对S弯扩压器内的分离流动进行了模拟,分析了不同湍流模型对计算结果的影响.将不同轴向位置的表面压力系数、分离区的大小、马赫数等三个量的实验值与计算值进行了比较,结果表明:分离区湍流模型对表面压力分布影响较大,未分离区不同湍流模型模拟精度相当;S-A模型对分离区的模拟较保守,RSM模型对分离的模拟精度较高,SST与RSM模拟结果很接近;各湍流模型计算的高速区域均比实验值大.扩压器的平均总压恢复系数受湍流模型影响不大,和实验结果吻合较好.      

旋转帽罩结冰相似准则的冰风洞试验研究

为了验证旋转帽罩冰风洞试验相似准则的适用性，开展了旋转帽罩表面结冰过程研究，推导了旋转帽罩结冰相似准则。采用底面直径100 mm、锥角80°的旋转锥作为原尺寸模型，底面直径50 mm、锥角80°的旋转锥作为缩比模型，进行不同工况冰风洞试验。通过图像法采集了帽罩表面冰形，并将冰形进行无量纲化处理，对比了原尺寸模型与缩比模型的结冰冰形。结果表明:对于明冰工况，帽罩前半部分冰形吻合较好，帽罩后半部分吻合欠佳；对于霜冰工况，全尺寸模型与缩比模型的冰形吻合良好。对造成这一现象的原因进行了分析，并提出改进的建议。     

机翼三维结冰数值模拟

采用欧拉-拉格朗日法计算空气、水滴两相流动,根据质量守恒和能量守恒计算结冰量和溢流水量,溢流水的流动方向及不同方向的流量分配由当地空气的速度矢量、网格边界的法向量计算得到。采用非稳态方法、非结构网格计算了NACA0012平直翼、截面为GLC-305的后掠翼的结冰,并与文献中的实验结果进行了对比,结果显示最大结冰厚度与实验吻合较好,预测的冰形与实验冰形发展趋势一致。计算结果表明,提出的三维结冰计算方法是有效、可行的,为三维结冰计算提供了一条可选的途径。     

人工制备冰雹的力学性能试验研究

人工制备的冰雹可用于多种航空附件冰撞测试。本文以ASTM F320—2010标准为基础,设计制作了多种冰雹,对冰雹的力学性能及影响冰雹力学性能的因素开展了研究。采用低温万能试验机对冰雹进行准静态试验,使用高速摄影技术深入分析冰雹破碎过程及特征,对比分析不同棉纤维种类、棉纤维特性和棉纤维含量下冰雹的力学性能。研究结果表明,含棉冰雹失效过程可分为3个阶段：线性增长阶段、屈服阶段和破碎阶段,失效过程中最大抗压强度出现在线性增长阶段末端。棉纤维的强度和弹性越好,冰雹破碎所需的加载力越大;棉纤维上附着的冰量越少,冰雹破碎所需的加载力越小。棉纤维含量增加可以加强棉纤维的连接作用,使含棉冰雹强度增加;但棉纤维含量到达15%时,继续增加棉纤维含量,冰雹抗压强度不会发生太大变化。亲水棉制成的冰雹抗压强度明显高于疏水棉制成的冰雹。研究结果可为中国民用大飞机适航符合性验证提供参考。     

改进的综合CFD与畸变合成法预测动态畸变(英文)

对综合CFD与畸变合成方法进行了改进。通过CFD求解,获得AIP上各稳态气动参数。采用BP人工神经网络直接求解高度非线性的湍流相关方程,该方程将CFD求解的稳态气动参数与试验测得湍流度进行关联。经过训练的网络可根据稳态CFD结果预测湍流度。之后根据CFD计算的稳态总压和神经网络预测的湍流度合成动态压力。最后根据合成的动态压力寻找最大瞬时畸变。文中使用了六个工况的飞行数据验证该方法的准确度并且与原方法的计算结果进行了比较。结果表明该方法预测的畸变值与试验畸变值吻合较好,且预测结果较原方法与试验值更接近,该方法能够在试验前用来预测动态畸变,具有很强的工程价值。     

气冷涡轮叶片三维换热问题计算

建立了三维气冷涡轮叶片换热计算的热-流耦合计算模型, 将辐射热量作为源项加入到方程中, 采用热-流耦合方法完成了某气冷涡轮叶片带辐射的三维换热计算.辐射热量的计算采用了P-1辐射模型, 比较了计算辐射和不计辐射时叶片表面的温度值及温度分布的差异.计算表明对于1600 K这样的高温燃气, 辐射热量不可忽略, 考虑辐射热量后, 叶片温度明显升高, 温度分布也产生了显著差异.      

航空发动机旋转帽罩结冰表面换热系数研究

旋转帽罩表面的对流换热系数是结冰、防冰研究的重要参数,旋转帽罩结冰模拟中对流换热量和蒸发量都需要已知对流换热系数。为了获取旋转帽罩表面的对流换热系数,采用数值模拟对旋转帽罩结冰表面的对流换热系数进行了研究。首先验证了网格与计算方法的合理性和可靠性,在此基础上,对影响表面换热的因素进行了分析,计算了不同转速、来流速度、来流温度对锥角40°,锥高176mm的旋转帽罩表面换热系数的影响。结果表明:锥尖区域,来流速度的影响与转速影响相比占主导;除锥尖外的其他区域,转速对换热系数的影响占据了主导地位;来流温度对整个表面的换热系数均有影响。采用不同工况下的数值模拟结果,在锥尖区域建立了来流雷诺数与努赛尔数之间的关联式;在除锥尖外的其他区域建立了旋转雷诺数与努赛尔数之间的关联式。