Abstract
Data are presented comparing upwind versus downwind operation of the National Renewable Energy Laboratory’s Phase VI wind turbine. Power was not reduced as expected for downwind, which may be attributed to inboard three-dimensional effects. Average flap bending loads were reduced with downwind coning and compared well with prediction. Fatigue loads were increased with downwind; however, fatigue was mitigated with a tower shroud. The shroud needs to align with the freestream, demonstrated by an increase in fatigue loads from a 10° shroud misalignment. Pressure data were acquired in the tower wake at the rotor location with and without the shroud. The bare tower wake data compared well with previous work. The shroud wake data at 10° and 20° misalignment showed velocity reduction and turbulence exceeding the bare tower values. Downwind operation, with an aligning tower shroud, should be considered for future designs given the load benefits of downwind coning.
Introduction
The predominant wind turbine configuration is horizontal-axis, three-bladed, upwind rotors (rotor upwind of tower). Several large downwind research turbines were constructed in the past; examples of these in the United States were the 100 kW MOD-0 (Glasgow et al., 1981) and the Hamilton Standard 4 MW WTS-4 (Spera, 2009), which were US government-sponsored research projects in response to the energy crises of the 1970s. Smaller downwind turbines available commercially in the 1980s were the 80-kW ESI-80 and the 25-kW Carter (Gipe, 1995). Downwind turbines have the potential for lower average blade and yawing loads and thus lower mass/cost. The design also effectively removes the constraint of blade tip clearance with the tower. However, problems such as fatigue and noise arise from the rotor interacting with the tower wake. Recent advances in offshore wind energy, including turbines on floating platforms, have renewed interest in large downwind designs that can be placed far from shore. Ichter et al. (2016) and Loth et al. (2017) present a large downwind concept, which was the inspiration for this writing.
This article reports on unpublished data from the National Renewable Energy Laboratory’s (NREL) Unsteady Aerodynamics Experiment (UAE) Phase VI (Hand et al., 2001) which was obtained at the National Full-Scale Aerodynamics Complex 80- by 120-foot wind tunnel at NASA Ames Research Center. The prior UAE literature was focused on model validations and not comparison between upwind and downwind operation. For example, Coton et al. (2002) report on modeling (prescribed wake) for the downwind cases in the UAE Phase VI data. They found that the modeling results were highly dependent on the wake model, which included a wake width and velocity deficit. They also report on a phase difference between the measurements and the modeling, which was suspected to be a rotor–tower wake interaction. This seemed to be confirmed in the modeling results (RANS) of Zahle et al. (2009).
Experimental comparisons between upwind and downwind operation have been rare. Glasgow et al. (1981) report on upwind and downwind operation of the MOD-0 (unconed). Mean blade bending moments were found to be the same, but cyclic moment showed an increasing trend with wind speed for downwind operation. Power measurements were not available for this experiment. Yoshida (2006) reports on a 100-kW turbine configured for upwind and downwind operation. Power was increased 7%–10% in downwind operation, which was attributed to a favorable combination of sloping terrain and shaft tilt. Water-channel experiments by Kress et al. (2015) showed higher power for downwind compared to upwind operation, which was attributed to favorable nacelle-flow confirmed in modeling by Frau et al. (2015).
This article addresses power and load differences between upwind and downwind operation, in addition to potential benefits from a tower aerodynamic shroud (fairing). The shroud was included in the experiment to study mitigating the tower wake. Tower wake measurements with and without the shroud will also be presented and compared to previous research. Previous research on tower wakes were from wind tunnel studies by Snyder and Wentz (1981) and Powles (1983). While Powles studied the wake of a 12-sided polygon, Snyder and Wentz also studied the polygon in addition to a cylinder and a cylinder with strakes. Wilmshurst et al. (1985) added an aerodynamic shroud to the tower in further work. O’Connor et al. (2013) reported on a self-aligning fairing design for wind turbine towers, with later experimental studies (O’Connor et al., 2015).
Methods
The NREL Unsteady Aerodynamics Experiment Phase VI wind turbine was two bladed with a 10-m diameter (Figure 1). The test was conducted at the NASA-Ames 24.4 by 36.6 m (80- by 120-foot) wind tunnel. The tunnel is the largest in the world and the error due to blockage from the turbine was determined to be less than 2% for all conditions (Hand et al., 2001). The turbine could be operated in upwind or downwind mode, with adjustable blade coning. For the results presented in this article, the coning was 0° for upwind operation and 3.4° (downwind) for downwind operation. These coning angles were maintained during the test over concerns of loading and stability (Simms et al., 1999). The teeter degree-of-freedom was fixed for upwind operation and free for downwind operation. There was an upwind sequence with teeter; however, the loads/power below rated were similar to fixed teeter. The fixed pitch, stall-regulated rotor was operated at a constant 72 r/min. Most of the comparisons in this article are made at wind speeds below rated to avoid the complication of blade stall. Turbine and tunnel data were acquired at 521 Hz and most data points represent 30 s of data (36 rotor revolutions). Blade loads were obtained with strain gage bridges which were calibrated in the tunnel. Blade surface pressures were also measured at various radial positions. One of the blades was equipped with five-hole pitot-tube probes at various radial positions to determined blade angles and velocities. Tower wake data were acquired by positioning the blade instrumented with the probes downstream of the tower (Figure 2).

UAE in upwind configuration. NASA image.

Tower shroud with instrumented blade positioned downwind for wake measurements. Photo by Lee Fingersh, NREL 36813.
Both blade loads and tower wake data were obtained with and without a tower aerodynamic shroud installed to investigate mitigating adverse effects of the tower wake. The tower diameter in the rotor swept area was a constant 0.4064 m. The maximum thickness of the symmetric shroud was 0.46 m with a 0.89-m chord length. The shroud coordinates are included in the Supplemental material. The shroud was installed with plywood ribs attached to the tower. Aluminum sheets were bent around the ribs with a finite seam at the nominally sharp trailing edge.
Note that a formal uncertainty analysis was not performed on the test; however, the expected accuracy of the data was within 1% given the calibration procedures and data review (Simms et al., 2001). For most plots in this article, the measurement uncertainty would be less than the marker size.
Table 1 lists the test sequences that were used in this study. Further details of the test campaign are available in Hand et al. (2001).
UAE test sequences.
UAE: Unsteady Aerodynamics Experiment.
Results and discussion
Power and loads
Figure 3(a) shows comparison of downwind versus upwind average power for the full range of wind speeds. With downwind operating at 3.4° coning, one would expect a lower power due to reduced rotor swept area (0.35%) and loss due to tower shadow. However, the annual energy production (AEP) for downwind operation is over 2% greater than upwind operation for all IEC wind classes (Rayleigh wind speed distribution) with this power curve. Figure 3(b) focuses on below-rated wind speeds, and downwind shows higher power except at the highest wind speed on the plot (9.3 m/s). In this range, the downwind AEP is 10% higher than upwind. The higher downwind power can be investigated with Figure 4, which shows a comparison of averaged spanwise aerodynamic performance between the upwind and downwind operation at 5 m/s, well below the rated speed of 10 m/s. Downwind torque coefficient is higher inboard compared to upwind, as shown in Figure 4(a). Downwind also shows higher inboard normal force (Figure 4(b)) and tangential force coefficient (Figure 4(c)). However, the local angle of attack is lower for downwind (Figure 4(d)); for example, alpha = 4.4° compared to 5.2° at 91% radius. The spanwise angle (angle between local chordline and relative wind speed vector) is higher for the downwind case (Figure 4(e)). The pressure coefficient (cp) at the 30% station and 0° azimuth (Figure 4(f)) shows a higher suction peak for downwind operation, which corresponds with the higher local torque coefficients in downwind operation. It is not clear why the downwind inboard rotor performance appears to be higher than upwind operations in fully attached flow conditions well away from stall. Potentially, the large rotating root region instrumentation package introduced additional three-dimensional effects and increased inboard radial flow in the upwind configuration, whereas the downwind rotor is affected by the wake of the stationary, streamlined nacelle (Figure 1). Frau et al. (2015) show similar power increases for downwind operation and similar inboard behavior.

Average power for upwind and downwind operation: (a) complete wind speed range and (b) below-rated wind speed. Downwind coning is 3.4°.

The 5-m/s data: (a) torque coefficient, (b) normal force coefficient, (c) tangential force coefficient, (d) local flow angle, (e) spanwise flow angle, and (f) pressure coefficient (30% radial station and 0° azimuth).
Figure 5 shows the difference in average flap bending moment for upwind and downwind operation. The downwind-coned rotor has a tendency to become unconed which produces a flap bending moment opposing the moment from aerodynamic loads. The predicted difference due to coning is 1040 N m, which compares very well with the measurements. With this setting of downwind coning, the moment becomes negative at low wind wind speeds.

Average Blade 1 root flap bending moment for upwind and downwind operation.
Blade fatigue is shown in Figure 6, with the damage equivalent load (DEL) for upwind and downwind operation, including downwind with shroud. The DEL was calculated using a rainflow-counting method with MLife (NWTC Information Portal, 2015). Results are shown in the region below stall operation of the turbine. In this region compared to upwind operation, blade fatigue is increased for downwind operation (e.g. 50% increase in DEL at 7 m/s). The use of a tower shroud results in significant reduction in fatigue; however, fatigue loads increase with shroud misalignment (right plot, e.g. 13.5% increase for 10° misalignment at 7 m/s). Here, misalignment is defined as the angle between the freestream velocity and mean chordline of the symmetrical shroud. Note that operating data were not acquired for −10° or −20° misalignment. Causes for these loads are examined below in the wake data study.

Blade 1 root flap bending moment damage equivalent load (DEL): (a) upwind, downwind, and downwind with shroud and (b) downwind with aligned and misaligned shroud.
Figures 7 and 8, respectively, show the azimuth averaged Blade 1 root flap and edge bending moments at 10 m/s and 0° misalignment. The moments have been averaged over 36 revolutions. The blade response is periodic for the loads and configurations shown. The difference in average flap moment due to coning is evident in Figure 7. The edge moments for downwind operation show higher harmonics in Figure 8 compared to the upwind trace, which primarily shows once-per-revolution variation due to the weight vector. The edge loads show periodic excitation from interaction with the tower/shroud wake. The shroud appears to reduce the magnitude of this wake induced mode.

Azimuth-averaged Blade 1 root flap moments at 10 m/s wind speed.

Azimuth-averaged Blade 1 root edge moments at 10 m/s wind speed.
Tower wake data
One of the turbine blades was instrumented with five-hole probes to measure upstream dynamic pressure and flow angle. This blade was fixed in position downwind of the tower to obtain wake measurements Figure 9 shows the wake velocity normalized to tunnel velocity at three tower diameters downwind and at 7 m/s wind speed. The Reynolds number for this condition is subcritical in the long-established experimental results for cylinder drag. Above this Reynolds number, the transition point moves aft and the drag coefficient lowers significantly. A typical cos2 tower wake model from the literature (Powles, 1983) is shown as well for comparison. The model wake velocity is given by the following
The comparison with the data is good; however, this model requires knowledge of the maximum wake deficit

Average wake velocity normalized by tunnel velocity (arrows with length according to grid), with cos2 wake model (dashes). 7 m/s,
Figure 10(a) shows the tower wake velocity normalized to tunnel velocity at 7, 15, and 20 m/s wind speeds. Figure 10(b) shows at the same speeds the turbulence intensity, defined as the standard deviation of the wake velocity divided by the average tunnel velocity. These speeds represent subcritical, transitional, and supercritical Reynolds numbers for the cylinder wake. The wake deficit and turbulence intensity reduces with Reynolds number and becomes asymmetric. The subcritical wake also shows a double peak in the turbulence intensity. These behaviors were similarly observed in experiments of cylinder wakes by Snyder and Wentz (1981) (this reference should be considered the comprehensive study of wind turbine tower wakes). The double peak was attributed to areas of high mixing in the subcritical wake. The asymmetry was attributed to turbulent separation on one side of the cylinder and laminar separation on the other side. Note that for modern MW-scale turbines, the tower Reynolds number would most likely be supercritical at all operational wind speeds.

(a) Tower wake normalized velocity at Reynolds numbers corresponding to 7, 15, and 20 m/s and (b) turbulence intensity at same Reynolds numbers.
Figure 11 shows the spectra of the wake pressure signal at 7 m/s, with a peak at 2.8 Hz. This corresponds to a Strouhal number

Spectra of pressure of 91% five-hole probe center port in tower wake at 7 m/s.
Figures 12 to 14 show wake measurements at 7 m/s with the shroud installed and at various misalignment angles. The cos2 wake model from Figure 9 is included for comparison. Figure 12 shows the wake with the shroud aligned with the freestream. The average velocity in the wake is actually higher than the freestream (arrow length greater than 1) where the measurements were taken. Figure 13 shows the shroud misaligned ±10° from the freestream. This ±10° shroud misalignment results in a larger wake deficit than the unshrouded tower configuration. Figure 14 shows the shroud misaligned ±20° from the freestream, with regions of significant reversed flow.

Measured normalized wake velocity at 7 m/s with tower shroud at 0° with cos2 (dashes) model from Figure 7.

Measured normalized wake velocity at 7 m/s with tower shroud at (a) 10° and (b) −10° misalignment with cos2 (dashes) model from Figure 7.

Measured normalized wake velocity at 7 m/s with tower shroud at (a) 20° and (b) −20° misalignment with cos2 (dashes) model from Figure 7.
Figure 15(a) shows the shroud wake velocity normalized to tunnel velocity at 7, 15, and 20 m/s wind speeds for 0° shroud misalignment. Figure 15(b) shows at the same speeds the turbulence intensity. The wake deficit increases slightly with Reynolds number and is more narrow than the bare tower wake. These same findings were found in experiments with a tower shroud by Wilmshurst et al. (1985). The turbulence intensity is much lower than the unshrouded tower wake, which is shown in Figure 10(b). No peaks were observed in the spectra of the shroud wake. In summary, the shroud wake is narrower and more steady than the bare tower which has a wider, more fluctuating wake. These results can explain the much lower DEL in Figure 6 with the shroud installed. Note that tower strakes, tower tapering, or lattice towers could also reduce the tower wake.

(a) tower wake normalized velocity at Reynolds numbers corresponding to 7, 15, and 20 m/s for 0° shroud misalignment and (b) turbulence intensity at same Reynolds numbers.
Figures 16 and 17 show wake normalized velocity and turbulence intensity for ±10° and ±20° shroud misalignment, respectively. With shroud misalignment, the wake velocity decreases and turbulence intensity increases dramatically which results in the increased DEL, as shown in Figure 6(b). These data also show a Reynolds number dependency, especially at the lowest value.

Shroud wake normalized velocity (left) and turbulence intensity (right) at chord Reynolds numbers corresponding to 7, 15, and 20 m/s for (a and b) ±10° and (c and d) −10° shroud misalignment.

Shroud wake normalized velocity (left) and turbulence intensity (right) at chord Reynolds numbers corresponding to 7, 15, and 20 m/s for (a and b) +20° and (c and d) −20° shroud misalignment.
Effect of tower wake on fatigue
As shown in Figures 7 and 8, the tower and shroud wake result in a periodic load on the blades. To study the effect of the wake on the blade fatigue, the blade was modeled as a single degree-of-freedom spring-mass-damper with periodic forcing from a pulse. The pulse is approximated by a 200-term (harmonics) Fourier series, and the blade response is a linear combination of the response of each harmonic. Balachandran and Magrab (2008) provide example response equations, which are shown in the Supplemental material. The response is linearly based on the pulse (load) magnitude, but is non-linearly related to the ratio of the pulse width (i.e. the wake width) to the period, the damping ratio, and the ratio of the pulse frequency to the blade natural frequency. Example 30 s time series were produced for a range of these parameters for values ranging from the UAE to an example 10 MW wind turbine (Stäblein et al., 2016). And example response is shown in Figure 18.

Example pulse (200 harmonics) response for 0.15 frequency ratio, 0.1 pulse ratio, and damping
The DEL was then determined using a rainflow-counting method with MLife (NWTC Information Portal, 2015) because there is no known closed-form DEL solutions for a periodic response. Figure 19 shows the normalized DEL for the pulse excitation as a function of the pulse ratio to the frequency ratio. Figure 19(a) is for a damping ratio of 0.05, which is representative of blade edge bending. Figure 19(b) is for a damping ratio of 0.7, which is representative of blade flap bending (Stäblein et al., 2016). The high damping plot shows increasing load with increasing pulse width and frequency ratio; however, this is not observed in the low damping case.

Damage equivalent load for pulse excitation as function of width of pulse to period (pulse ratio) and ratio of pulse frequency to natural frequency for damping ratio of (a) 0.05 and (b) 0.7.
Figure 20 shows the loading on the outboard portion of the blade resulting from the tower-alone wake and the shroud wake. The difference in width and depth of the wake disturbance is listed in Table 2. As mentioned above for the pulse, the difference in depth would result in a linear difference in DEL. The shroud peak is approximately 15% less then the tower-alone case for both load components. The shroud width is 35% less than the tower-alone case, which results in a different DEL for the pulse case depending on the damping ratio. The DEL is reduced by 11% for a damping ratio of 0.7, which would correspond to flap bending and normal coefficient loading. The DEL is increased 2% for a damping ratio of 0.05, which would correspond to edge bending and tangential coefficient loading. The actual DEL differences are listed at the end of Table 2, which show a decrease in flap bending DEL of 23% and an increase in edge DEL of 2%.

Azimuth-averaged and mean-removed (a) normal force coefficient and (b) tangential force coefficient at 95% station at 7 m/s nominal freestream wind speed.
Effect of differences in tower and shroud wake on damage equivalent load at 7 m/s nominal freestream wind speed.
DEL: damage equivalent load.
Other features to note about Figure 20 include an initial rise in the loads prior to the drop (minimum). The rise is more pronounced for the bare tower The recovery after the minimum is similar for the tower and the shroud. Also, the shroud minimum is offset from 180°, which may indicate a shroud misalignment or assymetric transition on the shroud. Also, the amount of minimum for the shroud is curious considering the lack of wake depth in Figure 15(a). This may indicate a coupling between the rotor and the tower/shroud wake, which is suggested for the tower wake in Zahle et al. (2009). The shroud design would be improved with further knowledge of this coupling. Also note that the time series (not averaged) for Figure 20 is repeatable and does not indicate unsteadiness from the wake itself.
The effect of the shroud on the loads indicate that the primary goal should be to reduce the magnitude of the disturbance. The effect of reducing the width of the wake is more complex; however, if the damping is increased in the edgewise direction (e.g. with tuned mass dampers), the decreased wake width should result in lower loads.
Conclusion
The comparison between upwind and downwind operation of NREL’s Phase VI Unsteady Aerodynamics Experiment shows that downwind operation should be reconsidered for future large wind turbines. Downwind coning offers a significant and predictable reduction in average loads for blade flap bending, in comparison with upwind operation and does not seem to result in lower energy capture.
The expected increased fatigue loads from downwind operation compared to upwind operation was confirmed in the experiment. The loads are the result of the tower wake, which was measured in this experiment and compare well with previous wake studies. The wake is dependent on the Reynolds numbers encountered in this study. However, the higher, supercritical Reynolds numbers and fully turbulent wakes for large modern turbines may remove this dependency. In regard to the tower wake measurements, the previous research (Snyder and Wentz, 1981) should be considered more comprehensive than the UAE measurements.
The fatigue loads in tower wake operation can be mitigated significantly with an aerodynamic shroud. However, the shroud must remain aligned with the wind direction, as the results show fatigue loads returning to tower-alone levels with a 10° misalignment. So a free-yawing wind-vane type shroud or a system tied to the rotor yaw would be required for practical implementation. Further research may demonstrate whether a shroud will lower the noise generation.
The Phase VI was machine was stall-regulated with fixed pitch. Future wind tunnel campaigns should consider modern rotor configurations with pitch regulation along with comprehensive rotor wake measurements.
Supplemental Material
WIE780387_Supplemental_material_CLN – Supplemental material for Reconsidering downwind operation by analysis of the National Renewable Energy Laboratory’s Phase VI data
Supplemental material, WIE780387_Supplemental_material_CLN for Reconsidering downwind operation by analysis of the National Renewable Energy Laboratory’s Phase VI data by Scott Larwood and Raymond Chow in Wind Engineering
Footnotes
Acknowledgements
The authors would like to thank the reviewers of this work for their time and helpful comments. Scott Schreck (formerly of NREL) handled several requests for data. Lee Jay Fingersh of NREL searched for and provided drawings of the tower shroud. The authors would also like to thank the entire UAE test team. Portions of this article were presented at TORQUE 2016: The Science of Making Torque from the Wind (Larwood and Chow, 2016).
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Supplemental material
Supplemental material is available for this article online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
