Integrating a new flow separation model and the effects of the vortex shedding for improved dynamic stall predictions using the Beddoes–Leishman method

This paper aims at improving dynamic stall predictions on the S809 aerofoil under 2D flow conditions. The method is based on the well‐known Beddoes–Leishman model; however, a new flow separation model and a noise generator are integrated to improve the predictions in the load fluctuations, including those induced by vortex shedding on the aerofoil upper surface. The flow separation model was derived from a unique approach based on the combined use of unsteady aerodynamic loads measurements, the Beddoes–Leishman model and a trial‐and‐error technique. The new flow separation model and random noise generator were integrated in the Beddoes–Leishman model through a new solution algorithm. The numerical predictions of the unsteady lift and drag coefficients were then compared with the Ohio State University measurements for the oscillating S809 aerofoil at several reduced frequencies and angles of attack. The results using the proposed models showed improved correlation with the experimental data. Hysteresis loops for the aerodynamic coefficients are in good agreement with measurements. Copyright © 2016 John Wiley & Sons, Ltd.     

风力机叶片翼型气动性能数值模拟

采用数值模拟方法对NACA23012,NACA4412,S809,S810等4种常用风力机叶片翼型进行了研究,分析了翼型静止与振荡时的气动性能.随着攻角的增加,静止翼型的升力系数先增大后减小,其阻力系数一直增大,显示出NACA4412翼型具有较好的低风速启动性能;振荡翼型的升力系数随着攻角的变化呈现一个闭合迟滞环曲线,显示出振荡翼型S809的动态失速迟滞效应最为明显.文章参照模拟结果和对比试验数据,验证了数值模拟的可靠性.     

A new stall delay algorithm for predicting the aerodynamics loads on wind turbine blades for axial and yawed conditions

Stall delay is a complicated phenomenon that has gained for many years the attention of industry and academics in the fields of helicopter and wind turbine aerodynamics. Since most of the potential flow theories still rely on the use of 2D aerofoil data for simulating loads on a rotating blade, less degree of accuracy is expected because of 3D rotational effects. In this work, a new model for correcting the 2D steady aerodynamic data for 3D effects is presented. The model can reduce the uncertainty in the blade design process and, subsequently, make wind turbines more cost‐effective. This model combines the stall delay model of Corrigan and Schillings, a modified version of an inviscid stall delay model, a new modification factor to account for the effect of the angle of attack changes and a new tip loss factor. Furthermore, the model applies the use of the separation factor of Du and Selig to evaluate the area on the rotor disc where stall delay is most prominent. The new stall delay model was embedded in a free‐wake vortex model to estimate the aerodynamic loads on the National Renewable Energy Laboratory Phase VI rotor blades consisting of the S809 aerofoil sections. The results in this study confirm the validity of the 3D corrections by the proposed new model under both axial and yawed flow conditions. Copyright © 2017 John Wiley & Sons, Ltd.     

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

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

风力机翼型等速上仰动态失速数值模拟

采用κ-ωSST模型,利用CFD软件模拟了NREL S809翼型正弦振荡动态失速,并将结果和俄亥俄州立大学（OSU）风洞试验值对比,显示出较好的一致性,验证了所用方法的有效性．在此基础上对该翼型在雷诺数Re=1．0×10^6时以攻角变化率α=34．54（＊）·s^-1等速上仰动态失速过程进行了数值模拟,详细描述了等速上仰动态失速过程涡的发展以及翼型周围流场的分布．结果表明,动态失速现象是由前缘主涡和尾缘逆向涡交替作用引起;其气动特性曲线的分析结果表明,其失速前气动性能较静态时有较大提升．     

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.     

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.     

A vortex panel model for the simulation of the wake flow past a vertical axis wind turbine in dynamic stall

A 2D vortex panel model with a viscous boundary layer formulation has been developed for the numerical simulation of a vertical axis wind turbine (VAWT), including the operation in dynamic stall. The model uses the ‘double wake’ concept to reproduce the main features of the unsteady separated flow, including the formation and shedding of strong vortical structures and the wake–blade interaction. The potential flow equations are solved together with the integral boundary layer equations by using a semi‐inverse iterative algorithm. A new criterion for the reattachment of the boundary layer during the downstroke of a dynamically stalled aerofoil is implemented. The model has been validated against experimental data of steady aerofoils and pitching aerofoils in dynamic stall at high and low Reynolds numbers (Re = 1.5 × 10⁶ and Re = 5 × 10⁴). For the low Reynolds number case, time‐resolved 2D particle image velocimetry (PIV) measurements have been performed on a pitching NACA 0012 aerofoil in dynamic stall. The PIV vorticity fields past the oscillating aerofoil are used to test the model capability of capturing the formation, growth and release of the strong leading edge vortex that characterizes the dynamic stall. Furthermore, the forces extracted from the PIV velocity fields are compared with the predicted ones for a quantitative validation of the model. Finally, the model is applied to the computation of the wake flow past a VAWT in dynamic stall; the predicted vorticity fields and forces are in good agreement with phase‐locked PIV data and CFD‐DES available in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Load control on a dynamically pitching finite span wind turbine blade using synthetic jets

The feasibility of active flow control, via arrays of synthetic jet actuators, to mitigate hysteresis was investigated experimentally on a dynamically pitching finite span S809 blade. In the present work, a six‐component load cell was used to measure the unsteady lift, drag and pitching moment. Stereoscopic Particle Image Velocimetry (SPIV) measurements were also performed to understand the effects of synthetic jets on flow separation during dynamic pitch and to correlate these effects with the forces and moment measurements. It was shown that active flow control could significantly reduce the hysteresis in lift, drag and pitching moment coefficients during dynamic pitching conditions. This effect was further enhanced when the synthetic jets were pulsed modulated. Furthermore, additional reduction in the unsteady load oscillations can be observed in post‐stall conditions during dynamic motions. This reduction in the unsteady aerodynamic loading can potentially lead to prolonged life of wind turbine blades. Copyright © 2014 John Wiley & Sons, Ltd.     

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

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

Modeling dynamic stall on wind turbine blades under rotationally augmented flow fields

This paper presents an investigation of two well‐known aerodynamic phenomena, rotational augmentation and dynamic stall, together in the inboard parts of wind turbine blades. This analysis is carried out using the following: (1) the National Renewable Energy Laboratory's Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation; (2) data from unsteady delayed detached eddy simulations (DDES) carried out using the Technical University of Denmark's in‐house flow solver Ellipsys3D; and (3) data from a reduced order dynamic stall model that uses rotationally augmented steady‐state polars obtained from steady Phase VI experimental sequences, instead of the traditional two‐dimensional, non‐rotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared with three select cases of the N‐sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared with those from the dynamic stall model. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in two‐dimensional flow to be investigated. Results indicated a good qualitative agreement between the model and the experimental data in many cases, which suggests that the current two‐dimensional dynamic stall model as used in blade element momentum‐based aeroelastic codes may provide a reasonably accurate representation of three‐dimensional rotor aerodynamics when used in combination with a robust rotational augmentation model. Copyright © 2015 John Wiley & Sons, Ltd.     

A new correlation on the MEXICO experiment using a 3D enhanced blade element momentum technique

The blade element momentum (BEM) theory is based on the actuator disc (AD) model, which is probably the oldest analytical tool for analysing rotor performance. The BEM codes have very short processing times and high reliability. The problems of the analytical codes are well known to the researchers: the impossibility of describing inside the one-dimensional code the three-dimensional (3D) radial flows along the span-wise direction. In this work, the authors show how the 3D centrifugal pumping affects the BEM calculations of a wind turbine rotor. Actually to ascertain the accuracy of the analytical codes, the results are compared with rotor performance, blade loads and particle image velocimetry measurements of the model experiment in controlled conditions. A reliable agreement with the measurement is obtained. A good improvement is gained when the blade stall state modified aerofoil data instead of the original aerofoil data are used in the calculations.     

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.     

An improved BEM model for the power curve prediction of stall‐regulated wind turbines

Blade element momentum (BEM) theory is the standard computational technique for the prediction of power curves of wind turbines; it is based on the two‐dimensional aerodynamic properties of aerofoil blade elements and some corrections accounting for three‐dimensional wing aerodynamics. Although most BEM models yield acceptable results for low‐wind and pitch‐controlled regimes where the local angles of attack are small, no generally accepted model exists up to date that consistently predicts the power curve in the stall regime for a variety of blade properties and operating conditions. In this article we present a modified BEM model which satisfactorily reproduces the power curves of four experimental wind turbines reported in the literature, using no free fit parameters. Since these four experimental cases comprehend a great variety of conditions (wind tunnel vs field experiments, different air densities) and blade parameters (no twist and no taper, no taper but twist, both twist and taper, different aerofoil families), it is believed that our model represents a useful working tool for the aerodynamic design of stall‐regulated wind turbines. Copyright © 2004 John Wiley & Sons, Ltd.     

Methods to control dynamic stall for wind turbine applications

Dynamic stall (DS) on a wind turbine is encountered when the sectional angles of attack of the blade rapidly exceeds the steady-state stall angle of attack due to in-flow turbulence, gusts and yaw-misalignment. The process is considered as a primary source of unsteady loads on wind turbine blades and negatively influences the performance and fatigue life of a turbine. In the present article, the control requirements for DS have been outlined for wind turbines based on an in-depth analysis of the process. Three passive control methodologies have been investigated for dynamic stall control: (1) streamwise vortices generated using vortex generators (VGs), (2) spanwise vortices generated using a novel concept of an elevated wire (EW), and (3) a cavity to act as a reservoir for the reverse flow accumulation. The methods were observed to delay the onset of DS by several degrees as well as reduce the increased lift and drag forces that are associated with the DSV. However, only the VG and the EW were observed to improve the post-stall characteristics of the airfoil.     

基于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种翼型动态失速预测结果与公开实验数据结论一致,证实所提出的动态失速气动模型计算结果准确可信,具有较强通用性。     

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

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

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.     

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

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

Performance analysis of the airfoil-slat arrangements for hydro and wind turbine applications

Standard airfoils historically used for wind and hydrokinetic turbines had maximum lift coefficients of around 1.3 at stall angles of attack, which is about 12°. At these conditions, the minimum flow velocities to generate electric power were about 7 m/s and 2 m/s for the wind turbine and the hydrokinetic turbine cases, respectively. In this study, NACA4412-NACA6411 slat–airfoil arrangement was chosen for these two cases in order to investigate the potential performance improvements. Aerodynamic performances of these cases were both numerically and experimentally investigated. The 2D and 3D numerical analysis software were used and the optimum geometric and flow conditions leading to the maximum power coefficient or the maximum lift to drag ratio were obtained. The maximum lift to drag ratio of 24.16 was obtained at the optimum geometric and flow parameters. The maximum power coefficient of 0.506 and the maximum torque were determined at the tip speed ratios of 5.5 and 4.0 respectively. The experimental work conducted in a towing tank gave the power coefficient to be 0.47 which is about %7 lower than the numerical results obtained. Hence, there is reasonable agreement between numerical end experimental values. It may be concluded that slat-hydrofoil or airfoil arrangements may be applied in the design of wind and hydrokinetic turbines for electrical power generation in lower wind velocities (3–4 m/s) and current velocities (about 1 m/s).