Blade loading on tidal turbines for uniform unsteady flow

An improved characterisation of the unsteady hydrodynamic loads on tidal turbine blades is necessary to enable more reliable predictions of their fatigue life and to avoid premature failures. To this end, this paper presents a set of blade-root bending moment responses for a scale-model tidal turbine subjected to an unsteady planar forcing in a towing tank. In cases where the boundary layer was believed to be attached to the outer sections of the blade, the out-of-plane bending moment amplitude for unsteady flow was up to 15% greater than the corresponding load measured in steady flow and exhibited a phase-lead of up to 4.5°. Both these observations are qualitatively consistent with the effects of dynamic inflow and non-circulatory forcing. The bending moment responses for a forcing time history that comprised multiple frequencies, as well as for a discrete half-sinusoidal perturbation, were able to be reconstructed reasonably well using the responses obtained from single-frequency oscillatory flows. This suggests that blade designers could utilise relatively low fidelity techniques and conduct potentially fewer experimental tests to acquire the fatigue load spectrum.     

Development of an aeroelastic code based on three‐dimensional viscous–inviscid method for wind turbine computations

Aerodynamic and structural dynamic performance analysis of modern wind turbines are routinely estimated in the wind energy field using computational tools known as aeroelastic codes. Most aeroelastic codes use the blade element momentum (BEM) technique to model the rotor aerodynamics and a modal, multi‐body or the finite‐element approach to model the turbine structural dynamics. The present work describes the development of a novel aeroelastic code that combines a three‐dimensional viscous–inviscid interactive method, method for interactive rotor aerodynamic simulations (MIRAS), with the structural dynamics model used in the aeroelastic code FLEX5. The new code, called MIRAS‐FLEX, is an improvement on standard aeroelastic codes because it uses a more advanced aerodynamic model than BEM. With the new aeroelastic code, more physical aerodynamic predictions than BEM can be obtained as BEM uses empirical relations, such as tip loss corrections, to determine the flow around a rotor. Although more costly than BEM, a small cluster is sufficient to run MIRAS‐FLEX in a fast and easy way. MIRAS‐FLEX is compared against the widely used FLEX5 and FAST, as well as the participant codes from the Offshore Code Comparison Collaboration Project. Simulation tests consist of steady wind inflow conditions with different combinations of yaw error, wind shear, tower shadow and turbine‐elastic modeling. Turbulent inflow created by using a Mann box is also considered. MIRAS‐FLEX results, such as blade tip deflections and root‐bending moments, are generally in good agreement with the other codes. Copyright © 2017 John Wiley & Sons, Ltd.     

Experimental measurements of the hydrodynamic performance and structural loading of the Transverse Horizontal Axis Water Turbine: Part 3

This paper is the third of three, which describe the procedures and results for a set of experiments on various configurations of the Transverse Horizontal Axis Water Turbine (THAWT), which is a horizontally orientated variant of the Darrieus cross-flow turbine. Tests were conducted in the combined wind, wave and current tank at Newcastle University on a 0.5 m diameter rotor, while the flow depth and velocity were varied over a range of realistic Froude numbers for tidal streams. Various configurations of the device were tested to assess the merits of varied blade pitch, rotor solidity, blockage ratio and truss oriented blades. Experiments were carried out using a speed-controlled motor/generator, allowing quasi-steady results to be taken over a range of tip speed ratios. Measurements of power, thrust, blade loading and free surface deformation provide extensive data for future validation of numerical codes and demonstrate the ability of the device to exceed the Lanchester–Betz limit for kinetic efficiency by using high blockage. This paper covers the instrumentation, hydrodynamic performance and loading of the truss bladed variant of the THAWT device. The first paper covers the experimental setup and hydrodynamic performance of the parallel bladed rotor and the second paper covers the instrumentation and hydrodynamic loading of the parallel bladed rotor.     

Experimental measurements of the hydrodynamic performance and structural loading of the Transverse Horizontal Axis Water Turbine: Part 1

This paper is the first of three, which outline the procedures and results for a set of experiments carried out on various configurations of the Transverse Horizontal Axis Water Turbine (THAWT), which is a horizontally orientated variant of the Darrieus cross-flow turbine. Tests were conducted in the combined wind, wave and current tank at Newcastle University on a 0.5 m diameter rotor, while the flow depth and velocity were varied over a range of realistic Froude numbers for tidal streams. Various configurations of the device were tested to assess the merits of varied blade pitch, rotor solidity, blockage ratio and truss oriented blades. Experiments were carried out using a speed controlled motor/generator, allowing quasi-steady results to be taken over a range of tip speed ratios. Measurements of power, thrust, blade loading and free surface deformation provide extensive data for future validation of numerical codes and demonstrate the ability of the device to exceed the Lanchester-Betz limit for kinetic efficiency by using high blockage. This paper covers the experimental procedures and results for the hydrodynamic performance for the parallel bladed variant of the THAWT device. The second paper covers the hydrodynamic loading of the parallel bladed rotor and the third covers both hydrodynamic performance and loading of the truss configured THAWT device.     

Experimental measurements of the hydrodynamic performance and structural loading of the transverse horizontal axis water turbine: Part 2

This paper is the second of three, which outline the procedures and results for a set of experiments on various configurations of the Transverse Horizontal Axis Water Turbine (THAWT), which is a horizontally orientated variant of the Darrieus cross-flow turbine. Tests were conducted in the combined wind, wave and current tank at Newcastle University on a 0.5 m diameter rotor, while the flow depth and velocity were varied over a range of realistic Froude numbers for tidal streams. Various configurations of the device were tested to assess the merits of varied blade pitch, rotor solidity, blockage ratio and truss oriented blades. Experiments were carried out using a speed controlled motor/generator, allowing quasi-steady results to be taken over a range of tip speed ratios. Measurements of power, thrust, blade loading and free surface deformation provide extensive data for future validation of numerical codes and demonstrate the ability of the device to exceed the Lanchester–Betz limit for kinetic efficiency, by exploiting high blockage. This second paper covers the instrumentation and analysis for the structural loading for the parallel bladed variant of the THAWT device. The first paper covers the experimental setup and hydrodynamic performance of the parallel bladed rotor, and the third paper covers both performance and loading of the truss configured THAWT device.     

Power output performance characteristics of a horizontal axis marine current turbine

A 0.4 m diameter horizontal axis marine current turbine was tested in a circulating water channel. The power output was measured over a range of flow speeds, blade pitch and rotor yaw angles. Experimental results were compared with the modelled output determined from a commercial blade element momentum computer package. The measured power output was found to be far in excess of predicted values for high blade inflow angles. This occurred where approximately half the blade was operating above the stall angle of attack. This represents 25% of the rotor disk area producing power under heavy stall. Values of overpower up to 140% were measured which are comparable to previous studies. The results show that power production and the optimum tip speed ratio reduced with yaw except for cases with high blade inflow angles.     

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

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

风剪切下塔影效应对水平轴风力机叶片及风轮气动载荷的影响

采用CFD方法,以NH1500三叶片大型水平轴风力机为研究对象,研究额定风速剪切来流下的塔影效应对水平轴风力机叶片和风轮非定常气动载荷的影响。结果表明：剪切来流下,叶片和风轮的气动载荷均呈余弦变化规律,塔影效应的主要影响叶片方位角范围为160°~210°,且该范围不随风剪切指数的变化而变化。相同风剪切指数下,塔影效应对叶片和风轮气动载荷的均方根影响较小,对其波动影响较大。当风剪切指数从0.12增至0.30时,塔影效应下,叶片气动载荷的均方根减小,推力和转矩的波动幅度增大,偏航力矩和倾覆力矩的波动幅度减小;风轮推力和转矩的均方根减小,波动幅度变化较小,而倾覆力矩和偏航力矩的均方根增大,且波动幅度也增大。     

Validation of the Beddoes–Leishman dynamic stall model for horizontal axis wind turbines using MEXICO data

The aim of this study is to assess the load predicting capability of a classical Beddoes–Leishman dynamic stall model in a horizontal axis wind turbine environment, in the presence of yaw misalignment. The dynamic stall model was tailored to the horizontal axis wind turbine environment and validated against unsteady thick airfoil data. Subsequently, the dynamic stall model was implemented in a blade element‐momentum code for yawed flow, and the results were compared with aerodynamic measurements obtained in the MEXICO (Model Rotor Experiments under Controlled Conditions) project on a wind turbine rotor placed in a large scale wind tunnel. In general, reasonable to good agreement was found between the blade element‐momentum model and MEXICO data. When large yaw misalignments were imposed, poor agreement was found in the downstroke of the movement between the model and the experiment. Still, over a revolution, the maximum normal force coefficient predicted was always within 8% of experimental data at the inboard stations, which is encouraging especially when blade fatigue calculations are being considered. Copyright © 2012 John Wiley & Sons, Ltd.     

Load alleviation of wind turbines by yaw misalignment

Vertical wind shear is one of the dominating causes of load variations on the blades of a horizontal axis wind turbine. To alleviate the varying loads, wind turbine control systems have been augmented with sensors and actuators for individual pitch control. However, the loads caused by a vertical wind shear can also be affected through yaw misalignment. Recent studies of yaw control have been focused on improving the yaw alignment to increase the power capture at below rated wind speeds. In this study, the potential of alleviating blade load variations induced by the wind shear through yaw misalignment is assessed. The study is performed through simulations of a reference turbine. The study shows that optimal yaw misalignment angles for minimizing the blade load variations can be identified for both deterministic and turbulent inflows. It is shown that the optimal yaw misalignment angles can be applied without power loss for wind speeds above rated wind speed. In deterministic inflow, it is shown that the range of the steady‐state blade load variations can be reduced by up to 70%. For turbulent inflows, it is shown that the potential blade fatigue load reductions depend on the turbulence level. In inflows with high levels of turbulence, the observed blade fatigue load reductions are small, whereas the blade fatigue loads are reduced by 20% at low turbulence levels. For both deterministic and turbulent inflows, it is seen that the blade load reductions are penalized by increased load variations on the non‐rotating turbine parts. Copyright © 2013 John Wiley & Sons, Ltd.     

Characterization of the unsteady flow in the nacelle region of a modern wind turbine

A three‐dimensional Navier–Stokes solver has been used to investigate the flow in the nacelle region of a wind turbine where anemometers are typically placed to measure the flow speed and the turbine yaw angle. A 500 kW turbine was modelled with rotor and nacelle geometry in order to capture the complex separated flow in the blade root region of the rotor. A number of steady state and unsteady simulations were carried out for wind speeds ranging from 6 m s^?1 to 16 m s^?1 as well as two yaw and tilt angles. The flow in the nacelle region was found to be highly unsteady, dominated by unsteady vortex shedding from the cylindrical part of the blades, which interacted with the root vortices from each blade, generating high‐velocity gradients. As a consequence, the nacelle wind speed and the nacelle flow angle were found to vary significantly with the height above the nacelle surface. The nacelle anemometry showed significant dependence on both yaw and tilt angles with yaw errors of up to 10 degrees when operating in a tilted inflow. Copyright © 2010 John Wiley & Sons, Ltd.     

The impact of yaw error on aeroelastic characteristics of a horizontal axis wind turbine blade

Horizontal axis wind turbines operate under yawed conditions for a considerable period of time due to the power control mechanism or sudden changes in the wind direction. This in turn can alter the dynamic characteristics of a turbine blade because the flow over the rotor plane may trigger complicated induced velocity patterns. In this study, an aeroelastic analysis under yawed flow conditions is carried out to investigate the effects of yaw error on the blade behaviors and dynamic stability. A beam model including geometric nonlinearity coupled with unsteady aerodynamics based on a free-vortex wake method with the blade element theory is employed in the present study. The aerodynamic approach for a horizontal axis wind turbine blade under yawed flow conditions is verified through comparison with measurements. It is also shown that the present method gives slightly better results at high yaw angles than does the method previously published in the literature. The dynamic instabilities of a National Renewable Energy Laboratory 5 MW reference wind turbine have subsequently been investigated for various wind speeds and yaw angles. Observations are made that yaw effects induce considerable changes in airloads and blade structural behavior. Also, the aeroelastic damping values for this particular blade under yawed flow conditions can be reduced by up to approximately 33% in the worst case. Therefore, it is concluded that the impacts of yaw misalignments adversely influenced the dynamic aeroelastic stability of the horizontal axis wind turbine blade.     

Wind turbine response to parameter variation of analytic inflow vortices

As larger wind turbines are placed on taller towers, rotors frequently operate in atmospheric conditions that support organized, coherent turbulent structures. It is hypothesized that these structures have a detrimental impact on the blade fatigue life experienced by the wind turbine. These structures are extremely difficult to identify with sophisticated anemometry such as ultrasonic anemometers. This study was performed to identify the vortex characteristics that contribute to high‐amplitude cyclic blade loads, assuming that these vortices exist under certain atmospheric conditions. This study does not attempt to demonstrate the existence of these coherent turbulent structures. In order to ascertain the idealized worst‐case scenario for vortical inflow structures impinging on a wind turbine rotor, we created a simple, analytic vortex model. The Rankine vortex model assumes that the vortex core undergoes solid body rotation to avoid a singularity at the vortex centre and is surrounded by a two‐dimensional potential flow field. Using the wind turbine as a sensor and the FAST wind turbine dynamics code with limited degrees of freedom, we determined the aerodynamic loads imparted to the wind turbine by the vortex structure. We varied the size, strength, rotational direction, plane of rotation, and location of the vortex over a wide range of operating parameters. We identified the vortex conformation with the most significant effect on the blade root bending moment cyclic amplitude. Vortices with radii on the scale of the rotor diameter or smaller caused blade root bending moment cyclic amplitudes that contribute to high damage density. The rotational orientation, clockwise or counter‐clockwise, produces little difference in the bending moment response. Vortices in the XZ plane produce bending moment amplitudes significantly greater than vortices in the YZ plane. Published in 2005 by John Wiley & Sons, Ltd.     

Prediction of aerodynamic performance of NREL offshore 5-MW baseline wind turbine considering power loss at varying wind speeds

In this study, a computational fluid dynamics (CFD) model was developed to simulate the aerodynamic performance of the National Renewable Energy Laboratory (NREL) offshore 5-MW baseline wind turbine with single rotor and full wind turbine. Using statistical methods, the relation between pitch angle and performance when the speed is above the rated wind speed was analyzed; furthermore, other published data were compiled to establish the functional equations of power, thrust with various inflow wind speeds, and pitch angles. In addition, according to shape optimization based on parametric modeling, the two-dimensional cross-section of the wind turbine blade can be defined through a metasurface approach, and the three-dimensional surface modeling of the wind turbine blade, nacelle, and tower is completed using the nonuniform rational B-splines (NURBS) interpolator. In terms of aerodynamic simulation, the unstructured polygon mesh was used herein to discretize the space and simulate the highly curved and twisted surfaces of the blade. In this study, the compact and accurate mesh form obtained through a grid independence test will be used to analyze the distribution of the pressure coefficient, shear stress coefficient, and limited streamline on the blade surface at various inflow wind speeds and pitch angles to understand the differences between different turbulence models and the causes of power and thrust attenuation at high inflow wind speeds. In addition, the phenomena of blade-tip vortices, dynamic stall, blade loading, and the interaction between nacelle and tower were collectively explored.     

Modeling turbine wakes and power losses within a wind farm using LES: An application to the Horns Rev offshore wind farm

A modeling framework is proposed and validated to simulate turbine wakes and associated power losses in wind farms. It combines the large-eddy simulation (LES) technique with blade element theory and a turbine-model-specific relationship between shaft torque and rotational speed. In the LES, the turbulent subgrid-scale stresses are parameterized with a tuning-free Lagrangian scale-dependent dynamic model. The turbine-induced forces and turbine-generated power are modeled using a recently developed actuator-disk model with rotation (ADM-R), which adopts blade element theory to calculate the lift and drag forces (that produce thrust, rotor shaft torque and power) based on the local simulated flow and the blade characteristics. In order to predict simultaneously the turbine angular velocity and the turbine-induced forces (and thus the power output), a new iterative dynamic procedure is developed to couple the ADM-R turbine model with a relationship between shaft torque and rotational speed. This relationship, which is unique for a given turbine model and independent of the inflow condition, is derived from simulations of a stand-alone wind turbine in conditions for which the thrust coefficient can be validated. Comparison with observed power data from the Horns Rev wind farm shows that better power predictions are obtained with the dynamic ADM-R than with the standard ADM, which assumes a uniform thrust distribution and ignores the torque effect on the turbine wakes and rotor power. The results are also compared with the power predictions obtained using two commercial wind-farm design tools (WindSim and WAsP). These models are found to underestimate the power output compared with the results from the proposed LES framework.     

Numerical investigation of unsteady aerodynamics of a Horizontal‐axis wind turbine under yawed flow conditions

In the present study, unsteady flow features and the blade aerodynamic loading of the National Renewable Energy Laboratory phase VI wind turbine rotor, under yawed flow conditions, were numerically investigated by using a three‐dimensional incompressible flow solver based on unstructured overset meshes. The effect of turbulence, including laminar‐turbulent transition, was accounted for by using a correlation‐based transition turbulence model. The calculations were made for an upwind configuration at wind speeds of 7, 10 and 15 m/sec when the turbine rotor was at 30° and 60° yaw angles. The results were compared with measurements in terms of the blade surface pressure and the normal and tangential forces at selected blade radial locations. It was found that under the yawed flow conditions, the blade aerodynamic loading is significantly reduced. Also, because of the wind velocity component aligned tangent to the rotor disk plane, the periodic fluctuation of blade loading is obtained with lower magnitudes at the advancing blade side and higher magnitudes at the retreating side. This tendency is further magnified as the yaw angle becomes larger. At 7 m/sec wind speed, the sectional angle of attack is relatively small, and the flow remains mostly attached to the blade surface. At 10 m/sec wind speed, leading‐edge flow separation and strong radial flow are observed at the inboard portion of the retreating blade. As the wind speed is further increased, the flow separation and the radial flow become more pronounced. It was demonstrated that these highly unsteady three‐dimensional aerodynamic features are well‐captured by the present method. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessing the influence of inflow turbulence on noise and performance of a tidal turbine using large eddy simulations

Large eddy simulations of a model scale tidal turbine encountering inflow turbulence have been performed. This has allowed both unsteady blade loading and hydrodynamic noise radiation to be predicted. The study is motivated by the need to assess environmental impact of tidal devices, in terms of their acoustic impact on marine species.Inflow turbulence was accounted for using a synthetic turbulence generator, with statistics chosen to represent the gross features of a typical tidal flow. The turbine is resolved in a fully unsteady manner using a sliding interface technique within the OpenFOAM^® libraries. Acoustic radiation is estimated using a compact source approximation of the Ffowcs Williams–Hawkings equation.It is observed that the long streamwise length scale of the inflow turbulence results in characteristic ‘humps’ in the turbine thrust and torque spectra. This effect is also evident in the far-field noise spectra. The acoustic sources on the blades are visualised in terms of sound pressure level and “Powell's source term”. These measures show that the dominant sources are concentrated at the blade leading edges towards the tip. This results from the high loading of the turbine blades, and causes the sound to radiate more akin to a monopole than a dipole.The full scale source level, obtained from scaling of the simulation results, is found to be lower than comparable measured data reported in the literature; this is attributed to additional sources not included in the present study. Based on the predicted source level, no physical impact on fish is expected.     

Determination of the angle of attack on rotor blades

Two simple methods for determining the angle of attack (AOA) on a section of a rotor blade are proposed. Both techniques consist of employing the Biot–Savart integral to determine the influence of the bound vorticity on the velocity field. In the first technique, the force distribution along the blade and the velocity at a monitor point in the vicinity of the blade are assumed to be known from experiments or CFD computations. The AOA is determined by subtracting the velocity induced by the bound circulation, determined from the loading, from the velocity at the monitor point. In the second method, the full pressure distribution on the blade is assumed to be known and used to determine the local distribution of circulation along the surface contour of the blade. Using the local distribution of circulation to determine the influence of the bound vorticity enables the velocity monitor points to be located closer to the blade, and thus to determine the AOA with higher accuracy. Data from CFD computations for flows past the Tellus 95 kW wind turbine at different wind speeds are used to test both techniques. Comparisons show that the proposed methods are in good agreement with existing techniques. The advantage of the proposed techniques, as compared with existing techniques, is that they can be used to determine the AOA on rotor blades under general flow conditions (e.g. operations in yaw or with dynamic inflow). Copyright © 2008 John Wiley & Sons, Ltd.     

Effect of rotor–tower interaction,tilt angle,and yaw misalignment on the aeroelasticity of a large horizontal axis wind turbine with composite blades

This paper presents a detailed analysis of the rotor–tower interaction and the effects of the rotor's tilt angle and yaw misalignment on a large horizontal axis wind turbine. A high‐fidelity aeroelastic model is employed, coupling computational fluid dynamics (CFD) and structural mechanics (CSM). The wind velocity stratification induced by the atmospheric boundary layer (ABL) is modeled. On the CSM side, the complex composite structure of each blade is accurately modeled using shell elements. The rotor–tower interaction is analyzed by comparing results of a rotor‐only simulation and a full‐machine simulation, observing a sudden drop in loads, deformations, and power production of each blade, when passing in front of the tower. Subsequently, a tilt angle is introduced on the rotor, and its effect on blade displacements, loads, and performance is studied, representing a novelty with respect to the available literature. The tilt angle leads to a different contribution of gravity to the blade deformations, sensibly affecting the stresses in the composite material. Lastly, a yaw misalignment is introduced with respect to the incoming wind, and the resulting changes in the blade solicitations are analyzed. In particular, a reduction of the blade axial displacement amplitude during each revolution is observed.     

基于滑移网格的垂直轴潮流水轮机的三维数值模拟

为了探究翼型对垂直轴水轮机水力效率的影响,基于叶素理论分析建立了垂直轴潮流水轮机在水槽中的物理模型,采用滑移网格技术在Fluent软件中对模型的流场进行了三维数值模拟。在保持转速一定、更改来流速度即改变叶尖速比的条件下,分析了两种不同NACA翼型直叶片的潮流水轮机内部流场以及水力性能。结果表明,翼型以及叶片安装角对垂直轴潮流水轮机的利用效率影响很大,其流动特性与来流速度、叶片布置形式有重要联系,为新型潮流水轮机的设计和翼型的选择应用提供了依据。