复杂工况下的大型风力机气动性能和尾迹研究Investigation of aerodynamic performance and wake of the large scale wind turbine in complicated conditions

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

复杂工况下的大型风力机气动性能和尾迹研究

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

基于涡尾迹方法的风力机气动性能分析

基于(势流)涡尾迹方法开发了水平轴风力机叶片气动性能分析程序,采用固定尾迹涡模型和自由尾迹模型分别对气动设计性能进行计算分析,得到风力机设计工况下的涡位置、诱导因子、功率系数及扭矩系数等气动性能参数,并与设计结果对比。结果表明,涡尾迹方法能够快速准确地计算风力机叶片气动性能参数,对风力机叶片气动分析,固定尾迹涡模型较自由尾迹模型计算时间短,具有较好的实用性。     

风力机气动特性的浸入边界法模拟

风力发电机的空气动力学性能是决定风力机安全与效率的最重要因素之一。但由于影响风力机气动性能的参数众多,更加高效精确地模拟风力机气动特性一直是风力机的重要发展方向。本文提出了基于浸入边界法的风力机建模,网格离散,以及数值模拟的统一性框架。利用同伦变形来生成光滑的叶片模型,并且使用仿射变换来处理叶片的渐缩与扭转问题。首先,针对二维翼型的升阻力,检验了算法的数值精度。表明此方法对于阻力的模拟具有非常严格的一阶精度,进而提出采用理查森外推法来精确高效修正升阻力模拟结果。同时,模拟研究了拱曲度以及厚度对二维翼型升阻力的影响。随后,模拟研究了单风力机(包含塔架)在不同尖速比下的功率系数,并对塔架与叶片间的相互气动作用进行了初步分析。最后,模拟研究了双风力机在风场中不同前后间隔距离下的气动干涉问题。本文主要意义在于验证建模,离散,与数值模拟的一体化框架的有效可行性,进而为后续研究(给定约束下风力机自动优化选型)提供坚实基础。     

基于Schmitz理论的风叶气动设计研究

应用Schmitz理论进行叶片气动设计,考虑了风力机叶片的气动损失,用Schmitz理论推导出风力机叶片的基本设计参数的计算公式,并考虑了风力机在启动和空载时风力机的实际工作点偏离了设计点,对叶片的气动性能参数进行了修正。通过对200kW风力机的算例表明:随着叶片半径的增大,入流角逐渐减小;叶片弦长先增大后减小,修正后得到的风力机在非设计点处的推力、驱动力矩、功率与实际风力机的特性规律相符。     

二维扩压叶栅非定常分离流控制途径探索

二维扩压叶栅非定常黏性数值模拟结果表明,在一定攻角范围内,叶片前缘点附近的周期性吹吸气激励能有效控制混乱的非定常分离流．详细研究了非定常激励频率、幅值、位置对流场的影响．满足一定条件的非定常激励能够使流动由无序变为有序,时均气动性能提高。     

基于CFD技术的垂直轴风力机动态尾流特性研究

利用CFD软件对麦克马斯特大学垂直轴风力机进行不同叶尖速比下的数值模拟,计算结果与风洞试验数据吻合良好。近场尾流中,与单叶片的风力机模拟结果比较,上游叶片产生并向下游延伸的旋涡影响下游运行轨道上叶片的升阻力特性,不仅使叶片扭矩输出峰值降低,而且峰值产生的时间延迟。对垂直轴风力机叶片叶梢进行修改,模拟结果显示,叶片扭矩输出峰值不变,但是谷值有所降低,修改后风力机沿风向推力幅值降低明显;远场尾流中,采用风速轮廓线原理,以瑞典的法尔肯贝里市200kW垂直轴风力机为原型,按照真实的空间排布进行数值模拟。模拟结果显示,上游风力机上下两端处产生较为集中的远场尾流,影响下游风力机叶片下半段的气动性能,下游风力机功率输出降低明显。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε（RNG）湍流模型,数值模拟美国可再生能源实验室（NREL）5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε(RNG)湍流模型,数值模拟美国可再生能源实验室(NREL)5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

风剪切对风力机叶片气动性能及尾迹形状的影响

风力机通常运行在非定常工况中,其气动性能及尾迹会随着工况的变化而变化. 风剪切是风力机长期所处的环境,它会影响到叶片气动载荷、尾迹形状、总体性能等,分析风剪切作用下的叶片气动性能对风力机的设计有重要意义.本文采用一种时间步进自由涡尾迹(free vortex wake, FVW)方法,耦合FVW方法与风剪切模型,计算不同风剪切因子作用下叶片的气动力系数、推力以及风轮后的尾迹形状变化,研究尾迹形状变化对风轮旋转平面诱导速度及风力机叶片气动性能的影响. 结果表明:在风剪切入流条件下, 随着风剪切因子的增大,风力机的气动力系数随时间做周期性波动的幅度加剧, 推力的平均值逐渐减小,尾迹倾斜程度增大, 尾迹在轮毂下方的倾斜程度更明显;尾迹形状的变化使风轮平面轴向诱导速度因子分布不均匀,同时使风力机的总体性能降低且偏离较大;倾斜尾迹相比于对称尾迹对风轮平面处的诱导影响有明显差别, 波动幅值增大,气动力系数在波谷处的偏差比波峰处大. 尾迹越倾斜,风轮旋转平面处的载荷不对称性越明显.      

旋转流场下模型大幅俯仰运动实验研究

在中航气动院FL-8低速风洞中,采用单自由度振荡机构进行了旋转流场下大幅俯仰运动的气动特性实验研究。模型在绕风速轴连续旋转的同时,进行给定频率和振幅绕体轴的俯仰振荡运动。测量了模型的动态气动特性,着重分析了不同运动参数对模型气动特性的影响。实验结果表明,旋转速度的存在使大振幅俯仰振荡实验中的滚转力矩和偏航力矩产生了明显的迟滞特性,但对俯仰力矩和法向力的迟滞特性影响不大。     

Investigation of horizontal-axis wind turbine (HAWT) blade threedimensional rotational effect based on field experiments

Field experiments are performed on a two-bladed 33 kW horizontal-axis wind turbine (HAWT). The pressures are measured with 191 pressure sensors positioned around the surfaces of seven spanwise section airfoils on one of the two blades. Three-dimensional (3D) and two-dimensional (2D) numerical simulations are performed, respectively, on the rotor and the seven airfoils of the blade. The results are compared with the experimental results of the pressure distribution on the seven airfoils and the lift coefficients. The 3D rotational effect on the blade aerodynamic characteristics is then studied with a numerical approach. Finally, some conclusions are drawn as follows. From the tip to the root of the blade, the experimental differential pressure of the blade section airfoil increases at first and then decreases gradually. The calculated 3D result of the pressure distribution on the blade surface is closer to that of the experiment than the 2D result. The 3D rotational effect has a significant impact on the blade surface flow and the aerodynamic load, leading to an increase of the differential pressure on the airfoils and their lift coefficient than that with the 2D one because of the stall delay. The influence of the 3D rotational effect on the wind turbine blade especially takes place on the sections with flow separation.     

风力机叶片翼型动态试验技术研究

风力机叶片动态振荡过程往往伴随着俯仰和横摆同时进行, 以前对许多动态问题不清楚的阶段, 工程上不惜以增加叶片重量为代价而采用偏安全的设计, 通常忽略横摆振荡的影响; 大型风力机设计对获取翼型更加全面、准确的动态载荷提出了更高要求, 研究横摆振荡对翼型动态气动特性的影响规律具有重要意义. 本文首次开展翼型横摆振荡动态风洞试验研究, 采用“电子凸轮”技术代替机械凸轮实现了振荡频率和振荡角度的无级变化, 基于设计的电子外触发装置实现了对动态流场的实时测量, 实现了风洞来流、模型角位移和动态压力数据的同步采集, 分别开展了翼型静态测压、俯仰/横摆动态测压、粒子图像测速和荧光丝线等试验研究, 试验结果准度较高、规律合理; 分析了动态试验洞壁干扰影响机制. 研究表明, 横摆振荡翼型的气动曲线也存在明显迟滞效应; 随着振荡频率升高, 翼型俯仰和横摆振荡下的气动迟滞性均增强; 翼型俯仰振荡正行程的动态失速涡破裂有所延迟; 洞壁与模型端部交界处的强三维效应对翼型压力分布影响较大; 建立的横摆振荡试验技术可为风力机动态掠效应的研究提供技术支撑.      

Dynamic stall model for wind turbine airfoils

A model is presented for aerodynamic lift of wind turbine profiles under dynamic stall. The model combines memory delay effects under attached flow with reduced lift due to flow separation under dynamic stall conditions. The model is based on a backbone curve in the form of the static lift as a function of the angle of attack. The static lift is described by two parameters, the lift at fully attached flow and the degree of attachment. A relationship between these parameters and the static lift is available from a thin plate approximation. Assuming the parameters to be known during static conditions, nonstationary effects are included by three mechanisms: a delay of the lift coefficient of fully attached flow via a second-order filter, a delay of the development of separation represented via a first-order filter, and a lift contribution due to leading edge separation also represented via a first-order filter. The latter is likely to occur during active pitch control of vibrations. It is shown that all included effects can be important when considering wind turbine blades. The proposed model is validated against test data from two load cases, one at fully attached flow conditions and one during dynamic stall conditions. The proposed model is compared with five other dynamic stall models including, among others, the Beddoes–Leishman model and the ONERA model. It is demonstrated that the proposed model performs equally well or even better than more complicated models and that the included nonstationary effects are essential for obtaining satisfactory results. Finally, the influence of camber and thickness distribution on the backbone curve are analysed. It is shown that both of these effects are adequately accounted for via the static input data.     

Dynamic study of a wind turbine blade with horizontal axis

The study of the dynamic behavior of a wind turbine with horizontal axis can be undertaken by various methods of analysis. The effects of the change of the aerodynamic flow (in the steady and unsteady cases), the variation of parameters of the cinematic movement (angle of attack, pitch angle and yaw angle) and the definition of subsystems characteristics that makes the wind turbine (blade, nacelle and pylon) allow one to characterize the structural dynamic behavior of the wind turbine. It is therefore necessary to develop these items. Once this is done, the structural dynamic behavior of the system can be improved. The term `improve' means the increase of the life duration by mastering the fatigue effects and the reduction of cost without sacrificing the aerodynamic output. The present study aims to examine the behavior of the blade, which is the main part of the wind turbine in that it that transmits forces to all other parts of the structure. The model is based on the theory of three-dimensional beams, under the assumption of variable sections of the type NACA 4415 airfoil, and takes into account membrane, transversal shear, flexion and free torsion effects. With regards to the aerodynamic loads (the lift, the drag and the pitching moment), a validation has been undertaken by considering experimental data and numerical results obtained by a CFD code (Fluent). The forces are obtained by means of a parametric CAD method interpolation of the aerodynamic poles by Bézier patch under geometrical constraints solved by a Simplex type algorithm. The emphasis is put on dynamic aspects by a complete processing of the dynamic equilibrium equation, applied to the wind turbine blade with horizontal axis.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

A brief review on wind turbine aerodynamics

This article briefly reviews wind turbine aerodynamics, which follows an explanation of the aerodynamic complexity. The aerodynamic models including blade momentum theory, vortex wake model, dynamic stall and rotational effect, and their applications in wind turbine aerodynamic performance prediction are discussed and documented. Recent progress in computational fluid dynamics for wind turbine is addressed. Wind turbine aerodynamic experimental studies are also selectively introduced.     

NUMERICAL STUDY ON DYNAMIC STALL FLOW CONTROL FOR WIND TURBINE AIRFOIL USING PLASMA ACTUATOR 1)

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.     

Lift enhancement of airfoil and tip flow control for wind turbine

Two techniques that improve the aerodynamic performance of wind turbine airfoils are described. The airfoil S809, designed specially for wind turbine blades, and the airfoil FX60-100, having a higher lift-drag ratio, are selected to verify the flow control techniques. The flow deflector, fixed at the leading edge, is employed to control the boundary layer separation on the airfoil at a high angle of attack. The multi-island genetic algorithm is used to optimize the parameters of the flow deflector. The results indicate that the flow deflector can suppress the flow separation, delay the stall, and enhance the lift. The characteristics of the blade tip vortex, the wake vortex, and the surface pressure distributions of the blades are analyzed. The vortex diffuser, set up at the blade tip, is employed to control the blade tip vortex. The results show that the vortex diffuser can increase the total pressure coefficient of the core of the vortex, decrease the strength of the blade tip vortex, lower the noise, and improve the efficiency of the blade.