Savonius Wind Turbine Numerical Parametric Analysis Using Space-Filling Design and Gaussian Stochastic Process

Wind energy is an alternative source of clean energy to address the growing energy demand and provide pollution-free electricity. With the rapid development of urban areas, high wind energy resources such as high-rise building rooftops are excellent locations for urban wind turbine installation. One of the practical and simple urban wind turbines is the Savonius design. It has a simple design, easy to maintain, and is very affordable. This work focuses on the design evaluation of a Savonius wind turbine (SWT) by varying the rotor diameter, rotor height, and twist angle for urban applications. A transient computational fluid dynamics (CFD) approach is applied to assess the various design treatments using a space-filling design of experiments. To address the spaces in the hypercube statistical design, a sphere packing design method was adopted which suited the evaluation of computational simulations results such as that of the CFD. The Gaussian stochastic process model was applied to establish the trend of the parametric performance of the optimized SWT design through the model fitting. The results have shown that optimized SWT performs well with its self-starting capability compared to the traditional Savonius design. In addition, the optimized SWT has shown a better peak power coefficient compared with the results of previous works on the design of SWT.


Introduction
The power generation industry has been one of the most important industries in the world especially as the demand of world energy consumption continues to increase. The IEA [1] projected the increase of energy demand by the year 2035 as twice that of the 1990s. Thus, the need to construct different power generation plants increases as well to meet the demand. However, most of these conventional power plants consume non-renewable raw materials such as coal, petroleum, and other fossil fuels. Accordingly, Calautit et al. [2] emphasized that the continuous consumption of fossil fuels would lead to more generation of greenhouse gas emissions. These greenhouse gas emissions play a significant role in worsening the global warming and climate change situation. To address these global concerns, the industry continues to seek renewable power generation methods that can sustainably generate the projected energy demand. Renewable energy sources for power generation have provided an alternative means of power generation for the past few decades. These vary between solar, hydroelectric, wind, and biomass, among others. Among these alternative energy sources, wind energy has gained interest from engineers and scientists to potentially generate energy through wind turbine design for the past Wind 2022, 2 115 further explored the expansion of the computational approach by the development of a fluid-structure interface (FSI) of an SWT. The FSI model was then validated with an experimental result and found to influence the centrifugal forces in the turbine operation. With the use of a computational fluid dynamic (CFD) approach, Marinić-Kragić et al. [25] introduced a scoop-based design of an SWT that enabled a 10% to 39% improvement in its power coefficient.
Various modelling techniques were developed to further enhance the performance of urban wind turbines. An adaptive neuro-fuzzy inference system was developed by Elsisi et al. [26] for the improved blade pitch control of urban wind turbines. A numerical approach coupled with an experimental test was performed by Zalhaf et al. [27] to understand the transient behavior of wind turbine blades that are hit by lightning. Recently, a hybrid numerical approach was developed by Sayed et al. [28] to incorporate dust loaded wire-duct precipitators in the computational domain of urban wind turbines. With the transmission of electricity from the wind turbine, Abouelatta et al. [29] developed a numerical method that integrated a full multigrid method to evaluate the fast corona discharge with varying wind speeds.
These studies are significant in enhancing the performance of the SWT. However, the design enhancement approach using the CFD method coupled with space-filling design with the sphere packing method has not been explored in previous works. CFD is a numerical approach that allows a design of an SWT to be evaluated computationally [30]. On the other hand, space-filling design is a type of statistical design of experiment tool that addresses the non-uniformity of the inside experimental space [31]. The sphere packing method is used in a hypercube which can be well represented by results of computational approaches such as the CFD method [32]. Moreover, a Gaussian stochastic process (GaSP) model was used to fit in a multivariate normal distribution the developed model results from the space-filling design [33]. Based on the available literature, no studies have applied CFD and space-filling design with the GaSP model to design an SWT. Hence, the novelty of this work lies in the methodology to redesign the SWT using a CFD approach with space-filling design of experiment and GaSP model. The designed SWT was evaluated based on the maximum power coefficient and assessed the influences of the considered design parameters.

Performance of the Savonius Wind Turbine
The performance of the SWT is best described by its power coefficients where the moment coefficient and the power coefficient are given as C m and C p , respectively. These coefficients were obtained from the ratios between the actual parameter over the theoretical parameter. Two of the important parameters are the moment and the power coefficients [34]. The moment coefficient is quantified using Equation (1) while the power coefficient is determined using Equation (2). (1) where A provides the swept area of the turbine in m 2 , M represents the rotor moment in N-mm, P w indicates the extracted power in W, V depicts the air velocity in m/s, r shows the radius of the rotor in m, ρ covers the density of air in kg/m 3 , and λ denotes the tip-speed ratio. Barlow et al. [35] described the model frontal area of the turbine concerning the cross-sectional area of the air domain as shown in Equation (3).

116
To establish a deterministic model validation for the simulation results, the Gaussian Stochastic process (GaSP) model was considered. The GaSP model consists of an n × 1 vector for its output response, y(x), which is given as N(µ1 n , σ 2 R(X, θ). The parameter R(X, θ) in the output response is a n × n correlation matrix and can take different forms depending on the situation [36]. Equation (4) provides the correlation function of the GaSP model while Equation (5) shows the fitted Gaussian prediction function. The GaSP equations shown in both Equations (4) and (5) were adopted from Jones and Johnson [37].
where R(X,θ) is a square correlation matrix wherein it is a function of the design space and points as well as some unknown thetas,μ provides the fitted mean,σ indicates the variance, and r(x,θ) represents the n × 1 vector with a condition if (θ k ≥ 0).

Methodology
Different analytical and simulation software was used to perform the study. The discussion on the methodology is divided into three portions such as the CFD model, the design of experiments through the space-filling design of experiment, and the GaSP model.

Computational Fluid Dynamics
The CFD approach was performed using the commercial software ANSYS Fluent. The K − ε turbulence model was used since this model provides accurate results with the SWT with relatively lesser computational time. These numerical simulations were conducted to obtain the various aerodynamic performance of each design treatment while determining its respective power and moment coefficients.

Simulation Set-Up
The geometries of the developed model were constructed using the Design Modeler of ANSYS. The geometry of the SWT rotor is shown in Figure 1 where the rotor diameter (D) was evaluated to change from 0.29 m to 0.40 m, the rotor height (H) was varied from 0.16 m to 0.35 m, and the twist angle was changed from 15 • to 65 • .
Stochastic process (GaSP) model was considered. The GaSP model consists of an vector for its output response, ( ), which is given as ( 1 , 2 ( , ). The para ( , ) in the output response is a × correlation matrix and can take different f depending on the situation [36]. Equation (4) provides the correlation function of the model while Equation (5) shows the fitted Gaussian prediction function. The GaSP tions shown in both Equations (4) and (5) were adopted from Jones and Johnson [37 where R(X,θ) is a square correlation matrix wherein it is a function of the design and points as well as some unknown thetas, ̂ provides the fitted mean, ̂ indicat variance, and ( ,̂) represents the × 1 vector with a condition if ( ≥ 0).

Methodology
Different analytical and simulation software was used to perform the study. Th cussion on the methodology is divided into three portions such as the CFD mode design of experiments through the space-filling design of experiment, and the model.

Computational Fluid Dynamics
The CFD approach was performed using the commercial software ANSYS F The − turbulence model was used since this model provides accurate results wi SWT with relatively lesser computational time. These numerical simulations were ducted to obtain the various aerodynamic performance of each design treatment determining its respective power and moment coefficients.

Simulation Set-Up
The geometries of the developed model were constructed using the Design Mo of ANSYS. The geometry of the SWT rotor is shown in Figure 1 where the rotor diam (D) was evaluated to change from 0.29 m to 0.40 m, the rotor height (H) was varied 0.16 m to 0.35 m, and the twist angle was changed from 15° to 65°.   In the development of the geometry, a 21.6% percent blockage ratio ( ) [38] was considered as shown in Equation (3). Meanwhile, the upwind distance of the rotor used was 0.30 m while a downstream distance of the rotor utilized was 1.20 m. In the development of the geometry of the SWT, two domains were identified, which are the stationary zone and the rotational zone as shown in Figure 2. The stationary zone is the wind tunnel domain which consists of a rectangular-shaped air domain; while the rotational zone is a cylindrical-shaped domain, as shown in the figure, that allows rotation of the SWT shaft axis.
Wind 2022, 2, FOR PEER REVIEW 5 In the development of the geometry, a 21.6% percent blockage ratio ( ) [38] was considered as shown in Equation (3). Meanwhile, the upwind distance of the rotor used was 0.30 m while a downstream distance of the rotor utilized was 1.20 m. In the development of the geometry of the SWT, two domains were identified, which are the stationary zone and the rotational zone as shown in Figure 2. The stationary zone is the wind tunnel domain which consists of a rectangular-shaped air domain; while the rotational zone is a cylindrical-shaped domain, as shown in the figure, that allows rotation of the SWT shaft axis. As shown in Figure 2, the inlet and the outlet zone were identified, while all the other sides of the rectangular zone were set to a walled domain. An air inlet velocity of 6.5 m/s was used together with an outlet pressure of 1 bar. A convergence criterion of 1 × 10 −6 was used with a maximum of 1000 iterations per simulation run. Figure 3 shows the separated cell zone of the simulation in a section view.

Evaluation of the Simulation Mesh
The quality of the mesh of the geometries for the simulation is an important step to consider to attain an accurate solution with low computational cost. The mesh setup of the developed SWT model is shown in Figure 4. The top view of the rotor that showcase the near wall mesh and the mesh surrounding the rotor is shown in Figure 5. As shown in Figure 2, the inlet and the outlet zone were identified, while all the other sides of the rectangular zone were set to a walled domain. An air inlet velocity of 6.5 m/s was used together with an outlet pressure of 1 bar. A convergence criterion of 1 × 10 −6 was used with a maximum of 1000 iterations per simulation run. Figure 3 shows the separated cell zone of the simulation in a section view.
Wind 2022, 2, FOR PEER REVIEW 5 In the development of the geometry, a 21.6% percent blockage ratio ( ) [38] was considered as shown in Equation (3). Meanwhile, the upwind distance of the rotor used was 0.30 m while a downstream distance of the rotor utilized was 1.20 m. In the development of the geometry of the SWT, two domains were identified, which are the stationary zone and the rotational zone as shown in Figure 2. The stationary zone is the wind tunnel domain which consists of a rectangular-shaped air domain; while the rotational zone is a cylindrical-shaped domain, as shown in the figure, that allows rotation of the SWT shaft axis. As shown in Figure 2, the inlet and the outlet zone were identified, while all the other sides of the rectangular zone were set to a walled domain. An air inlet velocity of 6.5 m/s was used together with an outlet pressure of 1 bar. A convergence criterion of 1 × 10 −6 was used with a maximum of 1000 iterations per simulation run. Figure 3 shows the separated cell zone of the simulation in a section view.

Evaluation of the Simulation Mesh
The quality of the mesh of the geometries for the simulation is an important step to consider to attain an accurate solution with low computational cost. The mesh setup of the developed SWT model is shown in Figure 4. The top view of the rotor that showcase the near wall mesh and the mesh surrounding the rotor is shown in Figure 5.

Evaluation of the Simulation Mesh
The quality of the mesh of the geometries for the simulation is an important step to consider to attain an accurate solution with low computational cost. The mesh setup of the developed SWT model is shown in Figure 4. The top view of the rotor that showcase the near wall mesh and the mesh surrounding the rotor is shown in Figure 5.  The mesh of the developed SWT model was performed to evaluate the skewness value, aspect ratio, and orthogonal quality. A tetrahedral grid shape was adopted as the mesh shape as it was able to capture the complex curvature shape of the SWT. With these factors of interest, a mesh independence test was performed by changing the mesh sizes of each geometry from −100 to 100 as shown in Table 1. Figure 6 illustrates the mesh independent test of the simulation setup in terms of skewness value and orthogonal quality.    The mesh of the developed SWT model was performed to evaluate the sk value, aspect ratio, and orthogonal quality. A tetrahedral grid shape was adopte mesh shape as it was able to capture the complex curvature shape of the SWT. Wi factors of interest, a mesh independence test was performed by changing the me of each geometry from −100 to 100 as shown in Table 1. Figure 6 illustrates the me pendent test of the simulation setup in terms of skewness value and orthogonal q  The mesh of the developed SWT model was performed to evaluate the skewness value, aspect ratio, and orthogonal quality. A tetrahedral grid shape was adopted as the mesh shape as it was able to capture the complex curvature shape of the SWT. With these factors of interest, a mesh independence test was performed by changing the mesh sizes of each geometry from −100 to 100 as shown in Table 1. Figure 6 illustrates the mesh independent test of the simulation setup in terms of skewness value and orthogonal quality.  The design of experiments was performed where each geometry model of the SWT was developed. The developed SWT designs were evaluated into two types of simulations, namely, the static condition and the dynamic condition. The static condition simulation was done using the moving reference frame (MRF) while the dynamic condition simulation was performed using the dynamic sliding mesh motion (SMM).

Design of Experiments
The JMP statistical analysis software was utilized to determine the simulation treatments for the study. The sphere-packing design method under the space-filling design is used to evaluate each design for the different parameters considered. The parametric factors considered are shown in Table 2 where the values were adopted from Lee et al. [38] which showed indicative promising results for the SWT.  Table 2, the total number of simulation treatments obtained from the spacefilling design of the experiment is 30. Hence, 30 unique SWT geometries were prepared and evaluated based on the moment coefficient and power coefficient performances.

Gaussian Stochastic Process (GaSP) Model
After the CFD simulation of the 30 treatments, the GaSP model was adopted to assess the effects of the height and the diameter of the rotor as well as its twist angle with regards to the moment coefficient and power coefficient performances. Afterward, model fitting was performed using the GaSP model to derive the model equations of the input and output factors. It was found that there were two responses to be obtained from this model, namely, the model coefficients and power coefficients. For the GaSP model validation, the model also considered the moment and the coefficient of moment as a factor on top of the considered 3 parametric factors defined in Table 2. These 5 parameters were evaluated based on the power coefficient response. The design of experiments was performed where each geometry model of the SWT was developed. The developed SWT designs were evaluated into two types of simulations, namely, the static condition and the dynamic condition. The static condition simulation was done using the moving reference frame (MRF) while the dynamic condition simulation was performed using the dynamic sliding mesh motion (SMM).

Design of Experiments
The JMP statistical analysis software was utilized to determine the simulation treatments for the study. The sphere-packing design method under the space-filling design is used to evaluate each design for the different parameters considered. The parametric factors considered are shown in Table 2 where the values were adopted from Lee et al. [38] which showed indicative promising results for the SWT. Table 2. Summary of the parameters used in the sphere-packing design method adopted from Lee et al. [38].

Factors Role Values
Rotor From Table 2, the total number of simulation treatments obtained from the spacefilling design of the experiment is 30. Hence, 30 unique SWT geometries were prepared and evaluated based on the moment coefficient and power coefficient performances.

Gaussian Stochastic Process (GaSP) Model
After the CFD simulation of the 30 treatments, the GaSP model was adopted to assess the effects of the height and the diameter of the rotor as well as its twist angle with regards to the moment coefficient and power coefficient performances. Afterward, model fitting was performed using the GaSP model to derive the model equations of the input and output factors. It was found that there were two responses to be obtained from this model, namely, the model coefficients and power coefficients. For the GaSP model validation, the model also considered the moment and the coefficient of moment as a factor on top of the considered 3 parametric factors defined in Table 2. These 5 parameters were evaluated based on the power coefficient response.

Results and Discussion
The results of the study are elaborated in this section along with its discussion.

Gaussian Stochastic Process Results
The results of the model fitting of the 30 simulation treatment results to the GaSP model are shown in Figure 7, which highlights the power coefficient prediction plot for the design treatments.

Results and Discussion
The results of the study are elaborated in this section along with its discussion.

Gaussian Stochastic Process Results
The results of the model fitting of the 30 simulation treatment results to the GaSP model are shown in Figure 7, which highlights the power coefficient prediction plot fo the design treatments. Based on this graph, the actual power coefficient is directly proportional to the pre dicted power coefficient. The 45° plot line shown in the figure signifies that the mode fitted the predicted model. Aside from this, the θ-value of each performance parameter i presented in Table 3. The scale of the θ-values ranges from 0 to 1 which indicates the relative correlation of each of the parameters. This means that the rotor height, the twist angle, and the mo ment have a relatively higher correlation with the respect to the power coefficient, whil the rotor diameter and the moment coefficient have a relatively lower correlation to th power coefficient. Moreover, Figure 8 shows the interaction plot of each of the considered parameters concerning the power coefficient. It was found that there is no interaction ef fect between the response and each performance parameters since the lines in the graph are all parallel. However, it is interesting to note how the range of values for each of th performance parameters varies in terms of the change in the power coefficient. The gap between the red line and the blue line represent the difference of the power coefficient a changing parameters. The larger the gap between the red and the blue lines indicates tha the power coefficient is sensitive to the changing values of the specific parameter. Hence Based on this graph, the actual power coefficient is directly proportional to the predicted power coefficient. The 45 • plot line shown in the figure signifies that the model fitted the predicted model. Aside from this, the θ-value of each performance parameter is presented in Table 3.
The scale of the θ-values ranges from 0 to 1 which indicates the relative correlation of each of the parameters. This means that the rotor height, the twist angle, and the moment have a relatively higher correlation with the respect to the power coefficient, while the rotor diameter and the moment coefficient have a relatively lower correlation to the power coefficient. Moreover, Figure 8 shows the interaction plot of each of the considered parameters concerning the power coefficient. It was found that there is no interaction effect between the response and each performance parameters since the lines in the graphs are all parallel. However, it is interesting to note how the range of values for each of the performance parameters varies in terms of the change in the power coefficient. The gaps between the red line and the blue line represent the difference of the power coefficient at changing parameters. The larger the gap between the red and the blue lines indicates that the power coefficient is sensitive to the changing values of the specific parameter. Hence, no gap between the red and the blue line signifies the insensitivity of the power coefficient on the specified parameter. R PEER REVIEW 9 no gap between the red and the blue line signifies the insensitivity of the power coefficient on the specified parameter. Subsequently, the profiler plots of each performance parameter in terms of the power coefficient are shown in Figure 9a-e. It can be seen in these graphs that the height and the diameter of the rotor are directly proportional to the power coefficient. It is to be noted that the power coefficient increases as the vacuum pressure and the moment length of the blade increases. Similarly, the power coefficient was found to be directly proportional to the moment coefficient as well as the moment. This can be attributed to the direct relationship of the moment or torque from the power calculation of the blades of the wind turbine. Among the performance parameters considered, it is only the twist angle that causes the power coefficient to decrease when the parameter increases. This is due to the decrease of suction pressure and moment arm as the twist angle is increased. Subsequently, the profiler plots of each performance parameter in terms of the power coefficient are shown in Figure 9a-e. It can be seen in these graphs that the height and the diameter of the rotor are directly proportional to the power coefficient. It is to be noted that the power coefficient increases as the vacuum pressure and the moment length of the blade increases. Similarly, the power coefficient was found to be directly proportional to the moment coefficient as well as the moment. This can be attributed to the direct relationship of the moment or torque from the power calculation of the blades of the wind turbine. Among the performance parameters considered, it is only the twist angle that causes the power coefficient to decrease when the parameter increases. This is due to the decrease of suction pressure and moment arm as the twist angle is increased. The resulting fitted model for the performance parameters with respect to the power coefficient is shown in Equation (6). The establishment of the coefficient is shown in the fitted equation as shown below.
where y is the projected power coefficient, D i indicates the diameter, H i shows the height, C mi denotes the moment coefficient, r is the stochastic process in a finite set of i. The values of the fitted mean and the variance for the GaSP model fitting are shown in Table 4 while the factor r' values are shown in Table 5.

Static Condition Analysis
The static condition of the simulation represents the starting position of the wind turbine. The moment coefficients of the new SWT were analyzed and compared with the traditional SWT with a semi-cylindrical shape SWT. The moment coefficients in the static condition of both the new and the old are provided in Figures 10 and 11. It can be seen in Figure 10 that the static moment coefficient of the new design remains positive. It is also shown here that the maximum static moment coefficient can be observed at the angle of attack values of 90 • to 160 • and 450 • to 540 • . Meanwhile, Figure 11 shows that the static moment coefficient of the traditional design crossed through the negative values. Additionally, it can be recognized that the minimum values of static moment coefficients are from 20 • to 80 • and 330 • to 380 • , respectively. With this said, it was found that the optimized design can start by itself since its moment coefficients remain positive while the SWT with a semi-cylindrical shape has negative moment coefficients which hinder its ability to self-start [39]. Therefore, the new SWT has a relatively superior self-starting capability compared to the traditional SWT.
Wind 2022, 2, FOR PEER REVIEW 12 ability to self-start [39]. Therefore, the new SWT has a relatively superior self-starting capability compared to the traditional SWT.   ability to self-start [39]. Therefore, the new SWT has a relatively superior self-starting capability compared to the traditional SWT.

Dynamic Condition Analysis
The parameters that were tested for the dynamic condition include the power coefficient and the dynamic moment coefficient of the wind turbines. Figure 12 showed the comparative study of the different SWT configurations with respect to the dynamic moment coefficient, while Figure 13 provided the comparative study of the different SWT configurations in terms of their power coefficients. It was recognized that the tip-speed ratio (TSR) also played a role in the performance of the SWT. The dynamic moment coefficient of each wind turbine decreases as the TSR is increased. However, the decrease of the new SWT is not as drastic compared to the other SWT configurations. A possible explanation for this is that the new design consists of smaller curvatures that help in maintaining the solution pressure and limit the wake-flow present at the rear of the rotor of the turbine. On the contrary, the power coefficient of each wind turbine followed a bell-curve relative to the TSR where the values increased then decreased until it reached a certain point that indicated its corresponding maximum power coefficient. The maximum power coefficient of the new design was observed at 0.2525 at a TSR of 0.80.

Dynamic Condition Analysis
The parameters that were tested for the dynamic condition include the power coefficient and the dynamic moment coefficient of the wind turbines. Figure 12 showed the comparative study of the different SWT configurations with respect to the dynamic moment coefficient, while Figure 13 provided the comparative study of the different SWT configurations in terms of their power coefficients. It was recognized that the tip-speed ratio (TSR) also played a role in the performance of the SWT. The dynamic moment coefficient of each wind turbine decreases as the TSR is increased. However, the decrease of the new SWT is not as drastic compared to the other SWT configurations. A possible explanation for this is that the new design consists of smaller curvatures that help in maintaining the solution pressure and limit the wake-flow present at the rear of the rotor of the turbine. On the contrary, the power coefficient of each wind turbine followed a bell-curve relative to the TSR where the values increased then decreased until it reached a certain point that indicated its corresponding maximum power coefficient. The maximum power coefficient of the new design was observed at 0.2525 at a TSR of 0.80.

Dynamic Condition Analysis
The parameters that were tested for the dynamic condition include the power coefficient and the dynamic moment coefficient of the wind turbines. Figure 12 showed the comparative study of the different SWT configurations with respect to the dynamic moment coefficient, while Figure 13 provided the comparative study of the different SWT configurations in terms of their power coefficients. It was recognized that the tip-speed ratio (TSR) also played a role in the performance of the SWT. The dynamic moment coefficient of each wind turbine decreases as the TSR is increased. However, the decrease of the new SWT is not as drastic compared to the other SWT configurations. A possible explanation for this is that the new design consists of smaller curvatures that help in maintaining the solution pressure and limit the wake-flow present at the rear of the rotor of the turbine. On the contrary, the power coefficient of each wind turbine followed a bell-curve relative to the TSR where the values increased then decreased until it reached a certain point that indicated its corresponding maximum power coefficient. The maximum power coefficient of the new design was observed at 0.2525 at a TSR of 0.80.

Discussion on the Comparison of the Results with Other Studies
In the dynamic condition analysis of the new SWT, the performance of the new design to the traditional SWT and other designs from past literature were compared. The maximum power coefficient of this new design is relatively higher than the previous designs presented by Wang and Zhan [40], the traditional semi-cylindrical SWT [40], the model of Wekesa et al. [41], and the prototype of Larin et al. [42] by 13.28%, 8.16%, 2.32%, and 1.38%, respectively. Hence, this work highlights the improvement of the maximum power coefficient compared with other studies.

Conclusions
A computational fluid dynamics (CFD) model was developed for a Savonius wind turbine (SWT) that considered various geometric parameters based on a space-filling design of experiments. The space-filling design considered the sphere packing method which worked well with the results of computational simulations such as CFD. A Gaussian stochastic process (GaSP) model was used to model fit the space-filling design results for the dependent variables and the response. With the application of the space-filling design, the study considered 30 distinct designs based on the combination of the parametric factors such as the diameter and height of the rotor as well as its angle of twist. The results have shown that the model fitted the GaSP model. The results revealed that the power coefficient is directly proportional to the performance parameters of the moment, moment coefficient, rotor height, and rotor diameter, while it was found that the power coefficient has shown a slight decrease while the twist angle increases. It was found that the maximum power coefficient of the optimized SWT design is 0.2525 at a tip speed ratio of 0.80. The maximum power coefficient was found to be superior with the SWT design of previous works. In addition, the results revealed that the optimized SWT design is capable of self-starting with relatively better performance compared to the other designs of previous works. Future work involves the simultaneous effects of the overlap ratio of the SWT and the influence of the inlet air velocity.  Data Availability Statement: Data available in a publicly accessible repository. The data presented in this study are openly available in [40,42].