Design of a vertical-axis wind turbine: how the aspect ratio affects the turbine’s performance

Abstract

This work analyses the link between the aspect ratio of a vertical-axis straight-bladed (H-Rotor) wind turbine and its performance (power coefficient). The aspect ratio of this particular wind turbine is defined as the ratio between blade length and rotor radius. Since the aspect ratio variations of a vertical-axis wind turbine cause Reynolds number variations, any changes in the power coefficient can also be studied to derive how aspect ratio variations affect turbine performance. Using a calculation code based on the Multiple Stream Tube Model, symmetrical straight-bladed wind turbine performance was evaluated as aspect ratio varied. This numerical analysis highlighted how turbine performance is strongly influenced by the Reynolds number of the rotor blade. From a geometrical point of view, as aspect ratio falls, the Reynolds number rises which improves wind turbine performance.


Highlights

Introduction

There are two types of wind turbine which produce electrical energy from the wind: they are horizontal-axis wind turbines (HAWTs) and vertical-axis wind turbines (VAWTs). The second, and in particular straight-bladed VAWTs, have a simplified geometry with no yaw mechanism or pitch regulation, and have neither twisted nor tapered blades [ 1 ]. VAWTs may be utilized to generate electricity and pump water, as well as in many other applications [ 1 ]. Furthermore, they can handle the wind from any direction regardless of orientation and are inexpensive and quiet [ 2 ]. Wind turbines have aroused the interest of both industry and the academic community [ 315 , 2931 ], which have developed different numerical codes for designing and evaluating wind rotor performance. Recent studies [ 1620 ] have highlighted that VAWTs can achieve significant improvements in efficiency.

VAWT can work even when the wind is very unstable making them suitable for urban and small-scale applications [ 21 ]. Their particular axial symmetry means they can obtain energy where there is high turbulence.

Their optimum operating conditions (maximum power coefficient) depend on rotor solidity and tip speed ratio [ 22 ]. For a VAWT rotor solidity depends on the number of blades, airfoil chord and rotor radius. Tip speed ratio is a function of angular velocity, undisturbed wind speed and rotor radius.

In the design process of a vertical-axis wind turbine it is crucial to maximize the aerodynamic performance [ 22 , 26 ]. The aim is to maximize the annual energy production by optimizing the curve of the power coefficient varying with the tip speed ratio [ 25 ]. For a fixed cross-sectional area of the turbine, to optimize the curve of the power coefficient it is possible to use different airfoil sections and/or rotors with different solidity [ 26 ].

To maximize energy extraction, other authors introduced guide vanes [ 27 ] and/or blade with a variable pitch angle [ 28 ] in vertical-axis wind turbines.

In the design process of a vertical-axis wind turbine, a wrong choice of the aspect ratio of the wind turbine may cause a low value of the power coefficient (wind turbine efficiency). This parameter (the aspect ratio) is often chosen empirically on the basis of the experience of the designer, and not on scientific considerations.

In this work, the link between the aspect ratio of a wind turbine and its performance has been studied, and a correlation between the aspect ratio and the turbine’s performance has been found.

Designing an H-Rotor

Designing a vertical-axis wind turbine with straight blades requires plotting power coefficient c p against tip speed ratio λ , as a function of rotor solidity σ (Fig.  1 ).

Fig. 1

Power coefficient for a VAWT, straight blades and symmetric airfoil

Figure  1 shows the behaviour of the power coefficient for a wind turbine with straight blades and a NACA 0018 airfoil.

Figure  1 curves were obtained using a calculation code based on MSTM theory [ 23 ].

From the graph in Fig.  1 , the solidity which maximizes power coefficient σ  = 0.3 can be identified, which has a c pmax  = 0.51 corresponding to λ  = 3.0.

Since solidity σ equals: chord c can be expressed as a function of solidity, rotor radius and blade number N b , as per Eq.  2 :

The power of a wind turbine with a vertical axis can be expressed as per Eq.  3 :

Having defined the turbine’s aspect ratio ( AR ) as the ratio between blade height and rotor radius ( AR  =  h/R ), rotor radius can be derived from Eq.  3 :

(in Eq.  4 power P and wind velocity V 0 are design data and ρ is air volume mass).

This design approach is iterative and from time to time it will be necessary to re-evaluate the blade’s Reynolds number and if necessary repeat the procedure with new power coefficient curves.

The local Reynolds number is: where c is the chord from Eq.  2 , ν is the kinematic air viscosity, and w is the air speed relative to the airfoil as Fig.  2 shows.

Fig. 2

Wind rotor rotational plane

Adopting a mathematical approximation, to evaluate the Reynolds number, w can be substituted by ωR with the advantage of having a mean Reynolds number independent of the angle of rotation (see Fig.  2 ).

To conclude the design cycle, simply calculate ωR directly from TSR relative to c pmax identifiable in Fig.  1 , and combining Eqs.  5 and 6 :

If the Reynolds number thus calculated is different to the one for the power coefficient curve adopted initially (Fig.  1 ), a new power coefficient curve should be plotted for a different Reynolds number (second attempt). Usually, the iterative design process needs only 2 or 3 iterations.

Figure  3 shows the power coefficient curves for the wind turbine with the NACA 0018 airfoil, at high Reynolds numbers.

Fig. 3

Characteristic curves for high Reynolds numbers

In conclusion, the Reynolds number strongly influences the power coefficient of a vertical-axis wind turbine. Furthermore, it changes as the main dimensions of the turbine rotor change. Increasing rotor diameter rises the Reynolds number of the blade.

The importance of aspect ratio

In Eq.  4 note how radius R increases as ratio AR decreases. In Eq.  2 if R increases, chord c increases too, and in Eq.  7 note how increasing the chord rises the Reynolds number. Finally, Fig.  3 shows how the power coefficient increases as the Reynolds number rises.

Moreover, the rotational velocity ω can be derived from Eq.  6 :

Equation ( 8 ) shows how ω is inversely proportional to R . From the Fig.  3 graph, note how λ cpmax decreases as Reynolds number rises.

So, to maximize the power coefficient, the rotor’s aspect ratio should be as small as possible. As aspect ratio diminishes there are two advantages: the local Reynolds number rises and simultaneously the rotational velocity diminishes. The effect of the rotor’s aspect ratio on Reynolds number and rotational velocity are shown in Fig.  4 for a twin-bladed 1 kW turbine in a wind velocity of 10 m/s.

Fig. 4

Effect of aspect ratio ( h/R ) on VAWT performance

Figure  5 shows two vertical-axis turbines with identical design power, blade number and aerodynamic profile (NACA 0018) but with two different aspect ratios ( AR 1  = 2; AR 2  = 0.4). As stated above, the turbine with the lowest AR will have the highest power coefficient and the lowest rotational velocity. This turbine will display two further advantages: firstly, a structural one (thicker blades are more stress resistant) and secondly, in-service stability (greater rotor inertia).

Fig. 5

Wind turbines with different aspect ratios

Figure  6 shows the same graphs as Fig.  4 but parameterized for design power variation. Note how the turbines with the highest power have higher Reynolds numbers and lower rotational velocity.

Fig. 6

How aspect ratio influences Reynolds number and rotational velocity, for different design powers

To conclude, VAWT performance is directly related to the blade’s Reynolds number. The higher this is, the better the turbine’s performance. This is because as the Reynolds number increases, the blade’s lift coefficient [ 24 ] rises and the drag coefficients decrease (Fig.  7 a, b) thus providing greater torque.

Fig. 7

a Lift coefficient for NACA 0018. b Drag coefficient for NACA 0018

Application for a case study

To better explain the design procedure investigated in this work, a case study is presented. From Fig.  4 , it is possible to notice that the lower the AR , the higher the Reynolds number; this can improve wind turbine performance. In the design procedure, a choice of a low value for the AR is therefore suitable. In this case study, a comparison between two wind rotors, designed with two different AR values ( AR  = 2 and AR  = 0.4), will be presented.

Suppose, for simplicity, we want to design a VAWT with a rated power of 1 kW, with only two straight blades with symmetrical aerodynamic profile (NACA 0018), and that the installation site is characterized from a mean wind speed of 10 m/s. The design procedure is showed in the Fig.  8 . In this figure a flowchart for designing the wind turbine is adopted. In this design procedure, first of all, a Reynolds number must be supposed. This Reynolds number will be called: “the first attempt Reynolds number”. Supposed a first attempt Reynolds number equal to 5 × 10 6 , from the graph in Fig.  3 a it is possible to read the highest performance ( c pmax  = 0.51) in correspondence of a particular rotor solidity and tip speed ratio: σ cpmax  = 0.3 e λ cpmax  = 3.0. For an aspect ratio of 2 ( h / R  = 2) (square VAWT cross section), Eq.  4 will provide the rotor radius R  = 0.904 m ( ρ  = 1.2 kg/m 3 ). From Eq.  2 , we obtain a chord c  = 0.136 m, and from Eq.  8 we obtain a rotational velocity of 317 rpm. Now a new Reynolds number (second attempt) must be evaluated from Eq.  7 : Re  = 2.8 × 10 5 .

Fig. 8

VAWT design flowchart

Now the design procedure should be repeated with the new value of the Reynolds number. This Reynolds number will be called: “the second attempt Reynolds number”. Taking into account Fig.  3 c, we will obtain the new geometric values for the wind turbine. These values are reported in Table  1 . The design procedure will end when the Reynolds number does not change in a meaningful manner. The design procedure goes to numerical convergence with only two iterations.

Table 1

VAWT––straight blade’s design procedure ( h/R  = 2)

Reference

1st attempt

2nd attempt

Power (kW)

1

1

Airfoil

NACA 0018

NACA 0018

Wind speed (m/s)

10

10

Air density (kg/m3)

1.2

1.2

Kinematic air viscosity (m2/s)

1.46 × 10−5

1.46 × 10−5

Rotor aspect ratio (h/R)

2

2

Number of blades (Nb)

2

2

First attempt Reynolds

5 × 106

2.8 × 105

cpmax; σcpmax; λcpmax

Fig. 3

0.51; 0.3; 3.0 (Fig. 3a)

0.464; 0.4; 2.96 (Fig. 3c)

Rotor radius (m)

Eq. 3

0.904

0.947

Airfoil Chord (m)

Eq. 4

0.136

0.189

Rotational speed (rpm)

Eq. 6

317

299

Second attempt Reynolds

Eq. 7

2.8 × 105

3.8 × 105

Next step

Go to the 2° attempt

END

From the data reported in Table  1 for h/R  = 2, the design procedure gives rise to a power coefficient of 0.464. Repeating the procedure described above for h/R  = 0.4, we will obtain an higher efficiency: c pmax  = 0.475 (see Table  2 )

Table 2

VAWT––straight blade’s design procedure ( h/R  = 0.4)

Reference

1st attempt

2nd attempt

Power (kW)

1

1

Airfoil

NACA 0018

NACA 0018

Wind speed (m/s)

10

10

Air density (kg/m3)

1.2

1.2

Kinematic air viscosity (m2/s)

1.46 × 10−5

1.46 × 10−5

Rotor aspect ratio (h/R)

0.4

0.4

Number of blades (Nb)

2

2

First attempt Reynolds

5 × 106

6.2 × 105

cpmax; σcpmax; λcpmax

Fig. 3

0.51; 0.3; 3.0 (Fig. 3a)

0.475; 0.3; 3.01 (Fig. 3c)

Rotor radius (m)

Eq. 3

2.021

2.094

Airfoil Chord (m)

Eq. 4

0.303

0.314

Rotational speed (rpm)

Eq. 6

142

137

Second attempt Reynolds

Eq. 7

6.2 × 105

6.5 × 105

Next step

Go to the 2° attempt

END

Conclusion

This work looks at designing a vertical-axis wind turbine to maximize its power coefficient. It has been seen that the power coefficient of a wind turbine increases as the blade’s Reynolds number rises. Using a calculation code based on the Multiple Stream Tube Model, it was highlighted that the power coefficient is influenced by both rotor solidity and Reynolds number.

By analysing the factors which influence the Reynolds number, it was found that the ratio between blade height and rotor radius (aspect ratio) influences the Reynolds number and as a consequence the power coefficient.

It has been highlighted that a turbine with a lower aspect ratio has several advantages over one with a higher value.

The advantages of a turbine with a lower aspect ratio are: higher power coefficients, a structural advantage by having a thicker blade (less height and greater chord), greater in-service stability from the greater inertia moment of the turbine rotor.


References

  1. Islam et al. (2008) Aerodynamic models for Darrieus-type straight-bladed vertical axis wind turbines (pp. 1087-1109) 10.1016/j.rser.2006.10.023
  2. Pope et al. (2011) Power correlation for vertical axis wind turbines with varying geometries (pp. 423-435) 10.1002/er.1703
  3. Eriksson et al. (2008) Evaluation of different turbine concepts for wind power (pp. 1419-1434) 10.1016/j.rser.2006.05.017
  4. Riegler (2003) (pp. 44-46) Elsevier Refocus
  5. Rohatgi and Barbezier (1999) Wind turbulence and atmospheric stability—their effect on wind turbine output (pp. 908-911) 10.1016/S0960-1481(98)00308-5
  6. Castro et al. (2006) Turbulence over urban-type roughness: deductions from wind-tunnel measurements (pp. 109-131) 10.1007/s10546-005-5747-7
  7. Gayev and Savory (1999) Influence of street obstructions on flow processes within urban canyons (pp. 89-103) 10.1016/S0167-6105(98)00212-8
  8. Akhmatov (2007) Influence of wind direction on intense power fluctuations in large offshore windfarms in the North Sea 31(1) (pp. 59-64) 10.1260/030952407780811384
  9. Sahin et al. (2006) Thermodynamic analysis of wind energy (pp. 553-566) 10.1002/er.1163
  10. Wright and Wood (2004) The starting and low wind speed behaviour of a small horizontal axis wind turbine (pp. 1265-1279) 10.1016/j.jweia.2004.08.003
  11. Unknown ()
  12. Unknown (2003) The CREATE Centre
  13. Biswas et al. (2007) Experimental investigation of overlap and blockage effects on three-bucket Savonius rotors 31(5) (pp. 313-368) 10.1260/030952407783418702
  14. Bontempo and Manna (2013) Solution of the flow over a non-uniform heavily loaded ducted actuator disk (pp. 1469-7645) 10.1017/jfm.2013.257
  15. Bontempo et al. (2014) Ducted propeller flow analysis by means of a generalized actuator disk model (pp. 1107-1115) 10.1016/j.egypro.2014.01.116
  16. Saha (2008) Optimum design configuration of Savonius rotor through wind tunnel experiments 96(8) (pp. 1359-1375) 10.1016/j.jweia.2008.03.005
  17. Hiraharaa et al. (2005) Testing basic performance of a very small wind turbine designed for multi-purposes (pp. 1279-1297) 10.1016/j.renene.2004.10.009
  18. Wakui et al. (2005) Hybrid configuration of Darrieus and Savonius rotors for stand-alone wind turbine-generator systems 150(4) (pp. 13-22) 10.1002/eej.20071
  19. Gupta et al. (2008) Comparative study of a three-bucket Savonius rotor with a combined three-bucket Savonius––three-bladed Darrieus rotor (pp. 1974-1981) 10.1016/j.renene.2007.12.008
  20. Unknown ()
  21. Armstrong et al. (2012) Flow separation on a high Reynolds number, high solidity vertical axis wind turbine with straight and canted blades and canted blades with fences (pp. 13-22) 10.1016/j.renene.2011.09.002
  22. Paraschivoiu (2002) Polytechnic International
  23. Unknown ()
  24. Unknown ()
  25. Bhutta et al. (2012) Vertical axis wind turbine––a review of various configurations and design techniques (pp. 1926-1939) 10.1016/j.rser.2011.12.004
  26. Hameed and Afaq (2013) Design and analysis of a straight bladed vertical axis wind turbine blade using analytical and numerical techniques (pp. 248-255) 10.1016/j.oceaneng.2012.09.007
  27. Chong et al. (2013) The design, simulation and testing of an urban vertical axis wind turbine with the omni-direction-guide-vane (pp. 601-609) 10.1016/j.apenergy.2012.12.064
  28. Unknown ()
  29. Ajayi et al. (2013) Wind profile characteristics and turbine performance analysis in Kano, north-western Nigeria 10.1186/2251-6832-4-27
  30. Morshed et al. (2013) Wind tunnel testing and numerical simulation on aerodynamic performance of a three-bladed Savonius wind turbine 10.1186/2251-6832-4-18
  31. Liu and Janajreh (2012) Development and application of an improved blade element momentum method model on horizontal axis wind turbines 10.1186/2251-6832-3-30