Blade load mitigation control design for a wind turbine operating in the path of vortices

As wind turbine rotor size continues to increase, load mitigation becomes an important control objective. Turbines with hub heights of nearly 100m operate in the stable, nocturnal boundary layer where coherent turbulence can be generated by atmospheric phenomena outside the surface layer. These coherent turbulent structures may contribute to blade fatigue loads that can be mitigated with advanced control algorithms. Disturbance accommodating control (DAC) methods were implemented in a wind turbine structural dynamics simulation code to mitigate transient blade load response induced by a simple, Rankine vortex in the inflow. As a best‐case scenario, a full‐state feedback controller (which included a very detailed disturbance model) showed that blade flap damage equivalent load caused by the vortex passing through the rotor could be reduced by 30% compared to one that resulted from simulation of a typical proportional‐integral (PI) controller. A realizable DAC controller that incorporates only the vertical shear component of the vortex reduced loads by 9% compared to that resulting from simulation of a PI controller. The load reduction was even greater when the vortex was superimposed over full‐field, homogeneous turbulence. DAC methods have the flexibility to incorporate properties of coherent turbulent inflow structures in the controller design to mitigate blade fatigue loads. Further work must be done to develop disturbance models as more details about the turbulent structures are identified. Copyright © 2007 John Wiley & Sons, Ltd.     

Large‐eddy simulation of stable boundary layer turbulence and estimation of associated wind turbine loads

Stochastic simulation of turbulent inflow fields commonly used in wind turbine load computations is unable to account for contrasting states of atmospheric stability. Flow fields in the stable boundary layer, for instance, have characteristics such as enhanced wind speed and directional shear; these effects can influence loads on utility‐scale wind turbines. To investigate these influences, we use large‐eddy simulation (LES) to generate an extensive database of high‐resolution ( ～ 10 m), four‐dimensional turbulent flow fields. Key atmospheric conditions (e.g., geostrophic wind) and surface conditions (e.g., aerodynamic roughness length) are systematically varied to generate a diverse range of physically realizable atmospheric stabilities. We show that turbine‐scale variables (e.g., hub height wind speed, standard deviation of the longitudinal wind speed, wind speed shear, wind directional shear and Richardson number) are strongly interrelated. Thus, we strongly advocate that these variables should not be prescribed as independent degrees of freedom in any synthetic turbulent inflow generator but rather that any turbulence generation procedure should be able to bring about realistic sets of such physically realizable sets of turbine‐scale flow variables. We demonstrate the utility of our LES‐generated database in estimation of loads on a 5‐MW wind turbine model. More importantly, we identify specific turbine‐scale flow variables that are responsible for large turbine loads—e.g., wind speed shear is found to have a greater influence on out‐of‐plane blade bending moments for the turbine studied compared with its influence on other loads such as the tower‐top yaw moment and the fore‐aft tower base moment. Overall, our study suggests that LES may be effectively used to model inflow fields, to study characteristics of flow fields under various atmospheric stability conditions and to assess turbine loads for conditions that are not typically examined in design standards. Copyright © 2013 John Wiley & Sons, Ltd.     

Using Wavelet Analysis to Assess Turbulence/Rotor Interactions

Large loading events on wind turbine rotor blades are often associated with transient bursts of coherent turbulent energy in the turbine inflow. These coherent turbulent structures are identified as peaks in the three‐dimensional, instantaneous, turbulent shearing stress field. Such organized inflow structures and the accompanying rotor aeroelastic responses typically have timescales of only a few seconds and therefore do not lend themselves to analysis by conventional Fourier spectral techniques. Wavelet analysis offers the ability to study more closely the spectral decomposition of short‐period events such as the interaction of coherent turbulence with a moving rotor blade. In this paper we discuss our initial progress in the application of wavelet analysis techniques to the decomposition and interpretation of turbulence/rotor interaction. We discuss the results of using both continuous and discrete wavelet transforms for our application. Several examples are given of the techniques applied to both observed turbulence and turbine responses and those generated using numerical simulations. We found that the presence of coherent turbulent structures, as revealed by the inflow Reynolds stress field, is a major contributor to large load excursions. These bursts of coherent turbulent energy induce a broadband aeroelastic response in the turbine rotor as it passes through them. Copyright © 2001 John Wiley & Sons, Ltd.     

HAWT dynamic stall response asymmetries under yawed flow conditions

Horizontal axis wind turbines can experience significant time‐varying aerodynamic loads, potentially causing adverse effects on structures, mechanical components and power production. As designers attempt lighter and more flexible wind energy machines, greater accuracy and robustness will become even more critical in future aerodynamics models. Aerodynamics modelling advances, in turn, will rely on more thorough comprehension of the three‐dimensional, unsteady, vortical flows that dominate wind turbine blade aerodynamics under high‐load conditions. To experimentally characterize these flows, turbine blade surface pressures were acquired at multiple span locations via the NREL Phase IV Unsteady Aerodynamics Experiment. Surface pressures and associated normal force histories were used to characterize dynamic stall vortex kinematics and normal force amplification. Dynamic stall vortices and normal force amplification were confirmed to occur in response to angle‐of‐attack excursions above the static stall threshold. Stall vortices occupied approximately one‐half of the blade span and persisted for nearly one‐fourth of the blade rotation cycle. Stall vortex convection varied along the blade, resulting in dramatic deformation of the vortex. Presence and deformation of the dynamic stall vortex produced corresponding amplification of normal forces. Analyses revealed consistent alterations to vortex kinematics in response to changes in reduced frequency, span location and yaw error. Finally, vortex structures and kinematics not previously documented for wind turbine blades were isolated. Published in 2000 by John Wiley & Sons, Ltd.     

Wind turbine extreme gust control

This paper focuses on the problem of extreme wind gust and direction change recognition (EG&DR) and control (EEC). An extreme wind gust with direction change can lead to large loads on the turbine (causing fatigue) and unnecessary turbine shutdowns by the supervisory system caused by rotor overspeed. The proposed EG&DR algorithm is based on a non‐linear observer (extended Kalman filter) that estimates the oblique wind inflow angle and the blade effective wind speed signals, which are then used by a detection algorithm (cumulative sum test) to recognize extreme events. The non‐linear observer requires that blade root bending moments measurements (in‐plane and out‐of‐plane) are available. Once an extreme event is detected, an EEC algorithm is activated that: (i) tries to prevent the rotor speed from exceeding the overspeed limit by fast collective blade pitching; and (ii) reduces 1p blade loads by means of individual pitch control algorithm, designed in an ?_∞ optimal control setting. The method is demonstrated on a complex non‐linear test turbine model. Copyright © 2009 John Wiley & Sons, Ltd.     

New vortex‐lift and tangential‐force models for HAWT aerodynamic load prediction

Horizontal axis wind turbines (HAWTs) experience three‐dimensional rotational and unsteady aerodynamic phenomena at the rotor blades sections. These highly unsteady three‐dimensional effects have a dramatic impact on the aerodynamic load distributions on the blades, in particular, when they occur at high angles of attack due to stall delay and dynamic stall. Unfortunately, there is no complete understanding of the flow physics yet at these unsteady 3D flow conditions, and hence, the existing published theoretical models are often incapable of modelling the impact on the turbine response realistically. The purpose of this paper is to provide an insight on the combined influence of the stall delay and dynamic stall on the blade load history of wind turbines in controlled and uncontrolled conditions. New dynamic stall vortex and nonlinear tangential force coefficient modules, which integrally take into account the three dimensional rotational effect, are also proposed in this paper. This module along with the unsteady influence of turbulent wind speed and tower shadow is implemented in a blade element momentum (BEM) model to estimate the aerodynamic loads on a rotating blade more accurately. This work presents an important step to help modelling the combined influence of the stall delay and dynamic stall on the load history of the rotating wind turbine blades which is vital to have lighter turbine blades and improved wind turbine design systems.     

Long‐term loads for a monopile‐supported offshore wind turbine

Accurate prediction of long‐term ‘characteristic’ loads associated with an ultimate limit state for design of a 5‐MW bottom‐supported offshore wind turbine is the focus of this study. Specifically, we focus on predicting the long‐term fore–aft tower bending moment at the mudline and the out‐of‐plane bending moment at the blade root of a monopile‐supported shallow‐water offshore wind turbine. We employ alternative probabilistic predictions of long‐term loads using inverse reliability procedures in establishing the characteristic loads for design. Because load variability depends on the environmental conditions (defining the wind speed and wave height), we show that long‐term predictions that explicitly account for such load variability are more accurate, especially for environmental states associated with above‐rated wind speeds and associated wave heights. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of light detection and ranging wind speed measurements for wind turbine control

Light detection and ranging (LIDAR) systems are able to measure the speed of incoming wind before it reaches a wind turbine rotor. These preview wind measurements can be used in feedforward control systems designed to reduce turbine structural loads. However, the degree to which such preview‐based control techniques can reduce loads by reacting to turbulence depends on how accurately the incoming wind field can be measured. This study examines the accuracy of different measurement scenarios that rely on coherent continuous‐wave or pulsed Doppler LIDAR systems, in terms of root‐mean‐square measurement error, to determine their applicability to feedforward control. In particular, the impacts of measurement range, angular offset of the LIDAR beam from the wind direction, and measurement noise are studied for various wind conditions. A realistic simulation case involving a scanning LIDAR unit mounted in the spinner of a MW‐scale wind turbine is studied in depth, with emphasis on preview distances that provide minimum measurement error for a specific scan radius. Measurement error is analyzed for LIDAR‐based estimates of point wind speeds at the rotor as well as spanwise averaged blade effective wind speeds. The impact of turbulence structures with high coherent turbulent kinetic energy on measurement error is discussed as well. Copyright © 2013 John Wiley & Sons, Ltd.     

Analytic analysis of load alignment for coning extreme‐scale rotors

Extreme‐scale wind turbines (rated powers greater than 10 MW) with large rotor diameters and conventional upwind designs must resist extreme downwind and gravity loads. This can lead to significant structural design challenges and high blade masses that can impede the reduction of levelized cost of wind energy. Herein, the theoretical basis for downwind load alignment is developed. This alignment can be addressed with active downwind coning to reduce/eliminate flapwise bending loads by balancing the transverse components of thrust, centrifugal, and gravitational force. Equations are developed herein that estimates the optimal coning angle that reduces flapwise loads by a specified amount. This analysis is then applied to a 13.2‐MW scale with 100‐m‐level wind turbine blades, where it is found that a load alignment coning schedule can substantially reduce the root flapwise bending moments. This moment reduction in this example can allow the rotor mass to be decreased significantly when compared with a conventional upwind three‐bladed rotor while maintaining structural performance and annual energy output.     

An experimental study on the effects of winglets on the tip vortex interaction in the near wake of a model wind turbine

An experimental study of the near wake up to four rotor diameters behind a model wind turbine rotor with two different wing tip configurations is performed. A straight‐cut wing tip and a downstream‐facing winglet shape are compared on the same two‐bladed rotor operated at its design tip speed ratio. Phase‐averaged measurements of the velocity vector are synchronized with the rotor position, visualizing the downstream location of tip vortex interaction for the two blade tip configurations. The mean streamwise velocity is found not to be strongly affected by the presence of winglet tip extensions, suggesting an insignificant effect of winglets on the time‐averaged inflow conditions of a possible downstream wind turbine. An analysis of the phase‐averaged vorticity, however, reveals a significantly earlier tip vortex interaction and breakup for the wingletted rotor. In contradistinction, the tip vortices formed behind the reference configuration are assessed to be more stable and start merging into larger turbulent structures significantly further downstream. These results indicate that an optimized winglet design can not only contribute to a higher energy extraction in a rotor's tip region but also can positively affect the wake's mean kinetic energy recovery by stimulating a faster tip vortex interaction.     

On wind turbine loads during the evening transition period

The late afternoon hours in the diurnal cycle precede the development of the nocturnal stable boundary layer. This “evening transition” (ET) period is often when energy demand peaks. This period also corresponds to the time of day that is a precursor to late‐afternoon downbursts, a subject of separate interest. To capture physical characteristics of wind fields in the atmospheric boundary layer (ABL) during this ET period, particularly the interplay of shear and turbulence, stochastic simulation approaches, although more tractable, are not suitable. Large‐eddy simulation (LES), on the other hand, may be used to generate high‐resolution ABL turbulent flow fields. We present a suite of idealized LES four‐dimensional flow fields that define a database representing different combinations of large‐scale atmospheric conditions (characterized by associated geostrophic winds) and surface boundary conditions (characterized by surface heat fluxes). Our objective is to evaluate the performance of wind turbines during the ET period. Accordingly, we conduct a statistical analysis of turbine‐scale wind field variables. We then employ the database of these LES‐based inflow wind fields in aeroelastic simulations of a 5‐MW wind turbine. We discuss how turbine loads change as the ET period evolves. We also discuss maximum and fatigue loads on the rotor and tower resulting from different ABL conditions. Results of this study suggest that, during the ET period, the prevailing geostrophic wind speed affects the mean and variance of longitudinal winds greatly and thus has significant influence on all loads except the yaw moment which is less sensitive to uniform and symmetric incoming flow. On the other hand, surface heat flux levels affect vertical turbulence and wind shear more and, as a result, only affect maximum blade flapwise bending and tower fore‐aft bending loads.     

Simulating the aerodynamic performance and wake dynamics of a vertical‐axis wind turbine

The accurate prediction of the aerodynamics and performance of vertical‐axis wind turbines is essential if their design is to be improved but poses a significant challenge to numerical simulation tools. The cyclic motion of the blades induces large variations in the angle of attack of the blades that can manifest as dynamic stall. In addition, predicting the interaction between the blades and the wake developed by the rotor requires a high‐fidelity representation of the vortical structures within the flow field in which the turbine operates. The aerodynamic performance and wake dynamics of a Darrieus‐type vertical‐axis wind turbine consisting of two straight blades is simulated using Brown's Vorticity Transport Model. The predicted variation with azimuth of the normal and tangential force on the turbine blades compares well with experimental measurements. The interaction between the blades and the vortices that are shed and trailed in previous revolutions of the turbine is shown to have a significant effect on the distribution of aerodynamic loading on the blades. Furthermore, it is suggested that the disagreement between experimental and numerical data that has been presented in previous studies arises because the blade–vortex interactions on the rotor were not modelled with sufficient fidelity. Copyright © 2010 John Wiley & Sons, Ltd.     

Changes in design driving load cases: Operating an upwind turbine with a downwind rotor configuration

This work considers the design driving load cases from a full design load basis analysis on an upwind turbine changed into a downwind configuration. The upwind turbine is a commercial class IIIA 2.1‐MW turbine, manufactured by Suzlon. The downwind turbine shows an increase in the normalized tower clearance by 6%, compared with the upwind concept. Removing the blade prebend increases the normalized minimum tower clearance by 17% in the downwind configuration compared with the upwind configuration. The extreme loads on the longitudinal tower bottom bending moment are seen to generally increase by 17% because of the overhanging gravity moment of the rotor‐nacelle assembly. The extreme blade root bending moments are reduced by 10% flapwise, because of the coning of the rotor in downwind direction. The fatigue loads suffer from the tower shadow, leading to an overall increase of the fatigue loads in the blades with up to 5% in flapwise direction in the downwind configuration. Because of blade deflection and coning direction, the downwind configuration shows a 0.75% lower annual energy production. Removing the prebend increases the annual energy production loss to 1.66%.     

低空急流条件下水平轴风力机风轮气动特性的研究

为阐明低空急流条件下风力机风轮的气动特性,基于工程化的边界层风速模型和Von Karman谱模型建立不同来流的脉动风场,对比研究低空急流条件下NREL 5 MW风力机风轮的输出功率和气动载荷的变化规律。结果表明:如果仅以轮毂高度处的风速作为风力机变桨控制的依据,与均匀来流和剪切来流相比较,低空急流条件下,虽然来流风功率明显增大,但风轮的输出功率在较高风速时反而减小;风轮所受的不平衡气动载荷,包括横向力、纵向力、偏航力矩和倾覆力矩在较高风速时小于剪切来流的结果;且仅以轮毂高度处的风速预测得到的风轮输出功率高于实际结果,其最大相对误差为89.4%。因此,低空急流条件下,为提高风能利用率和风轮输出功率的预测精度,应考虑不同高度位置处的风速大小对风力机进行变桨控制和功率预测。     

Flow characterization in the near‐wake region of a horizontal axis wind turbine

An experimental study is conducted to investigate the flow dynamics within the near‐wake region of a horizontal axis wind turbine using particle image velocimetry (PIV). Measurements were performed in the horizontal plane in a row of four radially distributed measurement windows (tiles), which are then patched together to obtain larger measurement field. The mean and turbulent components of the flow field were measured at various blade phase angles. The mean velocity and turbulence characteristics show high dependency on the blade phase angle in the near‐wake region closer to the blade tip and become phase independent further downstream at a distance of about one rotor diameter. In the near‐wake region, both the mean and turbulent characteristics show a systemic variation with the phase angle in the blade tip region, where the highest levels of turbulence are observed. The streamlines of the instantaneous velocity field at a given phase allowed to track a tip vortex which showed wandering trend. The tip vortices are mostly formed at r/R > 1, which indicates the wake expansion. Results also show the gradual movement of the vortex region in the axial direction, which can be attributed to the dynamics of the helical tip vortices which after being generated from the tip, rotate with respect to the blade and move in the axial direction because of the axial momentum of the flow. The axial velocity deficit was compared with other laboratory and field measurements. The comparison shows qualitative similarity. Copyright © 2015 John Wiley & Sons, Ltd.     

Further load reductions with individual pitch control

Previous work has demonstrated that significant reductions in fatigue loading on a wind turbine can be achieved by using individual pitch control, in which the pitch of each blade is adjusted individually, in response to measured loads. The asymmetrical out‐of‐plane rotor load is measured and an additional pitch action (dominated by the rotational frequency of the rotor) is calculated for each blade in order to minimize this load. This results in the near‐elimination of the dominant once‐per‐revolution (‘1P’) peak in the out‐of‐plane load spectrum seen by the rotating components, and fatigue loads can be reduced by 20%–40%. The load reduction is also transferred to the nacelle and tower, but here it is the low‐frequency loads which are removed, resulting in a load reduction of a few per cent at best, since the fatigue on the fixed components is dominated by the peak at the blade passing frequency (‘3P’ for a three‐bladed turbine), which is largely unaffected by the individual pitch control action. This article demonstrates a relatively straightforward addition to the individual pitch control algorithm which is capable of reducing the dominant load peak on the fixed components, resulting in significant fatigue load reductions on the whole structure. Copyright © 2005 John Wiley & Sons, Ltd.     

Load analysis of hydraulic‐pneumatic flywheel configurations integrated in a wind turbine rotor

In this paper, the impact on the mechanical loads of a wind turbine due to a previously proposed hydraulic‐pneumatic flywheel system is analysed. Load simulations are performed for the National Renewable Energy Laboratory (NREL) 5‐MW wind turbine using fatigue, aerodynamics, structures, and turbulence (FAST). It is discussed why FAST is applied although it cannot simulate variable rotor inertia. Several flywheel configurations, which increase the rotor inertia of the 5‐MW wind turbine by 15%, are implemented in the 61.5‐m rotor blade. Load simulations are performed twice for each configuration: Firstly, the flywheel system is discharged, and secondly, the flywheel is charged. The change in ultimate and fatigue loads on the tower, the low speed shaft, and the rotor blades is juxtaposed for all flywheel configurations. As the blades are mainly affected by the flywheel system, the increase in ultimate and fatigue loads of the blade is evaluated. Simulation results show that the initial design of the flywheel system causes the lowest impact on the mechanical loads of the rotor blades although this configuration is the heaviest.     

Performance enhancement and load reduction of a 5 MW wind turbine blade

A wind turbine rotor blade, based on the U.S. National Renewable Energy Laboratory (NREL) 5 MW reference turbine, is optimized for minimum cost of energy through simultaneous consideration of aerodynamics and bend-twist coupling. Eighty-three total design variables are considered, encompassing airfoil shapes, chord and twist distributions, and the degree of bend-twist coupling in the blade. A recently developed method requiring significantly less computation than finite element analysis is used for planning and predicting the bend-twist coupling behavior of the rotor. Airfoil performance is computed using XFOIL, while the wind turbine loads and performance are computed using the NREL FAST code. The objective function is annual cost of energy (COE), where reductions in flapwise bending loads and blade surface area are assumed to decrease rotor cost through reduced material requirements. The developed optimization process projects decreased blade loads while maintaining wind turbine performance.     

Influence of atmospheric stability on wind turbine loads

Simulations of wind turbine loads for the NREL 5 MW reference wind turbine under diabatic conditions are performed. The diabatic conditions are incorporated in the input wind field in the form of wind profile and turbulence. The simulations are carried out for mean wind speeds between 3 and 16 m s ^{? 1} at the turbine hub height. The loads are quantified as the cumulative sum of the damage equivalent load for different wind speeds that are weighted according to the wind speed and stability distribution. Four sites with a different wind speed and stability distribution are used for comparison. The turbulence and wind profile from only one site is used in the load calculations, which are then weighted according to wind speed and stability distributions at different sites. It is observed that atmospheric stability influences the tower and rotor loads. The difference in the calculated tower loads using diabatic wind conditions and those obtained assuming neutral conditions only is up to 17%, whereas the difference for the rotor loads is up to 13%. The blade loads are hardly influenced by atmospheric stability, where the difference between the calculated loads using diabatic and neutral input wind conditions is up to 3% only. The wind profiles and turbulence under diabatic conditions have contrasting influences on the loads; for example, under stable conditions, loads induced by the wind profile are larger because of increased wind shear, whereas those induced by turbulence are lower because of less turbulent energy. The tower base loads are mainly influenced by diabatic turbulence, whereas the rotor loads are influenced by diabatic wind profiles. The blade loads are influenced by both, diabatic wind profile and turbulence, that leads to nullifying the contrasting influences on the loads. The importance of using a detailed boundary‐layer wind profile model is also demonstrated. The difference in the calculated blade and rotor loads is up to 6% and 8%, respectively, when only the surface‐layer wind profile model is used in comparison with those obtained using a boundary‐layer wind profile model. Finally, a comparison of the calculated loads obtained using site‐specific and International Electrotechnical Commission (IEC) wind conditions is carried out. It is observed that the IEC loads are up to 96% larger than those obtained using site‐specific wind conditions.Copyright © 2012 John Wiley & Sons, Ltd.     

On the proof of concept of a ‘Smart’ wind turbine rotor blade for load alleviation

The trend with offshore wind turbines is to increase the rotor diameter as much as possible to decrease the costs per kilowatt‐hour. The increasing dimensions have led to the relative increase of the loads on the wind turbine structure. Because of the increasing rotor size and the spatial load variations along the blade, it is necessary to react to turbulence in a more detailed way; each blade separately and at several separate radial distances. In this paper, a proof of concept study is performed to show the feasibility of the load alleviation abilities of a ‘Smart’ blade, i.e. a blade equipped with a number of control devices that locally change the lift profile on the blade, combined with appropriate sensors and feedback controllers. Theoretical and experimental models are developed of a scaled non‐rotating rotor blade which is equipped with two trailing edge flaps and strain sensors to facilitate feedback control. A pitch actuator is used to induce disturbances with a similar character as a gust or turbulence. A feedback controller based on classical loop shaping is designed that minimizes the root bending moment in the flapping direction. We show that with appropriate control techniques, the loads for periodic disturbances and for turbulence generated disturbances can be reduced up to 90 and 55%, respectively. Copyright © 2008 John Wiley & Sons, Ltd.