An enthalpy-based model for the physics of ice crystal icing

Timothy Peters,Josh Shelton,Hui Tang,Philippe Trinh

2024-03-29

Abstract:Ice crystal icing (ICI) in aircraft engines is a major threat to flight safety. Due to the complex thermodynamic and phase-change conditions involved in ICI, rigorous modelling of the accretion process remains limited. The present study proposes a novel modelling approach based on the physically-observed mixed-phase nature of the accretion layers. The mathematical model, which is derived from the enthalpy change after accretion (the enthalpy model), is compared to an existing pure-phase layer model (the three-layer model). Scaling laws and asymptotic solutions are developed for both models. The onset of ice accretion, the icing layer thickness, and solid ice fraction within the layer are determined by a set of non-dimensional parameters including the Peclet number, the Stefan number, the Biot number, the Melt Ratio, and the evaporative rate. Thresholds for freezing and non-freezing conditions are developed. The asymptotic solutions presents good agreement with numerical solutions at low Peclet numbers. Both the asymptotic and numerical solutions show that, when compared to the three-layer model, the enthalpy model presents a thicker icing layer and a thicker water layer above the substrate due to mixed-phased features and modified Stefan conditions. Modelling in terms of the enthalpy poses significant advantages in the development of numerical methods to complex three-dimensional geometrical and flow configurations. These results improve understanding of the accretion process and provide a novel, rigorous mathematical framework for accurate modelling of ICI.

Fluid Dynamics

What problem does this paper attempt to address?

The main focus of this paper is the phenomenon of ice crystal icing (ICI) in aircraft engines, which poses a significant threat to flight safety. The current modeling methods are still limited in understanding the ice crystal icing process involving complex thermodynamics and phase change conditions. The study proposes a new modeling approach called the enthalpy model, based on the physically observed mixed-phase accretion layer, and compares it with the existing pure-phase layer model (three-layer model). By developing scaling laws and asymptotic solutions, the study determines non-dimensional parameters such as the Péclet number, Stefan number, Bi number, melting rate, and evaporation rate, which define the ice accretion initiation, ice layer thickness, and solid ice proportion within the layer. The asymptotic solutions of the model show good agreement with numerical solutions at low Péclet numbers, indicating that the enthalpy model exhibits thicker ice and water layers compared to the three-layer model, due to the characteristics of mixed-phase and modified Stefan conditions. The enthalpy model has significant advantages in the development of numerical methods for complex three-dimensional geometries and flow configurations, improving the understanding of the accretion process and providing a new mathematical framework for accurate modeling of ice crystal icing.