Survey of modelling methods for wind turbine wakes and wind farms

This article provides an overview and analysis of different wake‐modelling methods which may be used as prediction and design tools for both wind turbines and wind farms. We also survey the available data concerning the measurement of wind magnitudes in both single wakes and wind farms, and of loading effects on wind turbines under single‐ and multiple‐wake conditions. The relative merits of existing wake and wind farm models and their ability to reproduce experimental results are discussed. Conclusions are provided concerning the usefulness of the different modelling approaches examined, and difficult issues which have not yet been satisfactorily treated and which require further research are discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Self‐similarity and turbulence characteristics of wind turbine wakes via large‐eddy simulation

Mean and turbulent properties of the wake generated by a single wind turbine are studied in this paper with a new large eddy simulation (LES) code, the wind turbine and turbulence simulator (WiTTS hereafter). WiTTS uses a scale‐dependent Lagrangian dynamical model of the sub‐grid shear stress and actuator lines to simulate the effects of the rotating blades. WiTTS is first tested by simulating neutral boundary layers without and with a wind turbine and then used to study the common assumptions of self‐similarity and axisymmetry of the wake under neutral conditions for a variety of wind speeds and turbine properties. We find that the wind velocity deficit generally remains self similarity to a Gaussian distribution in the horizontal. In the vertical, the Gaussian self‐similarity is still valid in the upper part of the wake, but it breaks down in the region of the wake close to the ground. The horizontal expansion of the wake is always faster and greater than the vertical expansion under neutral stability due to wind shear and impact with the ground. Two modifications to existing equations for the mean velocity deficit and the maximum added turbulence intensity are proposed and successfully tested. The anisotropic wake expansion is taken into account in the modified model of the mean velocity deficit. Turbulent kinetic energy (TKE) budgets show that production and advection exceed dissipation and turbulent transport. The nacelle causes significant increase of every term in the TKE budget in the near wake. In conclusion, WiTTS performs satisfactorily in the rotor region of wind turbine wakes under neutral stability. Copyright © 2014 John Wiley & Sons, Ltd.     

Wind turbine boundary layer arrays for Cartesian and staggered configurations‐Part I,flow field and power measurements

Model wind turbine arrays were developed for the purpose of investigating the wake interaction and turbine canopy layer in a standard cartesian and row‐offset turbine array configurations. Stereographic particle image velocimetry was used to collect flow data upstream and downstream of entrance and exit row turbines in each configuration. Wakes for all cases were analyzed for energy content and recovery behavior including entrainment of high‐momentum flow from above the turbine canopy layer. The row‐offset arrangement of turbines within an array grants an increase in streamwise spacing of devices and allows for greater wake remediation between successive rows. These effects are seen in exit row turbine wakes as changes to statistical quantities including the in‐plane Reynolds stress, , and the production of turbulence. The recovery of wakes also strongly mitigates the perceived underperformance of wind turbines within an array. The flux of kinetic energy is demonstrated to be more localized in the entrance rows and in the offset arrangement. Extreme values for the flux of kinetic energy are about 7.5% less in the exit row of the cartesian arrangement than in the offset arrangement. Measurements of mechanical torque at entrance and exit row turbines lead to curves of power coefficient and demonstrate an increase in efficiency in row‐offset configurations. Copyright © 2014 John Wiley & Sons, Ltd.     

Application and validation of incrementally complex models for wind turbine aerodynamics,isolated wind turbine in uniform inflow conditions

A comparison of several incrementally complex methods for predicting wind turbine performance, aeroelastic behavior, and wakes is provided. Depending on a wind farm's design, wake interference can cause large power losses and increased turbulence levels within the farm. The goal is to employ modeling methods to reach an improved understanding of wake effects and to use this information to better optimize the layout of new wind farms. A critical decision faced by modelers is the fidelity of the model that is selected to perform simulations. The choice of model fidelity can affect the accuracy, but will also greatly impact the computational time and resource requirements for simulations. To help address this critical question, three modeling methods of varying fidelity have been developed side by side and are compared in this article. The models from low to high complexity are as follows: a blade element‐based method with a free‐vortex wake, an actuator disc‐based method, and a full rotor‐based method. Fluid/structure interfaces are developed for the aerodynamic modeling approaches that allow modeling of discrete blades and are then coupled with a multibody structural dynamics solver in order to perform an aeroelastic analysis. Similar methods have individually been tested by researchers, but we suggest that by developing a suite of models, they can be cross‐compared to grasp the subtleties of each method. The modeling methods are applied to the National Renewable Energy Laboratory Phase VI rotor to predict the turbine aerodynamic and structural loads and then also the wind velocities in the wake. The full rotor method provides the most accurate predictions at the turbine and the use of adaptive mesh refinement to capture the wake to 20 radii downstream is proven particularly successful. Though the full rotor method is unmatched by the lower fidelity methods in stalled conditions and detailed prediction of the downstream wake, there are other less complex conditions where these methods perform as accurately as the full rotor method. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical simulations of wake characteristics of a wind turbine in uniform inflow

The wake of a wind turbine operating in a uniform inflow at various tip speed ratios is simulated using a numerical method, which combines large eddy simulations with an actuator line technique. The computations are carried out in a numerical mesh with about 8.4·10⁶ grid points distributed to facilitate detailed studies of basic features of both the near and far wake, including distributions of interference factors, vortex structures and formation of instabilities. Copyright © 2009 John Wiley & Sons, Ltd.     

Actuator volumes and hr‐adaptive methods for three‐dimensional simulation of wind turbine wakes and performance

This paper presents a three‐dimensional numerical model for horizontal axis wind turbines, capable of simulating both the dynamic response of the turbine to changing flow conditions, and the full wake generated by the turbine. The turbine model is coupled to computational fluid dynamics software using Large Eddy Simulation to solve for unsteady flow conditions, with the Smagorinsky method handling sub‐grid turbulence. Wind tunnel simulations are compared with experimental data, and then a simulation of a real turbine in complex terrain is compared with LIDAR data from the same site, in both cases agreeing favourably. Copyright © 2011 John Wiley & Sons, Ltd.     

A stochastic wind turbine wake model based on new metrics for wake characterization

Understanding the detailed dynamics of wind turbine wakes is critical to predicting the performance and maximizing the efficiency of wind farms. This knowledge requires atmospheric data at a high spatial and temporal resolution, which are not easily obtained from direct measurements. Therefore, research is often based on numerical models, which vary in fidelity and computational cost. The simplest models produce axisymmetric wakes and are only valid beyond the near wake. Higher‐fidelity results can be obtained by solving the filtered Navier–Stokes equations at a resolution that is sufficient to resolve the relevant turbulence scales. This work addresses the gap between these two extremes by proposing a stochastic model that produces an unsteady asymmetric wake. The model is developed based on a large‐eddy simulation (LES) of an offshore wind farm. Because there are several ways of characterizing wakes, the first part of this work explores different approaches to defining global wake characteristics. From these, a model is developed that captures essential features of a LES‐generated wake at a small fraction of the cost. The synthetic wake successfully reproduces the mean characteristics of the original LES wake, including its area and stretching patterns, and statistics of the mean azimuthal radius. The mean and standard deviation of the wake width and height are also reproduced. This preliminary study focuses on reproducing the wake shape, while future work will incorporate velocity deficit and meandering, as well as different stability scenarios. Copyright © 2016 John Wiley & Sons, Ltd.     

Three‐dimensional viscous‐inviscid coupling method for wind turbine computations

In this paper, a computational model for predicting the aerodynamic behavior of wind turbine wakes and blades subjected to unsteady motions and viscous effects is presented. The model is based on a three‐dimensional panel method using a surface distribution of quadrilateral sources and doublets, which is coupled to a viscous boundary layer solver. Unlike Navier‐Stokes codes that need to solve the entire flow domain, the panel method solves the flow around a complex geometry by distributing singularity elements on the body surface, obtaining a faster solution and making this type of codes suitable for the design of wind turbines. A free‐wake model has been employed to simulate the wake behind a wind turbine by using vortex filaments that carry the vorticity shed by the trailing edge of the blades. Viscous and rotational effects inside the boundary layer are taken into account via the transpiration velocity concept, applied using strip theory with the cross sectional angle of attack as coupling parameter. The transpiration velocity is obtained from the solution of the integral boundary layer equations with extension for rotational effects. It is found that viscosity plays a very important role in the predictions of blade aerodynamics and wake dynamics, especially at high angles of attack just before and after boundary layer separation takes place. The present code is validated in detail against the well‐known MEXICO experiment and a set of non‐rotating cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of spatial and temporal resolution of the turbulent inflow on wind turbine performance estimation

The effects of spatial and temporal resolution of wind inflows generated using large eddy simulations (LES) on the scales of turbulence present in the wind inflow, and the resulting changes in wind turbine performance were investigated for neutral atmospheric boundary layer conditions. Wind inflows with four different spatial resolutions and five different temporal resolutions were used to produce different turbine responses. An aero‐elastic code assessed the dynamic response of two wind turbines to the different inflows. Auto‐spectral density functions (ASDF) of turbine responses, such as blade deflection and bending moment, that are representative of the turbine response were used to assess the effect of the inflow. The results indicated that, as additional turbulence scales were resolved, the wind turbines showed a similar increased response that was evident in both the ASDF and variance of the different wind turbine performance parameters. As a result, the amount to which turbulence is resolved in the inflow, particularly using tools such as LES, will be important to consider when using these inflows for wind turbine design and performance prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of a LES technique to characterize the wake deflection of a wind turbine in yaw

When a wind turbine works in yaw, the wake intensity and the power production of the turbine become slightly smaller and a deflection of the wake is induced. Therefore, a good understanding of this effect would allow an active control of the yaw angle of upstream turbines to steer the wake away from downstream machines, reducing its effect on them. In wind farms where interaction between turbines is significant, it is of interest to maximize the power output from the wind farm as a whole and to reduce fatigue loads on downstream turbines due to the increase of turbulence intensity in wakes. A large eddy simulation model with particular wind boundary conditions has been used recently to simulate and characterize the turbulence generated by the presence of a wind turbine and its evolution downstream the machine. The simplified turbine is placed within an environment in which relevant flow properties like wind speed profile, turbulence intensity and the anisotropy of turbulence are found to be similar to the ones of the neutral atmosphere. In this work, the model is used to characterize the wake deflection for a range of yaw angles and thrust coefficients of the turbine. The results are compared with experimental data obtained by other authors with a particle image velocimetry technique from wind tunnel experiments. Also, a comparison with simple analytical correlations is carried out. Copyright © 2009 John Wiley & Sons, Ltd.     

垂直轴风力机三维效应的数值模拟研究

文章针对二维和三维垂直轴风力机的数值模拟的差异,提出了风力机的三维效应是造成模拟差异的主要原因。运用计算流体力学方法对某直线翼垂直轴风力机模型进行了二维和三维的数值模研究。通过比对实验得到的风力机功率系数,发现三维模拟结果与实验值吻合。观察尖速比为1.5时二维和三维垂直轴风力机的速度型分布曲线、流向速度云图和涡量云图,研究了阻塞效应、叶梢涡、支撑结构和塔架对数值模拟结果的影响。研究发现:在二维的数值模拟中,风力机没有受阻塞效应影响,功率系数被严重高估;三维的数值模拟能够模拟出全部的流畅细节,受叶梢涡和支撑结构的影响,风力机的功率系数明显降低。     

Comparison of the near‐wake between actuator‐line simulations and a simplified vortex model of a horizontal‐axis wind turbine

The flow around an isolated horizontal‐axis wind turbine is estimated by means of a new vortex code based on the Biot–Savart law with constant circulation along the blades. The results have been compared with numerical simulations where the wind turbine blades are replaced with actuator lines. Two different wind turbines have been simulated: one with constant circulation along the blades, to replicate the vortex method approximations, and the other with a realistic circulation distribution, to compare the outcomes of the vortex model with real operative wind‐turbine conditions (Tjæreborg wind turbine). The vortex model matched the numerical simulation of the turbine with constant blade circulation in terms of the near‐wake structure and local forces along the blade. The results from the Tjæreborg turbine case showed some discrepancies between the two approaches, but overall, the agreement is qualitatively good, validating the analytical method for more general conditions. The present results show that a simple vortex code is able to provide an estimation of the flow around the wind turbine similar to the actuator‐line approach but with a negligible computational effort. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling and measurements of power losses and turbulence intensity in wind turbine wakes at Middelgrunden offshore wind farm

Understanding of power losses and turbulence increase due to wind turbine wake interactions in large offshore wind farms is crucial to optimizing wind farm design. Power losses and turbulence increase due to wakes are quantified based on observations from Middelgrunden and state‐of‐the‐art models. Observed power losses due solely to wakes are approximately 10% on average. These are relatively high for a single line of wind turbines due in part to the close spacing of the wind farm. The wind farm model Wind Analysis and Application Program (WAsP) is shown to capture wake losses despite operating beyond its specifications for turbine spacing. The paper describes two methods of estimating turbulence intensity: one based on the mean and standard deviation (SD) of wind speed from the nacelle anemometer, the other from mean power output and its SD. Observations from the nacelle anemometer indicate turbulence intensity which is around 9% higher in absolute terms than those derived from the power measurements. For comparison, turbulence intensity is also derived from wind speed and SD from a meteorological mast at the same site prior to wind farm construction. Despite differences in the measurement height and period, overall agreement is better between the turbulence intensity derived from power measurements and the meteorological mast than with those derived from data from the nacelle anemometers. The turbulence in wind farm model indicates turbulence increase of the order 20% in absolute terms for flow directly along the row which is in good agreement with the observations. Copyright © 2007 John Wiley & Sons, Ltd.     

An experimental investigation on the wake interferences among wind turbines sited in aligned and staggered wind farms

An experimental investigation was conducted for a better understanding of the wake interferences among wind turbines sited in wind farms with different turbine layout designs. Two different types of inflows were generated in an atmospheric boundary layer wind tunnel to simulate the different incoming surface winds over typical onshore and offshore wind farms. In addition to quantifying the power outputs and dynamic wind loads acting on the model turbines, the characteristics of the wake flows inside the wind farms were also examined quantitatively. After adding turbines staggered between the first 2 rows of an aligned wind farm to increase the turbine number density in the wind farm, the added staggered turbines did not show a significant effect on the aeromechanical performance of the downstream turbines for the offshore case. However, for the onshore case, while the upstream staggered turbines have a beneficial effect on the power outputs of the downstream turbines, the fatigue loads acting on the downstream turbines were also found to increase considerably due to the wake effects induced by the upstream turbines. With the same turbine number density and same inflow characteristics, the wind turbines were found to be able to generate much more power when they are arranged in a staggered layout than those in an aligned layout. In addition, the characteristics of the dynamic wind loads acting on the wind turbines sited in the aligned layout, including the fluctuation amplitudes and power spectrum, were found to be significantly different from those with staggered layout.     

On wind turbine loads during the evening transition period

The late afternoon hours in the diurnal cycle precede the development of the nocturnal stable boundary layer. This “evening transition” (ET) period is often when energy demand peaks. This period also corresponds to the time of day that is a precursor to late‐afternoon downbursts, a subject of separate interest. To capture physical characteristics of wind fields in the atmospheric boundary layer (ABL) during this ET period, particularly the interplay of shear and turbulence, stochastic simulation approaches, although more tractable, are not suitable. Large‐eddy simulation (LES), on the other hand, may be used to generate high‐resolution ABL turbulent flow fields. We present a suite of idealized LES four‐dimensional flow fields that define a database representing different combinations of large‐scale atmospheric conditions (characterized by associated geostrophic winds) and surface boundary conditions (characterized by surface heat fluxes). Our objective is to evaluate the performance of wind turbines during the ET period. Accordingly, we conduct a statistical analysis of turbine‐scale wind field variables. We then employ the database of these LES‐based inflow wind fields in aeroelastic simulations of a 5‐MW wind turbine. We discuss how turbine loads change as the ET period evolves. We also discuss maximum and fatigue loads on the rotor and tower resulting from different ABL conditions. Results of this study suggest that, during the ET period, the prevailing geostrophic wind speed affects the mean and variance of longitudinal winds greatly and thus has significant influence on all loads except the yaw moment which is less sensitive to uniform and symmetric incoming flow. On the other hand, surface heat flux levels affect vertical turbulence and wind shear more and, as a result, only affect maximum blade flapwise bending and tower fore‐aft bending loads.     

Characterizing the response of a wind turbine model under complex inflow conditions

A horizontal axis wind turbine model was tested in a closed‐circuit wind tunnel under various inflow conditions. Separate experiments placed the test turbine (i) in the wake of a three‐dimensional, sinusoidal hill, (ii) in the wake of another turbine and (iii) in the turbulent boundary layer, as a reference case. Simultaneous high‐frequency measurements of the turbine output voltage, rotor angular velocity along with streamwise and wall normal velocity components were collected at various locations through the turbine's miniature direct‐current (DC) generator, a high‐resolution laser tachometer and cross‐wire anemometer, respectively. Validation trials were conducted first in order to characterize the test turbine's output and response to the baseline turbulent boundary layer. Analysis was performed by comparing the cross‐wire anemometry measurements of the incoming flow with the turbine voltage output to investigate the unsteady rotor kinematics under different flow perturbations. Using spectral, auto‐correlation and cross‐correlation methods, it was found that the flow structures developing downwind of the hill leave a stronger signature on the fluctuations and spectrum of the rotor angular velocity, as compared with those flow structures filtered or deflected by placing a turbine upwind. In summary, we show that the effects on downwind turbines of complex terrain and multi‐turbine arrangements are consistent with the induced modifications by the hill or turbine on the large scale structures in the incoming flow. Copyright © 2014 John Wiley & Sons, Ltd.     

Reynolds number dependence of turbulence statistics in the wake of wind turbines

A wind tunnel experiment has been performed to quantify the Reynolds number dependence of turbulence statistics in the wake of a model wind turbine. A wind turbine was placed in a boundary layer flow developed over a smooth surface under thermally neutral conditions. Experiments considered Reynolds numbers on the basis of the turbine rotor diameter and the velocity at hub height, ranging from Re = 1.66 × 10⁴ to 1.73 × 10⁵. Results suggest that main flow statistics (mean velocity, turbulence intensity, kinematic shear stress and velocity skewness) become independent of Reynolds number starting from Re ≈ 9.3 × 10⁴. In general, stronger Reynolds number dependence was observed in the near wake region where the flow is strongly affected by the aerodynamics of the wind turbine blades. In contrast, in the far wake region, where the boundary layer flow starts to modulate the dynamics of the wake, main statistics showed weak Reynolds dependence. These results will allow us to extrapolate wind tunnel and computational fluid dynamic simulations, which often are conducted at lower Reynolds numbers, to full‐scale conditions. In particular, these findings motivates us to improve existing parameterizations for wind turbine wakes (e.g. velocity deficit, wake expansion, turbulence intensity) under neutral conditions and the predictive capabilities of atmospheric large eddy simulation models. Copyright © 2011 John Wiley & Sons, Ltd.     

A simple atmospheric boundary layer model applied to large eddy simulations of wind turbine wakes

A simple model for including the influence of the atmospheric boundary layer in connection with large eddy simulations of wind turbine wakes is presented and validated by comparing computed results with measurements as well as with direct numerical simulations. The model is based on an immersed boundary type technique where volume forces are used to introduce wind shear and atmospheric turbulence. The application of the model for wake studies is demonstrated by combining it with the actuator line method, and predictions are compared with field measurements. Copyright © 2013 John Wiley & Sons, Ltd.     

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.