A novel multiphase computational fluid dynamics model based on the coupled Lagrangian-Eulerian approach was developed to assess the performance of the S809 airfoil of a horizontal-axis wind turbine under rainy conditions. The model was utilized for simulations of water film formation over the airfoil. Our simulation showed significant performance loss due to an increase in drag and lift degradation.