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.     

Assessment of a comprehensive aeroelastic tool for horizontal‐axis wind turbine rotor analysis

This paper presents the development of a computational aeroelastic tool for the analysis of performance, response and stability of horizontal‐axis wind turbines. A nonlinear beam model for blades structural dynamics is coupled with a state‐space model for unsteady sectional aerodynamic loads, including dynamic stall effects. Several computational fluid dynamics structural dynamics coupling approaches are investigated to take into account rotor wake inflow influence on downwash, all based on a Boundary Element Method for the solution of incompressible, potential, attached flows. Sectional steady aerodynamic coefficients are extended to high angles of attack in order to characterize wind turbine operations in deep stall regimes. The Galerkin method is applied to the resulting aeroelastic differential system. In this context, a novel approach for the spatial integration of additional aerodynamic states, related to wake vorticity and dynamic stall, is introduced and assessed. Steady‐periodic blade responses are evaluated by a harmonic balance approach, whilst a standard eigenproblem is solved for aeroelastic stability analyses. Drawbacks and potentialities of the proposed model are investigated through numerical and experimental comparisons, with particular attention to rotor blades unsteady aerodynamic modelling issues. Copyright © 2016 John Wiley & Sons, Ltd.     

Predicting wind turbine blade loads and aeroelastic response using a coupled CFD–CSD method

The aeroelastic response and the airloads of horizontal-axis wind turbine rotor blades were numerically investigated using a coupled CFD–CSD method. The blade aerodynamic loads were obtained from a Navier–Stokes CFD flow solver based on unstructured meshes. The blade elastic deformation was calculated using a FEM-based CSD solver which employs a nonlinear coupled flap-lag-torsion beam theory. The coupling of the CFD and CSD solvers was accomplished in a loosely coupled manner by exchanging the information between the two solvers at infrequent intervals. At first, the present coupled CFD–CSD method was applied to the NREL 5MW reference wind turbine rotor under steady axial flow conditions, and the mean rotor loads and the static blade deformation were compared with other predicted results. Then, the unsteady blade aerodynamic loads and the dynamic blade response due to rotor shaft tilt and tower interference were investigated, along with the influence of the gravitational force. It was found that due to the aeroelastic blade deformation, the blade aerodynamic loads are significantly reduced, and the unsteady dynamic load behaviors are also changed, particularly by the torsional deformation. From the observation of the tower interference, it was also found that the aerodynamic loads are abruptly reduced as the blades pass by the tower, resulting in oscillatory blade deformation and vibratory loads, particularly in the flapwise direction.     

基于升力线理论的大型风力机气弹响应研究

针对风力机不断向大型化发展的趋势,导致结构柔度增加,气弹耦合特性和振动增强,研究了大型风力机高效精确的气弹响应分析方法。为了更准确模拟大型风力机气流沿叶片展向的三维流动现象,采用螺旋尾涡升力线模型代替传统叶素动量理论,建立了叶片气动载荷分析模型,进而结合风力机多体系统动力学模型,构建了机组的气弹耦合动力学方程和数值求解方法。以某10 MW风力机叶片为例,研究了稳态风况下不同风速的叶片气动性能,以及有效攻角、切向力等沿叶展方向的分布特点,并与采用修正叶素动量理论的气弹分析程序（HAWC）对比,结果表明,升力线理论无需引入经验修正模型即能获得叶素动量理论经修正后的分析精度。最后,通过非稳态风况下风力机的气弹响应分析,证明本文方法对大型风力机气弹耦合分析的有效性和准确性。     

Open‐loop frequency response analysis of a wind turbine using a high‐order linear aeroelastic model

Wind turbine controllers are commonly designed on the basis of low‐order linear models to capture the aeroelastic wind turbine response due to control actions and disturbances. This paper characterizes the aeroelastic wind turbine dynamics that influence the open‐loop frequency response from generator torque and collective pitch control actions of a modern non‐floating wind turbine based on a high‐order linear model. The model is a linearization of a geometrically non‐linear finite beam element model coupled with an unsteady blade element momentum model of aerodynamic forces including effects of shed vorticity and dynamic stall. The main findings are that the lowest collective flap modes have limited influence on the response from generator torque to generator speed, due to large aerodynamic damping. The transfer function from collective pitch to generator speed is affected by two non‐minimum phase zeros below the frequency of the first drivetrain mode. To correctly predict the non‐minimum phase zeros, it is essential to include lateral tower and blade flap degrees of freedom. Copyright © 2013 John Wiley & Sons, Ltd.     

Aeroelastic analysis of a wind turbine blade using the harmonic balance method

An aeroelastic model for wind turbine blades derived from the unsteady Navier‐Stokes equations and a mode shape–based structural dynamics model are presented. For turbulent flows, the system is closed with the Spalart‐Allmaras turbulence model. The computation times for the aerodynamic solution are significantly reduced using the harmonic balance method compared to a time‐accurate solution. This model is significantly more robust than standard aeroelastic codes that rely on blade element momentum theory to determine the aerodynamic forces. Comparisons with published results for the Caradonna‐Tung rotor in hover and the classical AGARD 445.6 flutter case are provided to validate the aerodynamic model and aeroelastic model, respectively. For wind turbines, flutter of the 1.5 MW WindPACT blade is considered. The results predict that the first flapwise and edgewise modes dominate flutter at the rotor speeds considered.     

Influence of wake dynamics on the performance and aeroelasticity of wind turbines

Wind turbines are currently a rapidly expanding form of renewable energy. However, there are numerous technological challenges that must be overcome before wind energy provides a significant amount of power in the United States. One of the primary challenges in wind turbine design and analysis is accurately accounting for the aerodynamic environment. This study is focused on a comprehensive verification and validation of the NREL FAST code, which is enhanced to include a free vortex wake model. The verification and validation is carried out through a comparison of blade lift distribution, wind turbine power and force and moment coefficients using a combination of CFD and experimental data. The results are also compared against Blade Element Momentum theory, and results from a 2001 double-blind NREL study on the prediction capabilities of wind turbine modeling tools. Results indicate that the enhanced aeroelastic code generally provides improved predictions. However, in several notable cases the predictions are only marginally improved, or even worse, than those generated using Blade Element Momentum theory aerodynamics. It is concluded that modeling of the aerodynamic environment remains incomplete, even after inclusion of wake effects. One important aspect identified is modeling of the unsteady aerodynamic lift characteristics of the rotor. Finally, the aeroelastic response in the combined presence of wake effects and inflow turbulence is examined. Significant differences are observed in loads, power, and structural response between results computed using the free wake model or simpler models, such as Blade Element Momentum theory.     

基于拟合气弹系数的气弹稳定性分析

鉴于多自由度大型风力机叶型气动系数计算以及气弹稳定性分析的复杂性,该文在一种四自由度叶型的基础上,采用拟合气弹系数简化了气动力计算过程及气弹稳定性分析过程,降低了非线性气动力模型的阶次。通过与相关文献以及Theodorsen方法进行气弹稳定性分析结果对比,验证了该方法的可靠性。拟合气动系数的方法可用于NACA的各类叶型。     

风力机叶片动态气弹变形及其对整机性能的影响

基于风力机整机刚柔耦合模型,文章提出了一种叶片动态气弹扭转变形分析的新方法。该方法采用SIMPACK和AeroDyn软件联合数值仿真对风力机在几种恶劣风况下进行动力学分析,通过对分析结果的变换处理,进而得到叶片在复杂工况下的动态气弹变形数据。采用该方法,重点分析了叶片气弹扭转变形对风力机气动功率及气弹稳定性的影响。该方法为大型风电叶片的气弹特性评价以及气弹剪裁设计提供了一种新的技术手段。     

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.     

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.     

Aeroelastic stability predictions for a MW‐sized blade

Classical aeroelastic flutter instability historically has not been a driving issue in wind turbine design. In fact, rarely has this issue even been addressed in the past. Commensurately, among the wind turbines that have been built, rarely has classical flutter ever been observed. However, with the advent of larger turbines fitted with relatively softer blades, classical flutter may become a more important design consideration. In addition, innovative blade designs involving the use of aeroelastic tailoring, wherein the blade twists as it bends under the action of aerodynamic loads to shed load resulting from wind turbulence, may increase the blade's proclivity for flutter. With these considerations in mind it is prudent to revisit aeroelastic stability issues for a MW‐sized blade with and without aeroelastic tailoring. Focusing on aeroelastic stability associated with the shed wake from an individual blade turning in still air, the frequency domain technique developed by Theodorsen for predicting classical flutter in fixed wing aircraft has been adapted for use with a rotor blade. Results indicate that the predicted flutter speed of a MW‐sized blade is slightly greater than twice the operational speed of the rotor. When a moderate amount of aeroelastic tailoring is added to the blade, a modest decrease (12%) in the flutter speed is predicted. By comparison, for a smaller rotor with relatively stiff blades the predicted flutter speed is approximately six times the operating speed. When frequently used approximations to Theodorsen's method are implemented, drastic underpredictions result, which, while conservative, may adversely impact blade design. These underpredictions are also evident when this MW‐sized blade is analysed using time domain methods. Published in 2004 by John Wiley & Sons, Ltd.     

Application of the actuator surface concept to wind turbine rotor aerodynamics

The structure of blade tip vortices is recognized as a key issue in wind turbine aerodynamic modelling by many researchers in the field. In the search for an intermediate model between full Navier–Stokes and blade‐element momentum simulations, this article presents a method using rotating actuator surfaces to model wind turbine aerodynamics. An actuator surface is a simple planar surface, porous to the flow, which is characterized by velocity and pressure discontinuities, whose action on the flow is achieved through an attached system of forces. These discontinuities and forces are determined from blade‐element analysis and the Kutta–Joukowski relation. After implementing this concept in a three‐dimensional CFD (Computational Fluid Dynamics) method, results are produced for the experimental rotors of NREL and TUDelft. The method is validated against both experimental measurements and the predictions of three other numerical models for wind turbine aerodynamic analysis. Qualitative and quantitative comparisons show that the actuator surface concept agrees well with the other numerical models. In addition to rotor aerodynamic analysis, the actuator surface concept can be used in the study of wake aerodynamics, or as the Eulerian flow solver in hybrid methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Application and validation of incrementally complex models for wind turbine aerodynamics,isolated wind turbine in uniform inflow conditions

A comparison of several incrementally complex methods for predicting wind turbine performance, aeroelastic behavior, and wakes is provided. Depending on a wind farm's design, wake interference can cause large power losses and increased turbulence levels within the farm. The goal is to employ modeling methods to reach an improved understanding of wake effects and to use this information to better optimize the layout of new wind farms. A critical decision faced by modelers is the fidelity of the model that is selected to perform simulations. The choice of model fidelity can affect the accuracy, but will also greatly impact the computational time and resource requirements for simulations. To help address this critical question, three modeling methods of varying fidelity have been developed side by side and are compared in this article. The models from low to high complexity are as follows: a blade element‐based method with a free‐vortex wake, an actuator disc‐based method, and a full rotor‐based method. Fluid/structure interfaces are developed for the aerodynamic modeling approaches that allow modeling of discrete blades and are then coupled with a multibody structural dynamics solver in order to perform an aeroelastic analysis. Similar methods have individually been tested by researchers, but we suggest that by developing a suite of models, they can be cross‐compared to grasp the subtleties of each method. The modeling methods are applied to the National Renewable Energy Laboratory Phase VI rotor to predict the turbine aerodynamic and structural loads and then also the wind velocities in the wake. The full rotor method provides the most accurate predictions at the turbine and the use of adaptive mesh refinement to capture the wake to 20 radii downstream is proven particularly successful. Though the full rotor method is unmatched by the lower fidelity methods in stalled conditions and detailed prediction of the downstream wake, there are other less complex conditions where these methods perform as accurately as the full rotor method. Copyright © 2014 John Wiley & Sons, Ltd.     

Flap/lead–lag aeroelastic stability of wind turbine blade sections

The scope of this article is to investigate the aeroelastic stability of wind turbine blade sections subjected to combined flap/lead–lag motion. The work is motivated by recent concern about destructive ‘edgewise' vibrations of modern, half‐megawatt‐scale, blades. The aeroelastic governing equations derive from the combination of a spring–mass–damper equivalent of the structure and a ‘non‐stationary' aerodynamic model. The aerodynamic model used in the present context is the differential dynamic stall model developed at ONERA. The resulting equations of motion are linearized and their stability characteristics are investigated in terms of the system entries, expressed through suitable, non‐dimensional, structural and aerodynamic parameters. Copyright © 1999 John Wiley & Sons, Ltd.     

A CFD study into the influence of unsteady aerodynamic interference on wind turbine surge motion

To improve knowledge of the unsteady aerodynamic characteristics and interference effects of a floating offshore wind turbine (FOWT), this article focuses on the platform surge motion of a full configuration wind turbine with the rotating blades, hub, nacelle, and tower shapes. Unsteady aerodynamic analyses considering the moving motion of an entire configuration wind turbine have been conducted using an advanced computational fluid dynamics (CFD) and a conventional blade element momentum (BEM) analyses. The present CFD simulation is based on an advanced overset moving grid method to accurately consider the local and global motion of a three-dimensional wind turbine. The effects of various oscillation frequencies and amplitudes of the platform surge motion have been widely investigated herein. Three-dimensional unsteady flow fields around the moving wind turbine with rotating blades are graphically presented in detail. Complex flow interactions among blade tip vortices, tower shedding vortices, and turbulent wakes are physically observed. Comparisons of different aerodynamic analyses under the periodic surge motions are summarized to show the potential distinction among applied numerical methods. The present result indicates that the unsteady aerodynamic thrust and power tend to vary considerably depending on the oscillation frequency and amplitude of the surge motion.     

Forced responses on a radial turbine with nozzle guide vanes

Radial turbines with nozzle guide vanes are widely used in various size turbochargers. However, due to the inter- ferences with guide vanes, the blades of impellers are exposed to intense unsteady aerodynamic excitations, which cause blade vibrations and lead to high cycle failures （HCF）.Moreover, the harmonic resonance in some frequency regions are unavoidable due to the wide operation conditions. Aiming to achieve a detail insight into vibration characteristics of radial flow turbine, a numerical method based on fluid .structure interaction （FSI） is presented. Firstly, the unsteady aerodynamic loads are determined by computational fluid dynamics （CFD）.And the fluctuating pressures are transformed from time domain to frequency domain by fast Fourier-transform （FFT）. Then, the entire rotor model is adopted to analyze frequencies and mode shapes considering mistuning in finite element （FE） method. Meanwhile, harmonic analyses, applying the pressure fluctuation from CFD, are conducted to investigate the impeller vibration behavior and blade forced response in frequency domain. The prediction of the vibration dynamic stress shows acceptable agreement to the blade actual damage in consistent tendency.     

An actuator line method with novel blade flow field coupling based on potential flow equivalence

A Reynolds‐averaged Navier–Stokes‐embedded actuator line model for wind and tidal turbine simulation has been implemented and validated using the National Renewable Energy Laboratory Phase VI wind tunnel experimental results. Actuator line models, first introduced by Sørensen and Shen, represent the blades virtually, enabling time‐resolved rotor simulations without requiring blade boundary layer discretization. This results in a lower computational cost than blade‐resolved simulations while preserving the predominant features of the rotor flow. The present method introduces a novel technique, based on potential flow equivalence, to determine the local flow velocity at the blade, and a method of projecting the resulting momentum sources to the flow field. These methods circumvent the requirement for smearing techniques used in other actuator line models. In addition, the model is adapted for use with an unstructured mesh, thereby enabling turbine components such as the tower and nacelle to be explicitly included in the domain. The model is validated through comparison of computed integrated loads and local force coefficients with the National Renewable Energy Laboratory Phase VI experimental results. Results for local force coefficients indicate general agreement with experiment, although discrepancies associated with three‐dimensional flow effects are observed at the tips. Copyright © 2014 John Wiley & Sons, Ltd.     

A stochastic model for the simulation of wind turbine blades in static stall

The aim of this work is to improve aeroelastic simulation codes by accounting for the unsteady aerodynamic forces that a blade experiences in static stall. A model based on a spectral representation of the aerodynamic lift force is defined. The drag and pitching moment are derived using a conditional simulation technique for stochastic processes. The input data for the model can be collected either from measurements or from numerical results from a Computational Fluid Dynamics code for airfoil sections at constant angles of attack. An analysis of such data is provided, which helps to determine the characteristics of stall. The model is applied to wind turbine rotor cases, including the stand still condition, and results are compared to experimental data. Copyright © 2009 John Wiley & Sons, Ltd.