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

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

作　　者：马乙楗 柴得林 王强 易贤[1,2] 孔满昭 MA Yijian;CHAI Delin;WANG Qiang;YI Xian;KONG Manzhao(Key Laboratory of Icing and Anti/De-icing,China Aerodynamics Research and Development Center,Mianyang 621000,China;State Key Laboratory of Aerodynamics,China Aerodynamics Research and Development Center,Mianyang 621000,China;AVIC The First Aircraft Design Institute,Xi’an 710089,China)

机构地区：[1]中国空气动力研究与发展中心结冰与防除冰重点实验室,绵阳621000 [2]中国空气动力研究与发展中心空气动力学国家重点实验室,绵阳621000 [3]航空工业第一飞机设计研究院,西安710089

出　　处：《航空学报》2023年第1期35-46,共12页Acta Aeronautica et Astronautica Sinica

基　　金：国家自然科学基金(12132019,12172372);国家科技重大专项(J2019-III-0010-0054)。

摘　　要：冰晶的运动、相变、黏附现象广泛存在于航空发动机内部结冰过程中,开展相应的计算方法研究,系统分析冰晶运动相变与黏附特性是研究发动机内部结冰、保证飞行安全的迫切需求。针对冰晶运动相变及黏附特性计算,建立了拉格朗日框架下冰晶运动-传热传质耦合的数值计算方法,分析了基于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.

关 键 词：混合相结冰 冰晶形状 融化率 黏附效率 收集系数 

分 类 号：V211.41[航空宇航科学与技术—航空宇航推进理论与工程]