For wing ice accretion problem,Euler method is used in calculating the air-water two-phase flow field. With the help of the Fluent’s user defined function, the water local collection coefficient and ice accretion on three dimensional wings are solved. During the computation process, water runback is divided along the chordwise and spanwise direction and the classical Messinger model is modified to fit the three dimensional airfoil; the analysis of three dimensional mass and heat transfer is also investigated according to this new model. Based on the traditional method of calculating heat transfer coefficient on two dimensional surfaces, an improved method is proposed to compute the heat transfer coefficient on three dimensional wings. Ice accretion is along the direction normal to the surface. According to the comparison of the calculated results with test results of icing tunnel abroad and those of software LEWICE, the computation method in this article is proved to be effective.