Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.