Sage Journals: Discover world-class research

Abstract

A numerical simulation study on the combination of a Darrieus and a Savonius wind turbine is conducted. Hybrid T-II, T-III, and T-IV turbines are suggested with the same Darrieus turbine T-I. In the T-II and T-III, the Savonius turbine is at the center of the Darrieus turbine, whereas in the T-IV, the Savonius turbine is above the Darrieus turbine. The T-III Savonius turbine has half the radius of that of the T-II turbine. Results reveal that variations of the power coefficients, C_p with the tip speed ratio, TSR for the hybrid turbines have different slopes. It is observed that C_p increases with increasing TSR for the T-II and T-IV and does not decrease for the range of TSRs considered, in contrast with the C_p behavior of the T-I. The proposed hybrid T-IV turbine has also a larger C_p than the T-I turbine at the highest TSR.

Keywords

Hybrid Darrieus-Savonius turbines vertical-axis wind turbine numerical simulation

Introduction

Wind energy plays an important role in global energy production from renewable energy sources. The overall world’s wind energy potential is estimated to be about 2700 trillion watts, which shows the potential impact of wind energy. Other important advantages of wind energy are its low environmental impact, creating local wealth and employment, and its lower price, compared to many other energy sources for electrical energy production, even some fossil fuels. Considering the above-mentioned points, wind energy can be claimed as one of the well-established clean energies, which has the lowest start-up cost, even in comparison with other clean energy sources.

Wind turbines convert wind kinetic energy to mechanical power and subsequently to electrical energy (Graham, 1999). Wind turbines are generally categorized into horizontal-axis wind turbines (HAWTs) and vertical-axis wind turbines (VAWTs). Traditionally, HAWTs dominated the wind power industry, but, nowadays the use of VAWTs has become more widespread in urban areas. One foible of HAWTs is that the blades use the aerodynamic lift to spin, hence they should always face toward the wind, which requires the use of expensive mechanical systems. Instead, VAWTs use lift, drag, or a mixture of the two so they can use wind coming from any direction. In addition, due to their vertical orientation, the required mechanical equipment can be placed closer to the ground, which substantially reduces maintenance costs, size, and unpleasant noise. Besides, VAWTs have a better self-starting ability at lower wind speeds. Due to these facts, VAWTs are considered more operational for urban areas with slower and more turbulent wind flows (Stathopoulos et al., 2018; Winslow, 2017). Many studies have been performed to understand the flow dynamics around VAWTs to improve their efficiency. Some recent examples of such studies are summarized in the following.

The effect of airfoil asymmetry on the efficiency of a VAWT was experimentally investigated by Islam (2008) and was found that the efficiency with symmetric airfoil sections was higher than with asymmetric ones. Wang et al. (2010) investigated the aerodynamic performance of a VAWT with straight blades at different wind speeds. Improving the self-starting capability of VAWTs was performed by Beri and Yao (2011), Howell et al. (2010) conducted numerical and experimental research on small VAWTs to examine the blade surface roughness effect on their performance. The effect of the number of blades on the performance of a straight-bladed VAWT was studied by Castelli et al. (2012). It was found that the efficiency of three-bladed VAWTs was higher than that of four or five-bladed ones. Gupta et al. (2008, 2012) examined the performance of two blades of a Savonius turbine with five overlaps of 16.2%, 20%, 25%, 30%, and 35%, among which, the case with 16.2% overlap had the highest output power. Stall dynamics at different azimuthal angles of a VAWT were studied by Buchner et al. (2015) through experiments and numerical simulations using unsteady RANS. A reasonable agreement between the two approaches was observed with differences in flow separation location at the highest tip-speed ratio (TSR). Ionescu et al. (2015) used computational fluid dynamics to conduct a study to increase the efficiency of wind turbines by attaching a flap to the end of VAWT blades. It was found that changing flap angles up to a certain degree could improve the performance of the VAWT. To determine optimal pitch angles, Abdalrahman et al. (2017) examined the aerodynamic performance of a straight-bladed Darrieus VAWT at different pitch angles and TSRs, using the finite-volume method. For each studied case, the power coefficient, $C_{P}$ was calculated and compared with experimental results. Siddiqui et al. (2018) examined the performance of a number of hybrid Darrieus-Savonius VAWTs in a laboratory set-up. It was found that at high TSRs, $C_{P}$ was abruptly reduced, which was attributed to the sudden increase of viscous forces at high rotation rates, compared to the production capacity. Hosseini and Goudarzi (2019) evaluated the performance of a hybrid VAWT with a two-bladed Savonius and a three-bladed Darrieus turbine to study the operating range and self-starting capabilities. The results showed that, while the Darrieus turbine had the highest $C_{P}$ of about 48% at TSR $\approx 2.5$ , it did not have the ability to self-start. Instead, the hybrid turbine achieved a maximum $C_{P}$ of about 41% at the operating TSR range of 2.5–4.5 and also demonstrated self-starting capabilities. Using computational fluid dynamics, Nuñez-Quispe et al. (2020) examined the performance of a hybrid Darrieus-Savonius turbine, considering wind characteristics in a rural area, and predicted a maximum electrical energy production of about 252 W in the wintertime. Ibadi et al. (2021) experimentally investigated the effect of added pitch angle on the performance of a hybrid Darrieus-Savonius turbine. The variables of this study were the wind speed and the pitch angle. A pitch angle in the range of 5°–30° was added with five variations. The results showed that increasing the pitch angle reduces the cut-in speed. The highest cutting speed of about 4.2 m/s, occurred at the pitch angle 5° and the lowest cut-in speed of about 3.8 m/s occurred at the pitch angle of 30°.

It is important to note that, stability is a crucial factor in the design, operation, and overall performance of wind turbines. The stability in wind turbines refers to the ability to maintain a steady and controlled operation under various environmental conditions. Structural stability, aeroelastic, and aerodynamic stabilities among other different aspects of stability are important to the turbine operation. There has been some interesting research on the stability of vertical axis wind turbines, such as Owens and Griffith (2014) and Wang et al. (2023).

There exist different extensive review works on the history of VAWTs, see for example the work by Möllerström et al. (2019). The most commonly used VAWTs, in general, are divided into Savonius and Darrieus types. The differences between the two geometries have caused differences in efficiency and ability to self-start. Savonius blades with concave and convex surfaces work with drag force due to the wind. Due to the proximity of the blades to the axis of rotation, Savonius turbines do not require high torque to start rotating and could even operate at low wind speeds. However, their efficiency is lower than the Darrieus turbines. Unlike the Savonius, Darrieus turbines operate with lift force and can have an efficiency of about 0.5. Despite this advantage, due to the blade location, away from the axis of rotation, these turbines require high torque to rotate and operate at higher wind speeds. One way to overcome this disadvantage is to use a control system to measure wind speed, but this system reduces power generation capacity and is costly. So, it is a necessity to design a Darrieus turbine with the ability to self-start (Kumar et al., 2018). Thus, hybrid configurations, employing Savonius turbines with self-start ability and helical Darrieus turbines with higher operating efficiency could potentially lead to the development of a hybrid wind turbine with both advantages.

In the present study, first, a single helical Darrieus wind turbine is numerically investigated. Then, by adding Savonius blades to the Darrieus turbine, three hybrid wind turbine structures are proposed and their performances are thoroughly examined and compared. Significant and novel results are achieved in terms of design and performance analysis of the hybrid VAWTs.

Geometry and mesh

Figure 1 shows the schematic of the reference helical Darrieus VAWT of Cheng et al. (2017), named T-I turbine, and the three hybrid Darrieus-Savonius geometries namely, T-II, T-III, and T-IV, designed in the present study with similar helical Darrieus blades as T-I, but with different layouts for the Savonius blades. The four-bladed T-I turbine, which serves as the validation test case, has a NACA0018 blade profile with a chord length of 0.1 m, rotor radius of 0.21 m, and blade height of 0.54 m. The Savonius geometry of the T-II, T-III, and T-IV turbines is located at the center of the hybrid wind turbine. The general idea of this arrangement is based on the laboratory structures presented by Siddiqui et al. (2018). With this in mind, the Savonius blades are at the center of the Darrieus blades in turbine T-II. The height of the Savonius blade is 0.2175 m and its radius is 0.13 m. In turbine T-III, the Savonius blades are also at the center of the Darrieus blades, but their radius is half of that of the T-II turbine. In the T-IV turbine, the Savonius blades are positioned above the Darrieus blades and the radius of the Savonius turbine is the same as the T-II turbine.

Figure 1.

Schematic of the helical wind turbine T-I (Cheng et al., 2017) and the three hybrid turbines, T-II, T-III, and T-IV, designed in the present study.

In the results and discussion section, performances of the T-I, T-II, T-III, and T-IV turbines are compared and the most efficient geometry is introduced, based on the $C_{p}$ values.

A hybrid tetrahedral and structured boundary-layer-type mesh is employed, see Figure 2(c). To properly resolve the boundary layer, 20 cells, in the form of a structured hexahedral mesh, span the boundary layer in the wall-normal direction. The first wall-adjacent cell has an average wall distance of less than $y^{+} \approx 1.8$ in wall units, which is suitable for the use of the transition-SST model 16. A representation of the surface mesh is shown in Figure 2(a) and the boundary layer mesh can be observed in Figure 2(b) for the Darrieus blades. In addition, Figure 3 shows the surface, volumetric, and boundary layer mesh for the Savonius blades of the T-II turbine. To perform the grid independence analysis, three grids denoted as cases G1, G2, and G3, were employed, see Table 1. A comparison of the static pressure predictions along the 50% blade span at 60° blade rotation angle is given in Figure 4, for the T-I turbine for the three grids. Only slight differences in the pressure predictions obtained with the G2 and G3 grids are observed. Hence, the G2 grid is considered the optimal mesh of the Darrieus turbine in the proceeding simulations. The number of elements in the hybrid geometries T-II to T-IV is not changed dramatically from that of the T-I turbine, due to the same number of meshes within the rotating region, and about 7,670,000 elements are created for these geometries. Figure 5 shows the number of elements with different quality levels for case G2. As can be seen, the number of elements below 0.5 is close to 0. Hence, the mesh has a proper quality, and the numerical predictions are expected to be reliable.

Figure 2.

Representation of the surface mesh on the Darrieus blade (a), boundary layer mesh at a cross section of the Darrieus blade (b) and a horizontal section of the volumetric mesh at the mid hight of the turbine throughout the domain (c) for the T-II turbine.

Figure 3.

Sections of the volumetric mesh used for the T-II turbine. Boundary layer mesh (a), surface mesh, (b) and volumetric mesh at the mid-section (c).

Table 1.

Number of mesh elements used in the simulations.

No.	Case	Number of elements
1	G1	5,850,000
2	G2	7,600,000
3	G3	9,880,000

Figure 4.

Comparison of the static pressure distribution along 50% span and 60° blade angle for grid : G1, : G2, and : G3.

Figure 5.

Grid quality diagram for the G2 grid.

Governing equations and numerical method

The governing equations are the unsteady Reynolds-averaged continuity and Navier–Stokes equations, which in the incompressible form, using the summation convention, read

\frac{\partial U_{i}}{\partial x_{i}} = 0,

(1)

\frac{\partial U_{i}}{\partial t} + U_{j} \frac{\partial U_{i}}{\partial x_{j}} = - \frac{1}{ρ} \frac{\partial P}{\partial x_{i}} + ν \frac{\partial^{2} U_{i}}{\partial x_{j} \partial x_{j}} - \frac{\partial}{\partial x_{j}} < u'_{i} u'_{j} >

(2)

where, $ρ$ is the fluid density, $ν$ is the kinematic viscosity, $U_{i}$ is the mean velocity, $p$ is the mean pressure, and $u'_{i}$ is the velocity fluctuation. The last term on the right-hand side, $< u'_{i} u'_{j} >$ is the Reynolds stress tensor, which needs to be modeled to close the equations. The Reynolds stress tensor is modeled via the Boussinesq hypothesis:

< u'_{i} u'_{j} > - \frac{1}{3} < u'_{k} u'_{k} > δ_{ij} = - 2 ν_{T} S_{ij} S_{ij} = \frac{1}{2} (\frac{\partial U_{i}}{\partial U_{j}} + \frac{\partial U_{j}}{\partial U_{i}}),

(3)

where, $S_{ij}$ is the mean strain-rate tensor, $δ_{ij}$ is the Kronecker delta, and $ν_{T}$ is the eddy viscosity, which has to be modeled using a suitable RANS turbulence model.

Choosing an appropriate turbulence model is crucial for the successful numerical simulation of a VAWT, which can significantly affect the accuracy of numerical predictions, such as the pressure coefficient. The transition-SST model with standard model coefficients is employed in the current simulations. The transition-SST model uses the SST $k - ω$ transport equations (Menter, 1994), along with a transport equation for the prediction of intermittency and transition onset, using the momentum thickness Reynolds number (Menter et al., 2006). Previous investigations, show that this model is suitable for the numerical simulation of VAWTs (Almohammadi et al., 2015; Lanzafame et al., 2013; Yao et al., 2012). Governing equations were discretized using the finite-volume method. A second-order upwind scheme was employed for the convective terms in the governing equations of momentum and turbulence models, whereas diffusion terms were discretized using a second-order central scheme. Time discretization was also carried out using a second-order implicit scheme and the SIMPLEC algorithm was employed for the pressure and velocity coupling. A time step of $Δ t = 1 e^{- 4} (s)$ is used and the Courant number is limited to 1 in the numerical simulations. ANSYS FLUENT 18.2 was used to solve the discretized governing equations.

Domain layout, size, and boundary conditions are similar to those of the study by Cheng et al. (2017), see Figure 6. The domain extent and position of the wind turbine are chosen such that the inlet and outlet boundary conditions do not affect the numerical predictions. The domain is divided into a stationary and a rotating part, which encloses the wind turbine. The sliding mesh method is used to connect the rotating and non-rotating parts of the domain, as suggested by Roychowdhury (2020). The matching method of the FLUENT software was used for the interface between the rotating and stationary parts of the domain. A constant 9 m/s velocity inlet boundary condition with a constant turbulence intensity of 1.5% is used at the inlet and a constant 1 atm pressure was employed at the outlet, see Figure 6.

Figure 6.

A horizontal cross-section of the schematic view of the domain, showing the layout, size, boundary conditions, and relative angle of the blades $θ$ . Turbine radius is denoted by $R$ .

Results and discussion

Tip-speed ratio (TSR) is used as the parameter in the analysis of the wind turbines. It is defined as the ratio of the blade tip speed to the free-stream velocity, written as (Ragheb and Ragheb, 2011)

T S R = \frac{R ω}{V_{\infty}},

(4)

where, $ω (r a d / s e c)$ is the angular velocity, $V_{\infty} (m / s e c)$ is the wind velocity, and $R (m)$ is the rotor radius. The six rotation rates used in this study and the corresponding TSRs are given in Table 2.

Table 2.

Numerical predictions of the $C_{P}$ for the T-I, T-II, T-III, and T-IV turbines of the present study at different TSRs.

$ω$	TSR	T-I	T-II	T-III	T-IV
38.57	0.9	0.0315	0.0286	0.0295	0.0274
48.86	1.1	0.0906	0.0717	0.0659	0.0683
53.57	1.2	0.0971	0.0774	0.0722	0.0798
98.57	1.4	0.1100	0.0890	0.0842	0.0997
77.14	1.8	0.1415	0.1004	0.0831	0.1343
98.57	2.3	0.0849	0.1155	0.0814	0.1625

Validation of the current numerical method

A comparison of the static pressure distribution over the T-I turbine blade, at 50% blade span and rotation angles $θ =$ 30°, 60°, 90°, and 120°, between the current predictions and those of Cheng et al. (2017) is given in Figure 7, for TSR = 1.46. A reasonable agreement between the current and reference predictions is observed with 2%–5% relative differences. Thus, this numerical solution process can be confidently used to simulate the similar geometries of the present study, namely wind turbines T-II, T-III, and T-IV.

Figure 7.

Static pressure distribution over the 50% blade span of the T-I turbine at angles $θ = 30 °$ (a) 60° (b) 90° (c), and 120° (d) and TSR = 1.46 for the present simulation: and numerical results of Cheng et al. (2017): .

Figure 8 shows the static pressure distribution over the mid-section of the blade for 90° rotations of the T-I turbine at TSR = 1.2. The figure demonstrates that the dynamic stall phenomenon, which has a complex physics, is captured in the current numerical simulations. It should be noted that the dynamic stall process is a phenomenon in wind turbines, during which the lift increases intensively, and performance improves for a short period of time. This continues with a rapid drop of the lift force and as a result, a drop in the overall turbine performance. Then, the boundary layer experiences a recovery process. The dynamic stall limits the operation of the VAWT and is known to be more intensive at lower TSRs. For brevity, pressure contours will not be examined for the other three cases.

Figure 8.

Static pressure distribution over the mid-section of the T-I turbine at TSR = 1.2 for 90° blade rotation, showing the dynamic stall phenomenon, different panels show the blade at angles $θ = 90 °$ (a), 108° (b), 126° (c), 144° (d), 162° (e), and 180° (f).

Static pressure predictions

To visualize the Savonius and Darrieus turbine flow interactions, instantaneous pressure contours at TSR = 0.9 and 1.2 are given at the mid-section of the T-I, T-II, T-III, and T-IV turbines in Figure 9. Pressure distribution over the T-I turbine shows the regular pressure distribution over a Darrieus turbine, see Figure 9(a) and (b). The figure shows small pressure differences between the lower and higher TSRs.

Figure 9.

Static pressure distribution at the mid-section of the Darreius turbine of the: (a) T-I @ $θ = 45^{o}$ , (b) T-I @ $θ = 15^{o}$ , (c) T-II @ $θ = 45^{o}$ , (d) T-II @ $θ = 35^{o}$ , (e) T-III @ $θ = 45^{o}$ , and (f) T-III @ $θ = 30^{o}$ turbines and at the mid-section of the Savonius turbine of (g) T-IV turbine @ $θ = 75^{o}$ and (h) T-IV turbine @ $θ = 90^{o}$ . In all figures: TSR = 1.2 (left) and TSR = 0.9 (right).

When the Savonius turbine is placed inside the Darrieus turbine (T-II and T-III turbines) pressure interactions are observed between the two turbines. The interaction is more intense in the T-II turbine, which has a larger radius, compared with the T-III turbine with the smaller Savonius turbine, compare Figure 9(c) and (d) with Figure 9(e) and (f).

Figure 9(g) and (h) show the pressure distribution at the mid-section of the Savonius blades of the T-IV turbine, where the Savonius turbine is placed on top of the Darrieus turbine. As is expected, an interaction between the two turbines is not observed.

Figure 10 shows a 3D view together with the top, bottom, and mid-section perspectives of the static pressure distribution for the hybrid Darrieus-Savonius wind turbine T-III. The 3D view of the pressure field is visualized to offer a more comprehensive view of the pressure field, throughout the entire turbine surface. In Figure 10(a), it is evident that the pressure level is comparatively lower in the inner Savonius turbine than in the outer helical Darrieus turbine. An intriguing observation reveals the presence of highly intense high- and low-pressure regions at the tip of the helical blades. Furthermore, as previously mentioned, notable pressure interactions can be observed between the Savonius and Darrieus turbines. Specifically, on the inner side of the Savonius turbine, discernible low-pressure regions are apparent, as depicted in Figure 10(c) and (d), which correspond to the mid-section view and the bottom-section perspective, respectively. Additionally, the pressure distribution at the bottom section of the hybrid turbine exhibits more pronounced pressure differences, as shown in Figure 10(d).

Figure 10.

Static pressure distribution of the hybrid Darrieus-Savonius wind turbine: (a) the 3D view, (b) the top view, (c) the mid-section view, and (d) the bottom view of the pressure field.

Blade surface static pressure is also given at 30%, 50%, and 70% blade span in Figure 11 for the helical Darrieus blades of the T-IV turbine at rotation angle $θ = 30 °$ and TSR = 1.14. Prediction at the 50% span is similar to that presented earlier in Figure 7(a) for the T-I turbine but at a higher TSR. At the higher blade span, that is, 70% blade span, the difference between the suction and pressure sides of the blade surface is reduced, hence, the net aerodynamic force is decreased.

Figure 11.

Static pressure distribution over blade sections at 30% span , 50% span , and 70% span of the helical blade at $θ = 30 °$ with TSR = 1.14 for the T-IV turbine.

Effect of the Savonius turbine position on the power production

To assess the power production by the hybrid turbines, the power coefficient, $C_{P}$ defined as is used. Here, $Γ_{t o t a l} (N . m)$ is the total torque exerted on the turbine by the aerodynamic forces, $ω (r a d / s)$ is the rotation speed, $ρ (kg / m^{3})$ is the air density, $A (m^{2})$ is the swept area by the turbine, and $V_{\infty} (m / s)$ is the inlet air velocity. The area $A$ is calculated by multiplying the blade height, $H$ by the revolution diameter $D$ . The torque is due to the pressure and viscous forces. The pressure force can also be divided into tangential, $F_{t}$ , and normal, $F_{n}$ forces. The former is parallel and the latter is perpendicular to the blade section chord. The total torque is the sum of the forces acting on the four blades, that is, $Γ_{t o t a l} = \sum_{i = 1}^{4} F_{ti} D / 2$ .

C_{P} = \frac{Extracted energy}{Wind kinetic} = \frac{Γ_{t o t a l} ω}{\frac{1}{2} ρ {A V}_{\infty}^{3}},

(5)

A comparison between the predicted $C_{P}$ for the T-I, T-II, T-III, and T-IV turbines at different TSRs is given in Table 2. A graphical comparison is also provided in Figure 12, between the current predictions and those of the reference experimental and numerical data of Cheng et al. (2017) for the T-I turbine. It is observed that $C_{P}$ for the T-I turbine increases with increasing TSR, reaching a maximum at TSR = 1.8, after which, it decreases by almost 50% at TSR = 2.3. The current $C_{P}$ predictions for this turbine at different TSRs mostly lie between those of the experimental and numerical predictions of the reference turbine with a better agreement with the latter at TSR = 0.9, 1.1, 1.8, and 2.3 and with the former at other TSRs.

Figure 12.

Comparison of the power coefficient, $C_{P}$ predicted for the T-I , T-II , T-III , and T-IV turbines with numerical simulations and experimental measurements of Cheng et al. (2017).

Adding the Savonius turbine in the T-II, T-III, and T-IV configurations leads to a different behavior in the $C_{P}$ predictions, where a decrease in the $C_{P}$ predictions is not observed in the range of TSRs considered, except for the T-III turbine, where a slight decrease is observed at the highest TSR. However, it is expected that $C_{P}$ eventually decreases with a further increase in the TSR in all cases. At TSR lower than 1.8, the hybrid turbines have lower $C_{P}$ , compared with the T-I turbine. However, the T-IV turbine, where the Savonius turbine is placed on the top of the configuration, has a higher $C_{P}$ than the T-I turbine at TSR greater than 1.8. Significant differences in the $C_{P}$ predictions are also observed for the hybrid T-II, T-III, and T-IV turbines starting at TSR = 1.4. One can observe about 20% and 15% increase in the $C_{P}$ at TSR = 2.3 for the T-IV and T-II turbines, respectively, relative to TSR = 2.1 predictions, indicating an increase in the effect of the viscous forces and vortex mixing with increasing TSR. Hence, the performance of the T-IV turbine at TSR higher than 1.2 is better than the other hybrid turbine geometries. Moreover, the T-II turbine has a better performance than the T-III turbine, with smaller Savonius blades, at all TSRs except at the lowest one.

It should also be noted that the hybrid T-II turbine has a slightly higher $C_{P}$ at TSR less than 1.2, compared with other hybrid turbines, suggesting that adding the Savonius geometry in the configuration could be suitable for the self-starting of the wind turbine at low TSRs. On the contrary, the presence of the Savonius blades inside the Darreius turbine at high TSR increases the viscous force, due to the interference of vortices resulting from the simultaneous rotation of the Darrieus-Savonius blades, thereby reducing the wind turbine power. Hence, in the T-IV turbine, where the Savonius geometry is located above the Darrieus turbine, the interference effect of the vortices is minimized.

Another important point in the experimental studies of hybrid Darrieus-Savonius wind turbines is related to the vibration effects of wind turbine blades, see the experimental work of Siddiqui et al. (2018). These vibrations could have a larger effect on the performance of the hybrid T-IV turbine, compared with the other geometries. Because in this geometry, the Savonius blade is positioned above the helical blade, and this structure at high rotation speeds could lead to a significant increase in blade vibrations.

Conclusions

The present study examined three different combinations of Savonius and helical Darrieus wind turbines to investigate the effect of the Savonius turbine position and radius on the overall aerodynamic performance of the hybrid wind turbine. The most optimal hybrid vertical wind turbine configuration was introduced.

According to the results, choosing the best case depends on the goal of adding the Savonius turbine. If it is merely due to performance, the addition of the Savonius turbine is only justified if the turbine is operating at tip speed ratios higher than 1.8. The best performance in this interval was observed when the Savonius turbine was positioned on top of the helical Darrieus turbine. The main reason for this difference in the $C_{P}$ predictions was associated with the interaction of the vortices due to the Savonius turbine located inside the helical Darrieus turbine, which would negatively affect the overall performance. At TSRs below 1.8, the single helical Darrieus wind turbine performed better than the hybrid ones, regardless of the radius or the position of the Savonius turbine. Thus, in this interval, the hybrid turbines do not have any advantages over the Darrieus turbine. However, at TSRs lower than about 0.9, one can speculate from the slopes of the $C_{P}$ curves that the hybrid turbines should have higher $C_{P}$ values than the Darrieus turbine. The smallest slope of the $C_{P}$ curves was observed when the Savonius turbine with the smallest radius was placed inside the helical turbine. Thus, in order to reach the optimum radius of the Savonius turbine, further investigations and numerical simulations are required.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Zeinab Pouransari

References

Abdalrahman

Melek

Lien

(2017) Pitch angle control for a small-scale Darrieus vertical axis wind turbine with straight blades (h-type VAWT). Renewable Energy 114: 1353–1362.

Almohammadi

Ingham

, et al. (2015) Modeling dynamic stall of a straight blade vertical axis wind turbine. Journal of Fluids and Structures 57: 144–158.

Beri

Yao

(2011) Effect of camber airfoil on self starting of vertical axis wind turbine. Environmental Science & Technology 4: 302–312.

Buchner

Lohry

Martinelli

, et al. (2015) Dynamic stall in vertical axis wind turbines: Comparing experiments and computations. Journal of Wind Engineering and Industrial Aerodynamics 146: 163–171.

Castelli

De Betta

Benini

(2012) Effect of blade number on a straight-bladed vertical-axis Darreius wind turbine. World Academy of Science, Engineering and Technology 61: 305–3011.

Cheng

Liu

, et al. (2017) Aerodynamic analysis of a helical vertical axis wind turbine. Energies 10: 575.

Graham

(1999) Wind Power. Chicago, IL: Heinemann/Raintree.

Gupta

Biswas

Sharma

(2008) Comparative study of a three-bucket Savonius rotor with a combined three-bucket Savonius–three-bladed Darrieus rotor. Renewable Energy 33(9): 1974–1981.

Gupta

Das

Gautam

, et al. (2012) CFD analysis of a two-bucket Savonius rotor for various overlap conditions. ISESCO Journal of Science and Technology 8(13): 67–74.

10.

Hosseini

Goudarzi

(2019) Design and CFD study of a hybrid vertical-axis wind turbine by employing a combined bach-type and h-Darrieus rotor systems. Energy Conversion and Management 189: 49–59.

11.

Howell

Qin

Edwards

, et al. (2010) Wind tunnel and numerical study of a small vertical axis wind turbine. Renewable Energy 35(2): 412–422.

12.

Ibadi

Sangidzun

Wijayanto

, et al. (2021) Effect of adding of pitch on the Darrieus blade against the cut in speed of the Savonius type s - Darrieus type h hybrid turbine. Journal of Physics Conference Series1808: 012003.

13.

Ionescu

Szava

Vlase

, et al. (2015) Innovative solution of vertical axis wind turbine, suitable for naval industry implementation (numerical methods and analytical calculus). Procedia Technology 19: 715–721.

14.

Islam

(2008) Analysis of fixed-pitch straight-bladed VAWT with asymmetric air-foils. PhD Thesis, University of Windsor, Canada.

15.

Kumar

Raahemifar

Fung

(2018) A critical review of vertical axis wind turbines for urban applications. Renewable and Sustainable Energy Reviews 89: 281–291.

16.

Lanzafame

Mauro

Messina

(2013) Wind turbine CFD modeling using a correlation-based transitional model. Renewable Energy 52: 31–39.

17.

Menter

(1994) Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal 32(8): 1598–1605.

18.

Menter

Langtry

Likki

, et al. (2006) A correlation-based transition model using local variables—part I: Model formulation. Journal of Turbomachinery 128: 413.

19.

Möllerström

Gipe

Beurskens

, et al. (2019) A historical review of vertical axis wind turbines rated 100 kw and above. Renewable and Sustainable Energy Reviews 105: 1–13.

20.

Nuñez-Quispe

Cedamanos-Cuenca

(2020) Performance analysis of a hybrid vertical axis wind turbine through CFD for distributed generation in rural areas of Peru. Journal of Physics Conference Series 1618: 042013.

21.

Owens

Griffith

(2014) Aeroelastic stability investigations for large-scale vertical axis wind turbines. Journal of Physics Conference Series 524: 012092.

22.

Ragheb

(2011) Wind turbines theory-the betz equation and optimal rotor tip speed ratio. In: Carriveau

(ed.) Fundamental and Advanced Topics in Wind Power, vol. 1, Croatia: IntechOpen, pp.19–38.

23.

Roychowdhury

(2020) Computational Fluid Dynamics for Incompressible Flows. Boca Raton: CRC Press.

24.

Siddiqui

Mian

Alam

, et al. (2018) Experimental study to assess the performance of combined savonius darrieus vertical axis wind turbine at different arrangements. In: 2018 IEEE 21st international multi-topic conference (INMIC), Karachi, Pakistan, 01–02 November, pp.1–8. New York, NY: IEEE.

25.

Stathopoulos

Alrawashdeh

Al-Quraan

, et al. (2018) Urban wind energy: Some views on potential and challenges. Journal of Wind Engineering and Industrial Aerodynamics 179: 146–157.

26.

Wang

Gao

, et al. (2023) Stability analysis of wind turbine blades based on different structural models. Journal of Marine Science and Engineering 11(6): 1106.

27.

Wang

Ingham

, et al. (2010) Numerical investigations on dynamic stall of low Reynolds number flow around oscillating airfoils. Computers & Fluids 39(9): 1529–1541.

28.

Winslow

(2017) Urban wind generation: Comparing horizontal and vertical axis wind turbines at Clark university in Worcester, Massachusetts. IDCE.

29.

Yao

Yuan

Wang

, et al. (2012) Numerical simulation of aerodynamic performance for two dimensional wind turbine airfoils. Procedia Engineering 31: 80–86.

Numerical investigation of the aerodynamic performance of a hybrid Darrieus-Savonius wind turbine

Abstract

Keywords

Introduction

Geometry and mesh

Governing equations and numerical method

Results and discussion

Validation of the current numerical method

Static pressure predictions

Effect of the Savonius turbine position on the power production

Conclusions

Footnotes

Declaration of conflicting interests

Funding

ORCID iD

References