航空学报 > 2020, Vol. 41 Issue (12): 124085-124085   doi: 10.7527/S1000-6893.2020.24085

二维机翼混合相结冰数值模拟

卜雪琴, 李皓, 黄平, 林贵平   

  1. 北京航空航天大学 航空科学与工程学院, 北京 100083
  • 收稿日期:2020-04-13 修回日期:2020-04-28 发布日期:2020-05-14
  • 通讯作者: 卜雪琴 E-mail:buxueqin@buaa.edu.cn
  • 基金资助:
    国家科技重大专项(2017-VIII-0003-0114);中央高校基本科研业务费专项资金(YWF-20-BJ-J-732)

Numerical simulation of mixed phase icing on two-dimensional airfoil

BU Xueqin, LI Hao, HUANG Ping, LIN Guiping   

  1. School of Aeronautic Science and Engineering, Beihang University, Beijing 100083, China
  • Received:2020-04-13 Revised:2020-04-28 Published:2020-05-14
  • Supported by:
    National Science and Technology Major Project of China (2017-VIII-0003-0114); the Fundamental Research Funds for the Central Universities (YWF-20-BJ-J-732)

摘要: 近年来冰晶引发的飞机结冰问题逐渐引起人们的重视。针对冰晶和过冷水滴同时存在的混合相结冰问题,通过数值模拟的方法实现结冰冰形预测。空气流场和对流换热的计算中采用了转捩剪切应力输运(SST)湍流模型,基于欧拉法计算冰晶和水滴收集系数。在Messinger结冰热力学模型基础上进行扩展,分析了二维结冰部件表面在混合相气象条件下的传热传质过程,建立了适用于混合相的结冰热力学模型,同时考虑冰晶的黏附效应,添加黏附效率经验公式。利用FLUENT的用户自定义函数(UDF)编程求解混合相热力学模型,计算了霜冰和明冰条件下NACA0012翼型表面结冰情况,与国外风洞试验结果进行比较,验证了计算模型和方法的有效性。结果表明,冰晶黏附效应对混合相结冰量及冰形有很大影响,明冰条件结冰形状偏向楔形。

关键词: 混合相, 结冰, 二维, 冰晶黏附, 机翼

Abstract: In recent years, the problem of aircraft icing caused by ice crystals has attracted wide attention. For the mixed phase (ice crystals and supercooled droplets coexisting) icing problem, the numerical simulation method is used to predict the ice shape. The transition Shear-Stress Transport (SST)turbulence model is adopted in the calculation of air flow field and convective heat transfer. The ice crystal and droplet collection coefficients are calculated based on the Eulerian method. Based on the classical Messinger icing model, the heat and mass transfer process on a two-dimensional surface under the mixed phase icing condition is analyzed, and the icing thermodynamic model suitable for the mixed phase icing is established. Meanwhile, considering the adhesion effect of ice crystals, an empirical formula of adhesion efficiency is added. The mixed phase icing thermodynamic model is solved based on FLUENT and its User Defined Function (UDF). The rime ice and glaze ice conditions are simulated on NACA0012 airfoil. The computational results are compared with foreign wind tunnel experimental results and show agreement. The results show that the effect of ice crystal adhesion has a significant influence on the amount and the shape of ice. Ice accumulation tends to be wedge-shaped under glaze ice condition.

Key words: mixed phase, icing, two-dimensional, ice crystal adhesion, airfoil

中图分类号: