CJ818客机机翼防冰范围的确定及防冰热载荷计算

通过数值计算的方法,对CJ818客机机翼周围的流场和水滴运动轨迹进行了计算;介绍了水滴撞击特性的计算方法,编写程序计算了机翼的水滴撞击特性,并据此确定出机翼的防冰范围;分析了防冰表面热流的组成,并给出了每种热流的计算方法,在此基础上编写程序计算得到了CJ818客机机翼的防冰热载荷。     

水滴撞击特性的数值计算方法研究

 对水滴撞击特性的数值计算方法进行了研究。在对防冰部件流场进行计算的基础上,通过对水滴所在单元的判断和单元内流场速度的插值求解以及水滴运动轨迹计算边界的判断,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数,为飞机防冰系统的设计奠定了基础;此外,以发动机进气道的水滴撞击特性的计算为算例,研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响。     

发动机进口支板防冰热负荷数值模拟

为研究航空发动机进口防冰问题,应用三维雷诺平均N—S方程和双方程k～ε模型,计算了某发动机进口支板外部空气流场,用欧拉法建立并求解了空气中过冷水滴控制方程;选取某些飞行工况,得到了支板防冰表面的水滴撞击特性,且分别在干、湿蒸发2种情况下得到了表面的防冰热负荷分布。可以此为基础进行支板表面加热形式和结构设计.     

发动机防冰支板水滴撞击特性的数值研究

在对发动机防冰支板附近的流场进行计算的基础上,使用离散相模型在拉氏坐标下模拟了该流场中过冷水滴的运动轨迹,对发动机防冰支板的水滴撞击特性进行了研究.研究表明,水滴撞击极限、总收集系数和局部收集系数随飞行高度、来流速度及水滴直径的增加而增大.     

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

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

水滴撞击特性的重力影响分析

以NACA0012翼型为研究对象,采用拉格朗日法分析了重力对水滴撞击特性的影响,获得了过冷大水滴（SLD）的撞击特性。计算结果表明：当水滴直径超过50μm时,重力对水滴轨迹以及极限轨迹的初始位置影响很大;SLD的水收集系数和收集区域都较直径为10～40μm的水滴（美国联邦航空条例FAR25部附录C规定的防冰系统设计用水滴直径范围）有较大幅度的提高,该结论可为飞机防冰系统的设计以及冰风洞中喷嘴安装位置的确定提供参考。     

蛇形进气道影响下发动机进口部件水撞击特性的数值研究

为了了解亚声速蛇形进气道影响下航空发动机进口段迎风表面处的水滴撞击特性,采用欧拉/欧拉法对蛇形进气道和发动机进口段在来流马赫数为0.3及环境温度为266.5K的大气结冰条件下进行了空气/过冷水滴两相流场的计算,再由已获得的水滴流场进行支板表面水撞击特性的计算。整流罩表面水滴的撞击计算时采用插值的方法将原两相流场的计算结果插值到新生成的整流罩的近壁网格处,之后采用新生成的整流罩面网格进行水撞击特性的计算。通过计算获得了发动机进口结冰参数、支板以及整流罩表面水收集系数的分布。计算结果表明整流罩和支板表面水滴撞击的主要区域与进气道的弯曲方向有关且位于其下表面一侧,最大水收集系数达到0.98。与之相对比的整流罩和支板其他区域的水收集系数则显得非常微小,在数值上约等于0。此外,对于与进气道对称面平行的竖直支板,其与机匣相交的区域存在水收集系数为0遮蔽区。     

基于CFD技术的飞机风挡加温系统热载荷计算

热载荷计算是飞机风挡加温系统设计的基础。选取不同的飞行条件和结冰气象,基于CFD技术,通过外流场和水滴撞击特性计算得到了风挡外表面的压力分布、热流密度、对流换热系数以及局部水收集系数,最终确定了风挡防冰和防雾热载荷,并对仿真结果进行了瞬态校核。经验证表明,本文使用的计算方法合理,可以应用于工程设计。     

防冰热载荷计算的一种新方法

提出了一种计算飞机防冰热载荷的新方法,该方法利用CFD软件FLUENT及其自带的用户自定义函数UDF进行数值模拟。首先用FLUENT的两相流欧拉模型模拟空气-水滴两相流动,进而使用FLUENT内部的数据结构,通过UDF二次开发计算水滴撞击量和防冰热载荷,并且利用FLUENT的后处理功能对计算结果进行可视化图形显示。以某3D迎角传感器为例,进行了三维空气-水滴两相流和防冰热载荷计算,防冰所需热功率计算结果和实验数据误差为7.5%,在工程计算允许范围内。     

考虑遮蔽区影响的旋翼三维水滴撞击特性计算新方法

针对直升机旋翼三维黏性流场特有的复杂环境,建立了一种基于欧拉法的旋翼三维水滴撞击特性计算的新方法。首先,在旋翼桨叶嵌套网格的基础上,发展了一套用于预测旋翼绕流流场的计算流体力学(CFD)模拟方法。然后,为克服传统直升机旋翼二维水滴撞击特性计算方法的不足,充分考虑旋翼流场的三维效应,在嵌套网格中基于欧拉法求解旋翼三维水滴撞击流场。其中,为解决尾流等区域的密度脉冲现象所引起的稳定性和收敛性问题,提出并建立了遮蔽区扩散模型。该模型通过判断遮蔽区变量,在计算域中动态生成遮蔽区域,并随迭代步数逐渐扩散。最后,通过与NACA0012翼型及国外UH-1H桨叶的试验和计算结果的对比,验证了旋翼三维水滴撞击特性计算新方法的可靠性,并进行了温度和水滴当量直径(MVD)对旋翼三维水滴撞击特性的影响分析。结果表明:遮蔽区扩散模型的加入,使二维情况的计算时间减少了22%,并增加了三维情况的计算稳定性,显著提高了旋翼三维水滴撞击特性的计算效率;沿着旋翼桨叶展向位置增大的方向,旋翼桨叶剖面水滴撞击范围有所增大,最大水滴局部收集系数呈先增加后减少再增加的变化趋势,其变化幅度接近50%;旋翼桨叶表面的水滴撞击区域和水滴局部收集系数随水滴当量直径的增加而增加。     

Effects of flow parameters on thermal performance of an inner-liner anti-icing system with jets impingement heat transfer

The multiple jets impingement heat transfer is widely applied in the wing anti-icing system. It is challenging to apply the similarity criterion to carry out the anti-icing experiments due to the complex flow and heat transfer behavior. In the present study, the full-scale slat model is used to carry out anti-icing experimental researches in a 2 m × 3 m icing wind tunnel of China Aerodynamics Research and Development Center. The effects of icing parameters Liquid Water Content (LWC) and Median Volume Diameter (MVD) and hot air parameters (mass flow rate and temperature) on the thermal performance of an inner-liner anti-icing system with jets impingement heat transfer are studied. The effects of the experimental parameters are analyzed in detail by combining impingement and evaporation heat transfer mechanisms. The impingement hot air mass flow rate dramatically affects the heat transfer performance of the impingement stagnation region within the range of the experimental parameters. The temperature of impingement hot air and that of wing skin are approximately linear correlated. The experimental results show the effects of LWC and MVD on water film formation and runback ice accretion. The formation of water film is analyzed by an analytical method based on the wing skin temperature difference of dry and wet air conditions.     

大型商用运输机机翼增升构型水滴撞击特性计算

基于欧拉-拉格朗日数值模拟框架,发展了一套适合求解三维多体模型水滴撞击特性的计算方法。采用粒子统计与面积比率计算相结合的方式,实现了水滴收集率的高效精确预测。在常规粒径条件下,该方法计算结果与国外标模试验数据吻合良好。应用上述算法,开展了某大型商用运输机增升构型的水滴撞击特性数值模拟研究,计算结果显示:水滴在缝翼背面缝道处出现聚集和轨迹交汇现象,聚集的水滴撞击在缝翼背面后缘,致使当地的水滴收集率出现了严重的激增,这一现象对飞机防冰系统设计提出了新的要求;另外,三维效应对增升构型的撞击特性影响显著,其中主翼下翼面水滴收集率沿展向分布存在较大差异。     

基于水膜流动与耦合传热的翼型防/除冰计算

针对飞行器防/除冰过程中翼面上空气-水膜-冰层-机翼之间的耦合传质传热现象,建立了一种基于水膜流动与耦合传热模型的翼型防/除冰数值模拟方法。基于Myers水膜流动模型建立了防/除冰热载荷作用下翼面溢流水流动、积冰及内部温度分布的数值计算理论。对于翼型及冰层内的传热现象,利用焓理论及有限体积法建立了复杂多层结构传热的数值模拟方法,对于冰层相变过程,提出了一种基于焓理论的相变修正方法以考虑相变潜热对温度变化的影响。最终实现了翼型防/除冰过程的耦合计算,结果表明:通过结合不同界面处的传热边界条件和考虑了相变潜热效应的焓理论对水膜流动与翼型/冰层传热模型进行耦合求解,能够对翼型/冰层内温度分布进行准确计算,可实现对翼型防/除冰过程中溢流水流动及积冰特性的有效预测与分析。      

基于内外传热耦合的热气防冰系统仿真计算

基于内外传热耦合原理,建立了热气防冰系统的AMESim仿真模型,研究热气防冰系统在不同引气状态下,管路流量、压力及蒙皮温度的变化规律.采用所建立的仿真方法,计算飞机机翼热气防冰系统的内部流动特性和内外耦合传热特性,将计算结果与热气防冰系统流量分配试验结果进行对比.结果表明:管路流量、压力和蒙皮温度的仿真计算结果与试验结果最大误差为9%,验证了该仿真模型的正确性.在此基础上,对飞机短舱热气防冰系统进行了仿真,分析了飞行包线下系统内外的换热、温度变化、加热效率等关键参数的瞬态特性.仿真结果为热气防冰系统的设计、分析与优化提供依据.      

航空发动机帽罩热气膜防冰的加热特性

针对复合材料帽罩采用的一种热气膜防冰系统,采用计算流体力学方法研究了其内部冲击换热和外部气膜加热效果。发展了二维轴对称计算方法并与全三维计算结果进行了对比,搭建了热气吹风实验台,利用红外测温仪和热电偶测量帽罩温度,验证了数值计算方法。全面研究了不同射流雷诺数、相对冲击距和气膜缝结构参数对热气加热特性的影响,结果表明:增大射流雷诺数有利于提高内部冲击换热效果和外部气膜加热效果;最佳相对冲击距随射流雷诺数的增大而增大,射流雷诺数为10000~40000范围内最佳相对冲击距在5~8内取得;气膜射流能大幅提高缝后外壁面温度,加热效率与气膜缝宽和位置有关,前缘开设的气膜缝还会提升前缘冲击换热效果。      

环境参数对飞机防冰热载荷的影响规律

为在飞机热气防冰系统设计过程中准确计算防护表面的防冰热载荷,以某飞机翼段模型为例,对防冰热载荷计算开展了研究,给出了计算流程,研究了来流速度、高度、温度、液态水含量和水滴直径等环境参数变化对翼段表面热载荷分布规律的影响。研究结果表明:①热载荷主要分布在机翼前缘,且越靠近翼尖部分热载荷数值越大,热载荷曲线呈现中间低两边高的分布;②表面热载荷分布范围和量级与来流速度和水滴直径呈正相关关系;飞行高度的变化对表面热载荷分布的影响较小;来流温度降低或液态水含量增加,热载荷数值增大但分布范围不变;水滴直径增加10 μm,上翼面热载荷分布范围增加约50%。      

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

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

运输机翼型结冰的计算和实验

采用数值计算和结冰风洞实验相结合的手段,对某运输机机翼剖面缩比翼型的结冰特性进行了研究.结冰实验在气动中心0.3m×0.2m结冰风洞中进行,实验模型为弦长0.18m的层流翼型.翼型的绕流流场通过求解低速黏流的时均Navier-Stokes(N-S)方程得到,采用拉格朗日法计算水滴撞击特性,在此基础上,求解结冰热力学模型...     

防冰部件表面流动换热与温度计算分析

在考虑防冰部件湿表面水膜流动与传热相互影响的基础上,提出了一种基于湿表面传热与流动耦合作用的防冰部件表面温度预测方法.应用此方法,计算了一种全蒸发防冰翼型工作时的表面流动换热和温度分布.分析了湿表面上各项热流的作用效果.结果表明:部件表面换热系数对壁面温度的计算最为敏感,影响很大:水膜流动能吸收并传递加热热量;全蒸发防...