In glaze ice conditions, beads on the surface usually grow to form roughness elements through coalescence, finally resulting in enhancement of local collection efficiency. However, the effects of roughness elements due to freezing of beads are not reflected on the local collection efficiency in CFD icing simulations. This is problematic for predicting the resultant ice shape, which may lead to inaccurate aerodynamic performance and load distribution. The aim of this study is to propose a macroscopic icing model which can reflect bead microscopic phenomena using the Eulerian approach. To this end, a correction was made for collection efficiency by introducing a novel parameter - the effective impinging angle- which is the angle to calculate the local collection efficiency depending on the physical state of surface. It is assumed that the parameter related to the contact angle represents the state of beads. The computational icing analysis of airfoil was performed using the proposed model both in the rime condition and glaze conditions. The results show that the icing characteristics in the feather region is captured with enhanced accuracy in both conditions.