多步法冰晶积冰数值模拟  

Multi⁃step method for numerical simulation of ice crystal icing

作　　者：卜雪琴[1] 黄平 林贵平[1,2] 楼燕霞 Xueqin BU;Ping HUANG;Guiping LIN;Yanxia LOU(School of Aeronautic Science and Engineering,Beihang University,Beijing 100191,China;Hangzhou International Innovation Institute of Beihang University,Hangzhou 311115,China)

机构地区：[1]北京航空航天大学航空科学与工程学院,北京100191 [2]北京航空航天大学杭州国际创新研究院,杭州311115

出　　处：《航空学报》2024年第14期144-157,共14页Acta Aeronautica et Astronautica Sinica

基　　金：国家自然科学基金(52272428);国家重点研发计划(2021YFB2601700);中央高校基本科研业务费(YWF-22-L-732)。

摘　　要：近年来,冰晶结冰问题逐渐成为飞机/发动机防除冰领域研究的热点。为了实现对冰晶积冰过程的准确预测,建立了冰晶结冰计算模型,包括考虑了不规则形状的冰晶运动轨迹模型、冰晶结冰质量动量守恒模型、黏附模型、侵蚀模型。采用时间多步法模拟积冰过程,在每个时间步长内更新空气流场、冰晶粒子运动轨迹、结冰热力学计算结果,并利用动网格技术更新冰形的几何边界。研究结果表明,利用时间多步法进行计算,能够实现稳定冰形的模拟。稳定冰形形成的原因在于:一方面,积冰减少了表面冰晶的黏附量;另一方面,积冰增大了对已形成积冰的侵蚀量,最终使黏附积冰过程与侵蚀过程达到动态平衡。本文模型和计算方法适用于低融化率情况,能够对三维结构冰晶结冰数值模拟提供指导,但对于高融化率情况还需要进一步考虑积冰向下游流走的情况。In recent years,ice crystal icing has emerged as a prominent and challenging research area in aircraft and engine anti-icing.This paper presents a computational model for accurate prediction of the ice crystal icing process.The proposed model encompasses a trajectory model that accounts for irregular ice crystal shapes,an adhesion model,an erosion model,and a thermodynamic model for icing.The multi-step method is employed to simulate the ice accumulation process,involving the updating of the airflow field,ice crystal particle trajectories,and thermodynamic calculations at each time step.Additionally,the dynamic mesh technology is utilized to update the geometric boundaries of the ice formation.The research findings demonstrate that the application of the multi-step method enables stable simulation of ice formation.The stability of the ice shapes is attributed to the combined effects of reduced adhesion on the surface of ice crystals and increased erosion on pre-existing ice,ultimately leading to a dynamic equilibrium between the adhesion and erosion processes.It should be noted that the proposed model and computational approach in this paper are specifically suited for low melting rate conditions,which can guide the precise simulation of three-dimensional structural ice crystal icing in future investigations.While high melting rate scenarios would necessitate consideration of the downstream flow of ice in the future.

关 键 词：冰晶结冰 黏附 侵蚀 多步法 动网格 

分 类 号：V244.1[航空宇航科学与技术—飞行器设计]