Extension to the Myers Model for Calculation of Three-Dimensional Glaze Icing

To simulate three-dimensional glaze ice accretion on a typical aircraft surface, the model of water film flow and ice accretion on a curved surface based on two coupled partial differential equations proposed by Myers et al. is extended in nonorthogonal curvilinear coordinates. Then, the governing equation for water film flow coupling with a Stefan model for glaze icing is built up. An efficient implicit-explicit numerical method is proposed for solving the governing equations. The glaze ice computations on two-dimensional airfoils and three-dimensional typical geometries in in-flight icing conditions are completely accomplished. Validation results of two-dimensional calculation cases with experimental ice shapes are presented. Three-dimensional simulations of glaze ice accretion on a sphere and a swept wing are successfully carried out and compared with ice profiles from icing tunnel tests and simulations based on the classical Messinger model.