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.     

尾缘气动弹片对翼型动态失速特性影响研究

通过研究尾缘气动弹片对翼型动态失速特性影响,提出一种基于气动弹片的主动控制策略,使其于大攻角时抬起,小攻角时闭合。并采用计算流体动力学方法对比分析主动式气动弹片对不同厚度翼型抑制流动分离作用的效果。结果表明：对于薄翼型,发生动态失速时,气动弹片可延缓翼型尾缘涡旋与前缘主流涡的相互作用,减小翼型升力系数骤降幅度;随翼型厚度增加,流动分离点从翼型前缘转向后缘,气动弹片可有效分割较大分离涡,减轻流动分离程度,限制分离涡发展,同时抑制尾缘伴随小涡产生,提高翼型升阻比。     

不同相对厚度前缘缝翼对S809翼型气动性能的影响

文章通过数值模拟方法研究了不同相对厚度的前缘缝翼对S809翼型气动性能的影响,并揭示了前缘缝翼相对厚度对流动控制产生影响的机理。研究结果表明:在大攻角下,空气流经过前缘缝翼会在其尾部产生涡旋,尾缘涡旋的形成有助于抑制S809翼型流动分离,进而改善翼型绕流场;不同相对厚度的前缘缝翼产生尾缘涡旋不同的流动轨迹,对翼型的流动控制作用效果不同;相同条件下,前缘安装最大相对厚度为35%的前缘缝翼能够将S809翼型最大升力系数提升至1.25,失速攻角推迟至17.21°;安装最大相对厚度为14%的前缘缝翼,能够使S809翼型最大升力系数提升至1.53,并使翼型在攻角为20.16°时仍未发生失速。     

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

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

Effect of acceleration on dynamic stall of airfoil in unsteady operating conditions

In this paper the effect of accelerated flow over a moving airfoil is considered and based on the flow field around the airfoil the dynamic stall is evaluated. In contrast to ordinary pitching motion, the dynamic stall evaluation in this study is performed with a different motion pattern, in which the airfoil has a heaving motion in one direction. This motion pattern is also similar to rotation of an element of blade in horizontal axis wind turbines (HAWTs). In present investigation, the Reynolds number is changed during simulation time and variations of this parameter from initial to final values are shown by acceleration parameter. The operating Reynolds number is more than 10⁶, and a S809 airfoil is selected to move with accelerations of 1, 4 and 6 m/s² in normal direction to free stream. To resolve accelerated flow filed in the two‐dimensional computational domain and to achieve results within a reasonable computation time, the unsteady Reynolds‐Averaged Navier–Stokes (URANS) equations are employed. The governing equations are discretized based on the finite volume approach and semi‐implicit method for pressure linked equations algorithm is used for pressure–velocity coupling. Furthermore, turbulence effect on flow field is accounted using shear stress transport (SST) k‐ω model. It is shown that the accelerated flow can significantly influence on the aerodynamic loads and dynamic stall trend. This study may introduce a new concept regarding dynamic stall and aerodynamic loads when the rotational acceleration is involved in HAWTs. Copyright © 2014 John Wiley & Sons, Ltd.     

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

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

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

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

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

前缘改形对翼型气动噪声特性影响

为得到高气动性能、低噪声的风力机专用翼型,基于参数化建模翼型,研究前缘外形对风力机翼型气动性能及气动噪声的影响规律。通过分离涡模拟方法和声学类比方程建立噪声预测方法。针对非对称翼型S809通过样条函数参数化处理前缘改形进行气动噪声计算。结果表明:翼型压力面前缘加厚,对翼型升阻力系数无明显影响,但大攻角时翼型周围压力分布均匀,流动相对稳定,且气动噪声声压级低于原始翼型,随压力面厚度增加气动噪声越大;吸力面加厚使得翼型升力系数增大,阻力系数减小,能抑制翼型失速时尾缘涡与前缘涡的生成,变形量越大气动噪声越小;翼型前缘上弯,翼型在失速区升力系数减小,阻力系数增大,流动越加不稳定,声压级随着攻角的增加呈递增趋势;翼型前缘下弯,翼型处于失速区升力系数增大,阻力系数减小,能抑制流动分离,未生成前缘涡和尾缘涡,当前缘下弯不变时,随加厚厚度增加翼型声压级呈减小趋势,且前缘下弯翼型声压级小于前缘上弯。     

Stall Behaviour in a Mixed-flow Compressor with Axial Slot Casing Treatment

Casing treatment is an effective technique in extending stall margin of axial and centrifugal compressor.However,its impacts on the stall behaviour of mixed-flow compressor are still not completely understood until now.To conquer this issue,unsteady full-annulus simulations were conducted to investigate the stall mechanism of a mixed-flow compressor with and without axial slot casing treatment(ASCT).The circumferential propagating speed of spike inception resolved by the numerical approach is 87.1%of the shaft speed,which is identical to the test data.The numerical results confirmed that the mixed-flow compressor fell into rotating stall via spike-type with and without ASCT.The flow structure of the spike inception was investigated at 50%design rotational speed.Instantaneous static pressure traces extracted upstream of the leading edge had shown a classic spiky wave.Furthermore,it was found that with and without ASCT,the mixed-flow compressor stalled through spike with the characteristic of tip leakage spillage at leading edge and tip leakage backflow from trailing edge,which is different from a fraction of the centrifugal compressor.The resultant phenomenon provides conoborating evidence for that unlike in axial-flow compressor,the addition of ASCT does not change the stall characteristics of the mixed-flow compressor.The flow structure that induced spike inception with ASCT is similar to the case with smooth casing.In the throttling process,tip leakage flow vortex had been involved in the formation of tornado vortices,with one end at the suction side,and the other end at the casing-side.The low-pressure region relevant to the downward spike is caused by leading-edge separation vortex or tornado vortex.The high-pressure region relevant to the upward spike is induced by blockage from the passage vortex.These results not only can provide guidance for the design of casing treatment in mixed-flow compressor,but also can pave the way for the stall waring in the highly-loaded compressors of next-generation aeroengines.     

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.     

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.     

Computational aerodynamic analysis of a blunt trailing‐edge airfoil modification to the NREL Phase VI rotor

The effects of modifying the inboard portion of the experimental NREL Phase VI rotor using a thickened, blunt trailing‐edge (or flatback) version of the S809 design airfoil are studied using a compressible, three‐dimensional, Reynolds‐averaged Navier–Stokes method. A motivation for using such a thicker airfoil design coupled with a blunt trailing edge is to alleviate structural constraints while reducing blade weight and maintaining the power performance of the rotor. The numerical results for the baseline Phase VI rotor are benchmarked against wind tunnel measurements obtained at freestream velocities of 5, 7 and 10ms^?1. The calculated results for the modified rotor are compared against those of the baseline rotor. The results of this study demonstrate that a thick, blunt trailing‐edge blade profile is viable as a bridge to connect structural requirements with aerodynamic performance in designing future wind turbine rotors. Copyright © 2007 John Wiley & Sons, Ltd.     

Wind tunnel quantification of dynamic stall on an S817 airfoil and its control using synthetic jet actuators

Wind tunnel experiments were performed to quantify the aerodynamic characteristics of the S817 airfoil in dynamic stall conditions, and the subsequent application of active flow control to modify the manner by which dynamic stall incepts. Both quasi‐2D and cantilevered finite span configurations were tested. Surface pressure, six‐component force‐torque sensor, and stereoscopic particle image velocimetry (SPIV) were used to quantify the baseline flow and the benefits of actuating synthetic jets (installed at x/c = 0.35, angled 45° into the flow, and at a momentum coefficient C_μ = 0.012). The airfoil was pitched at reduced frequencies of k_f = 0.025 and 0.05 and at shallow and deep stall. Vortex induced lift from dynamic stall was observed and was eliminated by the use of synthetic jets for nearly all conditions; pitching moment deviation was also observed to be significant, and was eliminated at shallow stall and significantly reduced during deep dynamic stall when the synthetic jets were actuated. Moreover, the activation of synthetic jets resulted in significant reduction in the hysteresis (area within the pitching up and pitching down load history) of the lift and pitching moment through all experimental conditions, as much as 41% and 85%, respectively. SPIV flow fields in shallow dynamic stall demonstrated that actuation of synthetic jets confined the separated region to the trailing edge, in both the instantaneous and time averaged sense. To further reduce the lift and pitching moment hysteresis at high angles of attack, a pulse modulation technique was used and showed a marked increase in synthetic jet performance compared with the continuously actuated case and achieved this result with approximately 65% less power consumption.     

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.     

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.     

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.     

Numerical investigation of the unsteady tip leakage flow and rotating stall inception in a transonic compressor

<正>It is well known that tip leakage flow has a strong effect on the compressor performance and stability. This paper reports on a numerical investigation of detailed flow structures in an isolated transonic compressor rotor-NASA Rotor 37 at near stall and stalled conditions aimed at improving understanding of changes in 3D tip leakage flow structures with rotating stall inception.Steady and unsteady 3D Navier-Stokes analyses were conducted to investigate flow structures in the same rotor.For steady analysis,the predicted results agree well with the experimental data for the estimation of compressor rotor global performance.For unsteady flow analysis, the unsteady flow nature caused by the breakdown of the tip leakage vortex in blade tip region in the transonic compressor rotor at near stall condition has been captured with a single blade passage.On the other hand, the time-accurate unsteady computations of multi-blade passage at near stall condition indicate that the unsteady breakdown of the tip leakage vortex triggered the short length-scale-spike type rotating stall inception at blade tip region.It was the forward spillage of the tip leakage flow at blade leading edge resulting in the spike stall inception. As the mass flow ratio is decreased,the rotating stall cell was further developed in the blade passage.     

翼型前缘对风力机翼型气动性能的影响

基于翼型参数化方法对翼型S809进行两类不同的前缘修改,采用翼型设计分析软件Xfoil对修改前、后的翼型进行气动性能计算分析,并采用计算流体力学（CFD）数值模拟方法进行流场特性分析。结果表明：翼型前缘下弯使得翼型在失速区升力系数增大,阻力系数减小,俯仰力矩系数减小,转捩现象延迟,翼型前缘上弯对气动性能的影响与之相反;翼型前缘上弯和下弯使得翼型表面压力系数分布均匀,吸力面及压力面压力系数增大;翼型前缘下弯能够抑制流动分离,抑制涡的形成,延迟翼型失速,翼型前缘上弯对翼型流场特性的影响则与之相反。     

非线性双方程湍流模型用于翼型分离流动模拟

采用非线性k-ε和q-ω双方程湍流模型,利用高收敛率、高精度和高分辨率的数值算法,对二维NRELS809翼型在0～18°攻角范围内的低速绕流进行了数值模拟,将翼型的升阻特性的计算值与实验值进行了比较。与采用线性k-ε和q-ω双方程湍流模型的计算结果对比表明,当翼型流动发生分离时,采用非线性双方程湍流模型可以明显地改进计算精度。本文采用的数理模型与数值方法可有效地用于风力机翼型大尺度分离流场的预测。