Abstract: Based on the parcel concept, a modified numerical method for simulating ice accretion on aircraft wing is presented in this study, which obeys the fundamental assumption of Lagrangian and Euler method. The method overcomes the shortages of Lagrangian method, such as large amount of computation work, low efficiency and the matching between droplet space distribution and grid density. NACA0012 rime ice accretion, local droplet collection and local impingement efficiency at 0° and 4° AOA (Angle of Attack) is simulated to validate the method. The computation results show that the modified method based on parcel concept is feasible and robust.