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

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

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

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

Gurney襟翼对水平轴风力机性能影响的实验研究

在小型低速风洞中对装有NACA4424翼型叶片的水平轴风力机及在其尾缘加装Gurney襟翼的风力机进行了一系列性能对比实验。Gurney襟翼的高度分别为2%b和4%b(b为翼型弦长),叶片安装角在6°~14°范围内,实验风速为6~15m/s。实验结果表明,Gurney襟翼对水平轴风力机性能有显著影响,特别是在大安装角(即大攻角和大升力)下;在小安装角(即小攻角和小升力)时,Gurney襟翼使风力机性能降低。同时,装2%b襟翼的风力机性能要高于装4%b襟翼的风力机;在12°安装角时,前者提高风力机功率最少有39%,而后者也可提高风力机功率在34%以上。对于风力机最常用的叶型FFA-W3-211加装2%b的Gurney襟翼后的风洞对比实验同样证明了上述结论。     

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

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

Indicial lift response function: an empirical relation for finite‐thickness airfoils,and effects on aeroelastic simulations

The aeroelastic response of wind turbines is often simulated in the time domain by using indicial response techniques. Unsteady aerodynamics in attached flow are usually based on Jones's approximation of the flat plate indicial response, although the response for finite‐thickness airfoils differs from the flat plate one. The indicial lift response of finite‐thickness airfoils is simulated with a panel code, and an empirical relation is outlined connecting the airfoil indicial response to its geometric characteristics. The effects of different indicial approximations are evaluated on a 2D profile undergoing harmonic pitching motion in the attached flow region; the resulting lift forces are compared with computational fluid dynamics (CFD) simulations. The relevance for aeroelastic simulations of a wind turbine is also evaluated, and the effects are quantified in terms of variations of equivalent fatigue loads, ultimate loads, and stability limits. The agreement with CFD computations of a 2D profile in harmonic motion is improved by the indicial function accounting for the finite‐thickness of the airfoil. Concerning the full wind turbine aeroelastic behavior, the differences between simulations on the basis of Jones's and finite‐thickness indicial response functions are rather small; Jones's flat‐plate approximation results in only slightly larger fatigue and ultimate loads, and lower stability limits. Copyright © 2012 John Wiley & Sons, Ltd.     

Fluid‐structure interaction of a morphing symmetrical wind turbine blade subjected to variable load

A brief summary of the main challenges of rotor design in wind energy conversion (WEC) systems, most notably the horizontal axis wind turbine (HAWT), are presented. One of the limiting factors in HAWT design is choosing the rated capacity to maximize power output and turbine longevity. One such strategy to accomplish this goal is to widen the operational range of the WEC system by using pitch or torque control, which can be costly and subject to mechanical failure. We present a morphing airfoil concept, which passively controls airfoil pitch through elastic deformation. As a justification of the concept, a two‐dimensional fluid‐structure interaction routine is used to simulate the aeroelastic response of a symmetric NACA 0012 blade subjected to variable loading. The results suggest that the morphing blade can be designed to offer superior average lift to drag ratios over a specified range of attack angles by up to 4.2%, and possibly even higher. This infers that the morphing blade design can increase the power production of WEC systems while conceivably reducing cost because the passive deformation of the morphing turbine does not require active control systems that come at an added upfront and maintenance cost. Copyright © 2011 John Wiley & Sons, Ltd.     

The Adaptive Camber Concept—A passive approach for gust load alleviation on wind turbines

A new passive sectional approach to alleviate gust loads on wind turbines—the Adaptive Camber Concept— is introduced. The concept entails fluid‐structure interaction, where flow conditions at the leading edge affect the airfoil shape and vice versa. A two‐dimensional airfoil equipped with the Adaptive Camber Concept is experimentally investigated under steady and unsteady conditions in a wind tunnel. Under steady conditions, the adaptive camber airfoil de‐cambers gradually with increasing angle of attack, yielding a lift curve with declined slope. Unsteady angle of attack fluctuations of various reduced frequencies are generated by means of an active grid. Under unsteady conditions, the adaptive camber airfoil is found to alleviate up to 60% of the fluctuating loads, while generating higher mean lift compared to a rigid reference airfoil.     

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

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

振动膜片的流动控制机理研究

采用计算流体力学方法,研究了主流风速为10 m/s,翼型弦长雷诺数为1.2×10~5条件下振动膜片对NACA0012翼型在18°攻角深失速下流动分离的影响。研究表明:振动膜片能明显提高翼型升力系数、降低阻力系数、改善流场状况;当无量纲频率处在1~1.5范围内时,翼型升阻比可大幅提升,最大可提高75.7%;无量纲振幅对翼型升阻比的影响也很显著,相对于原型存在一个最佳的振幅使得翼型升阻比能获得最大提升;不同振幅下,最佳升阻比对应的无量纲频率随振幅增大而减小。     

基于CFD技术与遗传算法的风力机大厚度翼型优化设计

针对目前风力机大厚度翼型设计参数空间有限、优化设计过程中气动力预测不准等问题,利用B样条函数表征通用翼型廓线,编制程序集成耦合翼型设计模块、任意翼型自适应网格模块、CFD流场计算模块、遗传算法优化模块,提出了基于CFD技术与遗传算法的风力机叶片大厚度翼型优化设计方法,并对比分析优化新翼型与DU97-W-300翼型的几何特性与气动性能。结果表明,优化方法设计的新翼型在主要攻角范围内具有较高的气动性能,在雷诺数为3.0×106的情况下,其升力系数、升阻比分别提高了13.555%、38.588%。该翼型优化设计方法为风力机大厚度通用翼型的设计与应用提供参考。     

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.     

Comparison of pneumatic jets and tabs for Active Aerodynamic Load Control

A fast, efficient way to control loads on utility scale wind turbines is important for the growth of the wind industry. Microtabs and microjets are two Active Aerodynamic Load Control devices, which address this need. Both act perpendicular to the surface of the airfoil, and these actively controlled devices are used to mitigate changes in aerodynamic loading experienced by wind turbine rotors due to wind gusts, wind shear, or other atmospheric phenomena. This work explores the aerodynamic effects of microjets and then compares them to those of microtabs. Flow around an airfoil with an activated microjet at the trailing edge has been simulated using the Reynolds‐averaged Navier–Stokes solver OVERFLOW‐2. Using a Chimera overset grid topology, a microjet has been placed near the trailing edge of the lower surface of a NACA 0012 airfoil. For a jet velocity about half of the freestream velocity, the microjet can change the lift up to ΔC_L = 0.2, but the amount of change varies with the momentum coefficient of the jet. The change in lift is not symmetric for positive and negative angles of attack due to changes in the boundary layer thickness with angle of attack. Increasing the Reynolds number reduces the effectiveness of the microjet only slightly. The effects of jet velocity, jet activation time, and airfoil angle of attack on airfoil lift, drag, and pitching moment are compared with previous work, which illustrates the deployment of a microtab at the 95% chord location of a NACA 0012 airfoil. This study shows that microjets and microtabs have very similar responses in lift and pitching moment, but the drag for the microjet is noticeably lower. Copyright © 2013 John Wiley & Sons, Ltd.     

Design of a new urban wind turbine airfoil using a pressure-load inverse method

This paper presents the design methodology of a new wind turbine airfoil that achieves high performance in urban environment by increasing the maximum lift. For this purpose, an inverse method was applied to obtain a new wind turbine blade section with constant pressure-load along the chord, at the design inlet angle. In comparison with conventional blade section designs, the new airfoil has increased maximum lift, reduced leading edge suction peak and controlled soft-stall behaviour, due to a reduction of the adverse pressure gradient on the suction side. Wind tunnel experimental results confirmed the computational results.     

Lift correction model for local shear flow effect on wind turbine airfoils

Wind turbines operate under various wind conditions in which turbulence virtually always exists. Therefore, unsteady wind turbine simulation methods to estimate wind loading in turbulent inflow conditions are very important for developing optimally designed wind turbines. Several methods have been developed for this purpose and are usually based on the blade element momentum theory (BEMT), which is used for calculation of the wind loading on turbine blades. The local shear flow effect induced by turbulence, however, is not explicitly considered in the popular BEMT-based simulations. Extreme situations can occur in a large-scale wind farm where the inflow field of a wind turbine may contain strong tip vortices generated from upstream turbines. In this study, the effects of idealized local shear flows around a two-dimensional airfoil, S809, on its aerodynamic characteristics were analyzed by CFD simulations. Various parameters including reference inflow velocity, shear rate, angle of attack, and cord length of the airfoil were examined. From the simulation results, several important characteristics were found. The shear rate in a flow causes some changes in the lift coefficient depending on its sign and magnitude, while the angle of attack does not have a distinguishable influence. The chord length and reference inflow also cause proportional and inversely proportional changes in the lift coefficient, respectively. Based on these observations, we adopted an analytic expression for the lift coefficient from the thin airfoil theory and proposed a lift correction model, which is easily applicable to the traditional load analysis procedure based on the BEMT.     

Numerical optimization of axial turbine with self‐pitch‐controlled blades used for wave energy conversion

Wells turbines are among the most practical wave energy converters despite their low aerodynamic efficiency and power produced. It is proposed to improve the performance of Wells turbines by optimizing the blade pitch angle. Optimization is implemented using a fully automated optimization algorithm. Two different airfoil geometries are numerically investigated: the standard NACA 0021 and an airfoil with an optimized profile. Numerical results show that each airfoil has its own optimum blade pitch angle. The present computational fluid dynamics optimization results show that the optimum blade pitch angle for NACA 0021 is +0.3° while that of the airfoil with an optimized profile equals +0.6°.The performance of the investigated airfoils is substantially improved by setting the blades at the optimum blade pitch angle. Both the turbine efficiency and tangential force coefficient are improved, especially at low flow rate and during turbine startup. Up to 4.3% average increase in turbine efficiency is achieved by optimizing the blade pitch angle. A slight improvement of the tangential force coefficient and decrease of the axial force coefficient are also obtained. A tangible increase of the stall‐free operating range is also achieved by optimizing the blade pitch angle. Copyright © 2013 John Wiley & Sons, Ltd.     

A Beddoes‐Leishman–type model with an optimization‐based methodology and airfoil shape parameters

Floating offshore wind turbines operate in a highly unsteady environment; thus, many flow transients occur at the blade cross‐sectional level, which affect the rotor aerodynamics. In every rotor aerodynamics modelling technique requiring the blade element theory, the blade cross‐sectional aerodynamics need to be predicted accurately on the basis of the flow conditions. At reduced frequencies of 0.01 and greater, the flow unsteadiness can be considered significant and cannot be treated as quasisteady. Floating offshore wind turbines can be expected to consistently operate in some degree of yaw or pitch, which may result in reduced frequencies greater than 0.01 over most of the blade when operating at rated wind speeds and rotor RPM. The Beddoes‐Leishman model is a comprehensive but complex model for predicting unsteady airfoil aerodynamics, containing 8 dimensionless time constants. In the present study, the Beddoes‐Leishman model was compared with experimental results of 10 different airfoil profiles, each performed under a range of Reynolds numbers, motion frequencies, mean, and amplitudes of angle of attack. An optimization was performed for all time constants in the model, the results of which were used to formulate a simplified model with fewer equations, without any reduction in accuracy. Further, optimizations were performed against the experimental results of each airfoil, and the optimized constants were compared with shape parameters of the airfoils, yielding possible correlations, which were then applied in the simplified Beddoes‐Leishman model to yield improved accuracy, measured as a 5% reduction in accumulated error between experimental and predicted coefficients of lift.     

前缘改形对风力机翼型气动性能影响的数值研究

基于翼型参数化方法对翼型S809进行4类不同的前缘修改,分别为前缘压力面加厚、前缘吸力面加厚、前缘上弯和前缘下弯,采用翼型设计分析软件Xfoil和商用CFD(Computational Fluid Dynamics)软件FLUENT分别对翼型气动参数和翼型周围流场进行计算。结果表明:翼型气动特性与流场特性受翼型压力面外形变化影响较小;在研究范围内,翼型吸力面加厚使得翼型在失速区升力系数增加,阻力系数减小;翼型前缘上弯使得翼型在大攻角工况下升力系数减小,阻力系数增大,且使翼型提前失速;在一定范围内翼型前缘下弯,使得翼型升力系数增大,阻力系数减小,且延迟失速。     

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.     

尾缘襟翼对风力机翼型气动特性影响研究

尾缘襟翼(TEF)因其对翼型气动特性的调控能力,被认为是降低叶片疲劳和局部载荷最具可行性的气动控制部件。对TEF进行建模,采用Xfoil和CFD软件分析了TEF对翼型气动特性的影响及其机理,并从叶素理论角度对变化来流下TEF的减载效果进行了验证,结果表明:TEF位于不同摆角时翼型升阻力系数均有不同程度的变化,TEF可有效实现对翼型气动特性的主动控制;TEF摆动改变了翼型表面的静压分布和流动状态,进而对翼型升阻力和失速攻角产生影响;TEF可快速有效降低风速突然增加后的叶素受力,进而控制并减小叶片载荷。     

USST高性能水平轴风力机翼型族

结合层流翼型与钝尾缘的特性,通过Hicks-Henne型函数对翼型参数化修型,基于多岛遗传算法及Xfoil气动分析,针对大型水平轴风力机翼型进行多目标函数、多设计工况、多约束条件下的优化设计,得到适用于大型风力机的高性能翼型族（USST翼型族）。其升阻比在大多数攻角下均高于同厚度的FFA、DU系列等现有风力机翼型族,且在同样的升力系数下具有更大的升阻比。最后为考核优化设计得到的翼型族,采用数值模拟方法对优化结果进行验证,证明设计得到的新型风力机翼型族具有优越的气动性能。