航空学报 > 2023, Vol. 44 Issue (1): 627817-627817   doi: 10.7527/S1000-6893.2022.27817

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

马乙楗1,2, 柴得林1, 王强1,2, 易贤1,2(), 孔满昭3   

  1. 1.中国空气动力研究与发展中心 结冰与防除冰重点实验室,绵阳 621000
    2.中国空气动力研究与发展中心 空气动力学国家重点实验室,绵阳 621000
    3.航空工业第一飞机设计研究院,西安 710089
  • 收稿日期:2022-07-15 修回日期:2022-08-01 接受日期:2022-08-11 出版日期:2023-01-15 发布日期:2022-08-17
  • 通讯作者: 易贤 E-mail:yixian_2000@163.com
  • 基金资助:
    国家自然科学基金(12132019);国家科技重大专项(J2019-III-0010-0054)

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

Yijian MA1,2, Delin CHAI1, Qiang WANG1,2, Xian YI1,2(), Manzhao KONG3   

  1. 1.Key Laboratory of Icing and Anti /De-icing,China Aerodynamics Research and Development Center,Mianyang 621000,China
    2.State Key Laboratory of Aerodynamics,China Aerodynamics Research and Development Center,Mianyang 621000,China
    3.AVIC The First Aircraft Design Institute,Xi ’an 710089,China
  • Received:2022-07-15 Revised:2022-08-01 Accepted:2022-08-11 Online:2023-01-15 Published:2022-08-17
  • Contact: Xian YI E-mail:yixian_2000@163.com
  • Supported by:
    National Natural Science Foundation of China(12132019);National Science and Technology Major Project(J2019-III-0010-0054)

摘要:

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

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

Abstract:

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.

Key words: mixed phase icing, ice crystal shape, melting ratio, sticking efficiency, collection coefficient

中图分类号: