An improved numerical method of simulating ice accretion on airfoil was developed based on parcel concept following the fundamental assumption of droplet movement.