Numerical simulation methods for airfoil anti-icing/deicing system based on water film flow and coupled heat transfer model were established to solve the coupled mass and heat transfer phenomena among air/water film/ice layer/airfoil on the airfoil surface in the aircraft anti-icing/deicing process. The numerical computation methods for runback flow water

ice accretion on the airfoil surface as well as inner temperature distribution under the effect of the thermal load of anti-icing/deicing system were established based on the Myers water film flow model. For the heat transfer process in the airfoil and ice layer

the enthalpy theory and the finite volume method were used to establish the numerical simulation method of heat transfer in the complex multilayer structure. For the phase change process in the ice layer

the phase change correction method based on enthalpy theory was proposed to consider the effect of phase change latent heat on temperature change. Finally

coupled calculation of the airfoil anti-icing/deicing process was realized. The results showed that through the coupled calculation of the water film flow and airfoil/ice layer heat transfer model with combination of the heat transfer boundary conditions at different interfaces and the enthalpy theory considering the latent heat effect of phase change

the temperature distribution can be accurately computed in the airfoil/ice layer

meanwhile the effective prediction and analysis of the flow and ice accretion characteristics of the runback water during the airfoil anti-icing/deicing process can be realized.