Numerical simulations of ice accretion on wind turbine blades: are performance losses due to ice shape or surface roughness?

Ice accretion on wind turbine blades causes both a change in the shape of its sections and an increase in surface roughness. These lead to degraded aerodynamic performances and lower power output. Here, a high-fidelity multi-step method is presented and applied to simulate a 3 h rime icing event on the National Renewable Energy Laboratory 5 MW wind turbine blade. Five sections belonging to the outer half of the blade were considered. Independent time steps were applied to each blade section to obtain detailed ice shapes. The roughness effect on airfoil performance was included in computational fluid dynamics simulations using an equivalent sand-grain approach. The aerodynamic coefficients of the iced sections were computed considering two different roughness heights and extensions along the blade surface. The power curve before and after the icing event was computed according to the Design Load Case 1.1 of the International Electrotechnical Commission. In the icing event under analysis, the decrease in power output strongly depended on wind speed and, in fact, tip speed ratio. Regarding the different roughness heights and extensions along the blade, power losses were qualitatively similar but significantly different in magnitude despite the well-developed ice shapes. It was found that extended roughness regions in the chordwise direction of the blade can become as detrimental as the ice shape itself.