Simulation of wind turbine wake interaction using the vorticity transport model

The aerodynamic interactions that can occur within a wind farm can result in the constituent turbines generating a lower power output than would be possible if each of the turbines were operated in isolation. Tightening of the constraints on the siting of wind farms is likely to increase the scale of the problem in the future. The aerodynamic performance of turbine rotors and the mechanisms that couple the fluid dynamics of multiple rotors can be most readily understood by simplifying the problem and considering the interaction between only two rotors. The aerodynamic interaction between two rotors in both co‐axial and offset configurations has been simulated using the Vorticity Transport Model. The aerodynamic interaction is a function of the tip speed ratio, and both the streamwise and crosswind separation between the rotors. The simulations show that the momentum deficit at a turbine operating within the wake developed by the rotor of a second turbine is governed by the development of instabilities within the wake of the upwind rotor, and the ensuing structure of the wake as it impinges on the downwind rotor. If the wind farm configuration or wind conditions are such that a turbine rotor is subject to partial impingement by the wake produced by an upstream turbine, then significant unsteadiness in the aerodynamic loading on the rotor blades of the downwind turbine can result, and this unsteadiness can have considerable implications for the fatigue life of the blade structure and rotor hub. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of the actuator surface concept to wind turbine rotor aerodynamics

The structure of blade tip vortices is recognized as a key issue in wind turbine aerodynamic modelling by many researchers in the field. In the search for an intermediate model between full Navier–Stokes and blade‐element momentum simulations, this article presents a method using rotating actuator surfaces to model wind turbine aerodynamics. An actuator surface is a simple planar surface, porous to the flow, which is characterized by velocity and pressure discontinuities, whose action on the flow is achieved through an attached system of forces. These discontinuities and forces are determined from blade‐element analysis and the Kutta–Joukowski relation. After implementing this concept in a three‐dimensional CFD (Computational Fluid Dynamics) method, results are produced for the experimental rotors of NREL and TUDelft. The method is validated against both experimental measurements and the predictions of three other numerical models for wind turbine aerodynamic analysis. Qualitative and quantitative comparisons show that the actuator surface concept agrees well with the other numerical models. In addition to rotor aerodynamic analysis, the actuator surface concept can be used in the study of wake aerodynamics, or as the Eulerian flow solver in hybrid methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Isogeometric based framework for aeroelastic wind turbine blade analysis

A framework based on isogeometric analysis is presented for parametrizing a wind turbine rotor blade and evaluating its response. The framework consists of a multi‐fidelity approach for wind turbine rotor analysis. The aeroelastic loads are determined using a low‐fidelity model. The model is based on isogeometric approach to model both the structural and aerodynamic properties. The structural deformations are solved using an isogeometric formulation of geometrically exact 3D beam theory. The aerodynamic loads are calculated using a standard Blade Element Momentum(BEM) theory. Moreover, the aerodynamic loads calculated using BEM theory are modified to account for the change in the blade shape due to blade deformation. The aeroelastic loads are applied in finite element solver Nastran, and both the stress response and buckling response are extracted. Furthermore, the capabilities of Nastran are extended such that design dependent loads can be applied, resulting in correct aeroelastic sensitivities of Nastran responses, making this framework suitable for optimization. The framework is verified against results from the commercial codes FAST and GH Bladed, using the NREL 61.5m rotor blade as a baseline for comparison, showing good agreement. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessment of a comprehensive aeroelastic tool for horizontal‐axis wind turbine rotor analysis

This paper presents the development of a computational aeroelastic tool for the analysis of performance, response and stability of horizontal‐axis wind turbines. A nonlinear beam model for blades structural dynamics is coupled with a state‐space model for unsteady sectional aerodynamic loads, including dynamic stall effects. Several computational fluid dynamics structural dynamics coupling approaches are investigated to take into account rotor wake inflow influence on downwash, all based on a Boundary Element Method for the solution of incompressible, potential, attached flows. Sectional steady aerodynamic coefficients are extended to high angles of attack in order to characterize wind turbine operations in deep stall regimes. The Galerkin method is applied to the resulting aeroelastic differential system. In this context, a novel approach for the spatial integration of additional aerodynamic states, related to wake vorticity and dynamic stall, is introduced and assessed. Steady‐periodic blade responses are evaluated by a harmonic balance approach, whilst a standard eigenproblem is solved for aeroelastic stability analyses. Drawbacks and potentialities of the proposed model are investigated through numerical and experimental comparisons, with particular attention to rotor blades unsteady aerodynamic modelling issues. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Numerical investigation of wind turbine wake development in directionally sheared inflow

Turbines in wind farms are subject to complex mutual aerodynamic interactions, which in detail depend upon the characteristics of the atmospheric boundary layer. Our two objectives with this paper were to investigate the impact of directionally sheared inflow on the wake development behind a single wind turbine and to analyse the impact of the wakes on the energy yield and loading of a downstream turbine, which is exposed to partial and full wake conditions. We performed simulations with a framework based on a coupled approach of large‐eddy simulation and an actuator line representation of an aeroelastic turbine model. Our results show that directionally sheared inflow leads to a non‐symmetrical wake development, which transfers to distinct differences in the energy yield and loading of downstream turbines of equal lateral offsets in opposite direction. Therefore, the assumption of wakes being axisymmetrical could lead to notable deviations in the prediction of wake behaviour and their impact on downstream turbines for atmospheric inflow conditions, which include directional shear. Copyright © 2016 John Wiley & Sons, Ltd.     

Validation of the dynamic wake meander model for loads and power production in the Egmond aan Zee wind farm

This paper investigates wake effects on load and power production by using the dynamic wake meander (DWM) model implemented in the aeroelastic code HAWC2. The instationary wind farm flow characteristics are modeled by treating the wind turbine wakes as passive tracers transported downstream using a meandering process driven by the low frequent cross‐wind turbulence components. The model complex is validated by comparing simulated and measured loads for the Dutch Egmond aan Zee wind farm consisting of 36 Vestas V90 turbine located outside the coast of the Netherlands. Loads and production are compared for two distinct wind directions—a free wind situation from the dominating southwest and a full wake situation from northwest, where the observed turbine is operating in wake from five turbines in a row with 7D spacing. The measurements have a very high quality, allowing for detailed comparison of both fatigue and min–mean–max loads for blade root flap, tower yaw and tower bottom bending moments, respectively. Since the observed turbine is located deep inside a row of turbines, a new method on how to handle multiple wakes interaction is proposed. The agreement between measurements and simulations is excellent regarding power production in both free and wake sector, and a very good agreement is seen for the load comparisons too. This enables the conclusion that wake meandering, caused by large scale ambient turbulence, is indeed an important contribution to wake loading in wind farms. Copyright © 2012 John Wiley & Sons, Ltd.     

风力机叶片动态气弹变形及其对整机性能的影响

基于风力机整机刚柔耦合模型,文章提出了一种叶片动态气弹扭转变形分析的新方法。该方法采用SIMPACK和AeroDyn软件联合数值仿真对风力机在几种恶劣风况下进行动力学分析,通过对分析结果的变换处理,进而得到叶片在复杂工况下的动态气弹变形数据。采用该方法,重点分析了叶片气弹扭转变形对风力机气动功率及气弹稳定性的影响。该方法为大型风电叶片的气弹特性评价以及气弹剪裁设计提供了一种新的技术手段。     

Simulating the aerodynamic performance and wake dynamics of a vertical‐axis wind turbine

The accurate prediction of the aerodynamics and performance of vertical‐axis wind turbines is essential if their design is to be improved but poses a significant challenge to numerical simulation tools. The cyclic motion of the blades induces large variations in the angle of attack of the blades that can manifest as dynamic stall. In addition, predicting the interaction between the blades and the wake developed by the rotor requires a high‐fidelity representation of the vortical structures within the flow field in which the turbine operates. The aerodynamic performance and wake dynamics of a Darrieus‐type vertical‐axis wind turbine consisting of two straight blades is simulated using Brown's Vorticity Transport Model. The predicted variation with azimuth of the normal and tangential force on the turbine blades compares well with experimental measurements. The interaction between the blades and the vortices that are shed and trailed in previous revolutions of the turbine is shown to have a significant effect on the distribution of aerodynamic loading on the blades. Furthermore, it is suggested that the disagreement between experimental and numerical data that has been presented in previous studies arises because the blade–vortex interactions on the rotor were not modelled with sufficient fidelity. Copyright © 2010 John Wiley & Sons, Ltd.     

Aeroelastic analysis of a wind turbine blade using the harmonic balance method

An aeroelastic model for wind turbine blades derived from the unsteady Navier‐Stokes equations and a mode shape–based structural dynamics model are presented. For turbulent flows, the system is closed with the Spalart‐Allmaras turbulence model. The computation times for the aerodynamic solution are significantly reduced using the harmonic balance method compared to a time‐accurate solution. This model is significantly more robust than standard aeroelastic codes that rely on blade element momentum theory to determine the aerodynamic forces. Comparisons with published results for the Caradonna‐Tung rotor in hover and the classical AGARD 445.6 flutter case are provided to validate the aerodynamic model and aeroelastic model, respectively. For wind turbines, flutter of the 1.5 MW WindPACT blade is considered. The results predict that the first flapwise and edgewise modes dominate flutter at the rotor speeds considered.     

Simulation comparison of wake mitigation control strategies for a two‐turbine case

Wind turbines arranged in a wind plant impact each other through their wakes. Wind plant control is an active research field that attempts to improve wind plant performance by coordinating control of individual turbines to take into account these turbine–wake interactions. In this paper, high‐fidelity simulations of a two‐turbine fully waked scenario are used to investigate several wake mitigation strategies, including modification of yaw and tilt angles of an upstream turbine to induce wake skew, as well as repositioning of the downstream turbine. The simulation results are compared through change relative to a baseline operation in terms of overall power capture and loading on the upstream and downstream turbine. Results demonstrated improved power production for all methods. Analysis of control options, including individual pitch control, shows potential to minimize the increase of, or even reduce, turbine loads.Copyright © 2014 John Wiley & Sons, Ltd.     

Wind tunnel investigation on the effect of the turbine tower on wind turbines wake symmetry

In wind farms, the wake of the upstream turbines becomes the inflow for the downstream machines. Ideally, the turbine wake is a stable vortex system. In reality, because of factors like background turbulence, mean flow shear, and tower‐wake interaction, the wake velocity deficit is not symmetric and is displaced away from its mean position. The irregular velocity profile leads to a decreased efficiency and increased blade stress levels for the downstream turbines. The object of this work is the experimental investigation of the effect of the wind turbine tower on the symmetry and displacement of the wake velocity deficit induced by one and two in‐line model wind turbines (,D= 0.9 m). The results of the experiments, performed in the closed‐loop wind tunnel of the Norwegian University of Science and Technology in Trondheim (Norway), showed that the wake of the single turbine expanded more in the horizontal direction (side‐wall normal) than in the vertical (floor normal) direction and that the center of the wake vortex had a tendency to move toward the wind tunnel floor as it was advected downstream from the rotor. The wake of the turbine tandem showed a similar behavior, with a larger degree of non‐symmetry. The analysis of the cross‐stream velocity profiles revealed that the non‐symmetries were caused by a different cross‐stream momentum transport in the top‐tip and bottom‐tip region, induced by the turbine tower wake. In fact, when a second additional turbine tower, mirroring the original one, was installed above the turbine nacelle, the wake recovered its symmetric structure. Copyright © 2017 John Wiley & Sons, Ltd.     

Wake impacts on downstream wind turbine performance and yaw alignment

Aerodynamic wake interaction between commercial scale wind turbines can be a significant source of power losses and increased fatigue loads across a wind farm. Significant research has been dedicated to the study of wind turbine wakes and wake model development. This paper profiles influential wake regions for an onshore wind farm using 6 months of recorded SCADA (supervisory control and data acquisition) data. An average wind velocity deficit of over 30% was observed corresponding to power coefficient losses of 0.2 in the wake region. Wind speed fluctuations are also quantified for an array of turbines, inferring an increase in turbulence within the wake region. A study of yaw data within the array showed turbine nacelle misalignment under a range of downstream wake angles, indicating a characteristic of wind turbine behaviour not generally considered in wake studies. The turbines yaw independently in order to capture the increased wind speeds present due to the lateral influx of turbulent wind, contrary to many experimental and simulation methods found in the literature. Improvements are suggested for wind farm control strategies that may improve farm‐wide power output. Additionally, possible causes for wind farm wake model overestimation of wake losses are proposed.Copyright © 2012 John Wiley & Sons, Ltd.     

LiDAR measurements for an onshore wind farm: Wake variability for different incoming wind speeds and atmospheric stability regimes

Wind measurements were performed with the UTD mobile LiDAR station for an onshore wind farm located in Texas with the aim of characterizing evolution of wind‐turbine wakes for different hub‐height wind speeds and regimes of the static atmospheric stability. The wind velocity field was measured by means of a scanning Doppler wind LiDAR, while atmospheric boundary layer and turbine parameters were monitored through a met‐tower and SCADA, respectively. The wake measurements are clustered and their ensemble statistics retrieved as functions of the hub‐height wind speed and the atmospheric stability regime, which is characterized either with the Bulk Richardson number or wind turbulence intensity at hub height. The cluster analysis of the LiDAR measurements has singled out that the turbine thrust coefficient is the main parameter driving the variability of the velocity deficit in the near wake. In contrast, atmospheric stability has negligible influence on the near‐wake velocity field, while it affects noticeably the far‐wake evolution and recovery. A secondary effect on wake‐recovery rate is observed as a function of the rotor thrust coefficient. For higher thrust coefficients, the enhanced wake‐generated turbulence fosters wake recovery. A semi‐empirical model is formulated to predict the maximum wake velocity deficit as a function of the downstream distance using the rotor thrust coefficient and the incoming turbulence intensity at hub height as input. The cluster analysis of the LiDAR measurements and the ensemble statistics calculated through the Barnes scheme have enabled to generate a valuable dataset for development and assessment of wind farm models.     

Parabolic RANS solver for low‐computational‐cost simulations of wind turbine wakes

A numerical framework for simulations of wake interactions associated with a wind turbine column is presented. A Reynolds‐averaged Navier‐Stokes (RANS) solver is developed for axisymmetric wake flows using parabolic and boundary‐layer approximations to reduce computational cost while capturing the essential wake physics. Turbulence effects on downstream evolution of the time‐averaged wake velocity field are taken into account through Boussinesq hypothesis and a mixing length model, which is only a function of the streamwise location. The calibration of the turbulence closure model is performed through wake turbulence statistics obtained from large‐eddy simulations of wind turbine wakes. This strategy ensures capturing the proper wake mixing level for a given incoming turbulence and turbine operating condition and, thus, accurately estimating the wake velocity field. The power capture from turbines is mimicked as a forcing in the RANS equations through the actuator disk model with rotation. The RANS simulations of the wake velocity field associated with an isolated 5‐MW NREL wind turbine operating with different tip speed ratios and turbulence intensity of the incoming wind agree well with the analogous velocity data obtained through high‐fidelity large‐eddy simulations. Furthermore, different cases of columns of wind turbines operating with different tip speed ratios and downstream spacing are also simulated with great accuracy. Therefore, the proposed RANS solver is a powerful tool for simulations of wind turbine wakes tailored for optimization problems, where a good trade‐off between accuracy and low‐computational cost is desirable.     

A pragmatic approach to wind farm simulations using the dynamic wake meandering model

This paper presents a computationally efficient method for using the dynamic wake meandering model to conduct simulations of wind farm power production. The method is based on creating a database, which contains the time and rotor‐averaged wake effect at any point downstream of a wake‐emitting turbine operating in arbitrary ambient conditions and at an arbitrary degree of wake influence. This database is later used as a look‐up table at runtime to estimate the operating conditions at all turbines in the wind farm, thus eliminating the need to run the dynamic wake meandering model at runtime. By using the proposed method, the time required to conduct wind farm simulations is reduced by three orders of magnitude compared with running the standalone dynamic wake meandering model at runtime. As a result, the wind farm production dynamics for a farm of 100 turbines at 10,000 different sets of ambient conditions run on a normal laptop in 1 h. The method is validated against full scale measurements from the Smøla and OWEZ wind farms, and fair agreement is achieved. Copyright © 2014 John Wiley & Sons, Ltd.     

Global sensitivity analysis of wind turbine power output

The dynamics of wind turbine behavior are complex and a critical area of study for the wind industry. Identification of factors that cause changes in turbine performance can sometimes prove to be challenging, whereas other times, it can be intuitive. The quantification of the effect that these factors have is valuable for making improvements to both power performance and turbine health. In commercial farms, large quantities of meteorological and performance data are commonly collected to monitor daily operations. These data can also be used to analyze the relationship between each parameter in order to better understand the interactions that occur and the information contained within these signals. In this global sensitivity analysis, a neural network is used to model select wind turbine supervisory control and data acquisition system parameters for an array of turbines from a commercial wind farm that exhibit signs of wake interaction. An extended Fourier amplitude sensitivity test is then performed for 2 years of 10‐min averaged data. The study examines the primary and combined sensitivities of power output to each selected parameter for two turbines in the array. The primary sensitivities correspond to single parameter interactions, whereas combined sensitivities account for interactions between multiple parameters simultaneously. Highly influential parameters such as wind speed and rotor rotation frequency produce expected results; the extended Fourier amplitude sensitivity test method proved effective at quantifying the sensitivity of a wide range of more subtle inputs. These include blade pitch, yaw position, main bearing and ambient temperatures as well as wind speed and yaw position standard deviation. The technique holds promise for application in full‐scale wake studies where it might be used to determine the benefits of emerging power optimization strategies such as active wake management. The field of structural health monitoring can also benefit from this method. Copyright © 2013 John Wiley & Sons, Ltd.     

Performance and near wake measurements of a model horizontal axis wind turbine

The performance characteristics and the near wake of a model wind turbine were investigated experimentally. The model tested is a three‐bladed horizontal axis type wind turbine with an upstream rotor of 0.90 m diameter. The performance measurements were conducted at various yaw angles, a freestream speed of about 10 m s ^?1, and the tip speed ratio was varied from 0.5 to 12. The time‐averaged streamwise velocity field in the near wake of the turbine was measured at different tip speed ratios and downstream locations. As expected, it was found that power and thrust coefficients decrease with increasing yaw angle. The power loss is about 3% when the yaw angle is less than 10° and increases to more than 30% when the yaw angle is greater than 30°. The velocity distribution in the near wake was found to be strongly influenced by the tip speed ratio and the yaw angle. At the optimum tip speed ratio, the axial velocity was almost uniform within the midsection of the rotor wake, whereas two strong peaks are observed for high tip speed ratios when the yaw angle is 0°. As the yaw angle increases, the wake width was found to be reduced and skewed towards the yawed direction. With increasing downstream distance, the wake velocity field was observed to depend on the tip speed ratio and more pronounced at high tip speed ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

Hybrid vortex simulations of wind turbines using a three‐dimensional viscous–inviscid panel method

A hybrid filament‐mesh vortex method is proposed and validated to predict the aerodynamic performance of wind turbine rotors and to simulate the resulting wake. Its novelty consists of using a hybrid method to accurately simulate the wake downstream of the wind turbine while reducing the computational time used by the method. The proposed method uses a hybrid approach, where the near wake is resolved by using vortex filaments, which carry the vorticity shed by the trailing edge of the blades. The interaction of the vortex filaments in the near vicinity of the wind turbine is evaluated using a direct calculation, whereas the contribution from the large downstream wake is calculated using a mesh‐based method. The hybrid method is first validated in detail against the well‐known MEXICO experiment, using the direct filament method as a comparison. The second part of the validation includes a study of the influence of the time‐integration scheme used for evolving the wake in time, aeroelastic simulations of the National Renewable Energy Laboratory 5 MW wind turbine and an analysis of the central processing unit time showing the gains of using the hybrid filament‐mesh method. Copyright © 2017 John Wiley & Sons, Ltd.