Unsteady aerodynamics of offshore floating wind turbines in platform pitching motion using vortex lattice method

As the flow states of an offshore floating wind turbine (OFWT) differ from those of an onshore fixed wind turbine, it is questionable as to whether the aerodynamic load prediction of a turbine using conventional blade element momentum theory (BEMT) is accurate. The aim of this paper is to show the characteristics of aerodynamic load predictions using the vortex lattice method (VLM). Washizu's experimental data, which was measured under a similar flow state of a floating wind turbine, is used for validation. The prediction shows good results compared to those of an experiment. To determine the unsteady aerodynamics of a floating wind turbine, the NREL 5 MW wind turbine model is used for the simulation of a floating wind turbine. These results show that a turbulent wake state (TWS), which is undesirable condition and cannot predicted in BEMT simulation, arises when a floating wind turbine is operated at a low-speed inflow condition. In addition, the rotor experiences a TWS when the floating platform undergoes upward pitching motion.     

1.5MW风力机整机定常及非定常数值模拟

本文运用计算流体力学方法,以1.5 MW风力机为例,对风力机整机三维模型的空气动力学特性开展研究。针对三翼型风力机叶片,利用改进的Wilson方法进行气动设计,并通过寻找各截面最佳雷诺数的方法进行优化修正。建立了整机三维模型,设计流域并划分网格,定义边界及区域。最后对上述模型进行额定工况下定常与非定常数值模拟,利用模拟结果开展有关压力、失速特性等空气动力学特性的分析。结果表明：非定常模拟在风轮背面上的平均压力比定常小,使风轮前后压差变大,输出功率加大,其主要原因是叶尖出力的增加;旋转使得风力机叶片发生流体分离延迟,且产生更高的升力系数。     

CFD simulation of a floating offshore wind turbine system using a variable-speed generator-torque controller

Prediction and control of rotor rotational velocity is critical for accurate aerodynamic loading and generator power predictions. A variable-speed generator-torque controller is combined with the two-phase CFD solver CFDShip-Iowa V4.5. The developed code is utilized in simulations of the 5 MW floating offshore wind turbine (FOWT) conceptualized by the National Renewable Energy Laboratory (NREL) for the Offshore Code Comparison Collaboration (OC3). Fixed platform simulations are first performed to determine baseline rotor velocity and developed torque. A prescribed platform motion simulation is completed to identify effects of platform motion on rotor torque. The OC3’s load case 5.1, with regular wave and steady wind excitation, is performed and results are compared to NREL’s OC3 results. The developed code is shown to functionally control generator speed and torque but requires controller calibration for maximum power extraction. Generator speed variance is observed to be a function of unsteady stream-wise platform motions. The increased mooring forces of the present model are shown to keep the turbine in a more favorable variable-speed control region. Lower overall platform velocity magnitudes and less rotor torque are predicted corresponding to lower rotor rotational velocities and a reduction in generated power. Potential improvements and modifications to the present method are considered.     

预弯对10MW级风力机叶片气动特性的影响

为分析预弯处理对10 MW级风力机叶片气动特性的影响,以DTU 10 MW风力机为例,采用CFD数值模拟方法,研究均匀来流不同风速下风力机的输出功率,并与BEM计算结果进行对比。同时,对比分析直叶片和预弯叶片风力机的功率特性、沿展向出力分布、沿展向不同截面翼型的流动特性。研究结果表明,直叶片各截面翼型的压力差较预弯叶片的大,做功能力较强。预弯通过对叶片的三维流动产生扰动,进而影响风力机的输出功率,且主要体现在叶片展向70%~90%的位置。研究成果可为风力机叶片气动性能的设计与优化提供参考。     

Wind turbine aerodynamic response under atmospheric icing conditions

This article deals with the atmospheric ice accumulation on wind turbine blades and its effect on the aerodynamic performance and structural response. The role of eight atmospheric and system parameters on the ice accretion profiles was estimated using the 2D ice accumulation software lewice Twenty‐four hours of icing, with time varying wind speed and atmospheric icing conditions, was simulated on a rotor. Computational fluid dynamics code, FLUENT, was used to estimate the aerodynamic coefficients of the blade after icing. The results were also validated against wind tunnel measurements performed at LM Wind Power using a NACA64618 airfoil. The effects of changes in geometry and surface roughness are considered in the simulation. A blade element momentum code WT‐Perf is then used to quantify the degradation in performance curves. The dynamic responses of the wind turbine under normal and iced conditions were simulated with the wind turbine aeroelastic code HAWC2. The results show different behaviors below and above rated wind speeds. In below rated wind speed, for a 5 MW virtual NREL wind turbine, power loss up to 35% is observed, and the rated power is shifted from wind speed of 11 to 19 m s^?1. However, the thrust of the iced rotor in below rated wind speed is smaller than the clean rotor up to 14%, but after rated wind speed, it is up to 40% bigger than the clean rotor. Finally, it is briefly indicated how the results of this paper can be used for condition monitoring and ice detection. Copyright © 2012 John Wiley & Sons, Ltd.     

Predicting wind turbine blade loads and aeroelastic response using a coupled CFD–CSD method

The aeroelastic response and the airloads of horizontal-axis wind turbine rotor blades were numerically investigated using a coupled CFD–CSD method. The blade aerodynamic loads were obtained from a Navier–Stokes CFD flow solver based on unstructured meshes. The blade elastic deformation was calculated using a FEM-based CSD solver which employs a nonlinear coupled flap-lag-torsion beam theory. The coupling of the CFD and CSD solvers was accomplished in a loosely coupled manner by exchanging the information between the two solvers at infrequent intervals. At first, the present coupled CFD–CSD method was applied to the NREL 5MW reference wind turbine rotor under steady axial flow conditions, and the mean rotor loads and the static blade deformation were compared with other predicted results. Then, the unsteady blade aerodynamic loads and the dynamic blade response due to rotor shaft tilt and tower interference were investigated, along with the influence of the gravitational force. It was found that due to the aeroelastic blade deformation, the blade aerodynamic loads are significantly reduced, and the unsteady dynamic load behaviors are also changed, particularly by the torsional deformation. From the observation of the tower interference, it was also found that the aerodynamic loads are abruptly reduced as the blades pass by the tower, resulting in oscillatory blade deformation and vibratory loads, particularly in the flapwise direction.     

海上风电结构横风向气动力阻尼特性研究

基于横风向气动力阻尼理论计算模型,以NREL-5 MW海上风电机组为例,对其运行过程中横风向气动力阻尼进行计算,并采用FAST软件对计算结果进行验证。之后,研究转速、叶片桨距角和运行方式对横风向气动力阻尼的影响。研究结果表明：NREL-5 MW海上风电机组结构运行状态下的横风向气动力阻尼在0%~0.8%范围内变化,其随风电机组运行转速及叶片桨距角的增大而增大;此外,海上风电机组不同运行方式对其横风向气动力阻尼也会产生较大影响。     

Shape optimization of wind turbine blades

This paper presents a design tool for optimizing wind turbine blades. The design model is based on an aerodynamic/aero‐elastic code that includes the structural dynamics of the blades and the Blade Element Momentum (BEM) theory. To model the main aero‐elastic behaviour of a real wind turbine, the code employs 11 basic degrees of freedom corresponding to 11 elastic structural equations. In the BEM theory, a refined tip loss correction model is used. The objective of the optimization model is to minimize the cost of energy which is calculated from the annual energy production and the cost of the rotor. The design variables used in the current study are the blade shape parameters, including chord, twist and relative thickness. To validate the implementation of the aerodynamic/aero‐elastic model, the computed aerodynamic results are compared to experimental data for the experimental rotor used in the European Commision‐sponsored project Model Experiments in Controlled Conditions, (MEXICO) and the computed aero‐elastic results are examined against the FLEX code for flow past the Tjæreborg 2 MW rotor. To illustrate the optimization technique, three wind turbine rotors of different sizes (the MEXICO 25 kW experimental rotor, the Tjæreborg 2 MW rotor and the NREL 5 MW virtual rotor) are applied. The results show that the optimization model can reduce the cost of energy of the original rotors, especially for the investigated 2 MW and 5 MW rotors. Copyright © 2009 John Wiley & Sons, Ltd.     

The effect of a leading edge erosion shield on the aerodynamic performance of a wind turbine blade

Armour EDGE is a novel shield developed to protect the leading edge of wind turbine blades from erosion. The aerodynamic impact on aerofoils of National Renewable Energy Laboratory (NREL) 5MW wind turbine has been investigated using 2D fully turbulent computational fluid dynamics (CFD), with three profiles at critical locations along the blade simulated both with and without the shield to compare aerodynamic performance. Two wind speeds were investigated that reflect regular operating conditions: at rated speed of 11.4 m/s and a below rated speed of 7 m/s. The results showed that the presence of the shield during rated wind speed reduced the drag by as much as 4.5%, where the lift‐to‐drag ratio increased by a maximum of 4%. At the below rated wind speeds, the shield had negligible impact on the performance of all but one National Advisory Committee for Aeronautics (NACA) 64‐618 profile, which resulted in an increase in the drag coefficient of 7%. It was also found that the suction side of the aerofoil is much more sensitive to leading edge protection placement than the pressure side. It was concluded that the erosion shield as a method of leading edge protection, with a gradual transition from shield to blade, will not have a major impact on the aerodynamic performance of a multi‐megawatt wind turbine blade and could slightly increase aerofoil efficiency at high wind speeds.     

Validation of a three-dimensional viscous–inviscid interactive solver for wind turbine rotors

MIRAS is a newly developed computational model that predicts the aerodynamic behavior of wind turbine blades and wakes subject to unsteady motions and viscous effects. The model is based on a three-dimensional panel method using a surface distribution of quadrilateral singularities with a Neumann no penetration condition. Viscous effects inside the boundary layer are taken into account through the coupling with the quasi-3D integral boundary layer solver Q³UIC. A free-wake model is employed to simulate the vorticity released by the blades in the wake. In this paper the new code is validated against measurements and/or CFD simulations for five wind turbine rotors, including three experimental model rotors [20–22], the 2.5 MW NM80 machine [23] and the NREL 5 MW virtual rotor [24]. Such a broad set of operational conditions and rotor sizes constitutes a very challenging validation matrix, with Reynolds numbers ranging from 5.0⋅10⁴ to 1.2⋅10⁷.     

垂直对称四叶片风力机气动特性仿真研究

采用CFD计算软件对垂直轴风力机气动性能进行计算.首先,使用ICEM软件对模型进行前处理,通过Fluent软件进行数值模拟,分析不同计算时间步长和湍流模型对风力机气动特性仿真结果的影响,确定符合该研究模型的计算方法.随后,对顺流垂旋型垂直轴风力机在不同叶尖速比下进行计算,发现该风力机在叶尖速比为0.42时获得最大功率系...     

The Darrieus wind turbine: Proposal for a new performance prediction model based on CFD

This paper presents a CFD model for the evaluation of energy performance and aerodynamic forces acting on a straight-bladed vertical-axis Darrieus wind turbine. The basic principles which are currently applied to BE-M theory for rotor performance prediction are transferred to the CFD code, allowing the correlation between flow geometric characteristics (such as blade angles of attack) and dynamic quantities (such as rotor torque and blade tangential and normal forces). The model is proposed as a powerful design and optimization tool for the development of new rotor architectures for which test data is not available.After describing and validating the computational model against experimental data, a full campaign of simulation is proposed for a classical NACA 0021 three-bladed rotor.Flow field characteristics are investigated for several values of tip speed ratio, allowing a comparison among rotor operation at optimum and lower C_p values, so that a better understanding of vertical-axis wind turbines basic physics is obtained.     

大型风力机风轮气动不平衡的特性研究与验证

为了研究和探索风轮气动不平衡的物理特性,以某2.0 MW三叶片水平轴风力机为研究对象,采用计算机仿真及试验相结合的方法,研究风轮气动不平衡对机组动力学特性、气动性能及气动载荷的影响研究。通过气动特性分析和动力学分析表明,随着风轮叶片安装角的不平衡度增大,其机组性能逐渐下降,塔顶的载荷波动逐渐增大,叶片的挥舞载荷出现明显差异,机舱振动加速度变大。对塔顶振动加速度进行快速傅里叶变换分析,出现明显特征变化。研究过程表明,监测机舱振动加速度和机组功率曲线能有效识别机组气动不平衡程度。     

Influence of wake dynamics on the performance and aeroelasticity of wind turbines

Wind turbines are currently a rapidly expanding form of renewable energy. However, there are numerous technological challenges that must be overcome before wind energy provides a significant amount of power in the United States. One of the primary challenges in wind turbine design and analysis is accurately accounting for the aerodynamic environment. This study is focused on a comprehensive verification and validation of the NREL FAST code, which is enhanced to include a free vortex wake model. The verification and validation is carried out through a comparison of blade lift distribution, wind turbine power and force and moment coefficients using a combination of CFD and experimental data. The results are also compared against Blade Element Momentum theory, and results from a 2001 double-blind NREL study on the prediction capabilities of wind turbine modeling tools. Results indicate that the enhanced aeroelastic code generally provides improved predictions. However, in several notable cases the predictions are only marginally improved, or even worse, than those generated using Blade Element Momentum theory aerodynamics. It is concluded that modeling of the aerodynamic environment remains incomplete, even after inclusion of wake effects. One important aspect identified is modeling of the unsteady aerodynamic lift characteristics of the rotor. Finally, the aeroelastic response in the combined presence of wake effects and inflow turbulence is examined. Significant differences are observed in loads, power, and structural response between results computed using the free wake model or simpler models, such as Blade Element Momentum theory.     

低空急流条件下水平轴风力机风轮气动特性的研究

为阐明低空急流条件下风力机风轮的气动特性,基于工程化的边界层风速模型和Von Karman谱模型建立不同来流的脉动风场,对比研究低空急流条件下NREL 5 MW风力机风轮的输出功率和气动载荷的变化规律。结果表明:如果仅以轮毂高度处的风速作为风力机变桨控制的依据,与均匀来流和剪切来流相比较,低空急流条件下,虽然来流风功率明显增大,但风轮的输出功率在较高风速时反而减小;风轮所受的不平衡气动载荷,包括横向力、纵向力、偏航力矩和倾覆力矩在较高风速时小于剪切来流的结果;且仅以轮毂高度处的风速预测得到的风轮输出功率高于实际结果,其最大相对误差为89.4%。因此,低空急流条件下,为提高风能利用率和风轮输出功率的预测精度,应考虑不同高度位置处的风速大小对风力机进行变桨控制和功率预测。     

Sensitivity of southern California wind energy to turbine characteristics

Using output from a high‐resolution meteorological simulation, we evaluate the sensitivity of southern California wind energy generation to variations in key characteristics of current wind turbines. These characteristics include hub height, rotor diameter and rated power, and depend on turbine make and model. They shape the turbine's power curve and thus have large implications for the energy generation capacity of wind farms. For each characteristic, we find complex and substantial geographical variations in the sensitivity of energy generation. However, the sensitivity associated with each characteristic can be predicted by a single corresponding climate statistic, greatly simplifying understanding of the relationship between climate and turbine optimization for energy production. In the case of the sensitivity to rotor diameter, the change in energy output per unit change in rotor diameter at any location is directly proportional to the weighted average wind speed between the cut‐in speed and the rated speed. The sensitivity to rated power variations is likewise captured by the percent of the wind speed distribution between the turbines rated and cut‐out speeds. Finally, the sensitivity to hub height is proportional to lower atmospheric wind shear. Using a wind turbine component cost model, we also evaluate energy output increase per dollar investment in each turbine characteristic. We find that rotor diameter increases typically provide a much larger wind energy boost per dollar invested, although there are some zones where investment in the other two characteristics is competitive. Our study underscores the need for joint analysis of regional climate, turbine engineering and economic modeling to optimize wind energy production. Copyright © 2012 John Wiley & Sons, Ltd.     

叶片转角对小型Savonius风机气动性能的影响

Savonius风机是一种典型的垂直轴风力发电机,通过对其进行流固耦合分析,研究叶片转角对风机气动性能的影响。利用ANSYS的CFX流体模块,流体湍流模型选择基于RANS的标准k-ε湍流模型,对风轮进行流固耦合分析,从而获得叶片产生的力矩情况,并计算了风机的功率特性。利用求解结果,得到了力矩系数与叶片转角之间的关系。分析了风机叶片在旋转一周中所产生的最大扭矩以及负扭矩所处的位置和范围。通过分析转角对风机性能的影响,可为今后的Savonius风机叶片形状优化和效率提升提供参考。     

Performance investigation of a power augmented vertical axis wind turbine for urban high-rise application

A shrouded wind turbine system has a number of potential advantages over the conventional wind turbine. A novel power-augmentation-guide-vane (PAGV) that surrounds a Sistan wind turbine was designed to improve the wind rotor performance by increasing the on-coming wind speed and guiding it to an optimum flow angle before it interacts with the rotor blades. The integration of the PAGV into the 3-in-1 wind, solar and rain water harvester on high-rise buildings has been illustrated. A particular concern related to public safety is minimized when the wind turbine is contained inside the PAGV and noise pollution can be reduced due to the enclosure. Besides, the design of the PAGV that blends into the building architecture can be aesthetic as well. Moreover, a mesh can be mounted around the PAGV to prevent the bird-strike problem. From the wind tunnel testing measurements where the wind turbine is under free-running condition (only rotor inertia and bearing friction were applied), the wind rotor rotational speed (with the PAGV) was increased by 75.16%. Meanwhile, a computational fluid dynamics (CFD) simulation shows that the rotor torque was increased by 2.88 times with the introduction of the PAGV. Through a semi-empirical method, the power output increment of the rotor with the PAGV was 5.8 times at the wind speed of 3 m/s. Also, the flow vector visualization (CFD) shows that a larger area of upstream flow was induced through the rotor with the PAGV.     

Multidisciplinary design optimization of offshore wind turbines for minimum levelized cost of energy

This paper presents a method for multidisciplinary design optimization of offshore wind turbines at system level. The formulation and implementation that enable the integrated aerodynamic and structural design of the rotor and tower simultaneously are detailed. The objective function to be minimized is the levelized cost of energy. The model includes various design constraints: stresses, deflections, modal frequencies and fatigue limits along different stations of the blade and tower. The rotor design variables are: chord and twist distribution, blade length, rated rotational speed and structural thicknesses along the span. The tower design variables are: tower thickness and diameter distribution, as well as the tower height. For the other wind turbine components, a representative mass model is used to include their dynamic interactions in the system. To calculate the system costs, representative cost models of a wind turbine located in an offshore wind farm are used. To show the potential of the method and to verify its usefulness, the 5 MW NREL wind turbine is used as a case study. The result of the design optimization process shows 2.3% decrease in the levelized cost of energy for a representative Dutch site, while satisfying all the design constraints.     

Numerical evaluation of aerodynamic and inertial contributions to Darrieus wind turbine blade deformation

This paper presents a model for the evaluation of aerodynamic and inertial contributions to a vertical-axis wind turbine (VAWT) blade deformation. Through the use of a specially designed coupling code, a solid modeling software, capable of generating the desired blade geometry depending on the design geometric parameters, is linked to a finite volume Computational Fluid Dynamic (CFD) code for the calculation of rotor performance and to a Finite Element Method (FEM) code for the structural design analysis of rotor blades. After describing the computational model and the relative validation procedure, a full RANS unsteady calculation is presented for a three-bladed rotor architecture, characterized by a NACA 0012 profile. Flow field characteristics are investigated for a constant unperturbed free-stream wind velocity of 9 m/s, determining the torque coefficient generated from the three blades as a function of rotor azimuthal coordinate. The emphasis is subsequently placed on obtaining an estimate for both pressure/tangential forces and centrifugal ones to blade structural loadings, thus assessing the influence of aerodynamic and inertial contributions to blade stresses and deformations.