Icing is a popular research area in wind energy, and the icing problem of the supercooled droplet–ice crystal mixed-phase condition is one of the new challenges. A numerical method for analyzing the icing characteristics of wind turbine airfoil under mixed-phase conditions is presented. The control equations for the dynamics of supercooled droplets and ice crystals are formulated using the Lagrangian method. Equations for the conservation of mass and energy during the icing process involving supercooled droplets and ice crystals are constructed. The impact of erosion phenomena on the mixed-phase icing process is examined, and methodologies for solving the control equations are introduced. The numerical method is utilized for modeling mixed-phase icing under a range of conditions. The results of these simulations are then compared with data obtained from icing wind tunnel tests to assess the validity of the method. The influence of various mixed-phase conditions on ice shapes is studied. It is found that higher icing temperatures correspond to a larger icing range and amount. The increase in supercooled droplet content, ice crystal content, and ice crystal diameter all contribute to enhanced ice accretion. However, the effects of ice crystal content and diameter are relatively minor.