Water film will be initiated and formed under the continuous impingement of incoming supercooled droplets, and directly affect the ice formation, which is always vital for the performance and safety of aero-engine. A hybrid algorithm integrating Discrete Phase Model (DPM) and Volume of Fluid (VOF) was established to simulate all the stages from droplet initiation to the film well-formed which can make the numerical simulation of icing more realistic and accurate. The transformation criterion from the particle to liquid volume fraction and the distribution of source terms has been improved compared with similar previous approaches. Three verification cases, namely single droplet train impact, staggered droplets impact and two glycerine-solution droplet impact, were conducted to verify the mass conservation and the applicability to both the sparse and dense two-phase flow. The results showed that the relative error of mass transformation is no more than 0.5%. The comparison between simulation and experiment result demonstrated that the hybrid algorithm has the ability to simulate the flow behaviour of droplets accurately. The thin water film formation on a typical wing of NACA0012 airfoil was then analysed, the results of which illustrated that the integrated algorithm could successfully capture the characteristics of liquid film formation and motion.