Near wake Reynolds‐averaged Navier–Stokes predictions of the wake behind the MEXICO rotor in axial and yawed flow conditions

In the present paper, Reynolds‐averaged Navier–Stokes predictions of the flow field around the MEXICO rotor in yawed conditions are compared with measurements. The paper illustrates the high degree of qualitative and quantitative agreement that can be obtained for this highly unsteady flow situation, by comparing measured and computed velocity profiles for all three Cartesian velocity components along four axial transects and several radial transects. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of different forest densities on atmospheric boundary‐layer development and wind‐turbine wake

The aim of this work is to investigate the atmospheric boundary‐layer (ABL) flow and the wind turbine wake over forests with varying leaf area densities (LAD). The forest LAD profile used in this study is based on a real forest site, Ryningsnäs, located in Sweden. The reference turbine used to model the wake is a well‐documented 5‐MW turbine, which is implemented in the simulations using an actuator line model (ALM). All simulations are carried out with openFOAM using the Reynolds averaged Navier‐Stokes (RANS) approach. Twelve forest cases with leaf area index (LAI) ranging from 0.42 to 8.5 are considered. Results show that the mean velocity decreases with increasing LAI within the forest canopy, but increases with LAI above the hub height. Meanwhile, the turbulent kinetic energy (TKE) varies nonmonotonically with forest density. The TKE increases with forest density and reaches to its maximum at an average LAI of 1.70, afterwards, it decreases gradually as the density increases. It is also observed that the forest density has a clear role in the wake development and recovery. Comparisons between no‐forest and forest cases show that the forest characteristics help in damping the added turbulence from the turbine. As a consequence, the forest with the highest upstream turbulence has the shortest wake downstream of the turbine.     

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.     

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.     

Tip‐vortex breakdown of wind turbines subject to shear

Sheared velocity profiles pervade all wind‐turbine applications, thus making it important to understand their effect on the wake. In this study, a single wind turbine is modeled using the actuator‐line method in the incompressible Navier–Stokes equations. The tip vortices are perturbed harmonically, and the growth rate of the response is evaluated under uniform inflow and a linear velocity profile. Whereas previous investigations of this kind were conducted in the rotating frame of reference, this study evaluates the excitation response in the fixed frame of reference, thus necessitating a frequency transformation. It is shown that increasing the shear decreases the spatial growth rate in the upper half of the wake while increasing it in the lower half. When scaled with the local tip vortex parameters, the growth rate along the entire azimuth collapses to a single value for the investigated wavenumbers. We conclude that even though the tip‐vortex breakdown is asymmetric in sheared flow, the scaled growth rates follow the behavior of axisymmetric helical vortices. An excitation amplitude reduction by an order of magnitude extends the linear growth region of the wake by one radius for uniform inflow. In the sheared setup, the linear growth region is extended further in the top half than in the bottom half because of the progressive distortion of the helical tip vortices. An existing model to determine the stable wake length was shown to be in close agreement with the observed numerical results when adjusted for shear.     

Large‐eddy simulations of wind‐farm wake characteristics associated with a low‐level jet

In this study, we performed a suite of flow simulations for a 12‐wind‐turbine array with varying inflow conditions and lateral spacings, and compared the impacts of the flow on velocity deficit and wake recovery. We imposed both laminar inflow and turbulent inflows, which contain turbulence for the Ekman layer and a low‐level jet (LLJ) in the stable boundary layer. To solve the flow through the wind turbines and their wakes, we used a large‐eddy simulation technique with an actuator‐line method. We compared the time series for the velocity deficit at the first and rear columns to observe the temporal change in velocity deficit for the entire wind farm. The velocity deficit at the first column for LLJ inflow was similar to that for laminar inflow. However, the magnitude of velocity deficit at the rear columns for the case with LLJ inflow was 11.9% greater because of strong wake recovery, which was enhanced by the vertical flux of kinetic energy associated with the LLJ. To observe the spatial transition and characteristics of wake recovery, we performed statistical analyses of the velocity at different locations for both the laminar and LLJ inflows. These studies indicated that strong wake recovery was present, and a kurtosis analysis showed that the probability density function for the streamwise velocity followed a Gaussian distribution. In a quadrant analysis of the Reynolds stress, we found that the ejection and sweep motions for the LLJ inflow case were greater than those for the laminar inflow case.     

A numerical study of tilt‐based wake steering using a hybrid free‐wake method

This study investigates the potential of using tilt‐based wake steering to alleviate wake shielding problems experienced by downwind turbines. Numerical simulations of turbine wakes have been conducted using a hybrid free‐wake analysis combining vortex lattice method (VLM) and an innovative free‐wake model called constant circulation contour method (CCCM). Simulation results indicate tilting a horizontal axis wind turbine's shaft upward causes its wake to ascend, carrying energy‐depleted air upward and pumping more energetic replacement air into downstream turbines, thereby having the potential to recover downstream turbine power generation. Wake cross section vorticity and velocity distributions reveal that the wake upward transport is caused by the formation of near‐wake streamwise vorticity components, and furthermore, the wake velocity deficit is weakened because of the skewed wake structure. Beyond the single turbine wake simulation, an inline two‐turbine case is performed as an assessment of the wake steering influence on the two‐turbine system and as an exploratory work of simulating turbine‐wake interactions using the hybrid free‐wake model. Individual and total turbine powers are calculated. A comparison between different tilting angles suggests turbine power enhancement may be achieved by tilting the upstream turbine and steering its wakes away from the downstream turbine.     

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.     

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.     

CFD study on the impact of yawed inflow on loads,power and near wake of a generic wind turbine

This article deals with the influence of yawed inflow conditions on the performance of a single generic 2.4MW wind turbine. It presents the results of studies performed at the Institute of Aerodynamics and Gas Dynamics by means of computational fluid dynamics , using a fully meshed wind turbine with all boundary layers being resolved. The block‐structured flow solver FLOWer is used; a dual‐time stepping method for temporal discretization and a second‐order Jameson–Schmidt–Turkel method for the calculation of the convective fluxes are applied. All simulations are carried out using a detached eddy simulation approach. In detail, two different wind speeds and a yaw angle range between ?50° and +50° are evaluated in the paper. Based on these data, it is shown that the reduction of power output follows a cosine to the power of X function of the yaw angle. Furthermore, the growing azimuthal non‐uniformity of the load distributions with increasing yaw angle magnitude is analysed by spanwise load distributions. As a central influence on the load distributions, the advancing and retreating blade effect is identified. Moreover, the deflection of the wake as a result of the inflow is investigated, and the deflection angles are compared with a modelling approach. A connection line between wake deflection and load asymmetry is drawn. The results are of particular importance for wind park situations with downstream turbines facing the distorted inflow created from upstream ones. Copyright © 2016 John Wiley & Sons, Ltd.     

Near‐wake flow simulation of a vertical axis turbine using an actuator line model

In the present work, the near‐wake generated for a vertical axis wind turbine (VAWT) was simulated using an actuator line model (ALM) in order to validate and evaluate its accuracy. The sensitivity of the model to the variation of the spatial and temporal discretization was studied and showed a bigger response to the variation in the mesh size as compared with the temporal discretization. The large eddy simulation (LES) approach was used to predict the turbulence effects. The performance of Smagorinsky, dynamic k‐equation, and dynamic Lagrangian turbulence models was tested, showing very little relevant differences between them. Generally, predicted results agree well with experimental data for velocity and vorticity fields in representative sections. The presented ALM was able to characterize the main phenomena involved in the flow pattern using a relatively low computational cost without stability concerns, identified the general wake structure (qualitatively and quantitatively), and the contribution from the blade tips and motion on it. Additionally, the effects of the tower and struts were investigated with respect to the overall structure of the wake, showing no significant modification. Similarities and discrepancies between numerical and experimental results are discussed. The obtained results from the various simulations carried out here can be used as a practical reference guideline for choosing parameters in VAWTs simulations using the ALM.     

Measurements on a wind turbine wake: 3D effects and bluff body vortex shedding

The velocity field in the wake of a two‐bladed wind turbine model (diameter 180 mm) has been studied under different conditions using a two‐component hot wire. All three velocity components were measured both for the turbine rotor normal to the oncoming flow as well as with the turbine inclined to the freestream direction (the yaw angle was varied from 0° to 20°). The measurements showed, as expected, a wake rotation in the opposite direction to that of the turbine. A yawed turbine is found to clearly deflect the wake flow to the side, showing the potential of controlling the wake by yawing the turbine. An unexpected feature of the flow was that spectra from the time signals showed the appearance of a low‐frequency fluctuation both in the wake and in the flow outside the wake. This fluctuation was found both with and without freestream turbulence and also with a yawed turbine. The frequency expressed as a Strouhal number was shown to be independent of the freestream velocity or turbulence level, but the low frequency was only observed when the tip speed ratio (or equivalently the drag coefficient) was high. The shedding frequency changed also with the yaw angle. This is in agreement with the idea that the turbine sheds structures as a bluff body. The phenomenon, noticeable in all the velocity components, was further investigated using two‐point cross‐correlations of the velocity signals. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigation of wind turbine wake superposition models using Reynolds‐averaged Navier‐Stokes simulations

It is well accepted that the wakes created by upstream turbines significantly impact on the power production and fatigue loading of downstream turbines and that this phenomenon affects wind farm performance. Improving the understanding of wake effects and overall efficiency is critical for the optimisation of layout and operation of increasingly large wind farms. In the present work, the NREL 5‐MW reference turbine was simulated using blade element embedded Reynolds‐averaged Navier‐Stokes computations in sheared onset flow at three spatial configurations of two turbines at and above rated flow speed to evaluate the effects of wakes on turbine performance and subsequent wake development. Wake recovery downstream of the rearward turbine was enhanced due to the increased turbulence intensity in the wake, although in cases where the downstream turbine was laterally offset from the upstream turbine this resulted in relatively slower recovery. Three widely used wake superposition models were evaluated and compared with the simulated flow‐field data. It was found that when the freestream hub‐height flow speed was at the rated flow speed, the best performing wake superposition model varied depending according to the turbine array layout. However, above rated flow speed where the wake recovery distance is reduced, it was found that linear superposition of single turbine velocity deficits was the best performing model for all three spatial layouts studied.     

Numerical study of the wake induced by a porous disk under deflected inflow

This study examines the wake of a porous disk that generates a velocity deficit equivalent to that of a wind turbine. Three-dimensional unsteady numerical simulations based on the finite volume method are performed. The URANS-SST (k–ω) model is applied for the turbulence closure. Two investigations are carried out in this study: (i) the influence of the disk porosity on the wake, for porosities values (p) ranging from 0 to 0.55 in the case of a perpendicular flow; and (ii) the influence of the yaw angle on the wake deviation, for yaw angles ranging from 5° to 30°. Good agreements with the available experimental data are obtained for the mean x-velocity component. The results confirm that wake length increases as porosity decreases. For nonporous disks, most part of the fluid is deflected toward the mast and above the disk. The y-velocity contours highlight two contra-rotating vortices in the vicinity of the disk. In both cases (nonporous and porous disks), a high turbulent kinetic energy is obtained near the disk area, with a higher maximum value for the nonporous disk.     

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.     

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.     

A generalized framework for reduced‐order modeling of a wind turbine wake

A reduced‐order model for a wind turbine wake is sought from large eddy simulation data. Fluctuating velocity fields are combined in the correlation tensor to form the kernel of the proper orthogonal decomposition (POD). Proper orthogonal decomposition modes resulting from the decomposition represent the spatially coherent turbulence structures in the wind turbine wake; eigenvalues delineate the relative amount of turbulent kinetic energy associated with each mode. Back‐projecting the POD modes onto the velocity snapshots produces dynamic coefficients that express the amplitude of each mode in time. A reduced‐order model of the wind turbine wake (wakeROM) is defined through a series of polynomial parameters that quantify mode interaction and the evolution of each POD mode coefficients. The resulting system of ordinary differential equations models the wind turbine wake composed only of the large‐scale turbulent dynamics identified by the POD. Tikhonov regularization is used to recalibrate the dynamical system by adding additional constraints to the minimization seeking polynomial parameters, reducing error in the modeled mode coefficients. The wakeROM is periodically reinitialized with new initial conditions found by relating the incoming turbulent velocity to the POD mode coefficients through a series of open‐loop transfer functions. The wakeROM reproduces mode coefficients to within 25.2%, quantified through the normalized root‐mean‐square error.  A high‐level view of the modeling approach is provided as a platform to discuss promising research directions, alternate processes that could benefit stability and efficiency, and desired extensions of the wakeROM.     

Flow characterization in the near‐wake region of a horizontal axis wind turbine

An experimental study is conducted to investigate the flow dynamics within the near‐wake region of a horizontal axis wind turbine using particle image velocimetry (PIV). Measurements were performed in the horizontal plane in a row of four radially distributed measurement windows (tiles), which are then patched together to obtain larger measurement field. The mean and turbulent components of the flow field were measured at various blade phase angles. The mean velocity and turbulence characteristics show high dependency on the blade phase angle in the near‐wake region closer to the blade tip and become phase independent further downstream at a distance of about one rotor diameter. In the near‐wake region, both the mean and turbulent characteristics show a systemic variation with the phase angle in the blade tip region, where the highest levels of turbulence are observed. The streamlines of the instantaneous velocity field at a given phase allowed to track a tip vortex which showed wandering trend. The tip vortices are mostly formed at r/R > 1, which indicates the wake expansion. Results also show the gradual movement of the vortex region in the axial direction, which can be attributed to the dynamics of the helical tip vortices which after being generated from the tip, rotate with respect to the blade and move in the axial direction because of the axial momentum of the flow. The axial velocity deficit was compared with other laboratory and field measurements. The comparison shows qualitative similarity. Copyright © 2015 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.     

A wind‐tunnel study of the wake development behind wind turbines over sinusoidal hills

In the present work, the wake development behind small‐scale wind turbines is studied when introducing local topography variations consisting of a series of sinusoidal hills. Additionally, wind‐tunnel tests with homogeneous and sheared turbulent inflows were performed to understand how shear and ambient turbulence influence the results. The scale of the wind‐turbine models was about 1000 times smaller than full‐size turbines, suggesting that the present results should only be qualitatively extrapolated to real‐field scenarios. Wind‐tunnel measurements were made by means of stereoscopic particle image velocimetry to characterize the flow velocity in planes perpendicular to the flow direction. Over flat terrain, the wind‐turbine wake was seen to slowly approach the ground while it propagated downstream. When introducing hilly terrain, the downward wake deflection was enhanced in response to flow variations induced by the hills, and the turbulent kinetic energy content in the wake increased because of the speed‐up seen over the hills. The combined wake observed behind 2 streamwise aligned turbines was more diffused and when introducing hills, it was more prone to deflect towards the ground compared to the wake behind an isolated turbine. Since wake interactions are common at sites with multiple turbines, this suggested that it is important to consider the local hill‐induced velocity variations when onshore wind farms are analysed. Differences in the flow fields were seen when introducing either homogeneous or sheared turbulent inflow conditions, emphasizing the importance of accounting for the prevailing turbulence conditions at a given wind‐farm site to accurately capture the downstream wake development.