涡扇发动机短舱结冰试验相似方法

通过分析短舱表面结冰缩比试验相似要求,在保证绕流流场、水滴撞击特性和结冰冰型相似的基础上,分别采用基于气流雷诺数,气流韦伯数和水滴韦伯数3种速度选取方法建立了试验相似准则。利用该准则获得了明冰和霜冰条件下涡扇发动机短舱1/2缩比模型的结冰试验参数,开展了短舱进气道表面结冰冰型预测,对参考模型和缩比模型表面过冷水滴撞击特性、结冰特性、溢流水分布特性进行了分析。结果表明:采用基于气流韦伯数和水滴韦伯数速度选取方法的相似准则,在明冰和霜冰条件下均能保证缩比模型与参考模型表面水滴撞击特性相似,且相似性不受结冰温度影响；溢流水变化趋势与参考模型相似,流动极限相近；缩比模型表面冰型均与参考模型相似,该方法能够为短舱结冰试验参数选取提供依据。      

过冷大水滴条件下结冰相似准则

虑到中国现有结冰风洞在过冷大水滴(SLD)结冰条件模拟能力方面的局限,发展SLD结冰风洞及其相应的缩比方法迫在眉睫。研究了大水滴运动机理,并依据绕流流场相似、水滴运动和撞击特性相似、结冰过程热力学特性相似原理,提炼出SLD结冰条件下的6种相似缩比准则,利用结冰数值模拟软件FENSAP-ICE对所提相似缩比准则进行筛选与比较,从而基于结冰相似百分比量化判断准则找出给定条件下的最优准则法。结果表明:考虑SLD运动特性的撞击参数与ONERA法相耦合的方法以及撞击参数与Ruff法相耦合的方法的计算结果与实际成冰构型吻合良好；此外,撞击参数与ONERA法相耦合的方法在计算直径为171 μm的SLD时优于其他。      

旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

MVD对翼型冰形及气动性能影响数值模拟

结冰问题严重威胁现代飞行器的飞行安全。当前,数值模拟结冰过程是一种高效的结冰问题研究方法。本文首先对OAXXX翼型生成六面体计算网格,然后对该翼型的流场、水滴撞击、结冰过程进行了数值模拟研究,重点研究了不同的MVD(平均水滴直径)对于水滴收集效率、冰形的影响,并将不同MVD情况下结冰后的翼型气动特性和结冰前进行了对比。     

一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程；四阶龙格一库塔法求解水滴运动轨迹方程；逐次二分法寻找水滴撞击极限，然后计算局部水滴收集率；分析翼型表面的传质、传热项，建立相应质量、能量方程并迭代求解，得到翼型表面结冰质量；根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好，结冰过程中的特征参数与冰型及相应冰形特征吻合。     

临界冰形确定方法及其对气动特性影响研究

临界冰形是指在适航规章结冰包线内,每个可用飞行构型下,对飞机操纵性和稳定性影响最严重的冰形,临界冰形分析是飞机适航取证中的重要工作。对临界冰形确定方法,及临界冰形对气动特性的影响进行了研究。发展了基于 CFD 方法计算临界冰形的一般方法,包括临界冰形分析状态、敏感性分析截面确定、结冰参数敏感性分析、临界结冰条件确定、临界冰形确定等。流场计算采用中航工业空气动力研究院气动力计算平台(UNSMB),基于 Jameson 中心格式的有限体积法求解 N-S 方程;水滴撞击特性计算采用 Eulerian 方法求解水滴轨迹运动方程;结冰计算采用经典的 Messinger 热力学模型。选取 CRM 飞机为研究对象,以机翼外翼50%展长处为敏感性分析截面,在典型飞行条件下,分析了结冰对环境温度、水滴直径、飞行速度、飞行迎角等参数的敏感性。利用“几何外形敏感性分析方法”,即通过对比冰形的上下冰角角度和冰角厚度等冰形几何参数来确定最严重冰形,得到了CRM 飞机的临界结冰条件和临界冰形,其中敏感性分析截面在水滴直径为30μm 时上冰角厚度和下冰角厚度最大,冰角最大厚度约41 mm。计算了结冰后的气动性能衰减规律,临界冰形对飞机气动性能影响严重,导致升力降低6.7%~23.8%,阻力增加17%~70.9%。发展了45min 待机临界冰形确定方法,基于几何外形敏感性分析方法进行环境温度、水滴直径、飞行条件等各类参数的结冰敏感性分析,得到飞机的临界结冰条件和临界冰形,对于民用飞机设计和适航取证具有一定的工程应用价值。     

三维翼面结冰过程及其影响数值研究

为得到三维机翼及翼身尾构型在结冰条件大气中的结冰过程和影响,采用结构化网格对M6机翼和翼身尾构结冰过程进行数值模拟,对飞机结冰过程的数值模拟采用欧拉两相流法求解水滴撞击特性,基于Messinger热力学模型模拟冰形增长过程,对其水滴撞击特性、结冰形状及结冰后对气动特性影响进行研究.分析结果可得:机翼展向,随着弦长变小,结冰量逐渐增多,尤其对大展弦比机翼更为明显.同时机身前端迎风面也是比较容易积冰的区域.     

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

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