Numerical simulation for in-cloud icing of three-dimensional wind turbine blades

Icing will occur when wind turbines are installed in cold and humid regions. Blade icing will reduce the efficiency of wind turbines. In this paper, a numerical simulation method is proposed for glaze, mixed, and rime ice accretion on three-dimensional wind turbine blades. The air flow field is calculated by the multiple reference frame model and droplet trajectories are obtained using the Eulerian two-phase flow method. In the ice accretion calculation, a three-dimensional icing calculation model is established and a new water flow distribution method based on shear force and centrifugal force is presented. The calculation results of ice shapes are compared with the experimental results and the good agreement proves the validity of the icing numerical simulation method. Based on the calculation results, the influence of centrifugal force on the ice shape is studied. Finally, the droplet impingement characteristics and ice shapes of the wind turbine blade are simulated. The results show that the in-cloud icing of wind turbine blades occurs mainly near the leading edge. The icing thickness is thicker and the icing range is wider with the increase of spanwise length.

[1]  Wei Dong,et al.  Calculation and Analysis of Runback Water Flow on Anti-Icing Airfoil Surface , 2016 .

[2]  Guy Fortin,et al.  Wind Turbine Icing and De-Icing , 2009 .

[3]  William Wright.  User's Manual for LEWICE Version 3.2 , 2008 .

[4]  Danmei Hu,et al.  A study on stall-delay for horizontal axis wind turbine , 2006 .

[5]  Robert P. Cassoni,et al.  Aircraft Anti-Icing and De-Icing Techniques and Modeling , 1996 .

[6]  Wagdi G. Habashi,et al.  FENSAP-ICE Simulation of Icing on Wind Turbine Blades, Part 1: Performance Degradation , 2013 .

[7]  Peter Frohboese,et al.  Effects of Icing on Wind Turbine Fatigue Loads , 2007 .

[8]  A. Rothmayer,et al.  Water film runback on an airfoil surface , 2000 .

[9]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[10]  L. Battisti,et al.  Influence of the Type and Size of Wind Turbines on Anti-Icing Thermal Power Requirements for Blades , 2007 .

[11]  Gary A. Ruff,et al.  Users Manual for the NASA Lewis Ice Accretion Prediction Code (LEWICE) , 1990 .

[12]  Sven Schmitz,et al.  Scaled ice accretion experiments on a rotating wind turbine blade , 2012 .

[13]  Wagdi G. Habashi,et al.  Three-Dimensional Eulerian Approach to Droplet Impingement Simulation Using FENSAP-ICE, Part 1: Model, Algorithm, and Validation , 2000 .

[14]  Qin Hu,et al.  Study on small wind turbine icing and its performance , 2017 .

[15]  Timo Laakso,et al.  Modelling and Prevention of Ice Accretion on Wind Turbines , 2001 .

[16]  P. A. Costa Rocha,et al.  k–ω SST (shear stress transport) turbulence model calibration: A case study on a small scale horizontal axis wind turbine , 2014 .

[17]  Xin Wang,et al.  Experimental Investigation of Energy Losses due to Icing of a Wind Turbine , 2007 .

[18]  Yonghui Xie,et al.  Effects of near-wall grid spacing on SST-K-ω model using NREL Phase VI horizontal axis wind turbine , 2012 .

[20]  Edward,et al.  Adverse Environment Rotor Test Stand Calibration Procedures and Ice Shape Correlation , 2010 .

[21]  E. Luke,et al.  Robust and Accurate Eulerian Multiphase Simulations of Icing Collection Efficiency Using Singularity Diffusion Model , 2010 .

[22]  L. Fingersh.  Unsteady Aerodynamics Experiment , 2001 .