CFD-Based Performance Analysis and Experimental Investigation of Design Factors of Vertical Axis Wind Turbines under Low Wind Speed Conditions in Thailand

This paper presents effects of design factors on mechanical performance of Vertical Axis Wind Turbines (VAWTs), and an experimental investigation of optimal VAWT performance under low wind speed conditions in Thailand. Design factors include types of wind turbines, number of blades, types of materials, height-to-radius ratios, and design modifications. Potential VAWT models with different design factors are numerically analyzed within a virtual wind tunnel at various wind speeds by utilizing Xflow Computational Fluid Dynamics (CFD) software. The performance curves of each VAWT are obtained as plots of power coefficients against tip speed ratios. It is found that the type of wind turbine, number of blades, and height-to-radius ratio have significant effects on mechanical performance whereas types of materials result in shifts of operating speeds of VAWTs. Accordingly, an optimal VAWT prototype is developed to operate under actual low speed wind conditions. The performance curve from experimental results agrees with the CFD results. The proposed methodology can be used in the computer design of VAWTs to improve mechanical performance before physical fabrication.


Introduction
Wind has been an important energy source due to its renewability and sustainability.There are two types of wind turbines: Vertical Axis Wind Turbines For Thailand, wind speeds are statistically recorded up to 6 m/s at a height of 40 m [1].VAWTs with wind boosters are highly suitable under the low wind speed conditions of Thailand [2].For standalone VAWTs, there were some sustainable developments of Savonius rotors for wind turbines to improve power generation in the past [3] [4] [5] [6].Recently, Schubel and Crossley reviewed aerodynamic designs of turbine blades for potentials of mechanical improvement [7].Mahmound et al. studied the improvement of Savonius VAWT performance [8].Different geometries of Savonius VAWT are investigated in order to determine the most effective conditions.D'Alessandro and et al. derived mathematical models for determining the mechanical performance of wind turbines [9].It was reported that the mathematical models led to efficient rotors for VAWTs.Salyers performed experimental investigations of aerodynamic improvement for VAWTs [10].
Bhutta et al. reviewed previous works, which are related to various configurations and design techniques for VAWTs [11].In wind engineering research, there are many attempts of improving the mechanical performance of Savonius rotors in various approaches.In an experimental investigation [12], the number of blades is increased for multiple blades in different operating conditions under ocean flows.With design modification of blades, a helical Savonius VAWT, twisted 90˚, was proposed to eliminate the excessive shaking of blades [13].The twisted angle of a Savonius VAWT blade is numerically optimized for maximum power generation [14] [15].The modified configuration of blades with an overlap at the center was introduced to improve mechanical torque [16] [17].Also, the number of stages has influences on the aerodynamic behavior of turbulent airflow around Savonius VAWTs where the power coefficient increases as the number of stages increases [18].Additionally, a cluster of two and three VAWTs was investigated to enhance mechanical power generation, compared to a standalone VAWT [19].The optimum gap distance is experimentally determined to avoid the wake of upstream VAWTs.
Those studies yield design knowledge for optimizing VAWTs with high efficiency.However, there is a deficiency of quantitative relationships between design factors and mechanical performance of VAWTs.Design factors of VAWTs are types (Savonius and Darrieus), number of turbine blades, types of materials, height to radius ratios, and design modifications.
In this study, VAWTs with various design factors are numerically simulated within virtual wind tunnels at different wind speeds by utilizing Xflow TM Computational Fluid Dynamics (CFD) software.Correspondingly, the power coefficients are obtained with respect to tip speed ratios.An optimal design of VAWTs is chosen according to those analytical experiences.For experimental investigation, the VAWT is fabricated accordingly, and it is operated under actual

Methodology
In this study, each VAWT model is developed for a given geometry by using computer-aided design software, such as AutoCAD TM and SolidWorks TM .Except for types of materials, design factors, such as types of wind turbines, number of turbine blades, height to radius ratios, and design modification, can be varied in this step.The types of materials are adjusted by setting the values of density in CFD simulation.Under a virtual wind tunnel, the VAWTs are simulated to investigate effects of the design factors on mechanical performance under low wind speed conditions from 1 m/s to 6 m/s.

Xflow TM CFD-Based Setup
For CFD-based case studies, Figure 1 illustrates a geometrical setup of rigidbody dynamics with fixed translation and y-axis rotational constraints.All types of VAWTs in this work have a swept area of 0.125 m 2 .The arrow heads indicate the flow direction of air within a virtual wind tunnel at a given wind speed.For the boundary conditions, a uniform flow of air arrives at the entrance through the virtual wind tunnel with a constant speed.A VAWT, which is installed within the virtual wind tunnel, has no rotation initially.Discretization parameters, such as number, grid resolution, and airflow conditions, are set up in the CFD software and listed in Table 1.The types of materials are defined from different densities of 1000 kg/m 3 , 2440 kg/m 3 , and 3500 kg/m 3 in this work.

Dimensionless Analysis of Mechanical Performance
In CFD-based analysis, the VAWT are simulated at various wind speeds with constant external torsional loads.The angular speeds of the VAWT are recorded under steady state conditions where the external torsional load is equal to the torque of the VAWT shaft.The mechanical power can be determined from multiplication of angular speed and torque.The performance curve is defined between the power coefficient and the tip speed ratio.The tip speed ratio can be determined by: where λ is the tip speed ratio, ω is the angular speed of wind turbine (rad/s), R is the radius of wind turbine (m) and V is the wind speed (m/s).
The power coefficient is defined as the ratio of mechanical power of wind turbine to wind power as expressed in Equation (2).
where p C is the power coefficient, t P is the mechanical power of turbine (W), and w P is the wind power (W).

Experimental Setup
In this section, an optimal VAWT design is chosen, according to analytical studies.The schematic diagram of the VAWT model is a 90˚ twisted two-layer Savonius-type VAWT with covers, as shown in Figure 2. The VAWT model is made of acrylic material, where the shape and dimension are identical to the optimal CFD model.Magnetic levitation is applied at the two supporting joints so as to minimize coulomb friction, as seen in Figure 3.The structure is placed under actual wind conditions.The wind speeds are measured by a hot-wire anemometer, while the angular speeds of the VAWT are detected by an encoder for the determination of power coefficient and tip speed ratio.The mechanical

Results and Discussion
In this section, CFD-based studies of mechanical performances are presented with respect to design factors: material types, number of turbine blades, height to radius ratios, design modification, and turbine patterns.

Effects of Material Types
A Savonius-type VAWT is set up within a virtual wind tunnel, as shown in Figure 1.The number of grids is an important factor for accuracy of numerical results.It is found that the greater the number of grids, the higher the accuracy.
However, the computational time increases exponentially when the number of grids increases linearly.The number of grids in this study is chosen in each case if an increase of the number of grids does not significantly change the simulated results.
The values of density are 1000 kg/m 3 (plastic), 2440 kg/m 3 (glass fiber), and 3500 kg/m 3 (metal) for the three case studies.Figure 4 shows plots of power coefficients against tip speed ratios of all three models.It is observed that a given geometry of VAWTs yields different curves of mechanical performance where the wind turbines are made of different materials.The operating range of angular speed is varied according to the inertia of the VAWTs.The wind turbine with high density rotates more slowly than low density while the same turbine torques, due to wind, are applied to the three VAWTs.In turn, the power coefficient of the wind turbine with low density is higher than high density.

Number of Turbine Blades
Savonius-type VAWTs models, with two, three, and four blades, are used in CFD simulations under various wind conditions.It is seen from Figure 5 that a VAWT with two blades yields the largest power coefficients among the three models.Also, a wide range of tip speed ratios is found from the VAWTs with two blades.The VAWTs can rotate fast with high efficiency.

Height to Radius Ratios
The radius of each model is adjusted with respect to a unity height such a way that the swept area remains constant.As shown in Figure 6, the power coefficient of VAWT increases rapidly as the radius of the VAWT blade increases.
The torque of a VAWT increasingly produces mechanical power.This reaches a maximum value at a ratio of height to radius of around 1:0.375.After that point, it decreases due to a high moment of inertia with slow rotation.

Design Modifications
There are three interesting design modifications: helix pattern, cover-added pattern, and double layer pattern, as shown in Figure 7. From Figure 8, the VAWTs with helix patterns in both 90˚ and 180˚ twisting angles yield a decrease of ability to harvest wind energy while rotating.The VAWTs with cover-added pattern are capable of taking up wind energy better than the conventional pattern.The power coefficients of the VAWTs with a cover-added pattern and double layers slightly decrease due to low height to radius ratio.However, the VAWTs with a cover-added pattern and perpendicular double layers can perform at maximum efficiency of mechanical work, among these design modification.than the wind speed.The value of the power coefficient reaches a maximum.

Darrieus-Type Turbine
However, if the torsional loading on the VAWT is less (or more), the power coefficient declines rapidly.A design for power generation must be considered for a suitable load torque in this case.Furthermore, the power coefficient of Darrieus VAWTs is higher than Savonius VAWTs, while they also run at a higher tip speed ratio.

CFD Analysis and Experiment of Optimal Model
The starting to rotate by itself, especially under low wind speed conditions.Figure 10 shows mechanical performance of an optimal VAWT, shown in For an experiment, the VAWT is installed under actual wind conditions.Angular speeds of the VAWT are measured by an electric encoder, while the wind speed is measured by a hot-wire anemometer.In Figure 11, plots of the angular speed and standard deviation are 5.32 m/s and 1.12 m/s, respectively.They are used to determine the tip speed ratios in Equation ( 1).Accordingly, the power of a VAWT is determined by the rate of change in kinetic energy.The power coefficient is calculated by using Equation (2).In Figure12, the performance curve of CFD results are fitted well to the power coefficients against tip speed ratios, which are obtained from the experiment.The averaged relative difference between experimental data and CFD results are 3.30%.The deviations may be caused by uncertainties in measurement, airflow, and/or VAWT structure.
However, it is confirmed that an optimal VAWT design can be obtained from the CFD analysis before physical fabrication.The power coefficient against tip speed ratio is used as a performance curve in power generation under various wind speed conditions.

Conclusion
Effects of design factors: material types, number of turbine blades, height to radius ratios, design modification, and turbine patterns, are analyzed by using CFD-based simulations within a virtual wind tunnel.The mechanical performances of VAWTs under low speed wind conditions are determined from plots of power coefficient against tip speed ratio.The Savonius-type VAWT can be selected for high torque at low angular speeds while the Darrieus-type VAWT can be chosen for power generation at high angular speeds.A number of blades and height to radius ratio have significant effects on operating conditions of rotation due to inertia.In this study, the Savonius-type VAWT, with a coveradded pattern and perpendicular double layers, yields great mechanical performance.The CFD results of the performance curve have good agreement with the experimental results of the optimal VAWT, which is operated under actual wind How to cite this paper: Unsakul, S., Sranpat, C., Chaisiriroj, P. and Leephakpreeda, T. (2017) CFD-Based Performance Analysis and Experimental Investigation of Design Factors of Vertical Axis Wind Turbines under Low Wind Speed Conditions in Thailand.Journal of Flow Control, Measurement & Visualization, 5, 86-98.https://doi.org/10.4236/jfcmv.2017.54007S. Unsakul et al.DOI: 10.4236/jfcmv.2017.5400787 Journal of Flow Control, Measurement & Visualization (VAWTs) and Horizontal Axis Wind Turbines (HAWTs), which are commonly used to convert the kinetic energy of wind into mechanical energy of wind turbines.In this work, VAWTs have efficiently omnidirectional capability to harvest wind energy at low wind speed conditions.

Figure 1 .
Figure 1.CFD setup: (a) wind turbine within virtual wind tunnel and (b) dimensions of VAWT.

Figure 3 .
Figure 3. Experiment of optimal VAWT under actual wind conditions.

Figure 4 .
Figure 4. Power coefficient curve of wind turbines made of different materials.

Figure 5 .
Figure 5. Power coefficient curves of Savonius-type VAWTs with different number of blades.

Figure 9 (Figure 6 .
Figure 9(a) shows a type of Darrieus VAWT, which is varied with a different height to radius ratio.From Figure 9(b), the mechanical performance of Darrieus VAWTs is better when the height is less than the width.It is observed that the power coefficient is low at low tip speed ratio since the Darriues VAWT has inherently a low torque for rotation.In fact, it requires a driving device for starting rotation.After tuning to high tip speed ratio, the Darrieus VAWT efficiently produces mechanical work while the tip of the VAWT blade runs faster

Figure 8 .
Figure 8. Mechanical performances of different design modifications.

Figure 2 .
The VAWT is simulated under various wind speeds and torsional loads within a virtual wind tunnel of X-Flow TM .The power coefficient is varied according to tip speed ratios.It is observed that the power coefficient increases as the tip speed ratio increases.Once it reaches a maximum value of 3.51% at around a tip speed ratio of 0.55, it decreases to zero, where the VAWT is run under no-load conditions.It should be noted that the power coefficient is lower than the commonly-known values of 0.15 to 0.2 since dimensions and geometry of wind turbines in this work are different from other turbines in literature.For example, influence of the end plates on both sides of the rotor causes different turbulent effects.

Figure 10 .Figure 11 .
Figure 10.Plots of power coefficient curve of VAWT from CFD simulation.

Figure 12 .
Figure 12.Comparison between experimental results and CFD results.

Table 1 .
CFD setup and air properties.