机翼混合相冰晶结冰现象的数值研究

To explore the characteristics of mixed-phase icing when ice crystals and supercooled water droplets coexist, a numerical study is conducted utilizing the FENSAPTCE simulation platform. The NACA0012 airfoil is selected as the subject of the study. The motion of non-spherical ice crystals is adjusted by modifying the drag coefficient, and both ice crystal adhesion and erosion effects arc taken into consideration. Within the Euler framework, a comprehensive numerical calculation method is developed to simulate the complete physical process of ice crystal motion, collision, adhesion, accretion, and erosion. The influence of ice water content (IWC), liquid water content (LWC), as well as ice crystal adhesion and erosion effects on the shape and thickness of icing is analyzed. The applicability conditions of two ice crystal adhesion models, namely NTI and NRC, arc compared. The research findings suggest that, when maintaining a constant total water content (TWO, an increase in the ratio of IWC to LWC leads to a gradual reduction in the ice coverage area. Additionally, the maximum ice thickness initially increases and then decreases, accompanied by a transition in ice shape from crown to angular. The maximum ice thickness at the stagnation point is achieved when IWC and LWC arc 0. 4 g • nC and 1. 0 g • rrC5, respectively. Under high IWC conditions, the NRC adhesion model demonstrates a stagnation point ice thickness that closely aligns with experimental results, with an error of 7. 2%. Conversely, under low IWC conditions, both the NTI and NRC models yield similar stagnation point ice thicknesses. The erosion effect primarily affects the region near the airfoil's stagnation point, and an increase in IWC or LWC intensifies this erosion effect. This study provides a theoretical foundation for the design of an efficient anti-icing system for aircrafts.