Numerical simulation of two dimensional ice accretion based on lubrication theory

Water film-ice layer model based on lubrication theory was extended to two-dimensional orthogonal curvilinear coordinate system on arbitrary sections. Two governing partial differential equations of water film flow and ice accretion on curved sections were obtained. The numerical methods for solving these unsteady governing equations were presented. An implicit-explicit discrete scheme was applied to obtain the algebraic form of governing equations. Examples of numerical icing calculation were verified on NACA0012 airfoil and circle cylinder section in typical aero icing and structural icing environment. The results of calculation using current model were compared with ice shapes of simulation using traditional Messinger model and results from ice wind tunnel tests. Ice shapes on airfoil obtained using this numerical method were close to icing curves calculated by Messinger model in low temperature glaze icing conditions. In relatively high temperature icing conditions, more accurate results were obtained than Messinger model compared with results of ice tunnel experiments. Ice accretion on transmission line wires which cannot be predicted properly by traditional Messinger model was also effectively predicted by current model in structural icing conditions.