﻿ Dependency of Torque on Aerofoilcamber Variation in Vertical Axis Wind Turbine

World Journal of Mechanics
Vol.06 No.11(2016), Article ID:72415,15 pages
10.4236/wjm.2016.611033

Dependency of Torque on Aerofoilcamber Variation in Vertical Axis Wind Turbine

1Department of Mechatronic Engineering, PAF Karachi Institute of Information Technology, Sahiwal, Pakistan

2Department of Mechanical Engineering, COMSATS Institute of Information Technology, Islamabad, Pakistan

3Department of Mechanical, Civil and Aerospace Engineering, Manchester, UK

Copyright © 2016 by authors and Scientific Research Publishing Inc.

Received: October 20, 2016; Accepted: November 27, 2016; Published: November 30, 2016

ABSTRACT

The advancements in the wind turbine technology specially associated with Vertical Axis Wind Turbines (VAWT), has been improved for last couple of years. This is due to extensive use of computational techniques. This paper investigates dependency of torque on aerofoil geometry by performing numerical simulation on Darrieustype VAWT with fix pitch blades. Coordinate points for aerofoil was generated using Java Foil software. Reynolds-Averaged Navier Stokes (RANS) turbulence modelling was used for predicting the flow and efficiency of the three blades VAWT. The unsteady flow condition was considered to make simulation as realistic as possible. In order to visualize high strain flow and separation, we used two equation models i.e. k-ε with RNG. NACA 0012 aerofoil was used and camber variations were carried out for developing samples of aerofoil to check the enhancement in performance of VAWT. Results demonstrate the torque and power along with its coefficients. It has been investigated that the performance efficiency was significantly improved by changing the aerofoil camber, demonstrating highest torque with camber (C3) aerofoil and the least performance was observed using camber (C0).

Keywords:

Co-Efficient of Moment, CFD, Torque, Power, Fluent, Java Foil

1. Introduction

2. Modeling and Simulation

2.1. Geometric Modeling

Geometry parameters were considered from Biadgo, Simonović [2] study. This study consist of wind turbine with three blades as it provides better efficiency [7] [12] [14] . VAWT model was made by importing Java Foil generated aerofoil to CATIA for generation of 2-D CAD model. Figure 1 represents aerofoil geometry of all samples whereas Figure 2 illustrates simulation domain and turbine rotor design.

Figure 1. Aerofoil models (a), (b), (c), (d) and (e) with coordinate.

Figure 2. Geometry Sketch and description.

After designing model of wind turbine an enclosure was made around the aerofoil. Figure 2 shows the geometry. It has three aerofoils with 200 mm chord length and they are part of VAWT with 2 m rotor radius. The design of Darrieus-type straight-bladed vertical axis turbines were studied using different models [15] [16] [17] . However, this research paper provides results based on six VAWT samples each with different cambered aerofoil; the design of aerofoil were varied from symmetric to asymmetric. Six wind turbines 2D models were generated using NACA 0012. Samples were generated by changing cambers with keeping other parameters of aerofoil fixed. After changing camber percentage values from 0 to 1, 3, 5, 7 aerofoils with code NACA 1412, NACA 3412, NACA 5412 and NACA 7412 respectively, were used for comparison of performance parameters. After generation of CAD geometry file, it was exported to mesh module of ANSYS.

2.2. CFD Methodology

Mesh

CFD analysis of VAWT was performed on ANSYS Fluent module. Triangular type mesh element with 108,384 number of nodes and 141,154 number of elements were used for meshing as it is considered more appropriate for 2D model [18] . It can be observed from the Figure 3 that mesh near the VAWT is fine for achieving higher accuracy solution near the wind turbine blades, for this reason minimum element size was used i.e. 0.0035 m. Large domain was used for smooth flow conditions and minimizing the disturbance created by back flow from wall.

2.3. Theoretical Formulation and Governing Equations Discretization

In CFD we want to compute velocities and pressure at different points in model. For

Figure 3. Mesh generated cells.

this we consider a control volume and by using conservation of mass and momentum principle we develop N-S equations as sown by following equations [19]

(1)

(2)

(3)

(4)

where represents density, first viscosity coefficient, second viscosity coefficient and are velocity and body force components. N-S equations cannot be solved as it have no analytical solution possible due to non-linear coupled partial differential equations. So we use numerical method in which we find solution at discrete points and use mathematical technique which in case of Fluent is Finite Volume Method (FVM) in this technique calculation is performed at the center of face in 2-D geometric case.

If we solve N-S equations numerically, method is known as DNS which is difficult to compute. To overcome this difficulty we compute average solution of actual solution and this is done by using RANS model in which value of any variable in turbulent flow can be found by sum of mean value plus fluctuating value. N-S equations solved after averaging operation is known as RANS equations as shown in Equation (5).

(5)

where is Reynolds stress tensor computed by transport equations. Standard k-ε model compute turbulence energy using the transport equation as in Equation (6).

(6)

where, is mean component of velocity, is dissipation rate, k is turbulence energy reference and is fluctuating component of velocity.

Figure 4 indicates a general overview of RANS approach through the family tree. In this study K-ε model RNG turbulence model was used for N-S equations normalization because of its capability of predicting flow separation, wind velocity and swirl [20] . For solution simple pressure-velocity coupling with Green-Gauss Cell Based gradient discretization was used. Furthermore, second order spatial discretization approach was used for pressure, momentum, turbulent kinetic energy and dissipation rate. Mesh motion was implemented for the blade region whereas outer and inner domain was kept stationary.

2.4. Boundary Conditions

Air is used as compressible fluid flowing around the wind turbine for power generation. Velocity of fluid at inlet was set to 5 m/s and symmetry wall was set to no slip condition. Internal rotor and outer domain was set stationary. Ring containing aerofoil was considered movable with 0.5 and 4 tip speed ratio. Monitors for moment, coefficients of drag (CD), lift (CL) and moments (CM) for individual as well as combined blades was established. Figure 5 provides a visual aid about the data collection procedure from geometry to contour results progression and then data compilation. These steps were followed throughout the study for all the wind turbine samples.

Figure 4. Scheme of general RANS approach.

Figure 5. Procedure of data collection.

2.5. Post Processing

Data for TSR 0.5 was compared with Biadgo, Simonović [2] as a standard and for TSR 4 data was collected directly from monitors. Values for CD, CL, CM and torque was collected for one complete revolution after periodic behavior of residuals was achieved.

2.6. Data Acquisition

Values of lift, drag and co-efficient of moment at different azimuth angles during one revolution collected from fluent was exported to excel. Performance of the wind turbine can be predicted using power, torque and there co-efficient which depends upon area of the wind turbine, velocity and density of the fluid. Torque was calculated using co-effi- cient of torque or moment formula shown in Equation (8) [21] , whereas power was calculated by multiplying torque with velocity. Co-efficient of torque was calculated directly through fluent and forco-efficient of power Equation (7) was used [21] . Angular velocity was calculated using TSR Equation (9) [21] .

(7)

(8)

(9)

where CP is the co-efficient of power, CM is the co-efficient of moment, A is area of the wind turbine, is the density of fluid and R is rotor radius. In fluent, coefficient of toque CT and coefficient of moment CM is same and can be calculated from same Equation (8) by replacing M with T.

3. Results

Results were collected right after seven revolutions were carried out. By doing this residuals were periodic and convergence criteria was met. Value of CM was obtained using the Fluent directly for all VAWT models for one revolution out of 7 revolutions. Comparison of aerofoils at different cambers were made and it can be observed from the Figure 6(a) between CM and azimuth angle that CM for camber 3 approached the highest value in the graph and CM for camber value 1 is in the lowest point in case of TSR 4.

Torque was calculated by utilizing CM values collected from fluent multiplying with available wind power as in Equation (4). It can be seen from Equation (4) that torque and CM is dependent upon each other so there graph trend will be same but values will

Figure 6. Fluctuation with rotation angles (a) Co-efficient moment and (b) Torque.

be different. Therefore, it can be visualized in Figure 6(b) that aerofoil NACA 0012 with camber value of 3% is more suitable and will provide best torque as compared to other cambered aerofoils.

Before attempting these simulations, initially the data of the torque of Camber 0 aerofoil with TSR 0.5 was compared with Biadgo, Simonović [2] results of similar camber and TSR values, this was performed as a control simulation. By doing a control simulation, it provides a validation that further prediction are reasonably correct. The comparison data has been flourished in Table 1.

Fluctuation of torque and CM with azimuth angle provide good approximations as viewed from Figure 6(a) and Figure 6(b), that fluctuation is less in case of CL and is maximum in case of C3 aerofoil.

Power values cannot be directly collected from software so it was calculated using torque value from Figure 6(b) numerical data and the angular velocity of the wind turbine. Graph for power shown in Figure 7(a) was generated using torque values of one revolution. It can be observed from the power graph that max power achieved is 130 watt for the case of TSR 4 by C3 cambered aerofoil VAWT model whereas for TSR value 0.5 maximum power is approximately 8 watt. This tells that TSR value changes the power output of wind turbine. CP graph provides a good observation of performance of aerofoil which will provides ease in selecting an aerofoil for VAWT case. Figure 7(b) represents co-efficient of power variations with azimuth angles which represent energy conversion efficiency by wind attacking on aerofoil.

4. Discussion

For each camber, CM data was collected directly from the Fluent and a mean value was used for comparison. Camber effect on CM and torque can be seen from Figure 8(a) and Figure 8(b) respectively. It can be observed that for camber 0 to 3, CM value shows increment however as the camber was increased beyond 3, the value of Cm decreased. Graphs indicate the similar trend for both Cm and torque that is due to dependency on each other both have identical reason for increase and decrease.

The bar chart in Figure 8(a) illustrates the variation in CM by changing the camber percentage from zero to seven. It can be observed that highest CM obtained is 0.089 for the case of 3% camber, whereas, the minimum CM value is 0.031 for 0% camber. Furthermore, chart illustrates that first trend of CM rise fast till 3% and then begin to fall. Chart in Figure 8(b) depicts that change in torque value with different camber value. It

Table 1. Percentage error of torque between the referenced result and this study.

Figure 7. Fluctuation with rotation angles; (a) Power and (b) Coefficient of power.

can be seen that value of torque was 4 Nm in case of symmetric aerofoil whereas by varying camber it reaches 10.5 Nm for 3% camber, which is considered optimum. The slight decrease in the trend was encountered by increasing camber percentage further. The torque average values from camber 0 to 7 for TSR 4.0 demonstrate to be within the range as predicted by Beri and Yao [6] which is between −20 Nm to 20 Nm. They are getting negative torque as they are using a non-self-startingaerofoil i.e. NACA 2415, this is also suggested by Micha Premkumar, Seralathan [22] .

The chart in Figure 9(a) depicts the power generation amount at different cambers. It can be visualized that the bar graph peaked from 38 W to 104 W by changing the camber percentage values from 0 to 3. After camber value 3% there is a gradual decline

Figure 8. Co-efficient variation by camber; (a) Change in moment and (b) Torque.

in the power generated which illustrates that by increasing camber further above 7%, the power generation rate will decrease below 65 W. Coefficient of power is a constant value that depicts conversion efficiency of wind energy to electricity. Figure 9(b) demonstrates the mean CP for different camber values, it can be determined that there is sudden jump in CP from camber C0 to C1 and it kept on increasing till C3 and after that a decline in performance value was observed. This provides result that at camber value 3% can provide maximum electricity generation as efficiency of conversion is at its peak.

It can be seen from Figure 7 and Figure 8 trends that cambered aerofoil C3 behavior is efficient. C5 model generate power near to C3 but difference between maximum and minimum periodic points is greater in C5 case resulting in more variation than C3. Variation from C0 to C7 is due to geometry. As camber increases aerofoil geometry

Figure 9. Camber variation effect due to changing; (a) Power and (b) Cp.

become more aerodynamic for moving in circular region but further increase in camber decrease performance as geometry become less aerodynamic. This probability occurs due to more drag and early flow separation from aerofoil. Thus, making flow more turbulent and causing interfering towards other aerofoil of VAWT model.

Bar graph in Figure 10 show VAWT performance variation by changing the camber values. It can be depicted that from C0 to C1 there is a sudden jump in performance by 92%, whereas, performance by changing camber from C1 to C3 is 41.6%. However, from C3 to C5 percentage begin to decrease resulting in demotion of performance.

Figure 10. Percentage increase and decrease in CP by changing camber value.

Table 1 depicts comparison of torque data of individual blade and combined, in case of literature [2] and this research. By initiating a formal planning process, this study involves reference geometry [2] for data comparison. Both have same geometric parameters but this study include more technically refined mesh and major contribution involves multiple sample for extending investigation to cambered aerofoil. Last column on the extreme right in table provides information about percentage error in case of blades. The cause of error is due to difference in number of nodes and elements between reference paper [2] and this study. Reference paper has 84,839 total nodes and 83,641 elements, however, our study has 108,384 no of nodes and 141,154 elements.

5. Conclusion

The camber increment has shown a relative increase in the performance of the VAWT, however, as the camber changes from 3% to 5% a downward trend is observed. This sudden drop in the performance shows that there is a certain limit for elevating the efficacy of the VAWT using camber. The initial increment of the performance from camber 0% to 3% can be due to the changing pressure differential profile around the aerofoil geometry. Furthermore, the drop after 3% camber is due to increased turbulence and flow separation. Future work may involve 3-D analysis of complex model of VAWT.

Cite this paper

Loya, A., Khan, M.Z.U., Bhutta, R.A. and Saeed, M. (2016) Dependency of Torque on Aerofoilcamber Variation in Vertical Axis Wind Turbine. World Journal of Mechanics, 6, 472-486. http://dx.doi.org/10.4236/wjm.2016.611033

References

1. 1. Singh, M., Biswas, A. and Misra, R. (2015) Investigation of Self-Starting and High Rotor Solidity on the Performance of a Three S1210 Blade H-Type Darrieus Rotor. Renewable Energy, 76, 381-387.
https://doi.org/10.1016/j.renene.2014.11.027

2. 2. Biadgo, M.A., et al. (2013) Numerical and Analytical Investigation of Vertical Axis Wind Turbine. FME Transactions, 41, 49-58.

3. 3. Klimas, P.C. and Worstell, M.H. (1981) Effects of Blade Preset Pitch/Offset on Curved-Blade Darrieus Vertical Axis Wind Turbine Performance. Sandia National Laboratories, Albuquerque.

4. 4. Paraschivoiu, I., Trifu, O. and Saeed, F. (2009) H-Darrieus Wind Turbine with Blade Pitch Control. International Journal of Rotating Machinery, 2009.

5. 5. Baker, J. (1983) Features to Aid or Enable Self Starting of Fixed Pitch Low Solidity Vertical Axis Wind Turbines. in Wind Engineering 1983 3C. Proceedings of the Sixth International Conference on Wind Engineering, Gold Coast, 21-25 March 1983, and Auckland, 6-7 April 1983.

6. 6. Beri, H. and Yao, Y. (2011) Effect of Camber Airfoil on Self Starting of Vertical Axis Wind Turbine. Environmental Science & Technology, 4, 302-312.
https://doi.org/10.3923/jest.2011.302.312

7. 7. Li, Q.A., et al. (2015) Effect of Number of Blades on Aerodynamic Forces on a Straight-Bladed Vertical Axis Wind Turbine. Energy, 90, 784-795.
https://doi.org/10.1016/j.energy.2015.07.115

8. 8. Nobile, R., et al. (2013) Unsteady Flow Simulation of a Vertical Axis Wind Turbine: A Two-Dimensional Study. EngD Conference, 2 July 2013.

9. 9. Nobile, R., et al. (2014) Unsteady Flow Simulation of a Vertical Axis Augmented Wind Turbine: A Two-Dimensional Study. Journal of Wind Engineering and Industrial Aerodynamics, 125, 168-179.
https://doi.org/10.1016/j.jweia.2013.12.005

10. 10. Lanzafame, R., Mauro, S. and Messina, M. (2014) 2D CFD Modeling of H-Darrieus Wind Turbines Using a Transition Turbulence Model. Energy Procedia, 45, 131-140.
https://doi.org/10.1016/j.egypro.2014.01.015

11. 11. Kaminsky, C., et al. (2012) A CFD Study of Wind Turbine Aerodynamics. A CFD Study of Wind Turbine Aerodynamics, Ohio, American Society for Engineering Education.

12. 12. Edwards, J., et al. (2007) Wind Tunnel and Numerical Study of a Small Vertical Axis Wind Turbine. 45th AIAA Aerospace Sciences Meeting and Exhibit, Reno.

13. 13. Roh, S.-C. and Kang, S.-H. (2013) Effects of a Blade Profile, the Reynolds Number, and the Solidity on the Performance of a Straight Bladed Vertical Axis Wind Turbine. Journal of Mechanical Science and Technology, 27, 3299-3307.
https://doi.org/10.1007/s12206-013-0852-x

14. 14. Abu-El-Yazied, T.G., et al. (2015) Effect of Number of Blades and Blade Chord Length on the Performance of Darrieus Wind Turbine. American Journal of Mechanical Engineering and Automation, 2, 16.

15. 15. Bhutta, M.M.A., et al. (2012) Vertical Axis Wind Turbine—A Review of Various Configurations and Design Techniques. Renewable and Sustainable Energy Reviews, 16, 1926-1939.
https://doi.org/10.1016/j.rser.2011.12.004

16. 16. Claessens, M. (2006) The Design and Testing of Airfoils for Application in Small Vertical Axis Wind Turbines. Master of Science Thesis.

17. 17. Islam, M., Ting, D.S.-K. and Fartaj, A. (2008) Aerodynamic Models for Darrieus-Type Straight-Bladed Vertical Axis Wind Turbines. Renewable and Sustainable Energy Reviews, 12, 1087-1109.
https://doi.org/10.1016/j.rser.2006.10.023

18. 18. Wu, J.-Y. and Lee, R. (1997) The Advantages of Triangular and Tetrahedral Edge Elements for Electromagnetic Modeling with the Finite-Element Method. IEEE Transactions on Antennas and Propagation, 45, 1431-1437.
https://doi.org/10.1109/8.623133

19. 19. Currie, I.G. (2002) Fundamental Mechanics of Fluids. 3rd Edition, Taylor & Francis, Milton Park.

20. 20. Li, L., Meng, H. and Wang, Y.M. (2013) The Applicability of Two-Equation Turbulence Models in Wind Velocity Prediction. Applied Mechanics and Materials, 291-294, 518-521.
https://doi.org/10.4028/www.scientific.net/AMM.291-294.518

21. 21. Cao, H. (2011) Aerodynamics Analysis of Small Horizontal Axis Wind Turbine Blades by Using 2D and 3D CFD Modelling.

22. 22. Micha Premkumar, T., Seralathan, S., Mohan, T. and Saran Reddy, N.N.P. (2015) Numerical Studies on the Effect of Cambered Airfoil Blades on Self-Starting of Vertical Axis Wind Turbine Part 1: NACA 0012 and NACA 4415. Applied Mechanics and Materials, 787, 250-254.
https://doi.org/10.4028/www.scientific.net/AMM.787.250

Nomenclature

Submit or recommend next manuscript to SCIRP and we will provide best service for you:

A wide selection of journals (inclusive of 9 subjects, more than 200 journals)

Providing 24-hour high-quality service

User-friendly online submission system

Fair and swift peer-review system

Efficient typesetting and proofreading procedure

Display of the result of downloads and visits, as well as the number of cited articles

Maximum dissemination of your research work

Submit your manuscript at: http://papersubmission.scirp.org/

Or contact wjm@scirp.org