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 acurved surface basedon two coupled partial differential equations proposed by Myerset 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 inflight 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.