Numerical and Computational Analysis of a New Vertical Axis Wind Turbine , Named KIONAS

Eleni Douvi 1,*, Dimitra Douvi 1,*, Dionissios Margaris 1,* and Ioannis Drosis 2 1 Fluid Mechanics Laboratory (FML), Mechanical Engineering and Aeronautics Department, University of Patras, GR-26500 Patras, Greece 2 Small Wind Turbines Development Manager, 97 Dimitros Street, GR-19200 Elefsina, Greece; giannisdrosis@hotmail.com * Correspondence: douvi@mech.upatras.gr (E.D.); dimdouvi@gmail.com (D.D.); margaris@mech.upatras.gr (D.M.)


Introduction
Since the reserves of fossil fuels are diminishing, there is increasing interest in renewable energy sources, among which wind energy is included.Wind energy is harnessed by wind turbines, which can be either Horizontal Axis Wind Turbines (HAWT) or Vertical Axis Wind Turbines (VAWT).A considerable amount of literature has been published on HAWTs, because they are more effective than the VAWTs.Recently, there has been growing interest in VAWTs, because of their easy installation, manufacture and maintenance.
A review of various configurations of vertical axis wind turbines, along with their advantages and disadvantages, was first conducted by Aslam Bhutta et al. [1].The techniques for VAWT design and the flow field over the blades were also reviewed.The results showed that the power coefficient is different for various configurations and can be optimized with tip speed ratio.The impact of structural parameters on the performance, particularly for Savonius wind turbine, is studied by Tang et al. [2].Their findings revealed that for different structures, the maximum power coefficient can be increased by more than 30%.
Numerous investigations have been conducted to show the effects of guide vanes on vertical axis wind turbine performance.The first systematic study of the effects of guide vanes was reported by Ejiri et al. [3], in 2006.A computational fluid dynamics (CFD) code was utilized to analyze the flow through the vertical axis cross-flow wind turbine and to propose an improved design configuration in order to achieve better performance, which had two guide vanes outside the turbine.
A year later, in 2007, Takao et al. [4] proposed a straight-bladed vertical axis wind turbine with a directed guide vane row, which was added to enhance its torque.The guide vane row improved the performance of the straight-bladed vertical axis turbine and the experimental results showed that the power coefficient of the proposed wind turbine was approximately 1.5 times higher than that of the original wind turbine which has no guide vane.A further study with more of a focus on the effect of guide vane geometry on the performance of wind turbine was conducted by the same researchers [5].An experimental investigation was carried out, in order to investigate the effects of distance between the guide vanes and the number of guide vanes on power and torque coefficients.Recently, Shahizare et al. [6] simulated a vertical axis wind turbine with omni-direction guide vanes of eight different shape ratios to determine the effects of guide vanes.It was concluded that all the shape ratios improved the power and torque coefficient.
Kim and Gharib [7] studied the effect on the power output of the vertical axis wind turbine configuration by adding an upstream deflector.They found out that the dimensions and the relative position of the deflector are related to the power output of the VAWT.The optimal design of the wind booster, a device which controls the airflow around a vertical axis wind turbine, was also analyzed lately by Korprasertsak and Leephakpreeda [8].This device consists of guide vanes mounted around the wind turbine, which not only accelerate but also direct the wind in order to impact the rotor blades at the most effective angle of attack.This results in an increase in the VAWT's angular speed and thus in power.
In order to improve a vertical axis wind turbine's performance, a stator is added around the rotor.The stator significantly increases the wind speed, and as a result it improves the self-starting behavior of the VAWT the power coefficient.The effect of a stator surrounding a vertical-axis wind turbine has been studied both experimental and numerical.
Simulations of a new power-augmented shroud integrated with a vertical axis wind turbine (VAWT) suitable for urban and suburban application were carried out by Wong et al. [9].The results of the numerical simulation indicated that the new design for the power-augmented shroud is able to increase the coefficient of power significantly for the VAWT, by about 147.1% compared to the bare VAWT.At the same time, Nobile et al. [10] proposed and analyzed a computational investigation of an augmented wind turbine.The results of their study also showed that the introduction of an omnidirectional stator around the wind turbine increases not only the power but also the torque coefficients by around 30%-35% when compared to the open case.They also concluded that for optimum performance, attention needs to be given to the orientation of the stator blades.
Chen et al. [11,12] developed a vortical stator assembly (VSA) surrounding a vertical-axis wind turbine to improve rotor performance.The VSA consisted of six guide vanes, in order to guide the air tangentially into the VSA.The experimental results showed that the VSA significantly increased rotor performance and the rotor starting speed was reduced to approximately 1 m/s.The optimal VSA diameter was found to be approximately 1.82 times the rotor diameter.
The first systematic numerical and experimental study of the operating performance and power output from a vertical axis wind turbine was reported by Pope et al. [13].They studied the effects of varying VAWT stator vane geometries of a Zephyr vertical axis wind turbine on the turbine's performance at constant and variable rotor velocities.Recently, Burlando et al. [14] have examined with a numerical wind turbine technique the effects of stator vanes on the flow around and inside a multi-stage vertical-axis wind turbine.They developed a numerical model to extend the results obtained by means of the physical model to more general conditions.
So far this type of vertical axis wind turbine has been applied to hybrid energy systems.Tong et al. [15] were the first to introduce a patented wind-solar hybrid renewable energy harvester for urban high rise application.The speed of the high-altitude free-stream wind is increased through fixed or yaw-able power-augmentation-guide-vane (PAGV) before entering the wind turbine at the center portion.The geometry of the PAGV was optimized by computational fluid dynamics (CFD) simulation.The generated power was 1.25 times higher after integrating the PAGV with the VAWT.Chong et al. [16] introduced the Eco-Greenergy™ hybrid wind-solar photovoltaic energy generation system to power LED lights or other appliances.The system consists of a novel omni-direction guide vane with a vertical axis wind turbine and a PV panel mounted on the top surface.The system was found to generate a total of 572.8 kWh of energy per year, and the ODGV increases the annual wind energy output by 438%.
This project was undertaken to describe the design of a new vertical axis wind turbine, named KIONAS, and evaluate its performance using mathematical modelling and its aerodynamic behavior and flow field through Computational Fluid Dynamics (CFD) simulations.
The wind turbine examined in this study is a vertical axis wind turbine, which is named KIONAS.It consists of a rotor which is surrounded by a stator.In the rotors center there is a shaft, where a number of flat blades are mounted on, at an angle.The stator is added in order to increase the wind speed that impacts the rotor blades.The above system is a one-stage wind turbine KIONAS.The main characteristic of this wind turbine is that a multi-stage turbine can be constructed by adding more of the same stages to the first one, directly above it.The basic geometry of a three-stage wind turbine KIONAS is presented in Figure 1.The current study contributes to our knowledge by drawing the following important conclusions.First of all, it was found that a multi-stage KIONAS could become advantageous in the long term because the output power per unit weight of the machine is higher as the number of KIONAS stages is increased and simultaneously the power of each individual stage is increased when an extra stage is added.Secondly, from the computational results it was obvious that part of the wind is directed both towards the upper and towards the lower stage for the two-stage KIONAS.

Mathematical Study
The first step in the mathematical study is the calculation of the resultant velocity of the air in each of the five sections of the front side of wind turbine KIONAS, Figure 2. The wind speed is considered equal to u and the average wind speed which enters each stator section, u , is The current study contributes to our knowledge by drawing the following important conclusions.First of all, it was found that a multi-stage KIONAS could become advantageous in the long term because the output power per unit weight of the machine is higher as the number of KIONAS stages is increased and simultaneously the power of each individual stage is increased when an extra stage is added.Secondly, from the computational results it was obvious that part of the wind is directed both towards the upper and towards the lower stage for the two-stage KIONAS.

Mathematical Study
The first step in the mathematical study is the calculation of the resultant velocity of the air in each of the five sections of the front side of wind turbine KIONAS, Figure 2. The wind speed is considered equal to u ∞ and the average wind speed which enters each stator section, u in , is where ϕ = (2ν−1)π 10 , for ν = 1, 2, 3, 4, 5. Thus, the average wind speed is: The wind speed that impacts the rotor blades, u, is calculated by the equation of continuity (3) where b 1 is the external and b 2 is the internal cross section length between the guide vanes.By substituting u in from Equation (2) in the above equation, the following equation emerges: for ν = 1, 2, 3, 4, 5.Moreover, each blade has a circumferential velocity, T, which is in the x-axis (ΩR) x = ΩR cos α, (6) and in the y-axis, namely the axis of wind direction, (ΩR) y = ΩR sin α ( In the above equations Ω is the rotational velocity of the rotor, R is the rotor's radius and α is the angle of attack.The resultant velocity is in the x-axis u x = (ΩR) x = ΩR cos α, (8) and in the y-axis u y = u ± (ΩR) y = u ± ΩR sin α, (9) where for the sections denoted by A the quantity ΩR sin α is positive while for the sections denoted by B is negative.The resultant velocity is given by Equation ( 10) By substituting in the above equation u x and u y from Equations ( 8) and ( 9) respectively, we get: u res = (ΩR) 2 sin 2 α + cos 2 α ± 2ΩR sin αu + u 2 ⇒ (13) The resultant velocity is given by the Equation (15), by substituting u from Equation (4) in Equation ( 14): where for the sections denoted by A the quantity 2ΩR sin α(b 1 /b 2 )u ∞ sin (2ν−1)π 10 is positive while for the sections denoted by B is negative.In the above equation Ω is the rotational velocity of the rotor, R is the rotor's radius, α is the angle of attack, b 1 is the external and b 2 is the internal cross section length between the guide vanes, u ∞ is the wind speed and ν is the number of each section, i.e., 1, 2, 3, 4 or 5.
Then, the aerodynamic forces of lift and drag that are exerted on the blades can be determined by Equations ( 16) and ( 17) where L is the lift, D is the Drag, c l and c d are the lift and drag coefficients respectively, ρ air is the density of the air and A is the blade's area.The resultant force, F, which is exerted on the blade and can be situated at any position is calculated by Equation ( 18) Then the aerodynamic torque can be estimated by Equation ( 19): And finally the power output of each blade at each position is given by Equation (20): Computation 2017, 5, 8 5 of 14 length between the guide vanes, ∞ is the wind speed and is the number of each section, i.e., 1, 2, 3, 4 or 5.
Then, the aerodynamic forces of lift and drag that are exerted on the blades can be determined by Equations ( 16) and (17) where L is the lift, D is the Drag, and are the lift and drag coefficients respectively, is the density of the air and A is the blade's area.
The resultant force, F, which is exerted on the blade and can be situated at any position is calculated by Equation ( 18) Then the aerodynamic torque can be estimated by Equation ( 19): And finally the power output of each blade at each position is given by Equation (20): The above equations are solved by the imperative programming language FORTRAN (FORmulae TRANslator), which was originally developed for scientific and engineering applications.The resulting relationships have enabled easy parametric sizing, and so the investigation of the performance of KIONAS on the dimensions of both blades and guide vanes was conducted.Mathematical study is necessary in order to calculate the multi-stage KIONAS power output in less computational time and with less required computational memory compared to the Computational Fluid Dynamics study.

Computational Fluid Dynamics Study
The commercial Computational Fluid Dynamics code ANSYS Fluent [17] was used in order to show the flow field around the vertical axis wind turbine KIONAS.The first step in the simulation process is the design of the geometries of the various configurations of KIONAS in DesignModeler.
A C-type mesh was constructed; in particular KIONAS is surrounded by a semicircle and a horizontal parallelepiped.The semicircle consists of five parts defined as velocity inlets, with The above equations are solved by the imperative programming language FORTRAN (FORmulae TRANslator), which was originally developed for scientific and engineering applications.The resulting relationships have enabled easy parametric sizing, and so the investigation of the performance of KIONAS on the dimensions of both blades and guide vanes was conducted.Mathematical study is necessary in order to calculate the multi-stage KIONAS power output in less computational time and with less required computational memory compared to the Computational Fluid Dynamics study.

Computational Fluid Dynamics Study
The commercial Computational Fluid Dynamics code ANSYS Fluent [17] was used in order to show the flow field around the vertical axis wind turbine KIONAS.The first step in the simulation process is the design of the geometries of the various configurations of KIONAS in DesignModeler.
A C-type mesh was constructed; in particular KIONAS is surrounded by a semicircle and a horizontal parallelepiped.The semicircle consists of five parts defined as velocity inlets, with different velocity magnitudes depending on the relative position of the area between stator vanes and the maximum wind speed location.The rotor and the stator are defined as wall boundary conditions, the rear side as a pressure outlet and the remaining sides as symmetry.
First of all, simulations of the flow over the simplest geometry-that is, KIONAS with one stage-were conducted, in which the stator guide vanes are flat.The rotor blades are also flat plates and their inclination is 25 • and 45 • .After that, simulations were performed for the one-stage KIONAS, with the difference that the stator guide vanes are curved.
Finally, simulations of the flow were performed over the two-stage wind turbine KIONAS with curved stator guide vanes.This type of stator guide vane was found to be more effective and for that reason it was selected for the simulations.The rotor blades inclination was 25 • , 30 • and 35 • .Figure 3 shows the rotor, the stator and the computational domain over the two-stage KIONAS.
Computation 2017, 5, 8 6 of 14 different velocity magnitudes depending on the relative position of the area between stator vanes and the maximum wind speed location.The rotor and the stator are defined as wall boundary conditions, the rear side as a pressure outlet and the remaining sides as symmetry.
First of all, simulations of the flow over the simplest geometry-that is, KIONAS with one stage-were conducted, in which the stator guide vanes are flat.The rotor blades are also flat plates and their inclination is 25° and 45°.After that, simulations were performed for the one-stage KIONAS, with the difference that the stator guide vanes are curved.
Finally, simulations of the flow were performed over the two-stage wind turbine KIONAS with curved stator guide vanes.This type of stator guide vane was found to be more effective and for that reason it was selected for the simulations.The rotor blades inclination was 25°, 30° and 35°. Figure 3 shows the rotor, the stator and the computational domain over the two-stage KIONAS.The grid independence study revealed that a mesh of 1,412,500 cells for the one-stage KIONAS and of 2,825,000 cells for the two-stage KIONAS would be sufficient to establish a grid-independent solution.The mesh was denser near the wind turbine, where changes in the flow occur and thus greater computational accuracy is necessary.
The rotation was simulated by the Moving Reference Frame Model (MRF) [18], a steady-state approximation in which individual cell zones move at different rotational and/or translational speeds, and the most appropriate turbulence model for such simulations is the k-ω shear-stress transport (SST) [19].In the MRF model, the rotor of KIONAS and the domain surrounding it were defined as a rotating reference frame, while the stator of KIONAS, in other words the flow outside the rotor region, was defined as a stationary frame.

Mathematical Results
This section summarizes the results obtained from the mathematical study.Figure 4 presents the total power at each blade's position for wind speeds of 5 m/s and 10 m/s and rotors with various numbers of blades.The results indicate that adding more blades results in higher values of total power for the wind speed that were tested.The results obtained from the analysis of rotors with 2 to 11 blades are summarized in Table 1.The results indicate that there is a significant positive correlation between the numbers of blades and the mean power output of KIONAS.Although most wind turbines achieve an optimized performance with a rotor which consists of three blades, there are plenty smaller scale turbines or at higher speeds, the efficiency of which increases by adding more blades.KIONAS is such a configuration of a vertical axis wind turbine, because not only it is a The grid independence study revealed that a mesh of 1,412,500 cells for the one-stage KIONAS and of 2,825,000 cells for the two-stage KIONAS would be sufficient to establish a grid-independent solution.The mesh was denser near the wind turbine, where changes in the flow occur and thus greater computational accuracy is necessary.
The rotation was simulated by the Moving Reference Frame Model (MRF) [18], a steady-state approximation in which individual cell zones move at different rotational and/or translational speeds, and the most appropriate turbulence model for such simulations is the k-ω shear-stress transport (SST) [19].In the MRF model, the rotor of KIONAS and the domain surrounding it were defined as a rotating reference frame, while the stator of KIONAS, in other words the flow outside the rotor region, was defined as a stationary frame.

Mathematical Results
This section summarizes the results obtained from the mathematical study.Figure 4 presents the total power at each blade's position for wind speeds of 5 m/s and 10 m/s and rotors with various numbers of blades.The results indicate that adding more blades results in higher values of total power for the wind speed that were tested.The results obtained from the analysis of rotors with 2 to 11 blades are summarized in Table 1.The results indicate that there is a significant positive correlation between the numbers of blades and the mean power output of KIONAS.Although most wind turbines achieve an optimized performance with a rotor which consists of three blades, there are plenty smaller scale turbines or at higher speeds, the efficiency of which increases by adding more blades.KIONAS is such a configuration of a vertical axis wind turbine, because not only it is a small scale wind turbine but also its rotor experiences high speeds as the wind velocity is accelerated when it passes through the stator that surrounds the rotor.
Computation 2017, 5, 8 7 of 14 small scale wind turbine but also its rotor experiences high speeds as the wind velocity is accelerated when it passes through the stator that surrounds the rotor.
(a) (b)  An estimation of the power output between KIONAS with blades made of galvanized sheet metal and KIONAS with blades made of aluminum is presented in Figure 5.The greater density of galvanized sheet metal, which means greater weight of the blades, results in lower rotational speed of the rotor and thus less power output.An estimation of the power output between KIONAS with blades made of galvanized sheet metal and KIONAS with blades made of aluminum is presented in Figure 5.The greater density of galvanized sheet metal, which means greater weight of the blades, results in lower rotational speed of the rotor and thus less power output.
Computation 2017, 5, 8 7 of 14 small scale wind turbine but also its rotor experiences high speeds as the wind velocity is accelerated when it passes through the stator that surrounds the rotor.
(a) (b)  An estimation of the power output between KIONAS with blades made of galvanized sheet metal and KIONAS with blades made of aluminum is presented in Figure 5.The greater density of galvanized sheet metal, which means greater weight of the blades, results in lower rotational speed of the rotor and thus less power output.Figure 6 demonstrates the total power at each blade's position for the same wind speeds and different rotor diameters.The rotor diameter is the diameter that occurs when the stator's diameter increases by 150% and 200%, while its width remains constant.From this data it is obvious that the power of KIONAS increases when increasing the rotor's number of blades and diameter.
Computation 2017, 5, 8 8 of 14 increases by 150% and 200%, while its width remains constant.From this data it is obvious that the power of KIONAS increases when increasing the rotor's number of blades and diameter.
(a) (b) Table 2 compares the mean power output in Watts for wind speed equal to 5 m/s and the initial rotor diameter of 0.72 m and the new rotor diameters of 1.37 m and 2.02 m.It can be seen from the data that the mean power output of the wind turbine increases as the rotor diameter and number of blades increase.Figure 7 shows the total power at each blade's position for eleven blades with and without inclination, rotor blade angle of 25 degrees and wind velocity 5 m/s and 10 m/s.It is obvious that the power output is reduced by the inclination of the blades, because it was assumed that the inclined blade has greater length and thus greater weight, and that leads to lower rotational speed.To determine the effects of adding more stages, the total power of different KIONAS types is compared and the results are presented in Figure 8.The results are referred in wind velocity of 5 m/s and diameter of the rotor equal to 0.72 m.There was a significant positive correlation between total Table 2 compares the mean power output in Watts for wind speed equal to 5 m/s and the initial rotor diameter of 0.72 m and the new rotor diameters of 1.37 m and 2.02 m.It can be seen from the data that the mean power output of the wind turbine increases as the rotor diameter and number of blades increase.Figure 7 shows the total power at each blade's position for eleven blades with and without inclination, rotor blade angle of 25 degrees and wind velocity 5 m/s and 10 m/s.It is obvious that the power output is reduced by the inclination of the blades, because it was assumed that the inclined blade has greater length and thus greater weight, and that leads to lower rotational speed.increases by 150% and 200%, while its width remains constant.From this data it is obvious that the power of KIONAS increases when increasing the rotor's number of blades and diameter.
(a) (b) Table 2 compares the mean power output in Watts for wind speed equal to 5 m/s and the initial rotor diameter of 0.72 m and the new rotor diameters of 1.37 m and 2.02 m.It can be seen from the data that the mean power output of the wind turbine increases as the rotor diameter and number of blades increase.Figure 7 shows the total power at each blade's position for eleven blades with and without inclination, rotor blade angle of 25 degrees and wind velocity 5 m/s and 10 m/s.It is obvious that the power output is reduced by the inclination of the blades, because it was assumed that the inclined blade has greater length and thus greater weight, and that leads to lower rotational speed.To determine the effects of adding more stages, the total power of different KIONAS types is compared and the results are presented in Figure 8.The results are referred in wind velocity of 5 m/s and diameter of the rotor equal to 0.72 m.There was a significant positive correlation between total To determine the effects of adding more stages, the total power of different KIONAS types is compared and the results are presented in Figure 8.The results are referred in wind velocity of 5 m/s and diameter of the rotor equal to 0.72 m.There was a significant positive correlation between total power and number of KIONAS stages.Further statistical tests revealed the exact percentage contribution of each stage in the total power of different KIONAS types and Figure 9 shows these results.The percentage increase in power of each stage when it is placed in KIONAS, which has an extra stage, was also investigated and is presented in Figure 10.The most striking result to emerge from the data is that, for example, when the first stage is alone, it has a power of 200.1 W, but when it is placed in a two-stage KIONAS, the first stage power is 249.4W; in other words, it increases 24.6%.Interestingly, if the first stage is placed in a ten-stage KIONAS, its increase is 91.4%.The percentage increase in power of each stage when it is placed in KIONAS, which has an extra stage, was also investigated and is presented in Figure 10.The most striking result to emerge from the data is that, for example, when the first stage is alone, it has a power of 200.1 W, but when it is placed in a two-stage KIONAS, the first stage power is 249.4W; in other words, it increases 24.6%.Interestingly, if the first stage is placed in a ten-stage KIONAS, its increase is 91.4%.The percentage increase in power of each stage when it is placed in KIONAS, which has an extra stage, was also investigated and is presented in Figure 10.The most striking result to emerge from the data is that, for example, when the first stage is alone, it has a power of 200.1 W, but when it is placed in a two-stage KIONAS, the first stage power is 249.4W; in other words, it increases 24.6%.Interestingly, if the first stage is placed in a ten-stage KIONAS, its increase is 91.4%. Figure 11 provides the power coefficient versus KIONAS weight and the power per weight unit for various KIONAS types.There is a clear trend of the power coefficient decreasing with KIONAS weight, due to the fact that the power coefficient is inversely proportional to the swept area, which increases by adding more stages because the surface between the stages is included too.Of particular interest is the correlation between power per KIONAS weight and KIONAS type, because despite the increase in the weight of the machine, which means an increase in the cost of construction, the power output is much higher, so a multi-stage KIONAS could become advantageous in the long term.

Computational Fluid Dynamics Results
In this section, the results obtained from the Computational Fluid Dynamics are summarized.In Figures 12 and 13 the static pressure distribution on the rotor with rotor blades angles of 45° and 25°, respectively, and flat stator guide vanes, which was calculated by ANSYS Fluent, is presented.It is obvious that the contested surface of the rotor has higher values of static pressure and the rotor rotates due to this pressure difference.The values of static pressure on the rotor blades with inclination of 25° are higher than the corresponding values of inclination of 45°.
Simulation of the flow over a one-stage KIONAS with rotor blade angle of 25° and stator with curved guide vanes was also conducted and the results of the static pressure distribution on the rotor are presented in Figure 14.From the comparison between Figures 13 and 14, it is concluded that in the case of curved stator vanes, the static pressure is greater.Figure 11 provides the power coefficient versus KIONAS weight and the power per weight unit for various KIONAS types.There is a clear trend of the power coefficient decreasing with KIONAS weight, due to the fact that the power coefficient is inversely proportional to the swept area, which increases by adding more stages because the surface between the stages is included too.Of particular interest is the correlation between power per KIONAS weight and KIONAS type, because despite the increase in the weight of the machine, which means an increase in the cost of construction, the power output is much higher, so a multi-stage KIONAS could become advantageous in the long term.Figure 11 provides the power coefficient versus KIONAS weight and the power per weight unit for various KIONAS types.There is a clear trend of the power coefficient decreasing with KIONAS weight, due to the fact that the power coefficient is inversely proportional to the swept area, which increases by adding more stages because the surface between the stages is included too.Of particular interest is the correlation between power per KIONAS weight and KIONAS type, because despite the increase in the weight of the machine, which means an increase in the cost of construction, the power output is much higher, so a multi-stage KIONAS could become advantageous in the long term.

Computational Fluid Dynamics Results
In this section, the results obtained from the Computational Fluid Dynamics are summarized.In Figures 12 and 13 the static pressure distribution on the rotor with rotor blades angles of 45° and 25°, respectively, and flat stator guide vanes, which was calculated by ANSYS Fluent, is presented.It is obvious that the contested surface of the rotor has higher values of static pressure and the rotor rotates due to this pressure difference.The values of static pressure on the rotor blades with inclination of 25° are higher than the corresponding values of inclination of 45°.
Simulation of the flow over a one-stage KIONAS with rotor blade angle of 25° and stator with curved guide vanes was also conducted and the results of the static pressure distribution on the rotor are presented in Figure 14.From the comparison between Figures 13 and 14, it is concluded that in the case of curved stator vanes, the static pressure is greater.

Computational Fluid Dynamics Results
In this section, the results obtained from the Computational Fluid Dynamics are summarized.In Figures 12 and 13 the static pressure distribution on the rotor with rotor blades angles of 45 • and 25 • , respectively, and flat stator guide vanes, which was calculated by ANSYS Fluent, is presented.It is obvious that the contested surface of the rotor has higher values of static pressure and the rotor rotates due to this pressure difference.The values of static pressure on the rotor blades with inclination of 25 • are higher than the corresponding values of inclination of 45 • .
Simulation of the flow over a one-stage KIONAS with rotor blade angle of 25 • and stator with curved guide vanes was also conducted and the results of the static pressure distribution on the rotor are presented in Figure 14.From the comparison between Figures 13 and 14, it is concluded that in the case of curved stator vanes, the static pressure is greater.Figure 15 shows the static pressure distribution on the rotor with a rotor blade angle of 25° and stator with curved guide vanes of two-stage KIONAS.The static pressure is greater on the contested surface of the rotors and on the upper rotor.This result implies that the rotational velocity of the rotor which is located on the upper stage is greater, and thus the power of the upper stage is greater as well.This finding is consistent with Figure 9, where the power of the second stage of a two-stage KIONAS is greater than the power of the first stage.Figure 15 shows the static pressure distribution on the rotor with a rotor blade angle of 25° and stator with curved guide vanes of two-stage KIONAS.The static pressure is greater on the contested surface of the rotors and on the upper rotor.This result implies that the rotational velocity of the rotor which is located on the upper stage is greater, and thus the power of the upper stage is greater as well.This finding is consistent with Figure 9, where the power of the second stage of a two-stage KIONAS is greater than the power of the first stage.Figure 15 shows the static pressure distribution on the rotor with a rotor blade angle of 25° and stator with curved guide vanes of two-stage KIONAS.The static pressure is greater on the contested surface of the rotors and on the upper rotor.This result implies that the rotational velocity of the rotor which is located on the upper stage is greater, and thus the power of the upper stage is greater as well.This finding is consistent with Figure 9, where the power of the second stage of a two-stage KIONAS is greater than the power of the first stage.Figure 15 shows the static pressure distribution on the rotor with a rotor blade angle of 25 • and stator with curved guide vanes of two-stage KIONAS.The static pressure is greater on the contested surface of the rotors and on the upper rotor.This result implies that the rotational velocity of the rotor which is located on the upper stage is greater, and thus the power of the upper stage is greater as well.This finding is consistent with Figure 9, where the power of the second stage of a two-stage KIONAS is greater than the power of the first stage.

Conclusions
This project was undertaken to describe the design and evaluate the aerodynamic performance of the wind turbine KIONAS.This study/research has shown that the power output of the wind turbine increases as the diameter and the number of rotor blades increase.It was also shown that the power output for KIONAS with inclined blades is reduced, due to the increased mass of each blade when it has any inclination.Another major finding was that the curved guide vanes in the stator of KIONAS result in higher power.
The most significant findings to emerge from this study are not only that the power of each individual stage is increased when an extra stage is added but also that the output power per unit weight of the machine is higher as the number of KIONAS stages is increased.That means that a multi-stage KIONAS could become advantageous in the long term.
The results of the computational investigation show that the static pressure is greater in the contested surface of KIONAS and the pressure difference between the front and the rear side of the rotor is responsible for the rotor rotation.The velocity vectors for the two-stage KIONAS indicate that part of the wind is directed both towards the upper and towards the lower stage.Further experimental investigations would be of great help in assessing the actual attitude of KIONAS and this information would help to include more factors affecting the flow in the mathematical analysis.
The evidence from this study suggests that, in general, the present wind turbine cannot compete with the large structures that produce power in the range of 1, 2 and 3 MW, but is a major competitor with the smallest structures of the 5, 10 and 100 kW range.This power may vary with larger diameter dimensions, lighter material, but also with higher wind speeds, in order to cover the needs of individual houses or small settlements.A great benefit of KIONAS is the simple construction of both the guide vanes and the rotor blades, along with its compact shape which can easily be mounted in spaces of small dimensions either on buildings or between them, without creating problems in terms of adaptability in space or public acceptance.

Conclusions
This project was undertaken to describe the design and evaluate the aerodynamic performance of the wind turbine KIONAS.This study/research has shown that the power output of the wind turbine increases as the diameter and the number of rotor blades increase.It was also shown that the power output for KIONAS with inclined blades is reduced, due to the increased mass of each blade when it has any inclination.Another major finding was that the curved guide vanes in the stator of KIONAS result in higher power.
The most significant findings to emerge from this study are not only that the power of each individual stage is increased when an extra stage is added but also that the output power per unit weight of the machine is higher as the number of KIONAS stages is increased.That means that a multi-stage KIONAS could become advantageous in the long term.
The results of the computational investigation show that the static pressure is greater in the contested surface of KIONAS and the pressure difference between the front and the rear side of the rotor is responsible for the rotor rotation.The velocity vectors for the two-stage KIONAS indicate that part of the wind is directed both towards the upper and towards the lower stage.Further experimental investigations would be of great help in assessing the actual attitude of KIONAS and this information would help to include more factors affecting the flow in the mathematical analysis.
The evidence from this study suggests that, in general, the present wind turbine cannot compete with the large structures that produce power in the range of 1, 2 and 3 MW, but is a major competitor with the smallest structures of the 5, 10 and 100 kW range.This power may vary with larger diameter dimensions, lighter material, but also with higher wind speeds, in order to cover the needs of individual houses or small settlements.A great benefit of KIONAS is the simple construction of both the guide vanes and the rotor blades, along with its compact shape which can easily be mounted in spaces of small dimensions either on buildings or between them, without creating problems in terms of adaptability in space or public acceptance.

Figure 1 .
Figure 1.(a) The front view; (b) the top view of the basic geometry of a three-stage wind turbine KIONAS.Units in mm.

Figure 1 .
Figure 1.(a) The front view; (b) the top view of the basic geometry of a three-stage wind turbine KIONAS.Units in mm.

Figure 2 .
Figure 2. Analysis of the half top view of the stator, in order to calculate the magnitude of velocity inlet and the resultant velocity at various stator sections.

Figure 2 .
Figure 2. Analysis of the half top view of the stator, in order to calculate the magnitude of velocity inlet and the resultant velocity at various stator sections.

Figure 4 .
Figure 4. Total power at each blade's position for wind velocity (a) 5 m/s; (b) 10 m/s for various numbers of blades.

Figure 5 .
Figure 5.Total power at each blade's position, for eleven rotor blades made of aluminum or galvanized sheet metal.

Figure 6
Figure6demonstrates the total power at each blade's position for the same wind speeds and different rotor diameters.The rotor diameter is the diameter that occurs when the stator's diameter

Figure 4 .
Figure 4. Total power at each blade's position for wind velocity (a) 5 m/s; (b) 10 m/s for various numbers of blades.

Figure 4 .
Figure 4. Total power at each blade's position for wind velocity (a) 5 m/s; (b) 10 m/s for various numbers of blades.

Figure 5 .
Figure 5.Total power at each blade's position, for eleven rotor blades made of aluminum or galvanized sheet metal.

Figure 6 Figure 5 .
Figure6demonstrates the total power at each blade's position for the same wind speeds and different rotor diameters.The rotor diameter is the diameter that occurs when the stator's diameter

Figure 6 .
Figure 6.Total power at each blade's position for wind velocity 5 m/s, various numbers of blades and rotor's diameter (a) 1.37 m; (b) 2.02 m.

Figure 7 .
Figure 7.Total power at each blade's position for eleven rotor blades with 25° inclination and without inclination and wind velocity (a) 5 m/s; (b) 10 m/s.

Figure 6 .
Figure 6.Total power at each blade's position for wind velocity 5 m/s, various numbers of blades and rotor's diameter (a) 1.37 m; (b) 2.02 m.

Figure 6 .
Figure 6.Total power at each blade's position for wind velocity 5 m/s, various numbers of blades and rotor's diameter (a) 1.37 m; (b) 2.02 m.

Figure 7 .
Figure 7.Total power at each blade's position for eleven rotor blades with 25° inclination and without inclination and wind velocity (a) 5 m/s; (b) 10 m/s.

Figure 7 .
Figure 7.Total power at each blade's position for eleven rotor blades with 25 • inclination and without inclination and wind velocity (a) 5 m/s; (b) 10 m/s.
of KIONAS stages.Further statistical tests revealed the exact percentage contribution of each stage in the total power of different KIONAS types and Figure9shows these results.

Figure 8 .
Figure 8.Total power output in Watts for various KIONAS types.

Figure 9 .
Figure 9.Each stage's contribution in the total power of each KIONAS type.

Figure 8 .
Figure 8.Total power output in Watts for various KIONAS types.

Figure 8 .
Figure 8.Total power output in Watts for various KIONAS types.

Figure 9 .
Figure 9.Each stage's contribution in the total power of each KIONAS type.

Figure 9 .
Figure 9.Each stage's contribution in the total power of each KIONAS type.

Figure 10 .
Figure 10.Each stage's increase in power by adding additional stages.

Figure 10 .
Figure 10.Each stage's increase in power by adding additional stages.

Figure 10 .
Figure 10.Each stage's increase in power by adding additional stages.

Figure 12 .Figure 13 .Figure 14 .
Figure 12.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of a one-stage KIONAS, rotor blade angle of 45 degrees and stator with flat guide vanes.

Figure 12 .
Figure 12.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of a one-stage KIONAS, rotor blade angle of 45 degrees and stator with flat guide vanes.

Figure 12 .Figure 13 .Figure 14 .
Figure 12.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of a one-stage KIONAS, rotor blade angle of 45 degrees and stator with flat guide vanes.

Figure 13 .Figure 12 .Figure 13 .Figure 14 .
Figure 13.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of one-stage KIONAS, rotor blade angle of 25 degrees and stator with flat guide vanes.

Figure 14 .
Figure 14.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of one-stage KIONAS, rotor blade angle of 25 degrees and stator with curved guide vanes.

Figure 15 .Figure 16 .Figure 17 .
Figure 15.Static pressure (Pa) distribution (a) on the contested surface; (b) on the non-contested surface of two-stage KIONAS, rotor blade angle of 25 degrees and stator with curved guide vanes.From Figures 16-18 the velocity vectors and the contours of velocity for two-stage KIONAS, stator with curved guide vanes and rotor blade angle of 25°, 30° and 35° are illustrated.The most interesting finding from these Figures is that part of the wind is directed towards the top and towards the lower stage for all the cases that examined, and this move of air mass into KIONAS structure affects the power of each stage when it is placed in KIONAS with an extra stage.

Figure 16 .Figure 15 .Figure 16 .Figure 17 .
Figure 16.(a) Velocity vectors (m/s); (b) contours of velocity (m/s) on xy plane and the center of the rotor for two-stage KIONAS, rotor blade angle of 25 • and stator with curved guide vanes.

Figure 18 .
Figure 18.(a) Velocity vectors (m/s); (b) contours of velocity (m/s) on xy plane and the center of the rotor for two-stage KIONAS, rotor blade angle of 35° and stator with curved guide vanes.

Figure 18 .
Figure 18.(a) Velocity vectors (m/s); (b) contours of velocity (m/s) on xy plane and the center of the rotor for two-stage KIONAS, rotor blade angle of 35 • and stator with curved guide vanes.

Table 1 .
Mean power output in Watts for wind speed of 5 m/s and 10 m/s for various numbers of blades.

Table 1 .
Mean power output in Watts for wind speed of 5 m/s and 10 m/s for various numbers of blades.

Table 1 .
Mean power output in Watts for wind speed of 5 m/s and 10 m/s for various numbers of blades.

Table 2 .
Mean power output in Watts for various rotor diameters and various numbers of blades.

Table 2 .
Mean power output in Watts for various rotor diameters and various numbers of blades.

Table 2 .
Mean power output in Watts for various rotor diameters and various numbers of blades.