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

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

旋转帽罩结冰相似准则的研究

为了研究航空发动机旋转帽罩结冰相似理论,在静止部件结冰相似准则的基础上,忽略次要因素的影响,建立了一种适用于旋转帽罩的结冰试验相似准则。利用建立的旋转结冰相似准则,通过数值模拟方法,对旋转帽罩缩比前后的水滴收集特性和结冰冰形进行了误差分析,并对此准则的可靠性进行评估。计算结果表明:用建立的相似准则所确定的参数对于全尺寸旋转帽罩转速分别为1kr/min和5kr/min时进行数值模拟。在全尺寸旋转帽罩及其1/2缩比模型上的对应区域得到的水滴收集特性和结冰冰形,误差约为2.7%和5.3%,表明选取缩比参数的方法是可靠的,并且提出的相似准则是有效的,对结冰风洞试验的理论研究和参数选取有一定的参考价值。     

离心力对于旋转叶片表面积冰影响的数值模拟研究

为了加深对旋转叶片表面积冰现象的认识,发展了一种适用于旋转表面积冰预测的数学模型,并采用该模型研究了离心力对于旋转叶片表面积冰的影响。通过对积冰控制体内的由于撞击水滴形成的薄水膜流动的质量、动量和能量守恒进行分析,发展了旋转贴体非正交曲线坐标系下的积冰模型并给出了相应的计算方法。所发展模型的动量方程中考虑了离心力对水膜流动的影响,并且在能量方程中考虑了净流出水带走的能量。采用文献中旋翼表面的积冰试验结果对所发展的模型和计算方法进行了验证,所计算得到的冰厚与试验结果吻合较好,计算得到的冰厚精度相比于LWICE软件有所提高。计算结果还表明,离心力增大会导致水膜速度在叶片展向上的分量变大,致使流出控制体的水质量增加,最终导致冰厚在驻点附近略微变小。     

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

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

单孔冲击式帽罩前缘流动换热特性数值分析

为了提高航空发动机帽罩冲击防冰结构的设计分析水平,对单孔冲击式帽罩前缘结构的流动换热特性进行数值研究,分析了不同冲击孔径与不同冲击雷诺数对帽罩前缘速度流场、换热系数与努塞尔数的分布规律。结果表明:在冲击雷诺数一定的条件下,冲击孔径越大,射流核心速度和前缘壁面附近的气流速度越小,前缘冲击区形成的涡流团越大,当孔径D=6 mm时,小孔径冲击下前缘区整体换热效果不如大孔径的,而在滞止区的换热效果则要优于大孔径的;当D>12 mm时,孔径大小对壁面换热基本没有影响;在冲击孔径相同时,增大冲击雷诺数使得冲击射流、前缘壁面附近及侧壁曲面通道内的气流流速增大,冲击区内的涡流团则逐渐减小;冲击雷诺数的增大也增强了前缘冲击区的换热特性。     

发动机旋转整流帽罩的水滴撞击特性分析

基于计算流体力学软件FLUENT,采用拉格朗日法,提出了一种求取旋转部件表面水滴收集系数的方法,并采用用户自定义函（UDF）予以实现。将该方法运用于发动机旋转整流帽罩的水滴撞击特性研究,结果表明：随着帽罩锥体锥角及转速的增大,帽罩表面的局部水收集系数沿其表面外形线下降较快,且撞击区域变小。     

发动机进气帽罩防冰热载荷的数值模拟研究

对某发动机进气帽罩防冰热载荷的计算方法进行研究。利用FLUENT软件对帽罩周围的流场进行模拟,同时采用FLUENT中的离散相模型计算帽罩周围的水滴运动轨迹,得到帽罩表面的水滴撞击效率;利用数值计算结果及流动换热经验关系式分析帽罩表面的各项热流,获得该发动机帽罩的防冰热载荷。计算中还比较了不同湍流模型对防冰热载荷的影响。     

锥角对旋转整流罩积冰影响的模拟实验

采用缩比模型对4种不同锥角的旋转整流罩进行了冰风洞积冰模拟实验.推导出旋转表面积冰相似参数,并根据重要参数的匹配确定缩比整流罩模型的积冰实验参数.在实验中采用高速摄像系统记录冰生长过程及最终冰形.结果表明:4种锥角的整流罩表面积冰均由初期生成、分布连续的明冰和后期增长迅速的白色霜冰构成;锥角小于等于74°的整流罩表面霜冰为针状或粒状,积冰厚度较小;锥角大于80°的整流罩表面霜冰为羽毛状,积冰厚度较大,并伴随冰脱落现象,锥角较大的整流罩冰脱落位置向下游移动.      

发动机旋转表面水滴撞击特性数值研究

介绍了发动机旋转机械内部过冷水滴轨迹的三维数值计算方法,阐述了旋转坐标系下气流及粒子的运动规律。采用ANSYS-CFX软件及其粒子输运模型,对某发动机风扇转子叶片外围空气及水滴流场进行了数值模拟。利用水滴的速度矢量、撞击区域等参数表征水滴的撞击特性,获得了转速、水滴直径等对风扇叶片表面水滴撞击特性的影响:水滴撞击区域集中在风扇叶片迎风面叶盆侧,且水滴撞击区域随着转速的增加而减小;水滴在叶片表面的撞击范围随着水滴直径的增大而减小。     

改进的翼型积冰数值模拟方法

基于微团概念,以雨滴运动的基本假设出发,对翼型积冰的拉格朗日方法进行了改进.新方法克服了原方法计算效率低、冰形同雨滴分布密度及网格密度相关的缺点,继承了原方法的鲁棒性,不存在欧拉方法的收敛性问题.为验证方法的正确性,计算了NACA0012翼型在0°和4°攻角下的霜冰积冰情况、表面水收集量和局部撞击系数.结果表明新方法减少了计算量,预测的冰形和实验结果一致,并分析了积冰对气动性能的影响.      

四自由度飞机结冰脱落模拟研究

飞机表面的冰粒脱落可能会对安装在尾部的发动机和其他飞机部件造成破坏,直接影响到飞机的安全飞行。在对实际结冰脱落过程进行合理简化的基础上,建立了四自由度结冰脱落轨迹模型,应用3阶Runge-Kutta方法求解冰粒运动方程,并嵌入了气动阻尼对冰粒轨迹的影响;同时将计算得到的冰粒轨迹及转角与文献中的计算和试验结果进行了对比验证;最后模拟并分析了不同形状、尺寸、初始俯仰角及气动阻尼对于冰粒运动轨迹的影响,为飞机表面结冰脱落的进一步研究工作提供了参考。