A permeable wall boundary condition is proposed to depict droplet impingement.
Critical ice thickness is derived as a criterion for determining existence of overflow.
Inner time step marching procedure immediately updates icing properties.
Direct 3-D simulation has reflected spanwise flow effect properly.
Effects of meterological parameters have been investigated and analyzed.