Open access peer-reviewed chapter

# A Review of Wind Turbine Yaw Aerodynamics

Written By

Daniel Micallef and Tonio Sant

Submitted: October 29th, 2015 Reviewed: April 4th, 2016 Published: July 27th, 2016

DOI: 10.5772/63445

From the Edited Volume

## Wind Turbines

Edited by Abdel Ghani Aissaoui and Ahmed Tahour

Chapter metrics overview

View Full Metrics

## Abstract

The fundamental physics of HAWT aerodynamics in yaw is reviewed with reference to some of the latest scientific research covering both measurements and numerical modelling. The purpose of this chapter is to enable a concise overview of this important subject in rotor aerodynamics. This will provide the student, researcher or industry professional a quick reference. Detailed references are included for those who need to delve deeper into the subject. The chapter is also restricted to the aerodynamics of single rotors and their wake characteristics. Far wake and wind turbine to turbine effects experienced in wind farms are excluded from this review. Finally, a future outlook is provided in order to inspire further research in yawed aerodynamics.

### Keywords

• horizontal axis wind turbine
• yaw
• wakes
• rotors

## 1. Introduction: The relevance of wind turbine yaw aerodynamics

### 1.1. Wind turbine yaw fundamentals

A wind turbine converts the kinetic energy from the wind to electrical energy. The power which can be extracted is proportional to the cube of the wind speed. There is, however, a theoretical limit to the power which can be extracted from the turbine. This is known as the Betz limit. This limit was derived by Betz [1] to correspond to 59% of the maximum available power which can be extracted by the turbine. The efficiency of a wind turbine is called the power coefficient and is defined by the following:

CP=P12ρU3AE1

where Pis the turbine power, ρis the density, Uis the wind velocity and Ais the rotor area. In practice, due to various practical reasons, only a fraction of the Betz limit can be exploited. Some reasons for this reduction in the overall efficiency of the turbines are mostly associated with the wind resource itself and include wind shear [2, 3], wind turbulence [4, 5], and yaw. In this chapter, we focus our attention on the latter for the case of a horizontal axis wind turbine (HAWT).

Yaw occurs when the wind direction is not perpendicular to the rotor plane. When the HAWT operates in yaw, the average power extracted by the turbine reduces as compared to the case when the wind is perpendicular to the rotor plane. This yawed flow situation is depicted in Figure 1 where γ is the yaw angle. The blade will experience a varying relative velocity and angle of attack with azimuthal blade position, leading to an unsteady aerodynamics problem.

### 1.2. Relevance of studying yawed flows

Although modern wind turbines allow for various yaw mechanisms in order to align the turbine with respect to the incoming wind, the response to variations in wind direction is very slow. Consequently, the turbine will be in yaw during most of its operational time. Eggleston and Starcher [6] and Madsen et al. [7] show time series results of wind direction as observed from a sensor installed on two different experimental turbines in the field. A time series for the wind yaw angle measured from the experimental turbines is shown in Figure 2. For the Tellus turbine (right figure), during a particular instant in time, the yaw angle increases and varies dramatically. The yaw angle can reach up to around 60°, depending on the wind flow characteristic on site.

While the yawed flow scenario presents an important and challenging problem, the basis of understanding the wind turbine aerodynamics is the axial flow case. The major challenges in the understanding of wind turbine flows are associated with the root and tip flows [810]. For an in‐depth review of wind turbine aerodynamics, particularly for wind turbines operating in axial flow, the reader is referred to [11] and lately in [12].

Unsteady aerodynamics inevitably results in unsteady loads which negatively affect the power quality and fatigue lifetime of wind turbines. Reference [13] provides a detailed overview of the various challenges in modelling wind turbines under unsteady conditions, including rotor yaw, given the lack of understanding of many aerodynamic phenomena associated with the wake structure developed behind the turbine and the dynamics of unsteady flows over the blade sections.

Wind turbines operating in wind farms are well known to suffer from significant aerodynamic interference losses, primarily due to the axial velocity deficit incurred by the wind flowing through the upstream turbines. The power output from full‐scale turbines operating in large wind farms may be as low as 40% of a stand‐alone turbine, depending on the turbine spacing as well as environmental flow conditions including the mean wind speed, turbulence intensity and atmospheric stability (see [14]). Yaw aerodynamics may eventually become more important for large offshore wind farm design, given the possibility of mitigating such wake losses by skewing the upstream turbine wake using smart rotor yaw control.

In this chapter, a review of the physics of HAWT aerodynamics in yaw is given by focusing on the aerodynamic phenomena of the near wake and the flow over the rotating blade sections. Far wake and turbine‐to‐turbine interference effects are excluded from this study. Furthermore, emphasis is made on single HAWT rotors operating in fixed yaw and under steady and uniform wind fronts. Reference is made to some of the latest scientific research in the field, including state‐of‐the‐art wind tunnel measurements undertaken on model turbines under controlled conditions.

## 2. The blade element momentum (BEM) theory

The momentum conservation principle can be applied in order to be able to determine the wake inductions, and hence, the turbine blade loading and performance. In the blade element momentum (BEM) method, the blade is divided into a number of 2D elements. The momentum principle is then applied on aerodynamically independent annular elements. In doing so, a number of limitations and assumptions arise:

1. Inviscid flow.

2. Annular air elements are independent (no radial flow from one element to the next).

3. Aerofoil polars are generally based on 2D wind tunnel measurements.

4. The loading of the blades on the air is assumed not to vary with azimuth meaning that an infinite number of blades assumed. This must be resolved using tip loss correction methods.

5. The theory alone cannot be used for highly loaded rotors, where turbulent mixing is prevalent.

The velocity diagram for a particular blade element is shown in Figure 3. The relative velocity is given by the following:

Vrel=(U+Vaxial,ind)2+(r+Vθ,ind)2E2

where Vaxial,indis the axial wake‐induced velocity (which would be opposing the freestream velocity), and Vθ,indis the tangential wake‐induced velocity. The inflow angle ϕ(indicated in Figure 3) can thus be found from:

tanφ=U+Vaxial,indr+Vθ,indE3

from which the angle of attack becomes the difference between the inflow angle and the blade's pitch and twist angles (the sum of which is here represented by θ).

=φθE4

From the aerofoil data, the normal and tangential force coefficients (with respect to the rotor plane) can be determined from the lift and drag coefficients (Cland Cd) and the inflow angle ϕas follows:

Cn=Clcosφ+CdsinφE5
Ct=ClsinφCdcosφE6

The induction factors defined as a1=Vaxial,ind/Uand a2 = Vθ,ind/rΩmay then be determined applying the law of conservation of momentum in the axial and azimuthal directions by assuming that the momentum expressions are equal to the blade element expressions:

4a1(1a1)U2=σVrel2CnE7
4a2(1a1)rUr=σVrel2CtrE8

where a1and a2are the axial and tangential induction factors. σis defined as the rotor local solidity:

σ=Bc2πrE9

where cis the chord and Bis the number of blades. This formulation does not include the effect of the finite number of blades. The most commonly used approach to deal with this problem is to use the Prandtl tip loss factor reported [1]. In this formulation, Prandtl uses a disc representation to represent the wake effect and the finite number of blades. This is shown in Figure 4.

Prandtl came up with an analytical and implementable solution from a tip loss factor Fof the following form:

F=2πcos1(exp(B(Rr)2R(R)2+U2(1a1)2U(1a1)))E10

One way to account for tip loss is to modify Eqs. (13) and (14) to include the tip loss factor and re‐arrange them for an iterative solution, as follows:

a1=14Fsin2φσCn+1E11
a2=14FsinφcosφσCt1E12

Various authors, such as in [15], assessed the validity of this model using a computational fluid dynamics (CFD) approach. While the accuracy of the Prandtl tip loss model is debatable, it nonetheless provides a simple solution to a complex problem and is still to date the most common method for accounting for a finite number of blades.

The algorithm for implementing the BEM theory can be found in various textbooks including [1518]. Here, the one given in [19] is adapted:

1. Set a1 = a2= 0.

2. Find the flow angle ϕand hence the angle of attack α.

3. Calculate the tip loss factor from Eq. (10).

4. Read Cland Cdfrom the 2D aerofoil polars.

5. The normal and tangential thrust coefficients Cnand Ctcan be calculated from Eqs. (5) and (6).

6. New values for a1and a2can be established from Eqs. (17) and (18).

7. If a1and a2have converged below a defined tolerance, then move to the next step. Otherwise, repeat starting from step 2.

8. Compute the thrust and torque for the wind turbine rotor using the following:

T=BRrR4πρU2a1(1a1)rdrE13
Q=BRrR4πρUra2(1a1)r2drE14

## 3. The BEM theory for yawed wind turbines

### 3.1. Modifications to the axial momentum theory

The linear momentum theory is modified for yawed wind turbines on the basis of Glauert's autogyro theory (see Glauert [20]). Only the normal wind flow component U cos γis assumed to be affected by the presence of the actuator disc. It is assumed that this component changes to U cos γ+uat the disc and the resultant flow velocity here becomes

u'=U2sin2γ+U2(cosγa)2E15

Glauert expresses the momentum equation for the axial thrust Ton the rotor as

T=2ρAaUu'E16

The thrust parallel to the plane of the disc is zero as it is assumed that the flow component Usinγis left unperturbed by the rotor.

In yawed conditions, the wind speed and inflow relative to the moving blades change with azimuth angle. The velocity diagram for the blade at the top‐most position is shown in Figure 5. The resultant velocities would hence be the following:

Vaxial=UcosγVaxial,indE17
Vθ=r+Vθ,indUsinγcosE18

Anderson showed, using a BEM approach that the power coefficient of a yawed turbine reduces with increasing yaw [16]. This hypothesis was later confirmed in [17, 18]. The results obtained by [18] for a 2‐ and 3‐bladed rotor are compared with those found in [16, 17] in Figure 6. CPis the power coefficient at a given yaw angle, while CPmaxis that under axial conditions. The rotors which were used in these studies are different, and hence, a qualitative agreement was only obtained.

The ratio CP/CPmaxrepresenting the loss in power due to rotor yaw may be approximated to cosx(γ). The exponent xis often thought to be equal to 3; however, this is only valid if the axial induction distribution is small compared to U cos γor the axial induction factor distribution is unchanged as a result of yaw. Recent wind tunnel measurements on the 0.9 m diameter 3‐bladed model at NTNU have shown that the CP/CPmaxfollowed well the cos3(γ) relation (see [21]). A different situation was observed in another measurement. In [22], it is reported that measurements involving both wind tunnel models and full‐scale turbines tested at Vattenfall/FFA and DNW have indicated that the exponent xmay vary between 1.8 and 5 for the NREL Phase VI rotor (see [23] pp. 181–183). Such studies have also consistently shown that the BEM theory over predicts the power at yawed conditions. This is possibly related to the limitations of the linear momentum equation modified for yawed flows [Eq. (16)].

Although yawed flow contradicts the fundamental assumption used in BEM of having radially independent annular elements, BEM models have still been used in yaw (refer to [24]). In order to correct for the asymmetric wake induction, a model was first developed in [25] in his work on the autogyro. Ignoring this correction will result in a zero restoring yawing moment, which has been shown experimentally to be unrealistic. Glauert proposed that the induction factor at the blades would take the following form:

a=a1(1+KrRsinϕ)E19

where a1is the average induction factor, Kis a function of yaw angle, ris the radial station in consideration, Ris the rotor radius and ϕ is the azimuthal position of the rotor blade. Much of the efforts of the early experiments in yawed flow HAWTs were aimed at understanding the wake expansion and deflection. This enabled various improvements in the classical BEM techniques by improving the parameter Kin the Glauert expression for the azimuthally varying induction [see Eq. (19)]. There are various models of the function Ksome of which are given here.

Coleman et al. [21]:

K=tan(2)E20

White and Blake [26]:

K=2sin()E21

Pitts and Peters [27]:

K=15π32tan(2)E22

Howlett [28]:

K=sin2()E23

where χis the skew angle; the angle between the wake and rotor axis. In the BEM formulation, this angle is defined as follows:

tanχ=(U sinγ)(U cosγua)E24

The Kfactors derived from the different models listed earlier are plotted in Figure 7 to compare the differences. As can be seen, there are dramatic discrepancies between the results of various authors. Since this function will ultimately affect the varying inductions on the plane, knowledge of the induction field is also necessary to create sound models.

Other engineering models were developed in the past years, some of which were examined against measurements in the EU JOULE Dynamic Inflow Projects (see [29, 30]). These models have been implemented in various BEM‐based aeroelastic codes. One such model is that developed by Øye [31] and is similar to Glauert's model but introduces a radial variation in induction. The model was derived through a curve fitting procedure with results from an actuator disc model. The induced velocity distribution was found to depend on the radial location, azimuth angle and the wake skew angle:

a=a1(1+f2,tudktan 2sin)E25
f2,tudk=rR+0.4(rR)3+0.4(rR)5E26

The Glauert correction model [Eq. (19)] gives acceptable results towards the tip region of the blade where the sinusoidal variation of the induction is prevalent. However, the root region of the blade does not exhibit sinusoidal behaviour as was found in the experiments in [32]. Based on this same experiment, Schepers used Fourier series expansions to establish a model for yaw which models better the radial dependency of the induction variation with azimuth (see [33]). This is shown here for a particular radial station along the blade:

a=a1(1A1cos(φ1)A2cos(2φ2))E27

where the amplitudes A1and A2and the phases ϕ1and ϕ2are curve fitting parameters, all functions of radial position. Also, Ref. [34] (pp. 403–405) recommends the inclusion of higher harmonics in BEM yaw modelling.

As discussed earlier, yaw results in an asymmetric flow field where the relative velocity magnitude and direction (and hence angle of attack) both vary with time. Yawed HAWTs are for this reason naturally susceptible unsteady aerodynamic phenomena. Such phenomena give rise to hysteresis effects, even at small angles of attack for which the flow is fully attached. At high angles of attack however, dynamic stall leads to more pronounced hysteresis on both the lift and drag, with the resulting loads often exceeding that static loads significantly. A description of different dynamic stall models applied to wind turbines is presented in [35].

Dynamic stall on a rotating blade is more complex than that experienced in 2D flow conditions. Schreck (see [36, 37]) investigated the 3D dynamic stall processes on the NREL UAE Phase VI rotor through the analysis of blade surface pressure data and the local inflow angle when operating under yawed conditions. In the EU project ‘Dynamic stall and 3D effects’ coordinated by FFA, it was attempted to use the IEA Annex XVIII measurement to understand dynamic stall effects. An example is given in Figure 8 for non‐yawed conditions in the field. They show measurements of normal force coefficients at 68% span as function of angle of attack (around a low and a high angle of attack) on the Risø experiments, carried out on a 3‐bladed, 19 m diameter turbine (see [38]). The angle of attack is derived from a local inflow angle as measured with a pitot probe and corrected for upwash. Hysteresis effects are visible at high angles of attack but in a very disorderly pattern.

The reduced frequency kcan be estimated as follows if the freestream velocity is assumed negligible compared to the rotational speed:

k=c2Vrelc2rc2rE28

where cis the chord length, ris the radial location along the blade, Ω is the rotational speed, and Vrelis the relative velocity to the blade section. In situations where k < 0.05 can be considered steady or quasi‐steady, values of k = 0.1 are typical for wind turbines in yaw. It is clear that this measure of unsteadiness is a function of the c/rratio.

In [7], various models are used including a BEM approach (HAWC), an Navier‐Stokes solver (EllipSys3D) and an actuator disc solver (HAWC3D) to find the local inflow angle (angle between local relative velocity and chord line) and angle of attack as well as the magnitude of the relative velocity and how it changes. Results are shown for the 91% station on the NREL UAE rotor in Figure 9. EllipSys3D results were very accurate as this model includes the effect of the wing geometry in its entirety, and hence, the local inflow angle was calculated. The other results, however, show quite a large difference since the models calculate the actual angle of attack. Despite the differences, the variation of angle of attack with azimuthal position is clear and leads to the unsteadiness that the yawed HAWT suffers from.

Lately, the loads on a yawed wind turbine have been established directly through stereo particle image (SPIV) measurements (see [39]) and reported in [40, 41], but these are limited to only one blade position, thus not allowing for establishing the time variation of the blade loads in yaw.

## 4. Wake characteristics

Experimental research has clearly demonstrated that in reality the wake of a wind turbine is by far more complex than that assumed in the BEM theory. It actually comprises helical vortex sheets, with one sheet originating from each individual rotating blade as illustrated in Figure 10. The circulation distribution in the vortex sheet originates from the bound circulation (ΓB) developed at the blades. The wake circulation is composed of two vector components: trailing circulation (ΓT) that is released from the blades in a direction normal to the trailing edges of the blades and is related to the spanwise variation of bound circulation (B/dr); and shed circulation (ΓS) that is released from the blades in a direction parallel to the blades’ trailing edges and is related to the variation of bound circulation with time (B/dt). The geometry of the vortex sheet emerging from each individual blade will change such that the edges roll‐up (similar to that observed on a wind in linear flight) to form a strong tip vortex and root vortex.

The helical wake of the turbine expands as a consequence of the slowdown of moving air. The higher the operating thrust coefficient, the higher the wake expansion. This has been shown by various researchers such as in [42, 43]. Tip vortices are easily noticeable in experiments on model rotors, even beyond one rotor diameter downstream. Smoke visualizations of tip vortices are shown in Figure 11. The root vortices are more difficult to track in experiments and have been only observed much closer to the rotor plane, as, for example, in the recent PIV measurements undertaken [39, 44]. It is most probable that the root vortex diffuses rapidly due to the interference effects of the turbine support structure.

In a yawed rotor, the blade sections have a geometric velocity component in a direction of the free wind speed resulting in the advancing and retreating effect. A variation in both the relative wind speed and angle of attack is experienced with the highest relative velocity reached by the advancing blade when it is vertical pointing upwards. Yawed conditions cause the wake to become skewed leading to an uneven distribution of induction at the rotor disc (see Figure 12). Measurements have shown that yawed wind turbines experience a restoring yawing moment. As noted in [23] (pp. 49–53), this moment is not generated by the advancing and retreating blade effect but by the skewed wake geometry which causes the trailing tip vorticity to be on average closer to the downwind side of the rotor plane.

Induced velocities on the upwind side or the rotor are hence lower, leading to large axial thrust loads on the blade sections here to develop a restoring yawing moment.

The influence of the skewed wake geometry on the induction at the rotor plane of yawed turbines was extensively studied in the EU JOULE Dynamic Inflow Projects (refer to Snel and Schepers [29, 30]). The implementation of advanced aerodynamic models based on vortex methods by Voutsinas et al. [45] and the acceleration potential technique by van Bussel [46, 47] has clearly indicated that, in addition to tip vorticity, root vorticity also has a considerable influence on the induction distribution at the rotor plane. However, root vorticity contributes to the generation of a destabilising yawing moment, hence opposing that induced by tip vorticity.

The axial thrust of a rotor with a finite yaw angle varies with azimuth and as a result, the expansion will vary along the azimuthal direction. This has been shown again by Grant et al. [48] and Haans et al. [49]. Some results are shown in Figure 13. As can be seen, there is a high expansion in the region where the rotor tip is pointing downwind and almost no expansion in the region where the rotor tip is pointing upwind. These results are coherent with the findings of [48]. The higher the yaw angle, the smaller the thrust coefficient, and hence, the smaller the expansion. Unfortunately, no empirical models have been proposed for the wake deformation in yaw. This would enable practical application in the design of wind turbine blades with applications in BEM or prescribed wake models. Haans et al. [49] have also derived the wake skew angle χdirectly from the smoke visualizations of tip vortices at yaw angles equal to 15°, 30°, and 45°. This angle was found to be larger than the yaw angle at all tested blade pitch angles and tip speed ratios. It was observed that for each tested yaw angle, χincreased linearly with the measured axial thrust coefficient CT(see Figure 14).

The analysis of wake development requires the use of velocity field measurements over a particular area. Due to this, techniques such as PIV must be used. Such measurements have been performed by [50] and in the MEXICO experiment (see [51]). More recent PIV measurements include those undertaken at TUDelft [39, 52], at NTNU [53] and at Monash University (see [44]). In the experiment by [50], PIV measurements were obtained at various blade azimuthal positions. With these measurements, the variation of the tip vortex circulation was calculated for yaw angles of +30° and -30° yaw as well as for the axial case. Figure 11 from [50] shows various vorticity contour plots of the tip vortex at various blade azimuth phase angles and at various yaw angles. Information on the magnitude of vorticity is not present so only insight on the iso‐vorticity lines can be gained. Soon after release, a ‘vorticity tail’ becomes apparent from the tip vortex for all cases particularly for yaw. The authors attribute the additional vorticity in these cases to the presence of additional shed vorticity in the wake of a yawed turbine due to the variation of bound vorticity with time. For the -30° yaw, the ‘vorticity tail’ is first located downstream since the measurement plane is in the upwind direction. For the +30° yaw, the ‘vorticity tail’ is pointing upwind since the measurement plane is in the downwind direction. As can be noticed, for the +30° yaw case, the vortex remains small. The trailing vorticity from the vortex sheet detaches itself and does not continue to feed into the tip vortex. This is the reason for the low circulation levels found in Figure 11 of [50]. On the other hand for the -30° yaw case, the interaction between the trailing vorticity and the tip vortex is more complex and is persistent up to 156° of azimuth blade position. The vortex becomes larger, and the circulation is similar to that found in the axial flow case.

In [54], the authors investigated the radial flows of the MEXICO rotor under both axial and yawed conditions. The investigation was based on both the experimental SPIV measurements and numerical analysis using a 3D unsteady free‐wake vortex model. The study found that the radial velocity increased by a factor of 2 when the turbine was yawed due to the in‐plane freestream component. Under axial conditions, it was observed that a slight inboard convection of the tip vortex occurs. This phenomenon was also confirmed later on the more detailed SPIV measurements by Micallef on the TUDelft rotor (see [39]). In these more recent measurements, the tip vortex evolution could be clearly traced as a function of rotor azimuth angle. It could be noted that under yawed conditions, outboard motion of the tip vortex was delayed to a blade azimuth of 10°, even though a radial component equal to Usinγwas present.

## 5. Summary of research on HAWT rotors in yaw

A comprehensive summary of various experiments on HAWT rotors under both controlled and open field conditions is given in Table 1. Table 2 includes references to various works covering the different aspects of HAWT aerodynamics in yaw discussed in this study.

Contributor Year Type Rotor diameter BladeVariablesInformation
Clayton and Filby, University College London1982HWA in the wake (for averages
velocities and turbulence
measurements) and power output
0.5 m (test section size not available)3‐bladed,
no details of the chord, twisted
Tip speed ratio, yaw angle, pitch angle1.4 ≤ λ≤ 9.8
0° ≤ γ≤ 40°
-15° ≤ θ≤ 15°
Constant chord, Twisted
N/AField experiment
Vermeer,
TUDelft
1987HWA2.24 m2‐bladed, constant chord, twisted, NACA0012Yaw angleMax Re = 2:4e5
Bruining,
TUDelft
measurements
untapered,untwisted
Tip speed ratio,
yaw angle
4 < λ < 10,
0:5e6 < Re < 1e6
Grant,
Heriot‐Watt
University
1997Smoke visualization
and power output
0.9 m in open
jet tunnel
2 and 3 blades, taperedYaw angleλ = 6:7
5:9e4 < Re < 1e5
Grant,
Heriot‐Watt
University
2000PIV0.9 m in open
jet tunnel
untwisted but pitched, NACA4613, NACA3712
and NACA4611
Yaw angleλ = 5:2
-30° ≤ γ ≤ 30°
NREL UAE
Phase VI Experiment
measurements
10 m in, 4.4 m
by 36.6 m
wind tunnel
twisted, S809
Tip speed ratio,
yaw angle
1.5 ≤ λ ≤ 7.6
0° ≤ γ ≤ 60°
Medici, KTH
Mechanics
measurement
0.18 m in 1.2 m by0.8 m wind
tunnel
untwisted,
Gottingen417A
Freestream,
tip speed
ratio
turbulence
and yaw angle
0° ≤ γ ≤ 60°
1.2e4 < Re < 3.4e4
Maeda (field
experiment),
Mie University
measurements and
tuft visualization
twisted, DU and NACA aerofoils
Wind speed
and yaw angle
-45° ≤ γ ≤ 45°
3.5e5 < Re < 5.5e5
Sant and Haans,
TUDelft
2005–2006HWA and smoke visualization1.2 m in 2.24 m
open‐jet tunnel
chord, twisted,
NACA0012
Yaw angle,
pitch, tip speed ratio
-45° ≤ γ ≤ 45°
Re ≈ 1.7e5
Maeda (tunnel
experiment), Mie
University
measurements
2.4 m in 3.6 m
open jet tunnel
twisted, DU91‐
W2‐250, DU93‐W‐
210 and NACA63‐215
Yaw angle
and tip speed
ratio
1.6 ≤ λ ≤ 6.5
0° ≤ γ≤ 45°
Re ≈ 2.1e5
MEXICO
Experiment
measurements
and PIV
4.5 m turbine
in 9.5 m by
9.5 m open
jet wind tunnel
twisted, DUW2‐
250, RisØ
A1‐21 and
NACA64‐418
Yaw angle,
Tip speed ratio, pitch,
Parked rotor
0° ≤ γ≤ 45°
andγ= ±45° (for PIV)
Re ≈ 6.5e5
Micallef,
TUDelft
2011–2012SPIV measurements2 m turbine in
octagonal open‐jet tunnel of
3 m equivalent diameter
tapered, twisted,
DU96‐W180
Yaw angle and
tip speed ratio
γ= 0° and 30°
λ = 5 (axial flow only) and 7
Re ≈ 3.6e5
measurements, wake
measurements using pitot‐static
probe and LDA
0.9 m turbine in a wind tunnel test‐section 1.9 m by 2.7 m3‐bladed, tapered,
twisted, S826
Freestream, yaw angle and tip speed ratio0° ≤ γ≤ 50°
0.2e5 < Re < 5.0e5

### Table 1.

Overview of experimental campaigns for HAWTs in axial and yawed flow.

Topic area Experimental campaigns Relevant literature
Gross performanceClayton and Filby (University College London)
Grant (Heriot‐Watt University)
Sant and Haans (TUDelft)
Glauert [25],
Clayton and Filby [17],
Grant et al. [18],
Haans et al. [49], Tongchitpakdee et al. [55], Krogstad et al. [53]
Wake geometry and developmentVermeer (TUDelft)
Grant (Heriot‐Watt University)
Sant and Haans (TUDelft)
MEXICO experiment
Micallef (TUDelft)
Øye [31], van Bussell [46, 47],
Snel and Schepers [29, 30],
Voutsinas et al. [45],
Grant et al. [18, 48],
Grant and Parkin [50], Hasegawa [56], Schepers [33], Vermeer [32, 43],
Vermeer et al. [11], Mikkelsen [56], Haans et al. [49, 57, 58], Sezer‐Uzol et al. [59], Schepers and Snel [51], Sant [34],
Schepers et al. [23, 60],
Sørensen et al. [61],
Micallef et al. [52, 54], Micallef [39], Krogstad et al. [53], Sherry et al. [44]
Rotor aerodynamicsBruining (TUDelft)
Vermeer (TUDelft)
Sant and Haans (TUDelft)
Micallef (TUDelft)
van Bussell [46, 47],
Snel and Schepers [29, 30],
Bruining [61], Schepers [23, 33, 60], Mikkelsen [62],
Tongchitpakdee et al. [55],
Haans et al. [49, 57, 63]
Sezer‐Uzol et al. [59], Sant [34, 64],
Maeda et al. [65], Micallef et al. [52, 53], Micallef [39],
Suzuki and Chattot [66],
Aerofoil aerodynamicsVestas Tjaereborg 2MW wind turbine
Bruining (TUDelft)
NREL UAE Phase
I–IV Experiment
NREL UAE Phase
VI Experiment
Maeda (Mie University)—field experiment
Maeda (Mie University)—wind
tunnel experiment
MEXICO experiment
Voutsinas et al. [45], Bruining [67],
Schreck et al. [36, 37, 68–70], Hand et al. [8],
Madsen et al. [7], Johensen et al. [71], Tangler [72], Tongchitpakdee et al. [55],
Maeda and Kawabuchi [73],
Sezer‐Uzol et al. [58], Sant [34],
Breton et al. [74],
Schmitz and Chattot [75],
Shen et al. [76, 77],
Schepers et al. [60], Schepers [23], Holierhoek et al. [35]

### Table 2.

Areas of research in yawed HAWT aerodynamics along with the experimental campaigns available and the relevant literature articles.

## 6. Conclusions and outlook

The fundamental physics of HAWT aerodynamics under yawed conditions has been discussed with reference to some of the latest works. Parallel research on helicopter rotor aerodynamics has been shown to consist of a wealth of information which can be transposed to wind turbine rotor research. While no Mach number dependence would be expected for a HAWT rotor, the Reynolds number effect on the performance of aerofoils becomes very relevant since the currently available controlled experiments reproduce the rotor model on a much smaller scale.

It is clear that the current trend is towards the understanding of detailed flow fields around the blades. While various researchers have successfully utilized vortex methods and CFD to provide fundamental understanding of the flow, there is still the need for simple physical models which can be used in industry with little computational overhead. While BEM approaches utilize a number of engineering models to refine the representation of the physics, the fundamental assumptions in this theory ignore flow three‐dimensionality. A simple approach such as BEM but which takes into account flow three‐dimensionality would be very useful especially when dealing with the issue of yaw. More advanced issues such as 3D aerofoil aerodynamics have also been discussed and an overview of the most fundamental research questions which still need to be addressed has been given.

Knowledge of the complex 3D flow phenomena over wind turbine blades in yaw is still limited. New and more reliable models are necessary for modelling 3D aerofoil data to account for the combined effects of stall delay and unsteady flows. Better models are also necessary for modelling aerofoil data at the tip. The undertaking of more experiments under both controlled and open field conditions and with the latest sophisticated measurements techniques remains crucial for deepening the level of knowledge about the flow physics of yawed turbines.

As larger wind turbines rotors are being developed, especially for offshore deployment, rotor aerodynamics will play a more critical role to ensure a successful future for the wind energy industry. The design of larger wind turbine blades presents new challenges for engineers, making it more important to integrate more sophisticated and efficient aerodynamic models in design codes allowing fully coupled multi‐physical simulation of entire wind turbine systems to solve complex fluid–structural interaction, control and cost optimisation problems. Further research in yaw aerodynamics should also cater for the latest technology development trends involving larger but more flexible blades and the implementation of smart control algorithms allowing to wind turbines in arrays to be controlled with a coordinated approach to minimise the negative impacts of aerodynamic turbine‐to‐turbine interactions on the energy yield from large wind farms.

## References

1. 1. Betz A (1919)Schraubenpropeller mit Geringstem Energieverlust(Dissertation, Gottingen Nachrichten, Gottingen)
2. 2. Dolan, D.S.L. and Lehn, P.W. (2006). Simulation model of wind turbine 3p torque oscillations due to wind shear and tower shadow. InPower Systems Conference and Exposition, 2006. PSCE ‘06. 2006 IEEE PES, Oct. 29 2006–Nov. 1 2006, pp. 2050–2057. doi:10.1109/PSCE.2006.296240
3. 3. Micallef, D., Simao Ferreira, CJ., Sant, T. and Van Bussel, G.J.W. (2010). An analytical model of wake deflection due to shear flow. InThe 3rd EWEA Conference‐Torque 2010: The Science of making Torque from Wind,Heraklion, Crete, Greece, 28–30 June 2010.
4. 4. Rosen, A. and Sheinman, Y. (1994). The average output power of a wind turbine in a turbulent wind.Journal of Wind Engineering and Industrial Aerodynamics, 51(3):287–302. ISSN 0167‐6105.
5. 5. Wagner, R., Antoniou, I., Pedersen, S.M., Courtney, M.S. and Jørgensen, H.E. (2009). The influence of the wind speed profile on wind turbine performance measurements.Wind Energy, 12:348–362. doi:10.1002/we.297
6. 6. Eggleston, D.M. and Starcher, K. (1990). A comparative study of the aerodynamics of several wind turbines using flow visualization.Journal of Solar Energy Engineering, 112:301–309.
7. 7. Madsen, H., Sørensen, N. and Schreck, S. (2003). Yaw aerodynamics analyzed with three codes in comparison with experiment. InAIAA Paper 2003‐519, 41st AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada.
8. 8. Hand, M., Simms, D., Fingersh, L., Jager, D., Cotrell, J., Schreck, S. and Larwood, S. (2001). Unsteady aerodynamics experiment phase VI: Wind tunnel test configurations and available data campaigns. Technical Report NREL/TP‐500‐29955, National Renewable Energy Laboratory, Colorado.
9. 9. Akay, B., Micallef, D., Simão Ferreira, C.J. and van Bussel, G.J.W. (2014). Effects of geometry and tip speed ratio on the HAWT blade's root flow. Journal of Physics: Conference Series 555 (2014) 012002; doi:10.1088/1742-6596/555/1/012002.
10. 10. Micallef, D., Ferreira, C.S., Sant, T., and Bussel, G.J.W. (2015). Experimental and numerical investigation of tip vortex generation and evolution on horizontal axis wind turbines.Wind Energy.doi:10.1002/we.1932
11. 11. Vermeer, L.J., Sørensen, J.N. and Crespo, A. (2003). Wind turbine wake aerodynamics.Progress in Aerospace Sciences, 98:467–510.
12. 12. Wang, T. (2012). A brief review on wind turbine aerodynamics.Theoretical and Applied Mechanics Letters, 2(6). Article 062001, ISSN 2095‐0349.
13. 13. Leishman, J.G. (2002). Challenges in modeling the unsteady aerodynamics of wind turbines.Wind Energy, 5(11):85–132.
14. 14. Barthelmie, R., Frandsen, S., Rathmann, O., Hansen, K., Politis, E., Prospathopoulos, J., Schepers, J., Rados, K., Cabezon, D., Schlez, W., Neubert, A. and Heath, M. (2011). Flow and wakes in large wind farms: Final report for UpWind WP8. Technical Report Risø‐R‐1765(EN), Risø DTU, National Laboratory for Sustainable Energy.
15. 15. Hansen, M.O.L. and Johansen, J. (2004). Tip studies using CFD and comparison with tip loss models.Wind Energy, 7(4), October/December 2004, pp. 343–356.
16. 16. Anderson, M. (1979). Horizontal wind turbine in yaw. InProc. 1st BWEA Wind Energy Conference, Cranfield University of Technology, pp. 68–77.
17. 17. Clayton, B. and Filby, P. (1982). Measured effects of oblique flows and change in blade pitch angle on performance and wake development of model wind turbines. InProc 4th BWEA Wind Energy Conference, Cranfield Institute of Technology, UK.
18. 18. Grant, I., Parkin, P. and Wang, X. (1997). Optical vortex tracking of a horizontal axis wind turbine in yaw using laser‐sheet, flow visualisation.Experiments in Fluids, 23(6):513–519.
19. 19. Hansen, M. (2008).Aerodynamics of Wind Turbines, second edition. Earthscan., UK pp. 45–62.
20. 20. Glauert, H. (1926a).The Analysis of Experimental Results in the Windmill Brake and Vortex Ring States of an Airscrew. ARCR R&M 1026, Aeronautical Research Council (ARS), UK.
21. 21. Coleman, R., Feingold, A. and Stempin C.W. (1945).Evaluation of the Induced Velocity Fields of an Idealized Helicopter Rotor.NACA ARR L5E10, National Advisory Committee for Aeronautics (NACA), USA.
22. 22. Schepers, J.G. and Schreck, S. (2013). The importance of aerodynamics and the role of aerodynamic measurements, ECN‐M‐13‐002, Energy Research Centre of the Netherlands.Presented at theScience of Making Torque from Wind Conference, Olderburg, Germany, 2012.
23. 23. Schepers, J.G. (2012). Engineering models in wind energy aerodynamics—Development, implementation and analysis using dedicated aerodynamic measurements. PhD thesis, Delft University of Technology. Available online fromhttp://repository.tudelft.nl
24. 24. Burton, T., Sharpe, D., Jenkins, N. and Bossanyi, E. (2011).Wind Energy Handbook.John Wiley & Sons Ltd., West Sussex, UK. pp. 139–141.
25. 25. Glauert, H. (1926b).A General Theory of the Autogyro. ARCR R and M 1111, Aeronautical Research Council (ARS), UK.
26. 26. White, F. and Blake, B. (1979). Improved method of predicting helicopter control response and gust sensitivity. InProceedings of the 35th Annual Forum of the American Helicopter Soc., Washington, DC.
27. 27. Pitts, D. and Peters, D. (1981).Theoretical prediction of dynamic inflow derivatives.Vertica, 5:21–34.
28. 28. Howlett, J. (1981).UH‐60A Blackhawk Engineering Simulation Program: Vol 1: Mathematical Model. NASA CR‐66309, National Aeronautics and Space Administration (NASA), USA.
29. 29. Snel, H. and Schepers, J.G. (Ed.) (1994).Joint Investigation of Dynamic Inflow Effects and Implementation of an Engineering Method. Energy Research Centre of the Netherlands, Petten, the Netherlands, ECN‐C‐94‐107.
30. 30. Snel H. and Schepers J.G. (Ed.) (1995).Dynamic Inflow: Yawed Conditions and Partial Span Pitch Control. Energy Research Centre of the Netherlands, Petten, the Netherlands, ECN‐C‐95‐056.
31. 31. Øye, S. (1992). Induced velocities for rotors in yaw. InProceedings of the Sixth IEA Symposium, ECN, Petten, The Netherlands.
32. 32. Vermeer, L. (1998). Wind tunnel experiments on a rotor model in yaw. InProceedings of the 12th symposium on Aerodynamics of Wind Turbines, Lyngby, Denmark.
33. 33. Schepers, J. (1999). An engineering model for yawed conditions, developed on basis of wind tunnel measurements. InProceedings of the 37th Aerospace Sciences Meeting and Exhibit, AIAA‐1999‐39, Reno, NV, USA.
34. 34. Sant, T. (2007). Improving BEM‐based aerodynamics models in wind turbines design codes. PhD thesis, Delft University of Technology. Available online fromhttp://www.tudelft.nl
35. 35. Holierhoek, de Vaal, J.B., van Zuijlen, A.H. and Bijl, H. (2013). Comparing different dynamic stall models.Wind Energy, 16(1):139–158.
36. 36. Schreck, S., Robinson, M., Hand, M. and Simms, D. (2000). HAWT dynamic stall response asymmetries under yawed flow conditions. InProceedings of the 36th ASME/AIAA Wind Energy Symposium, Reno, Nevada.
37. 37. Schreck, S., Robinson, M., Hand, M. and Simms, D. (2001). Blade dynamic stall vortex kinematics for a horizontal axis wind turbine in yawed conditions.Journal of Solar Energy Engineering, 123(4):272–281.
38. 38. Schepers, J., Brand, A., Madsen, H., Stefanatos, N., Simms, D., Hand, M., Bruining, A., van Rooij, R., Shimizu, Y., Maeda, T. and Graham, M. (2001). IEA annex xiv/xviii, field rotor aerodynamics. InProceedings of the EWEC Conference, Brussels.
39. 39. Micallef, D (2012c). 3D flows near a HAWT rotor: A dissection of blade and wake contributions. PhD Thesis, Delft University of Technology.
40. 40. del Campo, V., Ragni, D., Micallef, D., Akay, B., Diez, F.J. and Simão Ferreira, C.J. (2014). 3D load estimation on a horizontal axis wind turbine using SPIV.Wind Energy, 17:1645–1657. doi:10.1002/we.1658
41. 41. del Campo, V., Ragni, D., Micallef, D., Diez, F.J. and Ferreira, C.J.S. (2015). Estimation of loads on a horizontal axis wind turbine operating in yawed flow conditions.Wind Energy, 18:1875–1891. doi:10.1002/we.1794
42. 42. Alfredsson, P.‐H. and Dahlberg, J.‐A. (1981).A Preliminary Wind Tunnel Study of Windmill Wake Dispersion in Various Flow Conditions, The Aeronautical Research Institute of Sweden (FFA), Sweden. Technical Note AU‐1499, Part 5.
43. 43. Vermeer, L. (2001). A review of wind turbine wake research at TUDELFT. InA Collection of the 2001 ASME Wind Energy Symposium Technical Papers, USA pp. 103–113.
44. 44. Sherry, M., Sheridan, J. and Jacono, D.L. (2013). Characterisation of a horizontal axis wind Turbine's tip and root vortices.Journal Experiments in Fluids,54:1417.
45. 45. Voutsinas, S.G., Belessis, M.A. and Rados K.G. (1995). Investigation of the yawed operation of wind turbines by means of a vortex particle method.AGARD‐CP‐552 FDP Symposium on Aerodynamics and Aeroacoustics of Rotorcraft, Berlin, Germany.
46. 46. vanBussel, G. (1993).The Acceleration Potential Models PREDICHAT/PREDICDYN Applied for Calculation of Axisymmetric Dynamic Inflow Cases, IW‐93071R, Delft University of Technology, the Neterlands.
47. 47. vanBussel, G. (1995). The aerodynamics of horizontal axis wind turbine rotors explored with asymptotic expansion methods. PhD thesis, Technische Universiteit Delft.
48. 48. Grant, I., Mo, M., Pan, X., Parkin, P., Powell, J., Reinecke, H., Shuang, K., Coton, F. and Lee, D. (2000). An experimental and numerical study of the vortex filaments in the wake of an operational, horizontal‐axis, wind turbine.Journal of Wind Engineering and Industrial Aerodynamics, 23(6):513–519.
49. 49. Haans, W., Sant, T., van Kuik, G., and van Bussel, G. (2005). Measurement of tip vortex paths in the wake of a HAWT under yawed flow conditions.Journal of Solar Energy Engineering, 127(4):456–463.
50. 50. Grant, I. and Parkin, P. (1999). A DPIV study of the trailing vortex elements from the blades of a horizontal axis wind turbine in yaw.Experiments in Fluids, 28:368–376.
51. 51. Schepers, J.G. and Snel, H. (2007).Model Experiments in Controlled Conditions (MEXICO), Energy Research Centre of the Netherlands, Petten, the Netherlands. ECN‐E‐07‐042.
52. 52. Micallef, D., van Bussel, G.J.W., Ferriera, C.S. and Sant, T. (2012a). The origins of a wind turbine tip vortex. InProceedings from the Science of making Torque from Wind Conference, Olderburg, Germany.
53. 53. Krogstad, P.Å. and Adaramola, M.S. (2012). Performance and near wake measurements of a model horizontal axis wind turbine.Wind Energy, 15(5):743–756.
54. 54. Micallef, D., van Bussel, G.J.W., Ferriera, C.S. and Sant, T. (2012b). An investigation of radial velocities for a horizontal axis wind turbine in axial and yawed flows.Wind Energy, 16(4):529–544.
55. 55. Tongchitpakdee, C., Benjanirat, S. and Sankar, L.N. (2005). Numerical simulation of the aerodynamics of horizontal axis wind turbines under yawed flow conditions.Journal of Solar Energy Engineering,127(4):464–474.
56. 56. Hasegawa, Y. (1999). Numerical analysis of yawed inflow effects on a HAWT rotor. InProc. of 3rd ASME/JSME Joint Fluids Eng. Conf., FEDSM99‐7820.
57. 57. Haans, W., Sant, T., van Kuik, G. and van Bussel, G. (2008). HAWT near wake aerodynamics, part I: axial flow conditions.Wind Energy, 11(3):19.
58. 58. Haans, W. (2011). Wind turbine aerodynamics in yaw, unravelling the measured rotorwake. PhD Thesis, Delft University of Technology. Available online fromhttp://repository.tudelft.nl
59. 59. Sezer‐Uzol, N. and Long, L.N. (2006). AIAA Paper 2006‐0394.Presented at AIAA Aerospace Sciences Meeting, Reno, NV.
60. 60. Schepers, J.G., Boorsma, K., et al. (2011).Final Report of IEA Task 29, Mexnext (Phase 1), Analysis of Mexico Wind Tunnel Measurements. ECN‐E‐12‐004, Energy Research Center of the Netherlands, ECN.
61. 61. Sørensen, N.N., Bechmann, A., Réthoré, P.‐E. and Zahle, F. (2012). Near wake Reynolds‐averaged Navier–Stokes predictions of the wake behind the MEXICO rotor in axial and yawed flow conditions.Wind Energy., 17(1), 2014, p. 75–86. doi:10.1002/we.155
62. 62. Mikkelsen, R. (2003). Actuator disc methods applied to wind turbines. PhD thesis, Technical University of Denmark, Lyngby, Denmark.
63. 63. Haans, W., van Kuik, G. and van Bussel, G. (2007). Experimentally observed effects of yaw misalignment on the inflow in the rotor plane.Journal of Physics: Conference Series, 75:10.
64. 64. Sant, T., van Bussel, G. and van Kuik, G. (2009). Estimating the angle of attack from blade pressure measurements on the NREL Phase VI using a free wake vortex model: Yawed conditions.Wind Energy, 12:1–32.
65. 65. Maeda, T., Kamada, Y., Suzuki, J. and Fujioka, H. (2008). Rotor blade sectional performance under yawed inflow conditions.Journal of Solar Energy Engineering, 130:031018‐1–031018‐7.
66. 66. Suzuki, K. and Chattot, J.J. (2012). Unsteady hybrid Navier–Stokes/vortex model applied to wind turbine aerodynamics under yaw conditions. InProceedings from the Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii.
67. 67. Bruining, W. (1997). Aerodynamic characteristics of a 10 m diameter rotating wind turbine blade experimental results from the Delft University of Technology for the “Dynamic Stall and 3D Effects Project”. IW95‐084R, TU Delft, Institute of Wind Energy.
68. 68. Shreck, S. and Robinson, M. (2002). Rotational augmentation of horizontal axis wind turbine blade aerodynamic response.Wind Energy Special Issue: Analysis and Modeling of the NREL Full‐Scale Wind Tunnel Experiment, 5(2–3):133–150.
69. 69. Schreck, S. and Robinson, M. (2005). Blade three‐dimensional dynamic stall response to wind turbine operating condition.Journal of Solar Energy Engineering,127(4):488–495.
70. 70. Schreck, S., Sant, T. and Micallef, D. (2010). Rotational augmentation disparities in MEXICO and UAE phase VI experiments. InProceedings from the Conference Science of Making Torque from Wind, Organised by the European Academy of Wind Energy, Crete.
71. 71. Johansen, J. and Sørensen, N.N. (2004). Aerofoil characteristics from 3D CFD rotor computations.Wind Energy, 7(4):283–294.
72. 72. Tangler, J.L. (2004). Insight into wind turbine stall and post‐stall aerodynamics.Wind Energy, 7(3):247–260.
73. 73. Maeda, T. and Kawabuchi, H. (2005). Surface pressure measurement on a rotating blade of field horizontal axis wind turbine in yawed conditions.JSME International Journal, 48:156–163.
74. 74. Breton, S.P., Coton, F.N. and Moe, G. (2008). A study on rotational effects and different stall delay models using a prescribed wake vortex scheme and NREL phase VI experiment data.Wind Energy, 11(5):459–482.
75. 75. Schmitz, S. and Chattot, J. (2009). Flow physics and stokes theorem in wind turbine aerodynamics.Computational Fluid Dynamics,9:801–806.
76. 76. Shen W.Z., Hansen, M. and Sørensen, J. (2008). Determination of the angle of attack on rotor blades.Wind Energy,12(1):91–98.
77. 77. Shen, W.Z., Hansen, M.O.L. and Sørensen, J.N. (2009). Determination of the angle of attack on rotor blades.Wind Energy, 12(1):91–98.

Written By

Daniel Micallef and Tonio Sant

Submitted: October 29th, 2015 Reviewed: April 4th, 2016 Published: July 27th, 2016