To simulate the ice accretion on airfoil,a boundary moving technique is proposed to deal with the distortion of airfoil surface due to ice accretion on the leading edge.Incorporating the two-phase model of air-supercooled droplets in the Eulerian coordinates,this technique was applied to simulating the process of the rime ice accretion (the droplets freeze at the instant impinging on the airfoil) on the NACA 0012 airfoil,and the ice profile after ice accretion was achieved successfully.And the brief comparison between the results of this paper and the experiments indicates that the current method is applicable and effective.