Comparison of Simulations and Wind Tunnel Measurements For The Improvement of Design Tools For Vertical Axis Wind Turbines
Comparison of Simulations and Wind Tunnel Measurements For The Improvement of Design Tools For Vertical Axis Wind Turbines
Comparison of Simulations and Wind Tunnel Measurements For The Improvement of Design Tools For Vertical Axis Wind Turbines
Abstract. As wind turbines are getting larger and cost reductions through up-scaling are
reaching a limit, there is increasing pressure to reduce their installation and operating costs,
making the business case for Vertical Axis Wind Turbines (VAWTs) more interesting again.
However, in order to properly assess the economic potential of large-scale VAWTs, improved
design tools are required. As the research budget for VAWT projects is generally too small to
carry out expensive large-scale wind tunnel tests and field measurements, alternative methods
are required. In this work, the application of small-scale, lower budget methods for improving
VAWT design tools is assessed. It is shown that (a) currently available tools for VAWT design
have not been sufficiently validated, (b) lower budget, small-scale wind tunnel tests can be
effective for examining VAWT performance in terms of average power coefficient vs. tip speed
ratio as well as the dynamics of the torque, and (c) Detached Eddy Simulations on small-scale
VAWTs can be effective for examining their performance in terms of average power coefficient
vs. tip speed ratio as well as the forces on the blades. Based on these results, a process
for the transfer of the results of lower budget, small-scale measurements and simulations into
recommendations for improved design tools is being developed.
1. Introduction
Vertical Axis Wind Turbine (VAWT) technology has not yet managed to establish itself as
a reliable technology for wind energy projects. This is mainly because VAWTs have several
disadvantages over Horizontal Axis Wind Turbines (HAWTs), including a lower efficiency, higher
bearing loads and lower average height above the ground. However, advantages include the
positioning of the drivetrain at the base, leading to reduced installation and service costs, and
the independency of operation to wind direction, resulting in less fatigue. As wind turbines are
getting larger and cost reductions through up-scaling are reaching a limit, there is increasing
pressure to reduce their installation and operating costs. This is making the business case for
Vertical Axis Wind Turbines (VAWTs) more interesting again.
This increased interest is reflected in the recent research activities aimed at better
understanding the dynamics and aerodynamics of VAWTs - for example from the DeepWind
project ([1], [2] and [3]) and from Sandia National Laboratories ([4], [5] and [6]) as well as the
development of new offshore floating concepts such as [7], [8] and [9].
However, in order to properly assess the economic potential of large-scale VAWTs, improved
design tools are required. Usually, design tools are improved by comparing them to costly wind
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
tunnel and full-scale field tests in which the shaft torque and blade loads are measured at a
range of flow conditions, in order to better understand the physics. As the research budget for
VAWT development projects is generally small and because very few large-scale VAWTs are
actually in operation, it is currently not feasible to carry out such measurement campaigns.
In this work, the application of small-scale, lower budget methods for improving VAWT
design tools is assessed. This is done by firstly reviewing the currently available design tools,
then carrying out small-scale wind tunnel measurements and Detached Eddy Simulations (DES)
and comparing the results to simulations using an open-source design tool, QBlade. Finally, the
feasibility of these methods for improving design tools is assessed.
2
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
3
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
The key design parameter of a wind turbine is the shaft torque, allowing the rotor power
coefficient (CP ) to be calculated as shown in Equation 1, where τ = rotor torque (Nm), ω =
rotational speed (rad/s), ρ = air density (kg/m3 ), V0 = inflow wind speed (m/s) and A = rotor
cross-sectional area (m2 ).
τω
CP = 1 3
(1)
2 ρV0 A
A Lorenz Messtechnik in-line torque sensor type DR-2112R was mounted between the motor
and the rotor and measurements were carried out at a frequency of 2 kHz for at least 30 seconds
per measurement point. Blade loads are also key to assessing design tools; however, they were
not measured due to the expense and complexity of installing strain gauges in the rotating
system.
The results in terms of power coefficient vs. tip speed ratio are shown in Figure 3 for wind
speeds of 6.0 m/s, 6.5 m/s and 7.0 m/s, where the tip speed ratio (λ) is given by Equation
2 and R = rotor radius (m). They are compared to measurements on a similar-sized VAWT
model from the University of New Hampshire [6], showing that the approximate magnitude is
reasonable (the dimensions, aspect ratio, solidity and airfoil types are different). This shows that
the chosen set-up could be used successfully to measure the full range of the power coefficient
vs. tip speed ratio curve of the VAWT model.
Rω
λ= (2)
V0
Examination of the 2 kHz torque signal allows the periodical behaviour of the rotor to be
analysed. Figure 4 shows the frequency content of three typical measurements carried out for
60 seconds each at tip speed ratios of 0.82, 1.65 and 2.47 (200, 400 and 600 rpm), obtained by
carrying out Fast Fourier Transforms on the measurement data. As expected, clear peaks can
be seen at the rotational speed (1P ), at the blade passing frequency (three times the rotational
speed = 3P ), at twice the blade passing frequency (6P ) and at three times the blade passing
frequency (9P ). Large peaks at 2P can also be seen for every rotational speed, indicating a rotor
4
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
imbalance from an error in manufacture or pitch setting of one of the blades. Additionally, large
peaks at 4.3P can be seen. These are thought to be caused by the dynamics of the motor, which
turns 4.3 times faster than the rotor. A small peak can be seen on each plot at a frequency
of approximately 23 Hz, circled on the figure, indicating a dynamic that is independent of the
rotor speed. This could be the first torsional natural frequency of the structure.
Figure 4. Frequency content of a three typical torque signals at tip speed ratios of 0.82, 1.65
and 2.47 (200, 400 and 600 rpm).
As the frequency content at 4.3P is due to the motor and not to the wind turbine rotor, it
needs to be removed in order to examine the dynamic behaviour of the rotor. On the left-hand
graph of Figure 5, a typical torque signal at a tip speed ratio = 2 is shown with a low-pass
filter of 100 Hz applied. The periodic nature of the torque can be seen; however, the 4.3P
content is dominating. Applying a 30 Hz low pass filter removes the 4.3P frequency, but also
the interesting 6P and 9P frequencies (middle graph). Here only the blade passing frequency
can be seen, as well as the asymmetry due to the misalignment of one blade as discussed above.
The 4.3P peak was therefore removed manually from the Fourier transformed-data and then an
Inverse Fast Fourier Transform was performed on the resulting dataset, giving the graph shown
on the right-hand side of Figure 5.
This shows a distinctive periodical behaviour; the torque does not smoothly transition from
the maximum to the minimum value, but it increases again slightly in the middle. This is
thought to be caused by the interaction of the blades with the wake, and is examined further in
the next section.
Figure 5. Typical torque signal at tip speed ratio = 2 (rotational speed = 435 rpm = 7.25 Hz).
In summary, it has been found that lower budget small-scale wind tunnel tests can be used to
compare the performance in terms of power coefficient vs. tip speed ratio as well as to investigate
5
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
the natural frequencies and the dynamic behaviour of the total torque on the shaft within a rotor
revolution. This has the potential for helping to understand the physics of VAWTs and thus
improving design tools.
6
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
Figure 6. (a) Set-up of the Detached Eddy Simulations for one blade; (b) Results in terms of
lift coefficient vs. angle of attack.
The results in terms of power coefficient vs. tip speed ratio are shown in in Figure 8 compared
to the wind tunnel measurements. At and above the maximum CP , the results match fairly well.
At lower tip speed ratios, the simulations match less well, due to the increasing average angle
of attack (leading to more flow separation, which is difficult to simulate) and the decreasing
Reynolds Number (reducing the likelihood of a turbulent flow, which is easier to simulate).
Further differences are thought to be caused by the non-uniform wind input field in the wind
tunnel, which is not accounted for in the simulations.
In Figure 9, the variation of torque with time for a tip speed ratio = 2 is examined. In
part (a), the torque exerted on the shaft by one blade is shown. One large peak is visible per
revolution, when the blade passes the front of its swept area and is operating at its optimal angle
of attack and producing lift. A second, smaller peak is seen as the blade passes the back of its
swept area, producing a small amount of lift. The torque even becomes negative at two azimuth
locations, because here the drag dominates. The torque produced by the other two blades shows
a similar pattern, offset by one-third of a revolution. The total torque, shown in part (b), is
7
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
then given by the sum of these three torques at each azimuth position. Like the measurements
shown in the previous section, this signal does not smoothly transition from the maximum to the
minimum value, but it increases again slightly in the middle. This occurs due to the presence
of the second small peak as each blade passes behind the tower. Additionally, it can be seen
in part (b) that the magnitude of the torque fluctuations matches well with the measurements
(both approximately 0.15 Nm), confirming the ability of the simulations to predict the physics
correctly. Part (b) also shows a comparison between the positive torque produced by the blades
and the resistance caused by the structure of the wind turbine, comprising approximately 10%
of the total torque.
Figure 9. Variation of torque with time from DES: (a) For one blade; (b) For all blades (solid
line) and for the structure (dashed line).
In summary, it has been found that lower budget Detached Eddy Simulations on small-scale
VAWTs can be used to compare the performance in terms of power coefficient vs. tip speed
8
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
ratio as well as the dynamic behaviour of the wind turbine. This has the potential for helping
to understand the physics of VAWTs and thus improving design tools. Further analysis of these
results is on-going.
Figure 10. (a) Comparison of results from measurements, DES and QBlade in terms of power
coefficient vs. tip speed ratio; (b) Power coefficient vs. time for the total torque from the QBlade
LLFVW simulations for a tip speed ratio = 2.
9
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
and Detached Eddy Simulations were successfully used to examine the average torque as well
as its variation with time in order to examine the physics of VAWTs. This has the potential
for improving VAWT design tools and ultimately improving VAWT design and reducing wind
farm operating costs. Based on these results, a process for the transfer of the results of lower
budget, small-scale measurements and simulations into recommendations for improved design
tools is currently being developed.
4. Conclusions
This work has shown the following:
• Currently available aerodynamic tools for VAWT design have not been sufficiently validated
in terms of the shaft torque and blade loads, the most important design parameters.
• Lower budget, small-scale wind tunnel tests can be effectively used for examining VAWT
performance in terms of average power coefficient vs. tip speed ratio as well as the dynamics
of the torque.
• Lower budget Detached Eddy Simulations on small-scale VAWTs can be effectively used
for examining their performance in terms of average power coefficient vs. tip speed ratio as
well as the forces on the blades.
• A process for the transfer of the results of lower budget, small-scale measurements and
simulations into recommendations for improved design tools is currently being developed.
Acknowledgments
Thanks to Lukas Schnider, Simon Boller and Roland Peterer for help with the wind tunnel
measurements.
References
[1] Najafi N and Paulsen U S 2017 Energy 405–416
[2] Galinos C, Larsen T J, Madsen H A and Paulsen U S 2016 Energy Procedia 94 319–328
[3] Madsen H A, Larsen T J, Paulsen U S and Vita L 2013 Proceedings of 51st AIAA Aerospace Sciences Meeting
including the New Horizons Forum and Aerospace Exposition (Chapter AIAA 2013-0913) (AIAA)
[4] Owens B C and Griffiths D T 2014 Journal of Physics: Conferrence Series vol 012092 p 524
[5] Owens B C and Griffith T D 2014 32nd ASME Wind Energy Symposium, AIAA SciTech Forum vol 0518
(AIAA)
[6] Owens B C, Hurtado J E, Paquette J A, Griffith D T and Barone M 2013 54th
AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference vol
1552 (AIAA)
[7] Borg M, Shires A and Collu M 2013 Renewable and Sustainable Energy Reviews 39 1214–1225
[8] Dupont M, Pitance D, Kluczewska-Bordier J, Immas A and Silvert F 2017 WindEurope Conference and
Exhibition (Poster P.O. 262)
[9] Pedersen T F, Paulsen U S, Madesn H A, Nielsen P H, Enevoldsen K, Tesauro A, Kragh K A, Vita L, Ritchie
E, Leban K, Wedell-Heinen J and Larsen K H 2013 European Wind Energy Conference & Exhibition 2013
(European Wind Energy Association (EWEA))
[10] 2018 Dnv gl bladed URL https://www.dnvgl.com/energy/generation/software/bladed/index.html
[11] 2018 Hawc2 URL www.hawc.de
[12] 2018 Nwtc information portal (aerodyn) URL https://nwtc.nrel.gov/AeroDyn
[13] Glauert H 1935 Airplane propellers vol 4 (Springer)
[14] Glauert H 1926 A General Theory of the Autogyro (ARCR R&M No. 1111)
[15] Leishman J G and Beddoes T S 1989 Journal of the American Helicopter Society 34(3) 3–17
[16] Peters D A and He C J 1992 Journal of American Helicopter Society
[17] Ferreira C J S, van Zuijlen A, Bijl H, van Bussel G and van Kuik G 2010 Wind Energy 13 1–17
[18] Ferreira C S, van Kuik G, van Bussel G and Scarano F 2009 Experiments in Fluids 46 97–108
[19] Frounie P, Beguier C, Paraschivoiu I and Brochier G 1986 Journal of Propulsion and Power 5 445–449
[20] Almohammadi K M, Ingham D B, Ma L and Pourkashanian M 2015 Journal of Fluids and Structures 57
144–158
10
Global Wind Summit 2018 IOP Publishing
IOP Conf. Series: Journal of Physics: Conf. Series 1102 (2018)
1234567890 ‘’“” 012002 doi:10.1088/1742-6596/1102/1/012002
11