Predictions of ice formations on wind turbine blades and power production losses due to icing

Prediction of ice shapes on a wind turbine blade makes it possible to estimate the power production losses due to icing. Ice accretion on wind turbine blades is responsible for a significant increase in aerodynamic drag and decrease in aerodynamic lift and may even cause premature flow separation. All these events create power losses and the amount of power loss depends on the severity of icing and the turbine blade profile. The role of critical parameters such as wind speed, temperature, liquid water content on the ice shape, and size is analyzed using an ice accretion prediction methodology coupled with a blade element momentum tool. The predicted ice shapes on various airfoil profiles are validated against the available experimental and numerical data in the literature. The error in predicted rime and glime ice volumes and the maximum ice thicknesses varies between 3% and 25% in comparison with the experimental data depending on the ice type. The current study presents an efficient and accurate numerical methodology to perform an investigation for ice‐induced power losses under various icing conditions on horizontal axis wind turbines. The novelty of the present work resides in a unified and coupled approach that deals with the ice accretion prediction and performance analysis of iced wind turbines. Sectional ice profiles are first predicted along the blade span, where the concurrence of both rime and glaze ice formations may be observed. The power loss is then evaluated under the varying ice profiles along the blade. It is shown that the tool developed may effectively be used in the prediction of power production losses of wind turbines at representative atmospheric icing conditions.     

The relationship between chord length and rime icing on wind turbines

Numerical simulations of ice accumulation on four different wind turbine blade profiles, from 450 kW, 600 kW, 1 MW and 2 MW, fixed speed, stall controlled, wind turbines, were performed to determine how wind turbine size affects atmospheric icing. The simulations indicate that dry rime icing is less severe for larger wind turbines both in terms of local ice mass and in terms of relative ice thickness. Copyright © 2009 John Wiley & Sons, Ltd.     

风力机叶片覆冰状态监测基准值与分级诊断标准研究

以2 MW风力机为研究对象,基于实际风力机状态(SCADA)系统大数据,选取叶片正常状态和覆冰状态下的风速、功率、桨距角和偏航角数据,采用核密度-均值数据处理方法,得到叶片覆冰状态监测基准值及其定量表达式。同时,根据叶片不同覆冰时期桨距角和功率值随风速的变化情况,提出叶片覆冰状态分级诊断标准。应用结果表明,根据桨距角随风速的变化情况可判断在叶片覆冰过程中机组最大功率追踪情况以及气动性能损失情况,根据风速-功率值分布情况可较准确地判别叶片的覆冰状态。     

Performance losses due to ice accretion for a 5 MW wind turbine

A numerical study of power performance losses due to ice accretion on a large horizontal axis wind turbine blade has been carried out using computational fluid dynamics (CFD) and blade element momentum (BEM) calculations for rime ice conditions. The computed aerodynamic coefficients for the normal and iced blades from the CFD calculations were used together with the BEM method to calculate the torque, power and curves of the wind turbine for both normal and icing conditions. The results are compared with the published data. It is shown that icing results in a reduced power production from the turbine and that changing the turbine controller could improve the power production with iced blades. Copyright © 2011 John Wiley & Sons, Ltd.     

多参数模型风电机组叶片结冰监测与预警研究

详细分析叶片结冰对风电机组运行性能和运行参数的影响,采用功率、叶轮转速和环境温度作为监测叶片结冰的变量.采用高斯过程回归分别建立功率模型和叶轮转速模型实现2个参数的实时监测.引入序贯概率比检验方法分析功率和叶轮转速模型的预测残差以发现2个参数在叶片结冰时的异常变化.当风电机组功率异常、叶轮转速异常且环境温度在0℃附近这...     

Real-time rotor effective wind speed estimation based on actuator disc theory: Design and full-scale experimental validation

The use of state estimation techniques offers a means of inferring rotor effective wind speed from standard measurements of wind turbines. Typical wind speed estimators rely upon a pre-computed quasi-steady aerodynamic mapping, which describes the relationship between pitch angle and tip-speed ratio and the power coefficient. In practice, the static mapping does not capture the influence of turbine structural dynamics and atmospheric turbulence, inevitably resulting in poor performance of the wind speed estimation. In addition, the turbine aerodynamic properties might not be easily accessible. Thus, this paper presents a rotor effective wind speed estimation method that obviates the requirement for prior knowledge of turbine power coefficients. Specifically, the proposed method exploits a simple actuator disc model, where the aerodynamic power and thrust coefficients can be characterized in terms of axial induction factors. Based on this insight and standard turbine measurements, real-time estimation of rotor effective wind speed and axial induction factors can then be achieved using a simplified turbine drive-train model and an extended Kalman filter. In addition, the actuator disc model can be updated easily over time by calibrating solely two correction factors. Thus, the proposed algorithm presents an alternative for estimating the rotor effective wind speed, which is valuable for numerous applications, for example, LiDAR-assisted control and coherence studies.     

The effect of a leading edge erosion shield on the aerodynamic performance of a wind turbine blade

Armour EDGE is a novel shield developed to protect the leading edge of wind turbine blades from erosion. The aerodynamic impact on aerofoils of National Renewable Energy Laboratory (NREL) 5MW wind turbine has been investigated using 2D fully turbulent computational fluid dynamics (CFD), with three profiles at critical locations along the blade simulated both with and without the shield to compare aerodynamic performance. Two wind speeds were investigated that reflect regular operating conditions: at rated speed of 11.4 m/s and a below rated speed of 7 m/s. The results showed that the presence of the shield during rated wind speed reduced the drag by as much as 4.5%, where the lift‐to‐drag ratio increased by a maximum of 4%. At the below rated wind speeds, the shield had negligible impact on the performance of all but one National Advisory Committee for Aeronautics (NACA) 64‐618 profile, which resulted in an increase in the drag coefficient of 7%. It was also found that the suction side of the aerofoil is much more sensitive to leading edge protection placement than the pressure side. It was concluded that the erosion shield as a method of leading edge protection, with a gradual transition from shield to blade, will not have a major impact on the aerodynamic performance of a multi‐megawatt wind turbine blade and could slightly increase aerofoil efficiency at high wind speeds.     

Ice detection on wind turbines using the observed power curve

Icing on the blades of a wind turbine can lead to significant production losses during the winter months for wind parks in cold climate regions. However, there is no standard way of identifying ice‐induced power loss. This paper describes three methods for creating power threshold curves that can be used to separate iced production periods from non‐iced production periods. The first approach relies on a percentage deviation from the manufacturer's power curve. The other two approaches fit threshold curves based on the observed variance of non‐iced production data. These approaches are applied to turbines in four wind parks and compared with each other and to observations of icing on the nacelle of one of the turbines in each park. It is found that setting an ice threshold curve using 0.1 quantile of the observed power data during normal operation with a 2‐h minimum duration is the best approach for icing identification. The quantile should be fit based on at least 1year of data, and a smoothing function should be applied to the quantile results to remove any outliers caused by limited numbers of data points. Copyright © 2015 John Wiley & Sons, Ltd.     

Sensitivity of southern California wind energy to turbine characteristics

Using output from a high‐resolution meteorological simulation, we evaluate the sensitivity of southern California wind energy generation to variations in key characteristics of current wind turbines. These characteristics include hub height, rotor diameter and rated power, and depend on turbine make and model. They shape the turbine's power curve and thus have large implications for the energy generation capacity of wind farms. For each characteristic, we find complex and substantial geographical variations in the sensitivity of energy generation. However, the sensitivity associated with each characteristic can be predicted by a single corresponding climate statistic, greatly simplifying understanding of the relationship between climate and turbine optimization for energy production. In the case of the sensitivity to rotor diameter, the change in energy output per unit change in rotor diameter at any location is directly proportional to the weighted average wind speed between the cut‐in speed and the rated speed. The sensitivity to rated power variations is likewise captured by the percent of the wind speed distribution between the turbines rated and cut‐out speeds. Finally, the sensitivity to hub height is proportional to lower atmospheric wind shear. Using a wind turbine component cost model, we also evaluate energy output increase per dollar investment in each turbine characteristic. We find that rotor diameter increases typically provide a much larger wind energy boost per dollar invested, although there are some zones where investment in the other two characteristics is competitive. Our study underscores the need for joint analysis of regional climate, turbine engineering and economic modeling to optimize wind energy production. Copyright © 2012 John Wiley & Sons, Ltd.     

An investigation of variable power collective pitch control for load mitigation of floating offshore wind turbines

This paper investigates the loads on offshore floating wind turbines and a new control method that can be used to reduce these loads. In this variable power collective pitch control method, the rated generator speed, which is the set point that the collective pitch control attempts to drive the actual generator speed towards, is no longer a constant value but instead is a variable that depends on the platform pitch velocity. At a basic physical level, this controller achieves the following: as the rotor of a floating turbine pitches upwind, the controller adjusts so as to extract more energy from the wind by increasing the rated generator speed and thus damps the motion; as the rotor pitches downwind, less energy is extracted because the controller reduces the rated generator speed and again damps the motion. This method is applied to the NREL 5 MW wind turbine model, in above rated conditions where the platform motion is most problematic. The results indicate significant load reductions on key structural components, at the expense of minor increases in power and speed variability. The loads on the blades and tower are investigated more generally, and simple dynamic models are used to gain insight into the behavior of floating wind turbine systems. It is clear that for this particular design, aerodynamic methods for reducing platform motion and tower loads are likely inadequate to allow for a viable design, so new designs or possibly new control degrees of freedom are needed. Copyright © 2012 John Wiley & Sons, Ltd.     

Investigating the influence of dimensional scaling on aerodynamic characteristics of wind turbine using CFD simulation

Wind turbines are used in a variety of applications with different performance requirements. Investigating the influence of scaling on wind turbine characteristics can pave the way to utilize the experience gained from a smaller turbine for a larger one. In this paper, the effects of wind turbine size on aerodynamic characteristics of a rotor blade are examined using CFD simulation. NREL phase VI wind turbine rotor was simulated in order to validate the results and ensure the accuracy of the CFD model. A 2 MW wind turbine was then chosen as a large turbine and a scaled down model of its rotor was simulated numerically. The results of the simulation were introduced to Similarity Theory relations in order to predict the aerodynamic characteristics of the 2 MW wind turbine. The 2 MW turbine was also simulated and the results of the simulation were compared to predictions of Similarity Theory. It was observed that the results of the simulation completely follow the values predicted by Similarity Theory. Both Similarity Theory predictions and simulation results demonstrated that the torque increases with the cube of change in rotor diameter whereas the thrust value and aerodynamic forces grow with the square of change in diameter.     

Modeling,simulation and control of a wind turbine with a hydraulic transmission system

A complete mathematical model of a hydraulic transmission concept for use in wind turbines is presented. The hydraulic system transfers the power from the nacelle to ground level. The main focus has been to develop a model that takes into account the most important dynamics affecting the wind turbine and the hydraulic transmission system involved, such that the model can be used to analyze the dynamic feasibility of a hydraulic transmission concept. Further, dynamic analysis of a hydraulic transmission system for wind turbines is investigated. The nonlinear dynamic model is developed in MATLAB Simulink. Analytical calculation of natural periods of a linearized model corresponds well with simulations of the overall system. A valve control system is proposed to reduce pressure and power fluctuations at operation both below and above the rated wind speed for the wind turbine. Further, a blade pitch control system based on an aerodynamic power estimator is proposed for operation above the rated wind speed. System simulations for one case below and one case above the rated wind speed show that the dynamic response of the overall system is stable and that the wind turbine variables are within typical ranges for conventional variable speed wind turbines with mechanical transmission. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of wind turbine wake interaction using the vorticity transport model

The aerodynamic interactions that can occur within a wind farm can result in the constituent turbines generating a lower power output than would be possible if each of the turbines were operated in isolation. Tightening of the constraints on the siting of wind farms is likely to increase the scale of the problem in the future. The aerodynamic performance of turbine rotors and the mechanisms that couple the fluid dynamics of multiple rotors can be most readily understood by simplifying the problem and considering the interaction between only two rotors. The aerodynamic interaction between two rotors in both co‐axial and offset configurations has been simulated using the Vorticity Transport Model. The aerodynamic interaction is a function of the tip speed ratio, and both the streamwise and crosswind separation between the rotors. The simulations show that the momentum deficit at a turbine operating within the wake developed by the rotor of a second turbine is governed by the development of instabilities within the wake of the upwind rotor, and the ensuing structure of the wake as it impinges on the downwind rotor. If the wind farm configuration or wind conditions are such that a turbine rotor is subject to partial impingement by the wake produced by an upstream turbine, then significant unsteadiness in the aerodynamic loading on the rotor blades of the downwind turbine can result, and this unsteadiness can have considerable implications for the fatigue life of the blade structure and rotor hub. Copyright © 2009 John Wiley & Sons, Ltd.     

Phases of icing on wind turbine blades characterized by ice accumulation

Icing experiments on wind turbine blade profiles have been performed at the University of Manitoba Icing Tunnel Facility to facilitate a greater understanding of the mechanisms involved in the icing process for wind turbines exposed to cold climates. Blade icing results in the degradation of power performance and is a critical issue for the optimization of power performance and safe operation of wind turbines. Accumulation rate, the amount of ice that accumulates at the leading edge of the blade profile as a function of time, provides a characteristic measurement that can be used to classify the phases of icing in an icing event and further identify the severity of potential problems arising as a result of ice accumulation on wind turbine blades. To control this characteristic, the mitigation strategies that were employed involved coatings, heat treatments and the combination thereof, in both glaze and rime icing regimes. By understanding the icing process and its characteristic behavior to non-mitigated and mitigated scenarios, the phases of icing of both circumstances may be defined. This paper documents the data recorded from the experimental icing event and provides results of the comparative behavior of the icing mitigation strategies and extends this understanding to define the phases of icing on wind turbine blades.     

An investigation of variable power collective pitch control for load mitigation of floating offshore wind turbines

This paper investigates the loads on offshore floating wind turbines and a new control method that can be used to reduce these loads. In this variable power collective pitch control method, the rated generator speed, which is the set point that the collective pitch control attempts to drive the actual generator speed towards, is no longer a constant value but instead a variable that depends on the platform pitch velocity. At a basic physical level, this controller achieves the following: as the rotor of a floating turbine pitches upwind, the controller adjusts so as to extract more energy from the wind by increasing the rated generator speed and thus damps the motion; as the rotor pitches downwind, less energy is extracted because the controller reduces the rated generator speed and again damps the motion. This method is applied to the NREL 5 MW wind turbine model, in above‐rated conditions where the platform motion is most problematic. The results indicate significant load reductions on key structural components, at the expense of minor increases in power and speed variability. The loads on the blades and tower are investigated more generally, and simple dynamic models are used to gain insight into the behavior of floating wind turbine systems. It is clear that for this particular design, aerodynamic methods for reducing platform motion and tower loads are likely inadequate to allow for a viable design, and so new designs or possibly new control degrees of freedom are needed. Copyright © 2012 John Wiley & Sons, Ltd.     

Wind turbine performance under icing conditions

The wind energy market is in full growth in Quebec but technical difficulties due to cold climate conditions have occurred for most of the existing projects. Thus, icing simulations were carried out on a 0.2 m NACA 63 415 blade profile in the refrigerated wind tunnel of the Anti‐icing Materials International Laboratory (AMIL). The shapes and masses of the ice deposits were measured, as well as the lift and drag forces of the iced profiles. Scaling was carried out based on the 1.8 MW–Vestas V80 wind turbine technical data, for three different radial positions and two in‐fog icing conditions measured at the Murdochville wind farm in the Gaspé Peninsula. For both icing events, the mass of ice accumulated on the blade profile increased with an increase in the radial position. In wet regime testing (first icing event), glaze formed mostly near the leading edge and on the pressure side. It also accumulated by run‐off on the trailing edge of the outer half of the blade. In dry‐regime testing (second icing event), rime mostly accreted on the leading edge and formed horns. For both icing events, when glaze or rime accreted on the blade profile, lift decreased and drag increased. A load calculation using the blade element theory shows that drag force on the entire blade becomes too large compared to lift, leading to a negative torque and the stop of the wind turbine. Torque reduction is more significant on the outer third of the blade. Setting up a de‐icing system only on the outer part of the blade would enable significant decrease of heating energy costs. Copyright © 2007 John Wiley & Sons, Ltd.     

Effect of icing roughness on wind turbine power production

The objective of this work is a quantitative analysis of power loss of a representative 1.5‐MW wind turbine subject to various icing conditions. Aerodynamic performance data are measured using a combination of ice accretion experiments and wind tunnel tests. Atmospheric icing conditions varying in static temperature, droplet diameter and liquid water content are generated in an icing facility to simulate a 45‐min icing event on a DU 93‐W‐210 airfoil at flow conditions pertinent to 80% blade span on a 1.5‐MW wind turbine. Iced airfoil shapes are molded for preservation and casted for subsequent wind tunnel testing. In general, ice shapes are similar in 2D profile, but vary in 3D surface roughness elements and in the ice impingement length. Both roughness heights and roughness impingement zones are measured. A 16% loss of airfoil lift at operational angle of attack is observed for freezing fog conditions. Airfoil drag increases by 190% at temperatures near 0° C, 145% near 10° C and 80% near 20° C. For a freezing drizzle icing condition, lift loss and drag rise are more severe at 25% and 220%, respectively. An analysis of the wind turbine aerodynamic loads in Region II leads to power losses ranging from 16% to 22% for freezing fog conditions and 26% for a freezing drizzle condition. Differences in power loss between icing conditions are correlated to variance in temperature, ice surface roughness and ice impingement length. Some potential control strategies are discussed for wind turbine operators attempting to minimize revenue loss in cold‐climate regions. Copyright © 2016 John Wiley & Sons, Ltd.     

导流式翼缝对直线翼垂直轴风力机气动性能的影响

为减小流体从吸力面流至压力面的速度损失,基于小间距翼缝有助于减小气动损失的设计原理,针对NACA0021翼型,提出双侧导流式、内导流式和外导流式3种新型翼缝形式。通过数值模拟方法,分析不同翼缝对垂直轴风力机气动性能和流场结构的影响,并将其性能参数与原始翼型和非圆弧翼缝翼型进行对比。结果表明:内导流式翼缝风力机气动性能优于原始风力机,最佳尖速比减小8.06%,改善了叶片周围和整机流场结构,增强了风力机运行稳定性;在低尖速比下,双侧导流式翼缝风力机气动性能较高,而高尖速比时气动性能低于原始风力机;下游区叶片迎风速度较低,外导流式翼缝对流动分离现象改善效果不明显,导致气动性能较差;非圆弧翼缝的间距过大使最大风能利用系数降低了15.5%,不适用于直线翼垂直轴风力机。     

Shape optimization of wind turbine blades

This paper presents a design tool for optimizing wind turbine blades. The design model is based on an aerodynamic/aero‐elastic code that includes the structural dynamics of the blades and the Blade Element Momentum (BEM) theory. To model the main aero‐elastic behaviour of a real wind turbine, the code employs 11 basic degrees of freedom corresponding to 11 elastic structural equations. In the BEM theory, a refined tip loss correction model is used. The objective of the optimization model is to minimize the cost of energy which is calculated from the annual energy production and the cost of the rotor. The design variables used in the current study are the blade shape parameters, including chord, twist and relative thickness. To validate the implementation of the aerodynamic/aero‐elastic model, the computed aerodynamic results are compared to experimental data for the experimental rotor used in the European Commision‐sponsored project Model Experiments in Controlled Conditions, (MEXICO) and the computed aero‐elastic results are examined against the FLEX code for flow past the Tjæreborg 2 MW rotor. To illustrate the optimization technique, three wind turbine rotors of different sizes (the MEXICO 25 kW experimental rotor, the Tjæreborg 2 MW rotor and the NREL 5 MW virtual rotor) are applied. The results show that the optimization model can reduce the cost of energy of the original rotors, especially for the investigated 2 MW and 5 MW rotors. Copyright © 2009 John Wiley & Sons, Ltd.