翼面结冰过程中的冰晶运动相变与黏附特性

Phase change and adhesion characteristics of ice crystal movements in wing icing

冰晶的运动、相变、黏附现象广泛存在于航空发动机内部结冰过程中,开展相应的计算方法研究,系统分析冰晶运动相变与黏附特性是研究发动机内部结冰、保证飞行安全的迫切需求。针对冰晶运动相变及黏附特性计算,建立了拉格朗日框架下冰晶运动-传热传质耦合的数值计算方法,分析了基于Monte Carlo方法的冰晶撞击与黏附收集系数计算方法,基于NNWICE平台开发了相应的计算程序。以NACA0012翼型为对象,计算研究了扁椭球状、六角平板状和长椭球状3种非球状冰晶颗粒与球状冰晶的运动、传热传质过程,系统分析了来流温度等参数对冰晶黏附特性的影响,得到了冰晶形状对冰晶运动轨迹和融化过程的影响规律及来流总温与液态水含量/总水含量(LWC/TWC)对冰晶黏附特性的影响规律。相关工作为进一步分析混合相结冰数值模拟计算奠定了工作基础。

The movement,phase change and adhesion of ice crystals widely exist in the process of internal icing of aircraft engines.To study the internal icing of the engine and ensure flight safety,we explore the corresponding calculation methods,and systematically analyze the phase transition and adhesion characteristics of ice crystals.In view of the calculation of the phase transition and adhesion characteristics of ice crystals,a numerical calculation method of coupling the ice crystal motion-heat and the mass transfer under the Lagrange framework was established,the calculation method of the impact and adhesion collection coefficient of ice crystals based on the Monte Carlo method proposed,and the corresponding calculation procedure developed based on the NNWICE platform.Taking the NACA0012 airfoil as the research object,we calculated the motion and heat transfer process of three non-spherical ice crystal particles and spherical ice crystals in oblate,hexagonal flat and prolate,systematically analyzing the influence of incoming flow temperature on the adhesion characteristics of ice crystals.The influence of the ice crystal shape on the motion trajectory and melting process of ice crystals,and the effect of the total temperature of the incoming flow and Liquid Water Content/Total Water Content(LWC/TWC)on the adhesion characteristics of ice crystals were obtained.This work lays the foundation for further development of numerical simulation calculations for mixed-phase icing.