# Aerodynamic Design and Performance Prediction of a Micro Horizontal Axis Wind Turbine Blade

Text Only Version

#### Aerodynamic Design and Performance Prediction of a Micro Horizontal Axis Wind Turbine Blade

M. Shuwa*, G. M. Ngala, A. M. El- Jummah and F. M. Kajiama

Department of Mechanical Engineering, University of Maiduguri, Maiduguri, Nigeria

Keywords: Horizontal axis wind turbine, Blade design, Performance prediction.

NOMENCLATURE

– The Available Power in the wind

– Extractable or Mechanical Power

A – Swept Area of the Rotor

– Air Density

v – Wind Relative Velocity

C – Chord Length

– Angular Velocity of the Rotor

– Angle of Attack

– Wind Relative Angle

– Efficiency of the Blade on the Wind Turbine

d – Increase in Blade Span Length

– Rotor Torque

r – Tip Speed Ratio

CP – Power Coefficient

CL – Coefficient of Lift

CD – Coefficient of Drag

FL – Lift Force

FD – Drag Force

Re – Reynolds Number

N – Speed of the Rotor

TN – Thrust

Vrel – Wind Relative Velocity

dt – Change in Blade Thickness

– Meter

mm – Millimeter

– Watts

– Second

1.0 INTRODUCTION

Wind turbines are rotor dynamic machines that convert some of the kinetic energy in the wind to other useful forms of energy such as electricity. There are two types of wind turbines, horizontal and vertical axis wind turbines. Horizontal Axis Wind Turbines (HAWT) have the axis of rotation of their blades horizontal to the ground and almost parallel to the wind stream, while the Vertical Axis Wind Turbines (VAWT) have the axis of rotation and blades vertically oriented (Ngo and Natowitz, 2009). Most of the wind turbines fall under the category the Horizontal Axis Wind Turbines because they have some distinct advantages such as low cut-in wind speed and easy furling. In general, they show relatively high power coefficient (Rajakumar and Ravindran, 2010; Chandrala et al, 2012).

2.0 METHODOLOGY

The baseline for the aerodynamic design of the HAWT blade is the National Advisory Committee for Aeronautics (NACA) 4412 airfoil data. Base on the airfoil data, a modified blade geometry and orientation is aerodynamically design. The aerodynamic design is carried out using the Blade Element Momentum theory. In determining aerodynamic design features of the blade, the blade is divided in to 12 sections and each section is considered independently as two-dimensional airfoil. Aerodynamic design features of the blade airfoil are determined in each section and then summed to obtain the overall blade parameters which determine the blades ability to extract energy from the wind. The maximum available power in the wind increases with the cube of the wind speed, linear with density and swept area of the rotor. Theavailable energy in the wind based on a predefine rotor diameter is determined from Equation 1 (Polat and Tuncer, 2013).

= 0.53 (1)

The available wind speed at a given site is therefore often first measured over a period of time before a blade aerodynamic design is initiated. Aerodynamic design and performance prediction of the HAWT blade is based on correlated 30 years annual monthly mean wind speed data of the semi-arid region of Nigeria as obtained from Nigerian Meteorological Agency. However a physical limit exits to the quantity of energy that the blade can extract from the wind and it is a factor that depends on the blade airfoil design and its geometry. Complete extraction of the kinetic energy in the wind would imply zero velocity downwind, but this scenario is not possible. Thus all the kinetic energy in the wind cannot be extracted and this limitation is known as the Bezt limit which state that wind turbine efficiency cannot exceed 59.3%. Therefore, based on the aerodynamic design of the blade, the predicted extracted or mechanical power () to be generated by each of the blade section on a three blade rotor is determined from Equation 2 (Peter and Richard, 2012; Shekh and Salman, 2015).

= (2)

The efficiency of the blade depends on the ratio of the available power in the wind to the extracted or mechanical power generated by the blade and is called the power coefficient () determined from Equation 3. But all the kinetic energy in the wind cannot be extracted based on the limitation known as the Bezt limit, which state that wind turbine efficiency cannot exceed 59.3%. The coefficient of performance of the blade will only be 0.593 of its maximum efficiency (Kale and Varma, 2014; Shekh and Salman, 2015).

= = 0.593

(3)

0.53

The coefficient performance of a HAWT blade depends on its aerodynamic design and one of the important datum design parameter is the tip speed ratio. It is the relationship between rotor blade velocity and the velocity of the relative wind as in Equation 4. It is also a design parameter that is used in determining the relative wind angle () and angle of attack () at a predefined rotor radius (Peter and Richard, 2012; Bai et al, 2013).

=

=

2

(4)

60

The distribution of the angle of relative wind () which consists of angle of attack () and twist angle () along the blade sections are determined from Equation 5, 6 and 7 respectively (Peter and Richard, 2012; Shekh and Salman, 2015).

2 1

= 1 (

) (5)

3

2

= 1

3

()

(6)

= (7)

The angle of attack is the angle between the incoming flow stream and the chord line of the airfoil as in Figure 1. It is a required input in Design Foil Workshop software to find the values of the lift coefficient () and coefficient of drag () of the airfoil at a predetermined Reynolds number (). Coefficient of Lift to coefficient Drag ratio is calculated for angles of attack from 0o to 12o for wind speed from 0 to 10 m/s. Lift and drag coefficients are factors in determining the magnitude of lift () and drag () force acting on the airfoil section as shown in Figure1. Moreover other forces acting on the airfoil section as in Figure 1, where thrust ( ) defined by Equation 10 and torque ( ) define by Equation 11, both are generated by lift () and drag () forces (Shekh and Salman, 2015; Shuwa et al, 2016).

Figure 1 Airfoil Parameters

The lift force () which is the force acting on the blade airfoil perpendicular to the undisturbed wind flow, is determined from Equation 8. The drag force () which is the force acting on the blade airfoil in the direction of the undisturbed wind flow as given by Equation 9.

1

= 2

(8)

2

1

= 2

(9)

2

The lift and drag are available as functions of the angle of attack and the Reynolds number (Re). Also, the magnitude lift and drag force is influence by the blades angle of attack and the Reynolds number (Re). Reynolds number (Re) as given by Equation 10 is the ratio of inertial forces to viscous forces within a fluid which is subjected to relative internal movement due to different fluid velocity.

=

(10)

The chord length () of a blade airfoil is the length from the leading edge to the trailing edge of a blade cross section that is parallel to the vertical axis of symmetry as shown in Figure 1, and it is obtained from Equation 11.

=

8

(1 ) (11)

The blade thickness () along its span is determined as percentage function of the chord length of the airfoil. Base the selected airfoil NACA 4412 the change in thickness of the blade () along the blade span is 12 % of the chord length of each section or divided elements of the blade () as in Equation 12.

= 0.12 (12)

The thrust () and torque ( ) on each of the blade sections or elements are calculated from Equation 13 and 14 respectively and determines the magnitude of the power extracted by the blade. This based on the consideration of the mean axial and rotational velocity in the slipstream and the rate of decrease of momentum of the fluid.

= + (13)

= (14)

3.0 RESULTS AND DISCUSSION

Figure 2 Blade airfoil thickness and twist angle along the blade span

The blade twist angle is determined based on a wind relative angle of 50.05 degrees at the root of the blade to 9.36 degrees at the tip with a tip speed ratio 5. The aerodynamic design shows an optimum angle of attack of 7 degrees with an optimum coefficient of lift to drag ratio of 124 at Reynolds number () 3 x 106 and lift to drag force ratio of 118 as shown in Figure 3.

Figure 3 Lift to drag force ratio and coefficient of lift to coefficient of drag ratio at various angles of attack

Rajakumar and Ravindran, 2010 in their work analysis lift and drag forces at different section of a blade for angles of attack from 0 to 12 degrees. Their analysis showed that angle of attack of 5 degrees has high lift to drag force ratio. Biadgo and Aynekulu, 2017 developed an interface prgram to pridict the performance of a HAWT blade. Based on the program optimum angle of attack for small wind turbine blade lies between 4 to 9 degrees depending on the blade geometery and airfoil design, hence 7 degrees anagle of attack of this design is withein the predicted limit. Performance analysis of the blade determined interms of its ability to etract energy from the wind for wind speeds from 0 to 10 m/s based on a 3 m micro HAWT rotor. The result shows that the blade extract 32 percent (32%) of the power available in the wind. Maximum pridicted power available in the wind at 9.5 m/s is 320 Watts and the maximum power extracted by the blade considering a rotor diameter of 3 m is 102.4 Watts as in Figure 4.

Figure 4 Available power in the wind and extracted power by the bladesPeter and Richard, 2012 indicates that physical limit exist to the quantity of energy that can be extracted which is dependant on the blade design. The magnitude of energy harnessed is a function of the reduction in wind speed over the turbine. This priciple is also reported by Kale and Varma, 2014 and is known as the Betz limit which state that the efficiency of a wind turbine can not exceed 59.3%.

4.0 CONCLUSION

The micro HAWT blade has been aerodynamically design for a 3 m rotor and a wind speed range from 0 to 10 m/s. The blade twist angle is 50.05 degrees at the root and 9.36 degrees at the tip. The active blade span is 1,200 mm (1.2 m), maximum chord length is 180 mm (0.18 m) and the minimum chord length is 90 mm (0.09 m). The blade tip speed ratio is 5, optimum angle of attack is 7 degrees with an optimum coefficient of lift to drag ratio of 124 at Reynolds number () 3 x 106. The designed blade is analysis for performance, the analysis lift to drag force ratio is 118 and the blade is able to extract 32 % of the available power in the wind. However, the micro HAWT blade can further be modified and experimentally analysis for improved performance.

REFERENCES

1. Ahmed, M. R., Narayam, S., Zullah, M. A. and Lee, Y. H. (2011). Experimental and Numerical Studies on a Low Reynords Number Airfoil for Wind Turbine Blades. Journal of Fluid Science and Technology, 6: 357-371.

2. Bai, C., Hsiao, F., Li, M., Huang, G. and Chen, Y. (2013). Design of 10 kW Horizontal Axis Wind Turbine (HAWT) Blade and Aerodynamics Investigation Using Numerical Simulation. Elselvier Publishers Ltd., 5 (2): 34-40.

3. Biadgo, A. M. and Aynekulu, G. (2017). Aerodynamic Design of Horizontal Axis Wind Turbine Blades. FME Transactions , 45:647-660.

4. Chandrala, M., Choubeg, A. and Gupta, B. (2012). Aerodynamic Analysis of Horizontal Axis Wind Turbine Blade. International Journal of Engineering Research and Applications, 2(6):1244-1248.

5. Daminia, S. and Petal, H. (2013). Performance Prediction of Horizontal Axis Wind Turbine Blade. International Journal of Innovative Research in Science, Engineering and Technology, 2(5): 1401-1406.

6. Kale, A.S. and Varma, R. N. (2014). Aerodynamic Design of a Horizontal Axis Micro Wind Turbine Blade Using NACA 4412 Profile. International Journal of Renewable Energy Research, 4(1):69-72.

7. Ngo, C. and Natowitz, J. (2009). Resources, Alternatives and the Environment. New Jersey, USA: John Wiley and Sons Inc.

8. NIMET. (2017). Meteorological Data. Lagos: Nigeria Meteorological Agency.

9. Peter, S. and Richard, C. (2012). Wind Turbine Blade Design. Energies, 5: 3425-3449.

10. Polat, O. and Tuncer, H. I. (2013). Aerodynamic Shape Optimization of Wind Turbine Blades using a Parallel Genetic Algorithm. Elsevier Ltd, 61:28- 31.

11. Rajakumar, S. and Ravindran, D. (2010). Computational Fluid Dynamics of Wind Turbine Blade at Various Angle of Attack and Low Reynolds Number. International Journal of Engineering Science and Technology, 2: 11-17.

12. Shekh, M. and Salman, M. (2015). Design of Wind Energy Turbine Blades: A Review. International Advance Research Journal in Secience, Engineering and Technology, 2(1):193-197.

13. Shuwa, M., Ngala, G. M and Maina, B. M. (2016). Development and Performance Test of a Micro Horizontal Axis Wind Turbine Blade. International Journal of Engineering Resaerch and Applications, 6(2): 11-17.