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.     

On the S809 airfoil's unsteady aerodynamic characteristics

The S809 airfoil dynamic characteristics, which are based on the airfoil dynamic tests at the University of Glasgow, are presented in this paper. The airfoil tests include static, ramp‐type (ramp‐up and ramp‐down) and oscillatory motions at Reynolds numbers of 1.0 × 10⁶ and 1.5 × 10⁶ with and without the sand‐tripped leading edge. This study aimed to explore the unsteady aerodynamic features of the S809 airfoil, such as the progression of separation from leading edge to trailing edge, the large trailing edge separation before stalling, the stall onset inception and the re‐attachment convection, and to provide some useful data for tuning/refining the semi‐empirical dynamic stall (DS) models, such as the Leishman–Beddoes DS model or its variations. Experimental results show that the S809 airfoil has a complicated DS process that renders this airfoil a challenge for any modeller of the unsteady airloads. The leading‐edge roughness has small effect on the static features, but significantly invokes earlier stall onset inception under dynamic conditions, while has small influence on the convective phase of the re‐establishment of fully attached flow. 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.     

水平轴风力机翼型大攻角分离流动的数值模拟

翼型的失速特性是失调速节型水风力机的气动性能分析和颤振分析的基础，许多涉及这类问题的研究大多只给出了翼型刚开始失速时的计算结果，然而在正常运行工况下叶片端部翼型的深失速特性是风力机的最关键的一类问题，通过求解二维非常，可压的Ｎ－Ｓ方程计算了风力机常用翼型ＮＡＣＡ４４１８的绕流特性，Ｎ－Ｓ方程在贴体坐标系中给出用Ｐｏｉｓｓｏｎ方程法生成了Ｃ型网格，数值计算了采用了一种改进的ＬＵ－ＳＧＳ格式。将翼型的     

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.     

Dynamic stall analysis using harmonic balance and correlation‐based γ– transition models for wind turbine applications

In this paper, we use the harmonic balance method to study an oscillating S809 airfoil in dynamic stall. The periodic behavior of this problem makes it well suited for the harmonic balance method, which is able to model unsteady aerodynamics at greatly reduced computational costs when compared with time‐accurate unsteady‐flow solvers. A finite‐volume technique based on the lower–upper symmetric Gauss–Seidel scheme with Roe fluxes is used to solve the Reynolds‐averaged Navier–Stokes equations. The turbulent viscosity is computed with the one‐equation Spalart–Allmaras turbulence model. In addition, the laminar–turbulent transition is modeled using a correlation‐based approach originally developed by Langtry and Menter. Comparisons with experimental data for steady flows with the S809 airfoil highlight the necessity of the transition model to accurately predict the onset of static stall. For unsteady cases, the transition model provides improved agreement with experimental data, predicting dynamic stall when the fully turbulent model cannot. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamic stall control via adaptive blowing

An aerodynamic load control concept termed “adaptive blowing” was successfully tested on a NACA 0018 airfoil model at Reynolds numbers ranging from 1.5·10⁵ to 5·10⁵. The global objective was to eliminate lift oscillations typically encountered on wind turbine blade sections. Depending on the jet momentum flux, steady blowing from a control slot in the leading-edge region can be utilized to either enhance or reduce lift by suppressing or inducing boundary layer separation respectively. Furthermore, high momentum blowing effectively eliminated the dynamic stall vortex during deep dynamic stall conditions. Based on these previous findings, the present work explores the feasibility of controlling unsteady aerodynamic loads by dynamically varying the jet momentum flux to compensate for transient changes of the inflow. Various scenarios including high amplitude pitching, rapid freestream oscillations and combinations of both were investigated in a custom-built unsteady wind tunnel facility. An iterative control algorithm was implemented which successfully identified the momentum coefficient time profiles required to minimize the lift excursions. The combination of fully suppressing dynamic stall and dynamically adjusting the lift coefficient provided an unprecedented control authority, producing virtually constant phase averaged lift in all cases.     

Aerodynamic response of an airfoil section undergoing pitch motion and trailing edge flap deflection: a comparison of simulation methods

The study presents and compares aerodynamic simulations for an airfoil section with an adaptive trailing edge flap, which deflects following a smooth deformation shape. The simulations are carried out with three substantially different methods: a Reynolds‐averaged Navier–Stokes solver, a viscous–inviscid interaction method and an engineering dynamic stall model suitable for implementation in aeroelastic codes based on blade element momentum theory. The aerodynamic integral forces and pitching moment coefficients are first determined in steady conditions, at angles of attack spanning from attached flow to separated conditions and accounting for the effects of flap deflection; the steady results from the Navier–Stokes solver and the viscous–inviscid interaction method are used as input data for the simpler dynamic stall model. The paper characterizes then the dynamics of the unsteady forces and moments generated by the airfoil undergoing harmonic pitching motions and harmonic flap deflections. The unsteady aerodynamic coefficients exhibit significant variations over the corresponding steady‐state values. The dynamic characteristics of the unsteady response are predicted with an excellent agreement among the investigated methods at attached flow conditions, both for airfoil pitching and flap deflection. At high angles of attack, where flow separation is encountered, the methods still depict similar overall dynamics, but larger discrepancies are reported, especially for the simpler engineering method. 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.     

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.     

Computational fluid dynamics simulation of the aerodynamics of a high solidity,small‐scale vertical axis wind turbine

A computational fluid dynamics simulation was performed for a small‐scale, high solidity (σ = 0.48) H‐type Darrieus vertical axis wind turbine. Two‐dimensional unsteady Reynolds‐averaged Navier–Stokes equations were solved for the turbine numerical model, which has a large stationary domain and smaller rotating subdomain connected by a sliding mesh interface. The simulation results were first validated against steady‐state airfoil data. The model was then used to solve for three rotating blades with constant ambient flow velocity (Re = 360,000) over numerous blade speed ratios. The high solidity and the associated low blade speed ratio and rotational speed of the turbine result in complex flow–blade interaction mechanisms. These include dynamic stall resulting in vortex shedding, vortex impingement on the source blade and significant flow momentum extraction causing reduced power production from the downstream blade pass. Copyright © 2011 John Wiley & Sons, Ltd.     

基于Theodorsen理论的新型动态失速预测模型及其应用

风力机叶片动态失速时的非定常气动特性及严重的迟滞现象使得风力机功率实测值严重偏离其静态预测值。鉴于此,基于Theodorsen理论、基尔霍夫势流理论,在忽略低阶附加质量引起的下洗气流加速度项及状态变量转换后,提出一种包括翼型附着流和后缘动态分离流的新型动态失速模型。利用该模型分析NREL 5 MW海上风力机叶片6种翼型的非定常动态失速特性得出：通过翼型的气流在完全附着流与完全分离流之间不断转换,受附着流脱落尾诱导的动态下洗气流影响及边界层动态分离产生的压力滞后的双重作用,动态升力系数变化曲线和静态升力现象曲线偏差较大,6种翼型动态升力系数变化曲线均呈非常明显的迟滞环现象。DU40、DU35、DU30、DU25、DU21和NACA64这6种翼型动态升力系数增幅明显,分别达17.6%、60.9%、60.7%、55.1%、63.7%和40.8%。动态失速攻角极大地超过静态失速攻角,分别增大到36.53°、21.40°、20.20°、17.68°、16.97°和21.42°。6种翼型动态失速预测结果与公开实验数据结论一致,证实所提出的动态失速气动模型计算结果准确可信,具有较强通用性。     

Effects of nonuniform incident velocity on a dynamic wind turbine airfoil

For further insight into the performance of a horizontal axis wind turbine blade section under yaw loads, a 2D numerical simulation of a pitching S809 airfoil under dynamic stall with an unsteady incident velocity is presented. The streamwise incident velocity and pitch angle of incidence oscillated with the same frequency but with a range of phase differences, ? π ≤ Φ ≤ π. Changing Φ caused variation of the results, which can be highlighted as significantly augmented and dramatically damped dynamic stall loads, both increasing and decreasing trends for vortex growth time during Φ increase, a shifted location of the maximum loads and a change in the order of the vortex pair circulation in each cycle. The results showed strong dependency on the velocity and acceleration of the freestream during dynamic stall, which categorized the results in four individual subdomains with different behaviors. Copyright © 2014 John Wiley & Sons, Ltd.     

直接测力俯仰振荡翼型动态气动性能研究

在西北工业大学NF-3低速风洞二元实验段开展翼型俯仰振荡运动动态气动性能深入研究。实验模型为展向三段式测力模型,测力仅在模型中段进行以减小风洞侧壁干扰的影响。实验中采集模型的转动瞬态迎角、计算模型中段的惯性力和惯性力矩、并从天平采集数据中扣除以修正模型惯性对结果的影响。结果表明,迎角超过正向或负向静态失速迎角是升力系数和俯仰力矩系数产生大的迟滞环的必要条件。随着振荡缩减频率增大,动态失速会推迟,升力系数迟滞环增大,阻力系数增大,最大迎角附近的俯仰力矩系数减小。在迎角小于静态失速迎角或超过不大的迎角范围,随着缩减频率的增大,翼型振荡运动俯仰力矩系数上行时减小,下行时增大。随着振荡振幅的增大,翼型振荡运动动态升力系数和俯仰力矩系数的迟滞环增大。随着平均迎角的增大,翼型迎角更多地进入正向失速区,升力系数迟滞环增大,俯仰力矩系数最小值变小。雷诺数对升力系数、阻力系数和俯仰力矩系数迟滞环无明显影响;但是,在翼型模型下行过程,随着雷诺数的增大,升力恢复提前,同时迟滞环随雷诺数增大而减小。     

风力机叶片S832翼型的两种湍流模型数值分析方法及其比较

采用S-A和RNGK-ε湍流模型对风力机专用S832翼型的绕流流动建立了二维不可压缩湍流模型,利用计算流体力学软件Fluent,对两种模型进行数值模拟,得到了雷诺数为3×10^6时该翼型在-16°～30°攻角下的升力系数和阻力系数与来流攻角的关系以及压力分布图,并进一步分析了不同攻角下翼型表面压力分布特性,预测了大攻角（达30°）下翼型分离流动特性。结合NREL的试验数据,对两种湍流模型模拟的精度进行了分析比较,结果表明在小攻角范围内采用RNGK-~模型预测该翼型气动性,其结果更加有效。     

基于两种湍流模型的某风机专用翼型数值研究

针对风力机专用S832翼型绕流流动建立了二维不可压缩湍流模型,利用计算流体力学软件Fluent,分别选用S—A、RNGk-ε两种湍流模型对S832进行数值模拟,对比了两种湍流模型对气动模拟精度的影响,得出了雷诺数为3×10^6时,该翼型在-16°～30°攻角下的升力系数和阻力系数随来流攻角的变化关系及压力分布图,分析了不同攻角下翼型表面压力分布特性并进一步预测了大攻角（达30°）下翼型分离流动特性．并与NREL的试验数据进行比较,研究结果表明：RNGk—ε在预测该翼型小攻角范围气动性方面更加有效。     

Wind turbine blade vibration at standstill conditions—the effect of imposing lag on the aerodynamic response of an elastically mounted airfoil

The present study investigated physical phenomena related to stall‐induced vibrations potentially existing on wind turbine blades at standstill conditions. The study considered two‐dimensional airfoil sections while it omitted three‐dimensional effects. In the study, a new engineering‐type computational model for the aeroelastic response of an elastically mounted airfoil was used to investigate the influence of temporal lag in the aerodynamic response on the aeroelastic stability in deep stall. The study indicated that even a relatively low lag significantly increases the damping of the model. A comparison between the results from a model with lag imposed on all force components with the results from a model with lag imposed exclusively on the lift showed only marginal difference between the damping in the two cases. A parameter study involving positions of the elastic hinge point and the center of gravity indicated that the stability is relatively independent of these parameters. Another parameter study involving spring constants showed that the stability of each mode is dependent only on the spring constant acting in the direction of the leading motion of the mode. An investigation of the influence of the added mass terms showed that only the pitch‐rate and flapwise‐acceleration terms have any influence on the stability. An investigation of three different profiles showed that the stability is heavily dependent on the aerodynamic characteristics of the profiles—mainly on the lift. It was also shown that only the edgewise mode is unstable in deep stall. Moreover, independent of the amount of temporal lag in the aerodynamic response of the model, the inflow‐angle region in the vicinity of 180° remains unstable in the edgewise mode. Therefore, this inflow‐angle region may create stability problems in real life. The other type of vibrations potentially present at standstill conditions is vortex‐induced, being outside the scope of the present study. Copyright © 2014 John Wiley & Sons, Ltd.     

定常吸气对垂直轴风力机性能影响研究

定常吸气装置可有效提高垂直轴风力机气动性能,改善风轮流场结构及翼型动态失速特性.基于CFD方法对垂直轴风力机进行数值模拟,研究不同叶尖速比(TSR)下定常吸气对风力机气动及流场特性的影响,对比分析原始风力机及定常吸气作用下的风能利用率、整机转矩系数及涡量分布.结果表明:不同尖速比下定常吸气均可显著提高风力机气动性能,减...     

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

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