Rapid optimization of stall‐regulated wind turbine blades using a frequency‐domain method: Part 2, cost function selection and results

A fast and effective frequency‐domain optimization method was developed for stall‐regulated blades. It was found that when using linearized dynamics, typical cost functions employing damage‐equivalent root bending moments are not suitable for stall‐regulated wind turbines: when the cost function is minimized, the edgewise damping can be low, and the flapwise damping can approach zero during an extreme operating gust. A new cost function is proposed that leads to nicely balanced stall behavior and damping over the entire operating windspeed range. The method was used to design the blades of two multi‐MW, stall‐regulated, offshore wind turbines, comparable with the NREL 5 MW and NTNU 10 MW pitch‐regulated turbines. It is shown that the optimal stall‐regulated blade has a unique aerodynamic profile that gives high flapwise and edgewise damping and a uniform mean power output above the rated windspeed. The blades are described in sufficient detail that they can be used in further aeroelastic analyses, to compare large stall‐regulated and pitch‐regulated turbines. Copyright © 2014 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.     

凹槽-襟翼对翼型动态失速特性影响研究

风力机复杂运行环境使叶片常处于失速环境,导致翼型升力骤降,严重影响风力机气动性能.为改善翼型流动分离,延缓失速,对凹槽-襟翼对翼型动态失速特性作用效果开展研究,并利用计算流体力学方法分析不同折合频率与翼型厚度时凹槽-襟翼对翼型气动性能的影响.结果表明:俯仰振荡过程中,凹槽-襟翼可有效提升翼型吸力面流速,降低失速攻角下逆...     

Self‐induced vibrations of a DU96‐W‐180 airfoil in stall

This work presents an analysis of two‐dimensional (2D) and three‐dimensional (3D) non‐moving, prescribed motion and elastically mounted airfoil computational fluid dynamics (CFD) computations. The elastically mounted airfoil computations were performed by means of a 2D structural model with two degrees of freedom. The computations aimed at investigating the mechanisms of both vortex‐induced and stall‐induced vibrations related to a wind turbine blade at standstill conditions. In this work, a DU96‐W‐180 airfoil was used in the angle‐of‐attack region potentially corresponding to stall‐induced vibrations. The analysis showed significant differences between the aerodynamic stability limits predicted by 2D and 3D CFD computations. A general agreement was reached between the prescribed motion and elastically mounted airfoil computations. 3D computations indicated that vortex‐induced vibrations are likely to occur at modern wind turbine blades at standstill. In contrast, the predicted cut‐in wind speed necessary for the onset of stall‐induced vibrations appeared high enough for such vibrations to be unlikely. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of vortex‐induced and stall‐induced vibrations at standstill conditions using a free wake aerodynamic code

Vortex‐induced and stall‐induced vibrations of a 2D elastically mounted airfoil at high angles of attack in the vicinity of 90° are investigated using a vortex type model. Such conditions are encountered in parked or idling operation at extreme yaw angles provoked by control system failures. At very high angles of attack, massive flow separation takes place over the entire blade span, and vortex shedding evolves downstream of the blade giving rise to periodically varying loads at frequencies corresponding to the Strouhal number of the vortices shed in the wake. As a result, vortex‐induced vibrations may occur when the shedding frequency matches the natural frequency of the blade. A vortex type model formulated on the basis of the ‘double wake’ concept is employed for the modelling of the stalled flow past a 2D airfoil. By tuning the core size of the vortex particles in the wake, the model predictions are successfully validated against averaged 2D measurements on a DU‐96‐W‐180 airfoil at high angles of attack. In order to assess the energy fed to the airfoil by the aerodynamic loads, the behaviour under imposed sinusoidal edgewise motions is analysed for various oscillation frequencies and amplitudes. Moreover, stall‐induced and vortex‐induced vibrations of an elastically mounted airfoil section are assessed. The vortex model predicts higher aeroelastic damping as compared with that obtained using steady‐state aerodynamics. Excessive combined vortex‐induced and stall‐induced edgewise vibrations are obtained beyond the wind speed of 30 m s^?1. Copyright © 2014 John Wiley & Sons, Ltd.     

Vortex‐induced vibrations of a DU96‐W‐180 airfoil at 90° angle of attack

This work presents an analysis of vortex‐induced vibrations of a DU96‐W‐180 airfoil in deep stall at a 90° angle of attack, based on 2D and 3D Reynolds Averaged Navier Stokes and 3D Detached Eddy Simulation unsteady Computational Fluid Dynamics computations with non‐moving, prescribed motion and elastically mounted airfoil suspensions. Stationary vortex‐shedding frequencies computed in 2D and 3D Computational Fluid Dynamics differed. In the prescribed motion computations, the airfoil oscillated in the direction of the chord line. Negative aerodynamic damping, found in both 2D and 3D Computational Fluid Dynamics computations with moving airfoil, showed in the vicinity of the stationary vortex‐shedding frequency computed by 2D Computational Fluid Dynamics. A shorter time series was sufficient to verify the sign of the aerodynamic damping in the case of the elastic computations than the prescribed motion. Even though the 2D computations seemed to be capable of indicating the presence of vortex‐induced vibrations, the 3D computations seemed to reflect the involved physics more accurately. Copyright © 2013 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.     

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.     

Time‐accurate blade surface static pressure behaviour on a rotating H‐Darrieus wind turbine

Time‐accurate blade pressure distributions on a rotating H‐Darrieus wind turbine at representative tip speed ratios during start‐up are presented here, which allow blade dynamic stall and laminar separation bubbles to be observed clearly and which provide a rare experimental demonstration of the flow curvature effect inherent in H‐Darrieus turbine operation. The convection of a dynamic stall vortex along the blade surface at high reduced frequency has also been clearly identified. This study provides new information of the complex aerodynamics of the vertical axis wind turbines (VAWTs) and provides unique experimental data to validate the transient blade static surface pressure distribution predicted by CFD models. To the best of the authors' knowledge, this is the first time that the instantaneous pressure variation around the blade has been measured and recorded directly for an H‐Darrieus wind turbine.     

Dynamic stall modelling of the S809 aerofoil and comparison with experiments

An analysis of dynamic stall for the S809 aerofoil has been performed in conjunction with the Leishman–Beddoes dynamic stall model that was modified for wind turbine applications. Numerical predictions of the lift, drag and pitching moment coefficients were compared with measurements obtained for an oscillating S809 aerofoil at various reduced frequencies, mean angles of attack and angle of attack amplitudes. It was found that the results using the modified model were in good agreement with the experimental data. Hysteresis in the aerodynamic coefficients was captured well, although the drag coefficient was slightly underpredicted in the deep stall flow regime. Validation against the experimental data showed overall good agreement. The mathematical structure of the model is such that it can be readily incorporated into a comprehensive analysis code for a wind turbine. Copyright © 2006 John Wiley & Sons, Ltd.     

Insight into wind turbine stall and post‐stall aerodynamics

The objective of this study was to evaluate measured NASA Ames Unsteady Aerodynamic Experiment post‐stall blade element data and to provide guidelines for developing an empirical approach that predicts post‐stall aerofoil characteristics. Blade element data were analysed from the five radial stations of the baseline 5·03 m radius rotor. A lifting surface/prescribed wake performance prediction method was used to determine a reference angle of attack that corresponds to the measured blade element data. Using the measured normal and tangential force coefficients and estimated angle of attack, spanwise distributions of lift and drag performance characteristics were derived along with the circulation distributions. Guidelines for a new stall and post‐stall model based on the measured trends in the aerofoil performance characteristics, along with flat plate theory, are proposed for predicting the peak and post‐peak power. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Comparison of pitching and plunging effects on the surface pressure variation of a wind turbine blade section

Numerous experiments were conducted on an oscillating airfoil in a subsonic wind tunnel. The experiments involved measuring the surface pressure distribution when the model oscillated in two types of motion, pitch and plunge, at three different Reynolds numbers, 0.42, 0.63 and 0.84 million, and over a range of reduced frequencies, k = 0.03–0.09. The unsteady aerodynamic loads were calculated from the surface pressure measurements, 64 ports, along the chord for both upper and lower surfaces of the model. Particular emphasis was placed on the effects of different types of motion on the unsteady pressure distribution of the airfoil at pre‐stall, near‐stall and post‐stall conditions. It was found that variations of the pressure distribution and aerodynamic loads with angle of attack were strongly sensitive to the displacement, oscillation frequency and mean angle of attack. The width of the hysteresis loop, position of the ‘figure‐8 shape’ and slope of the pressure coefficient curve are influenced by both types of motion, pitch and plunge. The main difference between plunging and pitching motions is due to the presence of the pitch rate for the pitching motion case, which was absent in the plunging case. Pitch rate had the strongest influence on pressure data in the near‐stall and post‐stall conditions. The trend of increasing the width of the hysteresis loops of lift coefficients with changing reduced frequency was different in two motions in the pre‐stall and post‐stall regions. The aerodynamic damping was greater for the pitching case than for the plunging one at higher reduced frequencies due to the existence of the pitch rate in the pitching oscillation, which was reversed at lower reduced frequencies. In the near‐stall region, at higher reduced frequency, the dynamic stall angle for the pitching oscillation increased while for the plunging one the effect was minimal. Increasing the oscillation amplitude was more effective for the plunging motion than for the pitching one. The effects of surface grit roughness on the pressure signature for both types of motion were also investigated. Applying the surface roughness near the leading edge affected the performance of the airfoil significantly. Copyright © 2008 John Wiley & Sons, Ltd.     

水平轴失速型风力机主动非线性控制

讨论了大型主动失速型风力机在额定工况以上时的主动非线性控制问题。直接利用所推导的仿射性非线性模型．采用微分几何精确线性化理论,实现恒速风力机全局精确线性化控制,给出了反馈控制算法,并对闭环系统进行了数字仿真。     

Model improvements for evaluating the effect of tower tilting on the aerodynamics of a vertical axis wind turbine

If a vertical axis wind turbine is mounted offshore on a semi‐submersible, the pitch motion of the platform will dominate the static pitch and dynamic motion of the platform and wind turbine such that the effect of tower tilting on the aerodynamics of the vertical axis wind turbine should be investigated to more accurately predict the aerodynamic loads. This paper proposes certain modifications to the double multiple‐streamtube (DMS) model to include the component of wind speed parallel to the rotating shaft. The model is validated against experimental data collected on an H‐Darrieus wind turbine in skewed flow conditions. Three different dynamic stall models are also integrated into the DMS model: Gormont's model with the adaptation of Strickland, Gormont's model with the modification of Berg and the Beddoes–Leishman dynamic stall model. Both the small Sandia 17 m wind turbine and the large DeepWind 5 MW are modelled. According to the experimental data, the DMS model with the inclusion of the dynamic stall model is also well validated. On the basis of the assumption that the velocity component parallel to the rotor shaft is small in the downstream part of the rotor, the effect of tower tilting is quantified with respect to power, rotor torque, thrust force and the normal force and tangential force coefficients on the blades. Additionally, applications of Glauert momentum theory and pure axial momentum theory are compared to evaluate the effect of the velocity component parallel to the rotor shaft on the accuracy of the model. Copyright © 2013 John Wiley & Sons, Ltd.     

Plasma‐based feed‐forward dynamic stall control on a vertical axis wind turbine

Dynamic stall was controlled on a double‐bladed H‐Rotor vertical axis wind turbine model using pulsed dielectric barrier discharge plasma actuators in a feed‐forward control configuration. The azimuthal angles of plasma actuation initiation and termination, that produced the largest increases in power, were determined parametrically on the upstream half of the turbine azimuth in a low‐speed blow‐down wind tunnel at wind speeds of 7 m/s. A mathematical model, together with instantaneous turbine speed, was used to estimate transient torque and power developed by the turbine under the influence of plasma actuation. Overall performance improvements were based on changes between the final actuated and initial baseline results. A remarkable result of this investigation was that a net turbine power increase of 10% was measured. This was achieved by systematically reducing plasma pulsation duty cycles as well as the plasma initiation and termination angles. Nevertheless, it was determined that further performance increases could be achieved by changing the actuator's dielectric material, increasing the turbine radius and developing a method for control of dynamic stall on both the upwind (inboard of the blades) and downwind (outboard of the blades) halves of the turbine azimuth. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

考虑转捩的风力机翼型动态失速数值模拟

以风力机专用翼型的动态失速为对象,采用一种基于流场当地变量的Gamma-Theta转捩模型配合SSTk-ω湍流模型进行数值模拟,研究转捩对动态失速性能的影响和动态失速下的转捩规律。结果表明,使用考虑转捩效应,能够使动态失速过程中上仰段大迎角状态下失速和下俯段气流再附的模拟得到改善。在动态失速上仰段,上表面转捩由后缘分离泡向前缘分离泡的转变过程较快,导致转捩点迅速前移;而在下俯段,前缘分离泡向后缘分离泡的转变过程中经过了自然转捩和再层流化的过渡,因此转捩点的移动较上仰段平滑。     

对一种动态失速模型的改进及精度研究

基于现有的Hopf分岔法动态失速模型（Hopf bifurcation model）,引入Wagner函数计算其附着流下等效攻角,对原模型的边界层再附着项进行一定修改,使新模型可表示为状态空间的形式,并为原模型补齐了对于阻力和力矩系数的建模。相比于常见的ONERA和Leishman-Beddoes动态失速模型,新模型在附着流下拥有与解析理论最一致的幅值及相位特性;分离流下,新模型在大部分情况下的计算精度优于常见模型,且能更好地捕捉初级失速涡和深度失速下出现的多级失速涡现象。其中对轻度失速的分析表明,各动态失速模型在轻度失速下的环量项建模仍具有一定的提升空间。     

Digital tuft analysis of stall on operational wind turbines

Operational wind turbines are exposed to dynamic inflow conditions because of, for instance, atmospheric turbulence and wind shear. In order to understand the resulting three‐dimensional and transient aerodynamics effects at a site, a 10m stall‐regulated upwind two‐bladed wind turbine was instrumented for a novel digital tuft flow visualization study. High definition video of a tufted blade was acquired during wind turbine operation in the field, and a novel digital image processing algorithm calculated the blade stall directly from the video. After processing O(10⁵) sequential images, the algorithm achieved a ?5% bias error compared with previous manual analysis methods. With increasing wind speed (5m/s to 20m/s) the fraction of tufts exhibiting stalled flow increased from 5% to 40% on the outboard 40% of the blade. The independently measured instantaneous turbine power production correlates highly with the stall fraction. Some azimuthal variation in the stall fraction associated with dynamic stall induced by vertical wind shear was seen with a maximum in the 45–90° azimuthal location. The high detail, quantitative image processing method demonstrated good agreement with the expected behaviour for a stall‐regulated wind turbine and revealed azimuthal variation because of shear‐induced dynamic stall. The amount of reliable blade stall data to be obtained from digital tuft visualization has hereby been vastly increased. Copyright © 2015 John Wiley & Sons, Ltd.