An approach to numerically simulate ice accretion on 2D sections of a wind turbine blade is presented. The method uses standard commercial ANSYS-Fluent and Matlab tools. The Euler-Euler formulation is used to calculate the water impingement on the airfoil, and a UDF (Used Defined Function) has been devised to turn the airfoil's solid wall into a permeable boundary. Mayer's thermodynamic model is implemented in Matlab for computing ice thickness and for updating the airfoil contour. A journal file is executed to systematize the procedure: meshing, droplet trajectory calculation, thermodynamic model application for computing ice accretion, and the updating of airfoil contours. The proposed ice prediction strategy has been validated using iced airfoil contours obtained experimentally in the AMIL refrigerated wind tunnel (Anti-icing Materials International Laboratory). Finally, a numerical prediction method has been generated for anti-icing assessment, and its results compared with data obtained in this laboratory. (C) 2016 Elsevier Ltd. All rights reserved.