Observation of the Starting and Low Speed Behavior of Small Horizontal Axis Wind Turbine

1 Department of Mechanical Engineering, King Fahd University of Petroleum &Minerals, Dhahran 31261, Saudi Arabia 2 Institute of Mechatronics Engineering, University of Engineering & Technology, Peshawar 25000, Pakistan 3New York Institute of Technology, New York, NY 10001, USA 4Basic Engineering Department, College of Engineering, University of Dammam, Dammam 31400, Saudi Arabia 5 Center for Engineering Research, King Fahd University of Petroleum &Minerals, Dhahran 31261, Saudi Arabia 6Department of Chemical Engineering, King Fahd University of Petroleum &Minerals, Dhahran 31261, Saudi Arabia 7Department of Electrical Engineering, University of Engineering and Technology, Peshawar 25000, Pakistan


Introduction
From ancient times the kinetic energy of the wind is used for various household activities like milling of wheat and corns.In the 1980s blade element momentum theory was presented.Based on blade element momentum theory various wind turbine designs were proposed for various wind conditions.In areas of high wind speeds large wind turbines can be installed.For areas of low and medium wind speeds small wind turbines are installed.It is not possible to start wind turbine at low wind speeds, so changes should be made in the existing wind turbine models in order to increase their sensitivity [1,2].The wind turbine material should be such that it can withstand the environmental impacts and also it should have a density in a specific range [3].To increase the lift force in order to start the wind turbine at wind speeds various designs must be tested.In developing countries the practical experimentation is not possible so computer simulations are usually carried out to get the optimized design for a specific area [4].Clausen and Wood [5] also worked on improving the starting behavior of horizontal axis wind turbines and they concluded that the root region of the wind turbine blade is responsible for generating the initial torque.Singh and Ahmed [6] practically made various wind turbine models of wood with low Reynolds number and they suggest various models for low and medium wind speed areas.Singh and Ahmed [7] tested various low Reynolds profiles for various low wind speeds and at various pitch angles.Habali and Saleh [8] carried out the static proofload and field performance tests and then they designed a wind turbine model with 41.2% power coefficient.Song and Tan [9] designed 20 kW wind turbine blades.The optimized blade parameters were produced using MATLAB and the dynamic analysis was performed using ANSYS and Solid Works.Wichser and Klink [10] measured wind speed at a height between 70 and 75 meters above the ground level for three different sites of USA.They concluded that wind power production can be increased from 15 to 30 percent if small wind turbines are installed with improved starting behavior.
In this research work wind speed and direction data were taken from PMD (Pakistan Meteorological Department) for various sites of KP (Khyber Pakhtunkhwa).Optimized designs were then proposed for these locations.Optimized designs were made using MATLAB coding, ADAMS simulation, and PRO/E modeling.
The optimized designs were then installed at the selected location of KP (one model for each location made of wood) and the output parameters (torque and rpm) were measured which shows a huge improvement as compared to conventional wind turbine models installed in similar wind speed locations.

Simulation Setup
The main aim of this project was to propose optimized designs for the selected locations of KP.According to the wind speeds and work of Wood [10] sizes were selected for wind turbines in these locations of KP.The sizes and wind speeds of the selected locations are shown in Table 1.
MATLAB coding was done based on the blade element momentum theory.A MATLAB function with the name parameter was used to obtain the blade parameters (chord and twist distribution).The block diagram for this MATLAB function is given in Figure 1.
As shown in Figure 1 the parameter MATLAB function takes the station radii, number of blades, tip speed ratio, coefficient of lift, and angle of attack as inputs and at output we get the chord and twist angle distribution.
Wind turbine models were made in Pro/E environment based on the parameters obtained from MATLAB function.A wind turbine model made in Pro/E is shown in Figure 2.
A second MATLAB function calculates the aerodynamic forces on various stations of wind turbine blades.The block diagram for this function is shown in Figure 3.
The MATLAB function as shown in Figure 3 takes the blade parameters (chord and twist distribution) and wind speed as inputs and after a successive iteration process we get the aerodynamics forces at various stations on wind turbine blades.The Pro/E models were imported into ADAMS environment and the aerodynamic forces were applied on their blades.A wind turbine model with forces applied in ADAMS environment is shown in Figure 4.
Various models were simulated in ADAMS environment and the models with more starting torque were selected.The complete simulation process is explained in Figure 5.

Model Descriptions
The wind turbine model used for simulation has the following blade parameters.The distance between the leading and trailing ends of wind turbine blade is known as chord length ().The distance from the centre of hub to blade tip is termed as blade radius ().The distance from the hub to any station is termed as the station radii ().The width of any station is termed as station width (dr).Each station of the wind turbine blade is twisted at a certain angle with respect to the previous station; this angle is known as the twist angle ().For the same wind speed the aerodynamic effects on various stations of the blade will be different because each station has a different chord length and a different twist angle.According to blade element momentum theory the blade was divided into a number of segments known as stations (in our case 15).The overall effect on wind turbine blade was found out by adding the effects on the individual blade elements.The details of the blade parameters are given in Table 2.

Journal of Wind Energy
The torque output from the simulation was compared with the experimental results.The experiments were carried out in open environment with cross flow of wind on wind turbine rotor.The cross flow of wind caused a slight difference between the experimental and simulation results.The agreement between the output of simulation and experimental results shows us that we can analyze various wind turbine    models using this simulation setup.The comparison is shown in Figure 6 which shows a good agreement between the experimental and the simulation results.
In order to validate the simulation results from ADAMS software, the output torque is compared with the experimental results for the same wind conditions.For a wind turbine blade design modeled by [11,12], the torque comes out to be 5498 Nm.For this wind turbine design the blade diameter was 10 m and the wind speed was selected to be 20 m/sec.A MATLAB function within simulation gives the coefficient of performance for this wind turbine model to be 0.38.The coefficient of performance along with other parameters is used to find out the power from the equation:  =  1/2 2  3 .The next torque is found out from the equation  = /Ω.For the same model and for the same wind conditions ADAMS simulations give the output torque of 5567 Nm.

Results and Discussions
After simulation settings and forces application the initial Pro/e models for the selected locations were simulated.First of all the initial wind turbine designs were imported into ADAMS.After simulating, those initial models graphs from ADAMS postprocessor for each selected location were obtained.The output power of wind turbine depends on both the output rpm of wind turbine and the out torque of wind turbine.In order to study the starting behavior of wind turbines, at least one of these two parameters must be monitored.The torque at the output is monitored in this research work.The focus is to bring changes in wind turbine blade profile so that the output torque increases during the starting time.Various wind turbine models were simulated in ADAMS environment and the output torque was recorded.The output torque for the initial models is shown in Figures 7 and 8.
The wind turbine blade parameters were altered near the hub region and the new wind turbine models were then simulated in ADAMS environment.The output graphs for torque are shown in Figures 9,10,11,and 12. Figures 9-12 show that the starting torque increases by increasing the chord lengths and twist angles near the hub region of horizontal axis wind turbines.
The optimized wind turbine models for the selected locations of KP are shown in Table 3.

Conclusions
Wind turbine designs were proposed for low and medium wind speed locations of Pakistan based on the comparison between simulation and experimental results.Simulations were carried out by interfacing ADAMS and MATLAB software.For each low wind speed location the main aim was to improve the starting behavior of wind turbine models.Wind turbine models were modeled in Pro/E and were then simulated using ADAMS and MATLAB.
Various wind turbine models were modeled in Pro/E by varying chord length and blade angles.The simulation results show that increasing the chord lengths and blade angles near the hub region increases the output torque, angular velocity, and angular acceleration.
Increasing the chord lengths and blade angles near the hub decreases the idling time, so the wind turbine reaches its rated speed in minimum time.

Future Recommendations
Optimized wind turbine models were achieved based on simulation and experimental results.Actual wood models were made based on the chord and blade twist data of optimized designs.These wood models were tested and the output power was more as compared to previous wind turbine models designs for the same wind speeds.The next step of the research is to practically install wind turbine models in the proposed locations.Pitch control, if introduced in these wind turbine models, will help in staring and also help to minimize the rotational speed of wind turbine models during high speed wind.
next step.If no then successive values of a and a is <0.01.If yes, then move to the difference between return to step 1.

Figure 3 :
Figure 3: Logic diagram for the determination of aerodynamic forces.

Figure 4 :
Figure 4: ADAMS model with wind forces applied on its blades.

Figure 5 :
Figure 5: Flow chart for modeling and simulation of wind turbine blades.

Figure 6 :
Figure 6: Comparison between experimental and simulation results.

Figure 7 :
Figure 7: Output torque for the initial design of Peshawar region.

Figure 8 :
Figure 8: Output torque for the initial design of Cherat region.

Figure 9 :
Figure 9: Output torque for the increasing chord design of Cherat region.

Figure 10 :Figure 11 :
Figure 10: Output torque for the decreasing chord design of Cherat region.

Figure 12 :
Figure 12: Output torque for the decreasing blade angle design of Cherat region.

Table 1 :
Size of wind turbines at selected locations of KP.

Table 2 :
Blade parameters for the wind turbine model.

Table 3 :
Optimized wind turbine designs for the selected locations of KP.