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.     

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.     

带脊状结构风力机翼型的噪声特性研究

采用大涡模拟和声类比的方法,以NACA0018翼型为研究对象,研究脊状结构对翼型远场噪声的影响。分别模拟来流速度为12和24 m/s,在6°攻角下布置脊状结构的翼型流场,对应的基于弦长雷诺数约为8.0×104和1.6×10⁵。通过FW-H方程计算大涡模拟提取的声源项,得到Riblet-Q和Riblet-H翼型的流场和声场。非定常流场计算结果表明:6°攻角下Riblet-H翼型能够改善翼型边界层分离情况,抑制涡结构脱落,从而减小翼型表面压力脉动和接收点处声压波动,同时能提高翼型的气动性能。逆压梯度段脊状结构可有效减小频率在0～3000 Hz内的噪声。进一步研究表明,该状态下的噪声主要由边界层引起的涡脱落噪声主导。     

An experimental study on aerodynamic sound generated from wake interaction of circular cylinder and airfoil with attack angle in tandem

In our previous study, the effects of the interval between the cylinder and the airfoil on the aerodynamic sound were investigated and compared with the cases of single circular and single airfoil. In this study, the effects of the attack angle of the airfoil located downstream on the characteristics of aerodynamic sound and the wake structure are investigated at a given interval between the cylinder and the airfoil. It is found that the sound pressure level of DFN and the peak frequency decrease with increasing attack angle of airfoil because of the diffusive wake structure due to the increased back pressure of circular cylinder, which is caused by the blocking effect of airfoil. It is shown that the sound sources are corresponded to the attack points of shedding vortex form the upstream circular cylinder to the downstream airfoil. We conclude that the pressure fluctuation at the airfoil surface effects on the sound pressure level, from the flow visualizations and the exploration test of sound source.     

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.     

新型变桨立轴风力机涡效应分析

该文旨在通过变桨来改善升力型立轴风力机叶片气动特性,提高风力机最大运行效率。针对设计尖速比下风能利用系数较低的问题,提出减小叶片小攻角范围,增大叶片大攻角工作范围,以重点改善叶片低性能区域的气动特性为出发点,提高风能利用系数新变桨思路。以采用NACA0012翼型、2 m高和2 m旋转直径的两叶片H型风力机为研究对象,从涡理论来分析和比较在最佳尖速比为5的条件下,附着涡、尾随涡、脱体涡和桨距角对攻角、切向力和功率输出的影响规律。研究结果表明:变桨后,叶片的攻角、切向力和输出功率在原最大值两侧均有明显提高,拓宽了叶片高性能的工作区域;涡系中脱体涡对叶片气动特性影响最大,其中在上盘面影响较小,在下盘面影响较大;变桨前后涡系对上盘面的差异较小,对下盘面的影响差异较明显;变桨后,下盘面的叶片的涡尾迹弯曲程度在加大。     

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

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

Frequency Analyses for Mechanism Research of Unsteady Flow Inside Three Dimensional Compressor Cascade

Numerical simulations are carried out for unsteady flow field of certain kind of 3D compressor cascade. Emphasis is laid on vortex shedding and frequency analysis in a compressor cascade. Numerical simulations using unsteady Reynolds-averaged viscous turbulent equations are carried out. The results show that the flows in separated areas and wake areas are characterized as periodic or quasi-periodic vortex shedding and the frequencies of vortex shedding vary with incidences and Mach number. At the same Mach number, the frequency of vortex shedding will decrease as the incidence increases. Yet, the frequency will increase as the Mach number increases at the same incidence. In the same computation case, the frequency of vortex shedding will vary along the span of blade. The frequency is smaller at the middle of the span than that at the hub.     

A simple model for deep dynamic stall conditions

The present study is focused on modeling of dynamic stall behavior of a pitching airfoil. The deep stall regime is in particular considered. A model is proposed, which has a low implementation and computational complexity but yet is able to deal with different types of dynamic stall conditions, including those characterized by multiple vortex shedding at the airfoil leading edge. The proposed model is appraised against an extensive data set of experimental (α,C_L) curves for NACA0012. The results of an existing widely used model, having comparable complexity, are also shown for comparison. The proposed model is able to well reproduce not only the classic curves of deep dynamic stall but also the curves characterized by lift oscillations at high angles of attack due to the shedding of multiple vortices. Furthermore, the model appears to be robust to variations of its parameters from the optimal values and of the airfoil geometry. Finally, the model is successfully implemented in a commercial CFD software and applied to the simulation of a vertical axis wind turbine within the actuator cylinder approach. The accuracy of the prediction of the turbine power coefficient in the whole rotation cycle is very good for the optimal working condition of the turbine, for which the model parameters were calibrated. Fairly good accuracy is also obtained in significantly different working conditions without any further calibration.     

翼型表面随行波结构流动控制的数值研究

针对翼型表面的流动分离,采用数值模拟方法研究随行波结构流动控制的机制。为了验证计算方法的可靠性,将翼型表面静压曲线与实验测试结果进行对比,发现两者吻合程度较好。数值计算结果表明,适当增加随行波的相对弦长长度,有助于改善翼型的气动性能。当攻角小于3°时,随行波位于分离点之前,沟槽内形成顺时针旋转的二次涡,有助于加速边界层低速条带;在大攻角下,随行波位于分离点之后,顺时针卷起的分离涡在沟槽内形成逆时针旋转的二次涡,与分离涡互为反向涡对,减小了尾缘分离区范围。翼型表面随行波能有效地控制边界层流动。     

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

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

Vortex‐induced vibrations on a modern wind turbine blade

This article investigates the aero‐elastic response of the DTU 10‐MW RWT blade in deep stall conditions with angles of attack in the vicinity of 90 degrees. The simulations were conducted with the high‐fidelity fluid–structure interaction simulation tool HAWC2CFD employing the multi‐body‐based structural model of HAWC2 and the incompressible computational fluid dynamics solver EllipSys3D. The study utilizes detached eddy simulation computations and considers the three‐dimensional blade geometry including blade twist and taper. A preliminary frequency analysis of the load variations on a stiff blade showed that an inclined inflow with a velocity component along the blade axis can trigger a spanwise correlated vortex shedding over large parts of the blade. Moderate wind speeds were sufficient to generate vortex shedding with frequencies close to the first edgewise eigenfrequency of the blade. Aero‐elastic computations of the elastic blade confirmed the findings of the frequency analysis. Inflow conditions with inclination angles between Ψ = 20° and Ψ = 55° and relatively low to moderate wind speeds between V = 16 and V = 26ms^?1 were sufficient to trigger severe edgewise blade vibrations with blade tip amplitudes of several metres. The investigated inflow conditions are considered realistic and might occur when the wind turbine is idling or standing still and the yaw system is unable to align the wind turbine with the incoming wind. Copyright © 2016 John Wiley & Sons, Ltd.     

不同攻角下压气机叶栅涡流噪声辐射特性的研究

采用大涡模拟模型和边界元方法对压气机叶栅的非定常流场与声场进行了分析,研究了来流攻角对叶栅尾缘非定常涡脱落及其噪声辐射特性的影响.结果表明：压气机叶栅流场中存在着尾缘涡脱落现象,尾缘涡脱落的强度随着攻角的变化而变化.在较大的正攻角下,吸力面存在的分离涡与下游的尾涡相互影响,引起较大的升力脉动,使得噪声辐射增强;而在较大负攻角下,压力面的分离涡同样影响尾缘涡的脱落,噪声辐射水平也有增强趋势.来流攻角对噪声辐射的指向性影响不大.声压级的最小值出现在阻力系数最小的攻角下.     

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

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

尾随涡对风力机叶片气动性能影响机理研究

以Phase Ⅵ风力机叶片为研究对象,以r/R=30%、63%和95%处叶素为参考,建立与7、9、15 m/s试验风速下该风力机叶片附着涡环量沿展向分布相同的叶片模型,分析尾随涡对风力机当地翼型气动性能的影响机理。采用带转捩效应的SST k-ω湍流模型,对所建立的叶片模型和二维S809翼型的气动特性进行研究和对比分析。结果表明：旋转叶片尾随涡对分离现象产生抑制作用且随攻角的增大减弱;尾随涡的影响表现出多重效应,除了减小当地翼型的攻角,还降低其吸力面负压系数和压力面正压系数。     

Extraction of airfoil data using PIV and pressure measurements

A newly developed technique for determining the angle of attack (AOA) on a rotating blade is used to extract AOAs and airfoil data from measurements obtained during the MEXICO (Model rotor EXperiments in COntrolled conditions) rotor experiment. Detailed surface pressure and Particle Image Velocimetry (PIV) flow fields at different rotor azimuth positions are examined for determining sectional airfoil data. The AOA is derived locally by determining the local circulation on the blade from pressure data and subtracting the induction of the bound circulation from the local velocity. The derived airfoil data are compared to 2D data from wind tunnel experiments and XFOIL computations. The comparison suggests that the rotor is subject to severe 3D effects originating from the geometry of the rotor, and explains why the Blade Element Momentum technique with 2D airfoil data over‐predicts the loading of the rotor. The extraction technique is verified by employing the derived airfoil characteristics as input to computations using the BEM technique and comparing the calculated axial and tangential forces to the measured data. The comparison also demonstrates that the used technique of determining the AOA is a reliable tool to extract airfoil data from experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

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

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

Prediction of the wind turbine performance by using BEM with airfoil data extracted from CFD

Blade element momentum (BEM) theory with airfoil data is a widely used technique for prediction of wind turbine aerodynamic performance, but the reliability of the airfoil data is an important factor for the prediction accuracy of aerodynamic loads and power. The airfoil characteristics used in BEM codes are mostly based on 2D wind tunnel measurements of airfoils with constant span. Due to 3D effects, a BEM code using airfoil data obtained directly from 2D wind tunnel measurements will not yield the correct loading and power. As a consequence, 2D airfoil characteristics have to be corrected before they can be used in a BEM code. In this article, we consider the MEXICO (Model EXperiments In Controlled cOnditions) rotor where airfoil data are extracted from CFD (Computational Fluid Dynamics) results. The azimuthally averaged velocity is used as the sectional velocity to define the angle of attack and the coefficient of lift and drag is determined by the forces on the blade. The extracted airfoil data are put into a BEM code without further corrections, and the calculated axial and tangential forces are compared to both computations using BEM with Shen's tip loss correction model and experimental data. The comparisons show that the recalculated forces by using airfoil data extracted from CFD have good agreements with the experiment.     

An integral boundary layer engineering model for vortex generators implemented in XFOIL

To assess and optimize vortex generators (VGs) for flow separation control, the effect of these devices should be modelled in a cost and time efficient way. Therefore, it is of interest to extend integral boundary layer models to analyse the effect of VGs on airfoil performance. In this work, the turbulent boundary layer formulation is modified using a source term approach. An additional term is added to the shear‐lag equation, to account for the increased dissipation due to streamwise vortex action in the boundary layer, forcing transition at the VG leading edge where applicable. The source term is calibrated and a semi‐empirical relation is set up and implemented in XFOIL . The modified code is capable of addressing the effect of the VG height, length, inflow angle, and chordwise position on the airfoil's aerodynamic properties. The predicted polars for airfoils with VGs show a good agreement with reference data, and the code robustness is demonstrated by assessing different airfoil families at a wide range of Reynolds numbers.