Numerical Investigation of the Aeroelastic Behavior of a Wind Turbine with Iced Blades

Wind turbines installed in cold-climate regions are prone to the risks of ice accumulation which affects their aeroelastic behavior. The studies carried out on this topic so far considered icing in a few sections of the blade, mostly located in the outer part of the blade, and their influence on the loads and power production of the turbine are only analyzed. The knowledge about the influence of icing in different locations of the blade and asymmetrical icing of the blades on loads, power, and vibration behavior of the turbine is still not matured. To improve this knowledge, multiple simulation cases are needed to run with different ice accumulations on the blade considering structural and aerodynamic property changes due to ice. Such simulations can be easily run by automating the ice shape creation on aerofoil sections and two-dimensional (2-D) Computational Fluid Dynamics (CFD) analysis of those sections. The current work proposes such methodology and it is illustrated on the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine model. The influence of symmetrical icing in different locations of the blade and asymmetrical icing of the blade assembly is analyzed on the turbine’s dynamic behavior using the aeroelastic computer-aided engineering tool FAST. The outer third of the blade produces about 50% of the turbine’s total power and severe icing in this part of the blade reduces power output and aeroelastic damping of the blade’s flapwise vibration modes. The increase in blade mass due to ice reduces its natural frequencies which can be extracted from the vibration responses of the turbine operating under turbulent wind conditions. Symmetrical icing of the blades reduces loads acting on the turbine components, whereas asymmetrical icing of the blades induces loads and vibrations in the tower, hub, and nacelle assembly at a frequency synchronous to rotational speed of the turbine.

[1]  Ma Dongli,et al.  Effects of relative thickness on aerodynamic characteristics of airfoil at a low Reynolds number , 2015 .

[2]  William B. Wright,et al.  User Manual for the NASA Glenn Ice Accretion Code LEWICE: Version 2.0 , 1999 .

[3]  Jason Jonkman,et al.  FAST User's Guide - Updated August 2005 , 2005 .

[4]  Torgeir Moan,et al.  Wind turbine aerodynamic response under atmospheric icing conditions , 2014 .

[5]  J. Astolfi,et al.  Computational and experimental investigation of flow over a transient pitching hydrofoil , 2009 .

[6]  David Kennedy,et al.  Numerical aerodynamic simulations of a NACA airfoil using CFD with block-iterative coupling and turbulence modelling , 2012 .

[7]  J. Jonkman,et al.  Definition of a 5-MW Reference Wind Turbine for Offshore System Development , 2009 .

[8]  Corrado Groth,et al.  Reliable mesh morphing approach to handle icing simulations on complex models , 2014 .

[9]  Brouwers,et al.  The Experimental Investigation of a Rotor Hover Icing Model with Shedding , 2010 .

[10]  C. Lindenburg,et al.  Aero-elastic modelling of the DOWEC 6 MW pre-design in PHATAS , 2003 .

[11]  Per Johan Nicklasson,et al.  Ice sensors for wind turbines , 2006 .

[12]  Per Johan Nicklasson,et al.  Performance losses due to ice accretion for a 5 MW wind turbine , 2012 .

[13]  W. A. Timmer,et al.  Roughness Sensitivity Considerations for Thick Rotor Blade Airfoils , 2003 .

[14]  Ira H. Abbott,et al.  Summary of Airfoil Data , 1945 .

[15]  X Munduate,et al.  Wind Tunnel Tests of Wind Turbine Airfoils at High Reynolds Numbers , 2014 .

[16]  Ismail H. Tuncer,et al.  Ice Accretion Prediction on Wind Turbines and Consequent Power Losses , 2016 .

[17]  Jonghwa Kim,et al.  Study on correlation between wind turbine performance and ice accretion along a blade tip airfoil using CFD , 2018 .

[18]  Ville Lehtomäki,et al.  Overview of cold climate wind energy: challenges, solutions, and future needs , 2016 .

[20]  G. S. Bir,et al.  User's Guide to MBC3: Multi-Blade Coordinate Transformation Code for 3-Bladed Wind Turbine , 2010 .

[21]  H. E. Addy,et al.  CFD Analysis of the Aerodynamics of a Business-Jet Airfoil with Leading-Edge Ice Accretion , 2004 .

[22]  Wagdi G. Habashi,et al.  Development of a Second Generation In-Flight Icing Simulation Code , 2006 .

[23]  Devi Prasad Mishra,et al.  CFD Simulations for the Selection of an Appropriate Blade Profile for Improving Energy Efficiency in Axial Flow Mine Ventilation Fans , 2014, Journal of Sustainable Mining.

[24]  Georgios Alexandros Skrimpas,et al.  Detection of icing on wind turbine blades by means of vibration and power curve analysis: Icing detection in wind turbines , 2016 .

[25]  Jan-Olov Aidanpää,et al.  Identification of ice mass accumulated on wind turbine blades using its natural frequencies , 2018 .

[26]  M. Reggio,et al.  Numerical Study of Flow Around Iced Wind Turbine Airfoil , 2012 .

[27]  Jean Perron,et al.  Wind turbine performance under icing conditions , 2008 .

[28]  W. A. Timmer,et al.  Summary of the Delft University Wind Turbine Dedicated Airfoils , 2003 .

[29]  Jan-Olov Aidanpää,et al.  Influence of Icing on the Modal Behavior of Wind Turbine Blades , 2016 .

[30]  Rune Brincker,et al.  Modal identification of output-only systems using frequency domain decomposition , 2001 .

[31]  Qin Hu,et al.  Study of ice accretion feature and power characteristics of wind turbines at natural icing environment , 2018 .

[32]  Carlo L. Bottasso,et al.  Model-Independent Periodic Stability Analysis of Wind Turbines ∗ , 2015 .

[33]  Petr Straka,et al.  Comparison of several models of the laminar/turbulent transition , 2013 .

[34]  Xingliang Jiang,et al.  3D numerical simulation of aerodynamic performance of iced contaminated wind turbine rotors , 2018 .

[35]  H. E. Addy,et al.  Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation , 2002 .

[36]  Michael S. Selig,et al.  Wind Turbine Performance Under Icing Conditions , 1998 .

[37]  Douvi C. Eleni,et al.  Evaluation of the turbulence models for the simulation of the flow over a National Advisory Committee for Aeronautics (NACA) 0012 airfoil , 2012 .

[39]  W. Kieffer,et al.  CFD study of section characteristics of Formula Mazda race car wings , 2006, Math. Comput. Model..

[40]  Morten Hartvig Hansen,et al.  On the similarity of the Coleman and Lyapunov–Floquet transformations for modal analysis of bladed rotor structures , 2009 .

[41]  Christian Masson,et al.  Atmospheric icing impact on wind turbine production , 2014 .

[42]  Chenxing Hu,et al.  Wind turbines ice distribution and load response under icing conditions , 2017 .

[43]  Simo Rissanen,et al.  Modelling load and vibrations due to iced turbine operation , 2016 .

[44]  F. Menter.  Improved two-equation k-omega turbulence models for aerodynamic flows , 1992 .

[45]  Holger Koss,et al.  Ice Accretion on Wind Turbine Blades , 2013 .

[46]  Tomas Wallenius,et al.  Method for Estimating Wind Turbine Production Losses Due to Icing , 2012 .

[47]  Richard Crossley,et al.  Wind Turbine Blade Design , 2012 .

[48]  Michele De Gennaro,et al.  Wind energy harnessing of the NREL 5 MW reference wind turbine in icing conditions under different operational strategies , 2018 .