Numerical investigation of unsteady aerodynamics of a Horizontal‐axis wind turbine under yawed flow conditions

In the present study, unsteady flow features and the blade aerodynamic loading of the National Renewable Energy Laboratory phase VI wind turbine rotor, under yawed flow conditions, were numerically investigated by using a three‐dimensional incompressible flow solver based on unstructured overset meshes. The effect of turbulence, including laminar‐turbulent transition, was accounted for by using a correlation‐based transition turbulence model. The calculations were made for an upwind configuration at wind speeds of 7, 10 and 15 m/sec when the turbine rotor was at 30° and 60° yaw angles. The results were compared with measurements in terms of the blade surface pressure and the normal and tangential forces at selected blade radial locations. It was found that under the yawed flow conditions, the blade aerodynamic loading is significantly reduced. Also, because of the wind velocity component aligned tangent to the rotor disk plane, the periodic fluctuation of blade loading is obtained with lower magnitudes at the advancing blade side and higher magnitudes at the retreating side. This tendency is further magnified as the yaw angle becomes larger. At 7 m/sec wind speed, the sectional angle of attack is relatively small, and the flow remains mostly attached to the blade surface. At 10 m/sec wind speed, leading‐edge flow separation and strong radial flow are observed at the inboard portion of the retreating blade. As the wind speed is further increased, the flow separation and the radial flow become more pronounced. It was demonstrated that these highly unsteady three‐dimensional aerodynamic features are well‐captured by the present method. Copyright © 2012 John Wiley & Sons, Ltd.     

Performance and near wake measurements of a model horizontal axis wind turbine

The performance characteristics and the near wake of a model wind turbine were investigated experimentally. The model tested is a three‐bladed horizontal axis type wind turbine with an upstream rotor of 0.90 m diameter. The performance measurements were conducted at various yaw angles, a freestream speed of about 10 m s ^?1, and the tip speed ratio was varied from 0.5 to 12. The time‐averaged streamwise velocity field in the near wake of the turbine was measured at different tip speed ratios and downstream locations. As expected, it was found that power and thrust coefficients decrease with increasing yaw angle. The power loss is about 3% when the yaw angle is less than 10° and increases to more than 30% when the yaw angle is greater than 30°. The velocity distribution in the near wake was found to be strongly influenced by the tip speed ratio and the yaw angle. At the optimum tip speed ratio, the axial velocity was almost uniform within the midsection of the rotor wake, whereas two strong peaks are observed for high tip speed ratios when the yaw angle is 0°. As the yaw angle increases, the wake width was found to be reduced and skewed towards the yawed direction. With increasing downstream distance, the wake velocity field was observed to depend on the tip speed ratio and more pronounced at high tip speed ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

Wind turbine aerodynamics and loads control in wind shear flow

Wind turbine is subjected to some asymmetrical effects like wind shear, which will lead to unsteady blade airloads and performance. Fatigue loads can lead to damage of turbine components and eventually to failures. It is evident that the variation of the velocity over the rotor disc has an influence on the blade and introduces both flap-wise and edge-wise fatigue damage on the blade as a result of moment fluctuations in the two directions. The flap-wise moments on the blade are the origin of the rotor yaw and tilt moments which transmit to the turbine structure through the drive train to the yaw system and the tower. A lifting surface method with time marching free wake model is used to investigate the periodic unsteady nature in the wind shear. Individual pitch control (IPC) that is applied nowadays is the most advanced active control to reduce the fatigue. The blade airloads and performance of the turbine are also predicted under IPC control. It is found that IPC of the fluctuating blade root flap-wise moment can reduce the flap-wise fatigue damage remarkably while the blade root edge-wise moments are less sensitive to the varying blade pitch than the blade root flap-wise moments.     

Potential of power gain with improved yaw alignment

One of the primary criteria for extracting energy from the wind using horizontal axis upwind wind turbines is the ability to align the rotor axis with the dominating wind direction. The conventional way of estimating the direction of the incoming flow is by using transducers placed atop the nacelle and downwind of the rotor. Recent studies have suggested methods based on advanced upwind measurement technologies for estimating the inflow direction and improving the yaw alignment. In this study, the potential of increased power output with improved yaw alignment is investigated by assessing the performance of a current measurement and yaw control system. The performance is assessed by analyzing data containing upwind wind speed and direction measurements from a met mast, and yaw angle and power production measurements from an operating offshore wind turbine. The results of the analysis indicate that the turbine is operating with a wind speed‐dependent yaw error distribution. The theoretical annual energy production loss due to the yaw error distribution of the existing system is estimated to approximately 0.2%. Copyright © 2014 John Wiley & Sons, Ltd.     

Calibration of a spinner anemometer for wind speed measurements

The power curve of a wind turbine can be measured, according to IEC 61400‐12‐2 with a nacelle‐mounted anemometer. Typically, a sonic anemometer or a cup anemometer and a wind vane are mounted on the back of the nacelle roof. Another option is to use a spinner anemometer. The measurement principle of the spinner anemometer is based on the flow distortion caused by the wind turbine spinner. The flow on the spinner surface is measured by means of three 1D sonic sensors mounted on the spinner and a conversion algorithm to convert the wind velocity components measured by the three sonic sensors to horizontal wind speed, yaw misalignment and flow inclination angle. The algorithm utilizes two calibration constants that are specific to the spinner shape, blade root design and to the mounting positions of the sonic sensors on the spinner. The present analysis describes methods to determine the calibration constant related to wind speed measurements. The first and preferred method is based on the definition of the calibration constant and uses wind speed measurements during the stopped condition of the wind turbine. Two alternative methods that did not require the turbine to be stopped were investigated: one used relatively high wind speed measurements during normal operation of the wind turbine, while the other one used a CFD simulation of the flow over the spinner. The method that entails stopping the turbine in good wind conditions showed the best results and is recommended. The evaluation of uncertainty was not included in the present analysis. Copyright © 2016 The Authors Wind Energy Published by John Wiley & Sons Ltd.     

Numerical study of turbulent flow around a wind turbine nacelle

This article presents a numerical method for the simulation of turbulent flow around the nacelle of a horizontal axis wind turbine. The flow field around the turbine and nacelle is described by the Reynolds‐averaged Navier–Stokes equations. The k–? model has been chosen for closure of the time‐averaged turbulent flow equations. The rotor disc is modelled using the actuator disc concept. The main purpose of this article is to assess the impacts of the variation of some operational parameters (e.g. blade pitch angle changes) and atmospheric turbulence upon the relationship between wind speed measured near the nacelle and freestream wind speed established for an isolated turbine. Simulation results were compared with experimental data (from a typical stall‐controlled, commercially available wind turbine rated higher than 600 kW). In general, good qualitative agreements have been found that validate the proposed method. It has been shown that a level of accuracy sufficient for use in power performance testing can be obtained only when a proper aerodynamic analysis of the inboard non‐lifting cylindrical sections of the blade is included. Furthermore, the numerical method has proven to be a useful tool for locating nacelle anemometers. Copyright © 2005 John Wiley & Sons, Ltd.     

The impact of yaw error on aeroelastic characteristics of a horizontal axis wind turbine blade

Horizontal axis wind turbines operate under yawed conditions for a considerable period of time due to the power control mechanism or sudden changes in the wind direction. This in turn can alter the dynamic characteristics of a turbine blade because the flow over the rotor plane may trigger complicated induced velocity patterns. In this study, an aeroelastic analysis under yawed flow conditions is carried out to investigate the effects of yaw error on the blade behaviors and dynamic stability. A beam model including geometric nonlinearity coupled with unsteady aerodynamics based on a free-vortex wake method with the blade element theory is employed in the present study. The aerodynamic approach for a horizontal axis wind turbine blade under yawed flow conditions is verified through comparison with measurements. It is also shown that the present method gives slightly better results at high yaw angles than does the method previously published in the literature. The dynamic instabilities of a National Renewable Energy Laboratory 5 MW reference wind turbine have subsequently been investigated for various wind speeds and yaw angles. Observations are made that yaw effects induce considerable changes in airloads and blade structural behavior. Also, the aeroelastic damping values for this particular blade under yawed flow conditions can be reduced by up to approximately 33% in the worst case. Therefore, it is concluded that the impacts of yaw misalignments adversely influenced the dynamic aeroelastic stability of the horizontal axis wind turbine blade.     

Effect of rotor–tower interaction,tilt angle,and yaw misalignment on the aeroelasticity of a large horizontal axis wind turbine with composite blades

This paper presents a detailed analysis of the rotor–tower interaction and the effects of the rotor's tilt angle and yaw misalignment on a large horizontal axis wind turbine. A high‐fidelity aeroelastic model is employed, coupling computational fluid dynamics (CFD) and structural mechanics (CSM). The wind velocity stratification induced by the atmospheric boundary layer (ABL) is modeled. On the CSM side, the complex composite structure of each blade is accurately modeled using shell elements. The rotor–tower interaction is analyzed by comparing results of a rotor‐only simulation and a full‐machine simulation, observing a sudden drop in loads, deformations, and power production of each blade, when passing in front of the tower. Subsequently, a tilt angle is introduced on the rotor, and its effect on blade displacements, loads, and performance is studied, representing a novelty with respect to the available literature. The tilt angle leads to a different contribution of gravity to the blade deformations, sensibly affecting the stresses in the composite material. Lastly, a yaw misalignment is introduced with respect to the incoming wind, and the resulting changes in the blade solicitations are analyzed. In particular, a reduction of the blade axial displacement amplitude during each revolution is observed.     

Calibration of a spinner anemometer for yaw misalignment measurements

The spinner anemometer is an instrument for yaw misalignment measurements without the drawbacks of instruments mounted on the nacelle top. The spinner anemometer uses a non‐linear conversion algorithm that converts the measured wind speeds by three sonic sensors on the spinner to horizontal wind speed, yaw misalignment and flow inclination angle. The conversion algorithm utilizes two constants that are specific to the spinner and blade root design and to the mounting positions of the sonic sensors on the spinner. One constant, k₂, mainly affects the measurement of flow angles, while the other constant, k₁, mainly affects the measurement of wind speed. The ratio between the two constants, k_α=k₂/k₁, however, only affects the measurement of flow angles. The calibration of k_α is thus a basic calibration of the spinner anemometer. Theoretical background for the non‐linear calibration is derived from the generic spinner anemometer conversion algorithm. Five different methods were evaluated for calibration of a spinner anemometer on a 500 kW wind turbine. The first three methods used rotor yaw direction as reference angular, while the wind turbine, was yawed in and out of the wind. The fourth method used a hub height met‐mast wind vane as reference. The fifth method used computational fluid dynamics simulations. Method 1 utilizing yawing of the wind turbine in and out of the wind in stopped condition was the preferred method for calibration of k_α. The uncertainty of the yaw misalignment calibration was found to be 10%, giving an uncertainty of 1° at a yaw misalignment of 10°. © 2014 The Authors. Wind Energy published by John Wiley & Sons, Ltd.     

Power output performance characteristics of a horizontal axis marine current turbine

A 0.4 m diameter horizontal axis marine current turbine was tested in a circulating water channel. The power output was measured over a range of flow speeds, blade pitch and rotor yaw angles. Experimental results were compared with the modelled output determined from a commercial blade element momentum computer package. The measured power output was found to be far in excess of predicted values for high blade inflow angles. This occurred where approximately half the blade was operating above the stall angle of attack. This represents 25% of the rotor disk area producing power under heavy stall. Values of overpower up to 140% were measured which are comparable to previous studies. The results show that power production and the optimum tip speed ratio reduced with yaw except for cases with high blade inflow angles.     

Uncertainty of power curve measurement with a two‐beam nacelle‐mounted lidar

Nacelle lidars are attractive for offshore measurements since they can provide measurements of the free wind speed in front of the turbine rotor without erecting a met mast, which significantly reduces the cost of the measurements. Nacelle‐mounted pulsed lidars with two lines of sight (LOS) have already been demonstrated to be suitable for use in power performance measurements. To be considered as a professional tool, however, power curve measurements performed using these instruments require traceable calibrated measurements and the quantification of the wind speed measurement uncertainty. Here we present and demonstrate a procedure fulfilling these needs. A nacelle lidar went through a comprehensive calibration procedure. This calibration took place in two stages. First with the lidar on the ground, the tilt and roll readings of the inclinometers in the nacelle lidar were calibrated. Then the lidar was installed on a 9m high platform in order to calibrate the wind speed measurement. The lidar's radial wind speed measurement along each LOS was compared with the wind speed measured by a calibrated cup anemometer, projected along the LOS direction. The various sources of uncertainty in the lidar wind speed measurement have been thoroughly determined: uncertainty of the reference anemometer, the horizontal and vertical positioning of the beam, the lack of homogeneity of the flow within the probe volume, lidar measurement mean deviation and standard uncertainty. The resulting uncertainty lies between 1 and 2% for the wind speed range between cut‐in and rated wind speed. Finally, the lidar was mounted on the nacelle of a wind turbine in order to perform a power curve measurement. The wind speed was simultaneously measured with a mast‐top mounted cup anemometer placed two rotor diameters upwind of the turbine. The wind speed uncertainty related to the lidar tilting was calculated based on the tilt angle uncertainty derived from the inclinometer calibration and the deviation of the measurement height from hub height. The resulting combined uncertainty in the power curve using the nacelle lidar was less than 10% larger on average than that obtained with the mast mounted cup anemometer. Copyright © 2015 John Wiley & Sons, Ltd.     

Calculation of hoisting forces of the wind turbine rotor based on wind conditions

The assembly and hoisting process of the wind turbine rotor in an open wind environment are regarded to improve the hoisting safety, efficiency and quality. The wind turbine rotor model of a 1.5 MW wind turbine are given, and the hoisting forces of the wind turbine rotor in different poses with various azimuth angles, yaw angles and pitch angles in 3D coordinate system are calculated based on the defined wind conditions model. The maximum and minimum hoisting forces of the wind turbine rotor are acquired and the corresponding azimuth angle, yaw angle and pitch angle of the wind turbine rotor are obtained with respect to the wind conditions in the hoisting process. For four specific poses with particular azimuth angles, yaw angles and pitch angles of the wind turbine rotor, the hoisting forces of the wind turbine rotor are calculated along its hoisting height increment. The change processes of the hoisting forces of the wind turbine rotor in the hoisting process are analyzed and the conclusions are drawn.     

基于流固耦合偏航下风力机尾迹湍流特征的研究

偏航状态下风力机叶片与流场之间相互作用会导致风力机近尾迹流场的湍流特征变化,采用双向流固耦合对不同偏航工况下水平轴风力机近尾迹流场进行数值模拟研究,获得不同偏航角下尾迹湍流特征演化规律。结果表明：随着偏航角的增大,正偏航侧会出现“速度亏损圆环”,且此圆环的范围呈扩大趋势;偏航角的增大对叶根处速度亏损影响最大,对叶尖处速度亏损影响最小,与正偏航侧相比,负偏航侧的速度亏损值减为约1/2;随着偏航角的增大,正负偏航侧的湍流强度变化呈不对称性,正偏航侧对湍流耗散的影响程度较负偏航侧大;涡流黏度越来越小,且在偏航10°涡流黏度相对于偏航5°减小约1/2,沿着轴向叶尖涡的管状环涡结构变得不稳定,出现明显耗散,且在偏航15°之后涡结构的耗散破裂程度越来越剧烈,进而对风力机气动噪声产生较大影响。     

Prediction of aerodynamic performance of NREL offshore 5-MW baseline wind turbine considering power loss at varying wind speeds

In this study, a computational fluid dynamics (CFD) model was developed to simulate the aerodynamic performance of the National Renewable Energy Laboratory (NREL) offshore 5-MW baseline wind turbine with single rotor and full wind turbine. Using statistical methods, the relation between pitch angle and performance when the speed is above the rated wind speed was analyzed; furthermore, other published data were compiled to establish the functional equations of power, thrust with various inflow wind speeds, and pitch angles. In addition, according to shape optimization based on parametric modeling, the two-dimensional cross-section of the wind turbine blade can be defined through a metasurface approach, and the three-dimensional surface modeling of the wind turbine blade, nacelle, and tower is completed using the nonuniform rational B-splines (NURBS) interpolator. In terms of aerodynamic simulation, the unstructured polygon mesh was used herein to discretize the space and simulate the highly curved and twisted surfaces of the blade. In this study, the compact and accurate mesh form obtained through a grid independence test will be used to analyze the distribution of the pressure coefficient, shear stress coefficient, and limited streamline on the blade surface at various inflow wind speeds and pitch angles to understand the differences between different turbulence models and the causes of power and thrust attenuation at high inflow wind speeds. In addition, the phenomena of blade-tip vortices, dynamic stall, blade loading, and the interaction between nacelle and tower were collectively explored.     

The influence of the wind speed profile on wind turbine performance measurements

To identify the influence of wind shear and turbulence on wind turbine performance, flat terrain wind profiles are analysed up to a height of 160 m. The profiles' shapes are found to extend from no shear to high wind shear, and on many occasions, local maxima within the profiles are also observed. Assuming a certain turbine hub height, the profiles with hub‐height wind speeds between 6 m s^?1 and 8 m s^?1 are normalized at 7 m s^?1 and grouped to a number of mean shear profiles. The energy in the profiles varies considerably for the same hub‐height wind speed. These profiles are then used as input to a Blade Element Momentum model that simulates the Siemens 3.6 MW wind turbine. The analysis is carried out as time series simulations where the electrical power is the primary characterization parameter. The results of the simulations indicate that wind speed measurements at different heights over the swept rotor area would allow the determination of the electrical power as a function of an ‘equivalent wind speed’ where wind shear and turbulence intensity are taken into account. Electrical power is found to correlate significantly better to the equivalent wind speed than to the single point hub‐height wind speed. Copyright © 2008 John Wiley & Sons, Ltd.     

基于叶片变形偏航下风力机叶片绕流流场的研究

针对偏航工况下风力机叶片与流场之间的相互作用而产生的变形影响叶片绕流流场问题,基于叶片变形对不同偏航工况下水平轴风力机叶片绕流流场进行双向流固耦合数值计算,分析偏航工况对风力机叶片变形和表面应力的影响,在此基础上研究不同偏航工况对叶片绕流流场的影响.结果表明,不同叶片上的变形和应力呈现不均匀性,且随偏航角增大,不均匀性...     

Computational fluid dynamics analysis of the wake behind the MEXICO rotor in axial flow conditions

This paper presents a computational investigation of the wake of the MEXICO rotor. The compressible multi‐block solver of Liverpool University was employed, using a low‐Mach scheme to account for the low‐speed flow near the blade and in the wake. In this study, computations at wind speeds of 10, 15 and 24 m s ^{? 1} were performed, and the three components of the velocity were compared against experimental data around the rotor blade up to one and a half rotor diameters downstream. Overall, fair agreement was obtained with the computational fluid dynamics showing good vortex conservation near the blade. Vorticity values revealed discontinuities in the wake at approximately 70%R, where two different aerofoils with different zero‐lift angles are blended. The results suggest that all‐Mach schemes for compressible computational fluid dynamics methods can deliver good performance and accuracy over all wind speeds for flows around wind turbines, without the need to switch between incompressible and compressible flow methods. Copyright © 2014 John Wiley & Sons, Ltd.     

An experimental and numerical study of the atmospheric stability impact on wind turbine wakes

In this paper, the impact of atmospheric stability on a wind turbine wake is studied experimentally and numerically. The experimental approach is based on full‐scale (nacelle based) pulsed lidar measurements of the wake flow field of a stall‐regulated 500 kW turbine at the DTU Wind Energy, Risø campus test site. Wake measurements are averaged within a mean wind speed bin of 1 m s^?1 and classified according to atmospheric stability using three different metrics: the Obukhov length, the Bulk–Richardson number and the Froude number. Three test cases are subsequently defined covering various atmospheric conditions. Simulations are carried out using large eddy simulation and actuator disk rotor modeling. The turbulence properties of the incoming wind are adapted to the thermal stratification using a newly developed spectral tensor model that includes buoyancy effects. Discrepancies are discussed, as basis for future model development and improvement. Finally, the impact of atmospheric stability on large‐scale and small‐scale wake flow characteristics is presently investigated. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimating the angle of attack from blade pressure measurements on the National Renewable Energy Laboratory phase VI rotor using a free wake vortex model: yawed conditions

Wind turbine design codes for calculating blade loads are usually based on a blade element momentum (BEM) approach. Since wind turbine rotors often operate in off‐design conditions, such as yawed flow, several engineering methods have been developed to take into account such conditions. An essential feature of a BEM code is the coupling of local blade element loads with an external (induced) velocity field determined with momentum theory through the angle of attack. Local blade loads follow directly from blade pressure measurements as performed in the National Renewable Energy Laboratory (NREL) phase IV campaign, but corresponding angles of attack cannot (on principle) be measured. By developing a free wake vortex method using measured local blade loads, time‐dependent angle of attack and induced velocity distributions are reconstructed. In a previous paper, a method was described for deriving such distributions in conjunction with blade pressure measurements for the NREL phase VI wind turbine in axial (non‐yawed) conditions. In this paper, the same method is applied to investigate yawed conditions on the same turbine. The study considered different operating conditions in yaw in both attached and separated flows over the blades. The derived free wake geometry solutions are used to determine induced velocity distributions at the rotor blade. These are then used to determine the local (azimuth time dependent) angle of attack, as well as the corresponding lift and drag for each blade section. The derived results are helpful to develop better engineering models for wind turbine design codes. Copyright © 2008 John Wiley & Sons, Ltd.     

Validation of the Beddoes–Leishman dynamic stall model for horizontal axis wind turbines using MEXICO data

The aim of this study is to assess the load predicting capability of a classical Beddoes–Leishman dynamic stall model in a horizontal axis wind turbine environment, in the presence of yaw misalignment. The dynamic stall model was tailored to the horizontal axis wind turbine environment and validated against unsteady thick airfoil data. Subsequently, the dynamic stall model was implemented in a blade element‐momentum code for yawed flow, and the results were compared with aerodynamic measurements obtained in the MEXICO (Model Rotor Experiments under Controlled Conditions) project on a wind turbine rotor placed in a large scale wind tunnel. In general, reasonable to good agreement was found between the blade element‐momentum model and MEXICO data. When large yaw misalignments were imposed, poor agreement was found in the downstroke of the movement between the model and the experiment. Still, over a revolution, the maximum normal force coefficient predicted was always within 8% of experimental data at the inboard stations, which is encouraging especially when blade fatigue calculations are being considered. Copyright © 2012 John Wiley & Sons, Ltd.