A numerical simulation approach is presented to simulate the ice accretion on airfoil which can deal with the ice accretion under various atmospheric conditions.Incorporating the two-phase model of air supercooled droplets in the Eulerian coordinate,this approach was carried out to simulate the ice accretion on the NACA 0012 airfoil under different atmospheric conditions,and the final ice shapes(rime,mixed and glaze ice shapes) were predicted.The comparison between the results of this paper and the experimental data as well as other icing codes indicates that the current approach is applicable and effective.Also the effect of different ice shapes on lift and drag coefficient were assessed.The results indicate that the glaze ice has the most adverse impact on airfoil aerodynamic performance.