Catchment Based Aerodynamic Performance Analysis of Small Wind Turbine Using a Single Blade Concept for a Low Cost of Energy

: For low and medium wind conditions, there is a possibility to harness maximum wind potential reducing the cost of energy by employing catchment-based wind turbine designs. This paper aims to study catchment-based small wind turbine aerodynamic performance for improved e ﬃ ciency and reduced cost of energy. Hence, design parameters are considered based on speciﬁc conditions within a catchment area. The bins and statistical methods implemented with Weibull distribution of wind data for selected sites to characterize the wind conditions and a weighted average method proposed to create representative wind conditions implementing a single blade concept. The blade element method was applied using Matlab code (version R2017a, MathWorks Inc., Natick, MA, US) for aerodynamic design and analysis, and computational ﬂuid dynamics employed using ANSYS—Fluent (version 18.1, ANSYS Inc., Canonsburg, PA, USA) for validation. The performance of the designed blade is evaluated based on annual energy production, capacity factor and power coe ﬃ cient. Then, for site-speciﬁc wind conditions, yearly energy production, and relative cost of energy are examined against rated power. Appropriate rated power for a low cost of energy identiﬁed and performance measures evaluated for each site. As a result, a maximum power coe ﬃ cient of around 51.8% achieved at a design wind speed of 10 m / s, and higher capacity factors of 28% and 50.9% respectively attained for the low and high wind conditions at the proposed rated powers. Therefore, for di ﬀ erent wind condition sites, enhanced performance at a low cost of energy could be achieved using a single blade concept at properly selected rated powers employing suitable design conditions and procedures.


Introduction
With a growing installation of high wind speed turbines globally, low and medium wind speed turbines are getting more attention [1]. For low and medium wind speed conditions, specific wind turbines shall be designed to utilize the maximum possible wind potential of the site, to reduce loading and costs of energy [2]. Fuglsang et al. showed the potential for a site-specific design that design loads and annual energy production depend on-site with considerably large maximum variations comparing two typical wind turbines at representative wind conditions. The study showed that site-specific turbines generated a smaller cost of energy than a fixed design turbine in every wind climate and appeared to offer great benefits for low-wind-speed sites. In addition to this, the study indicated that the redesign of only blades and towers results in significant reductions in the cost of energy and On the other hand, controllers that can track maximum power are commonly used to ensure maximum energy yield and safety of the wind turbine. A review of the available maximum power point tracking algorithms shows that there are different options to consider depending on the parameters to control. These are tip speed ratio, optimal torque, power signal feedback, and perturbation and observation controls [17]. Along with, for monitoring and maintenance purposes, the state of the wind turbine can be predicted using sensor measurements with repeat-accumulate coded communication methods remotely [18]. A method of a particle filter, which is a nonlinear state estimator, can also be used to estimate the wind turbine parameters like rotor speed, tower top velocity, and tower top displacement [19] to enhance the reliability of the wind turbines.
Even though a lot is done, on the aerodynamic design and analysis of small wind turbines, focusing on the power performance, site-specific conditions, the cost of energy, and related issues do not get due attention. Hence, considering these issues in the design and analysis process are expected to enhance the performance and reduce associated costs. Besides, for off-grid application with one or few installed turbines per site, designing and producing different wind turbine blades is costly for the reason that design and mold making costs are highly expensive when utilizing different blades than a single blade. On the other hand, for sites with various wind conditions, using a single blade at the design rated power might result in a higher cost of energy and low capacity factor than better possible ones.
Thus, the paper considers a single blade concept for use in all the sites in the catchment area. The bins and statistical methods are employed to the wind data of selected sites to characterize the wind condition.A proposed weighted average method is applied to create representative wind characteristics for designing the aerodynamics of the blade using a single blade concept. Catchment-based horizontalaxis wind turbine aerodynamics is analyzed implementing BEM using Matlab code validating it with CFD using Ansys-Fluent. The performance of the designed blade is evaluated based on the annual energy production, capacity factor, and power coefficient applying a maximum power point tracking mechanism. In the end, it examines the relationships between yearly energy production and rated power in terms of the site-specific wind conditions. It also assesses the performance of the blade at each site with the single blade concept and recommending rated-powers for a low cost of energy.

Study Area and Wind Data Sites
A study on wind energy resource potential shows that Ethiopia has total wind energy reserves of 3030 GW with exploitable potential of 1599 GW and installable potential of 1350 GW [20,21]. The sites considered in this paper are in the northern part of Ethiopia called Geba catchment located in the surrounding of the city of Mekelle, found at a latitude of 13.48 • N, 39.49 • E, and around 2208 m elevation. Wind data has been recorded at four wind measurement sites around Mekelle for at least two years. The summary of the selected study sites in the Geba catchment is as shown in Table 1.

Wind Resource Characteristics
The bins method will be implemented based on the available records and best practice to summarize wind data. In addition to this, statistical methods are considered as per Manwell [8] to define the wind energy potential of the sites and evaluated the expected turbine productivity of the designed wind turbine. Considering the wind data categorized into N b bins of width w i , with midpoint wind speeds of v i , and frequency of occurrence f i , in each bin, the total number of data N is defined as: Then the average speed V and standard deviation σ V of the wind are defined as: The average power density P/A, is expressed as: where A is the swept area of the rotor, and ρ is the density of air. The average machine power P w , and the annual energy production by the rotor are then: where P w (v i ) is the power produced defined by the power curve of the rotor, and ∆t is the time of operation in a year.
In the statistical analysis, the probability distribution is used to describe the likelihood of the occurrence of certain wind speeds. The probability density function defines the probability of occurrence of wind speed between V i and V i+1 . There are a few probability distributions in the literature to represent wind behavior. Among these, the Rayleigh and the Weibull are the most often used ones [8]. The Weibull distribution, which needs information of two factors, a shape parameter k, and shape factor c, better represents the wind behavior as it incorporates both the exponential and Rayleigh distributions [8,22].
The Weibull probability density function f (V) and the cumulative distribution function F(V) are defined as: where the shape parameter k, for 1 ≤ k < 10, and the shape factor c, are defined as: where Γ(1 + 1 /k) is gamma function defined: Then, using the statistical methods, an equivalent value of the average machine power Pw in Equation (5), is evaluated as: which, in turn, is used to calculate capacity factor CF that is a performance parameter of the wind turbine at a specific site, defined as the ratio of the average machine power to the rated power P r over a given time. Thus: CF = P w P r (10)

Evaluating Representative Wind Parameters
Design and production of wind turbine blades are expensive for each different nearby sites with different wind resources or wind power production potentials within the same catchment area. Therefore, it is important to design a single blade based on a weighted average parameter of the wind potential of the sites. The weighted average parameter proposed in this paper wasthe wind power density per unit area of the selected sites within the catchment area. Therefore, the proposed weighting factor is: where W j is the weighting factor for site j = 1, 2, 3 . . . n, and n is the number of sites within the catchment area with available data.
Then, the representative parameters are the sum of the weighted parameters of each site. For instance, consider the midpoint wind speed of v i , for the bins, the weighted average midpoint wind speeds v ia , is: Based on these, weighted average values of the parameter wind bins, Weibull probability density function, and cumulative distributive function parameters will be available to design and analyze the representative small wind turbine blade.

Cost of Energy
The cost of energy (CoE) is one of the key parameters required in the design and study of wind turbines, defined as the production cost of unit energy of the rotor [8,20]. Thus: where C ce is the capital cost of the wind turbine, FCR is a fixed charge rate, and C O&M is the operational cost. In this paper, a variation of lifetime cost is estimated in terms of changes in power rating to evaluate relative-CoE assuming elements of the cost fraction that vary with changes of rated power. Considering this assumption is reasonable for limited scopes related to cost effects. Hence, it is assumed that cost of power electronics varies proportional to the rated power. Blade cost also increases with power rating levels due to more cyclic loads around the rated-power and is likely to vary as to the square root of power [6]. The gearbox is not considered in the scaling analysis as direct-drive small wind turbines are assumed. The costs applied are assuming the component level CoE contribution reported by National Renewable Energy Laboratory (NREL) [23]. Thus, relative-CoE is defined evaluating Equation (13) Energies 2020, 13, 5838 6 of 20 accordingly as a ratio of nominal CoE value of the respective sites. The rated power corresponding to a relative-CoE of one offers the operation of the blade with minimum CoE at the respective site.

Selection of Design Parameters
The paper aimed to design a 5 kW small wind turbine with appropriate parameters that suit the selected site conditions. Hence, the use of the best design parameters is very important to maximize the performance of the rotor. Among these, the tip-speed ratio and the number of blades are the most common ones that need careful selection. Hau and Barton et al. illustrates the variances in the efficiencies for rotors of different configurations qualitatively [24,25]. Accordingly, a three-bladed horizontal axis wind turbine at a tip-speed ratio of seven demonstrates maximum efficiency in comparison to the others. Therefore, a three-bladed wind turbine at a tip speed ratio of seven is the selected most efficient combination.
The design wind speed is another parameter that requires careful determination. There are different methods to select the design wind speed. For design considerations based on the wind classes, the International Electro-technical Commission (IEC) standard provides the average wind speeds for each wind class, and the design wind speed recommended is 40% more than the average wind speed [26]. On the other hand, with available wind data for a specific site, the recommended practice is selecting the wind speed at around the highest power density. The practice based on the power density distribution allows us to design the blade at around wind speed corresponding to the highest wind power distribution. It will enable the wind turbine to harness the power at the wind speed with the highest power density at the optimized design. Usually, the mean wind speed of a specific site is slightly higher than the most repeated one. Thus, the design wind speed should be slightly higher than the one that corresponds to the densest wind power. This method of selecting design wind speed is a common practice in the wind-turbine design process [27].
In addition to these density and the kinematic viscosity of air affects the design of wind turbines. For the selected catchment area, these values are taken based on the elevation, annual average temperature, and atmospheric pressure. Therefore, considering the data for environmental conditions of the Geba catchment, an average temperature of 23 • C and an atmospheric pressure of 760 hPa are considered. Accordingly, using the ideal gas law, an average air density of 0.92 kg/m 3 and kinematic viscosity of 1.84 × 10 −5 were found.
Airfoils are the major parts of the section requirements of the rotor blade design. These days it is easier to find and use existing airfoils for small wind turbine applications. National Renewable Energy Laboratory NREL (wind.nrel.gov/airfoils) and other similar institutions have been developing some groups of airfoils for wind rotors. The airfoils with decent performance features for small rotors are listed in the literature [15,28,29]. Hence, a set of airfoils of NREL, namely S823 and S822, developed for small wind turbines (2-20 kW) at various conditions, were found appropriate to the sites considered. A rectangular section is selected as the third airfoil to connect the rotor blades with the hub plate at the root part. Lift coefficient C l , and drag coefficient C d , at a specific angle of attack α that corresponds to around the maximum lift to drag ratio was taken for the NREL airfoils while the recommended force coefficient by the IEC 61400-2 2006 [26] was considered for the rectangular section.

Blade Element Method and Iterative Process
The BEM is the most typically used method for the aerodynamic and aero-elastic study. Different modifications and improvements have been made, in the design of rotor blades, for better evaluation of forces and velocities. Madson and Dossing explained the BEM [30,31] that is also commonly available in the literature [7,8]. Therefore the modified BEM outlined by Madson and Dossing will be implemented [30,31], in combination with the tip loss factor modified by Shen [32], which is a modification of the Glauert's tip loss factor that gives an improvement in predicting the aerodynamic forces.
An iterative optimization process has been undertaken in two steps using Matlab code. The first step was the optimization process in the BEM. Here the tangential and axial induction factors have been Energies 2020, 13, 5838 7 of 20 the constraints required to converge to an ideal value within an acceptable tolerance level. The flow angle or the angle of attack, force coefficients, and element forces were variables to be determined and the tip speed ratio was considered as constant when designing for selected values. Taking the results of the first step as input, the second step considered the design power as the constraint. In the second step, the number of blades, wind velocity, density, and viscosity were the constants while the radius, angular velocity, and chord length were considered as variables.

CFD Modeling and Simulation
Ansys-Fluent, a widely used CFD modeling software, is used for the CFD simulation of the rotor. The CFD modeling considers one-third of a sector of the full cylindrical domain using one blade only with periodic boundary conditions to reduce computation time considerably. Figure 1 shows the boundary conditions and computational-domain of the simulation. The cylindrical computational-domain at the upstream inlet waslocated in front of the blade at three times the blade radius (3R), and its diameter wasfive times the diameter of the blade (5D). On the downstream side, it waslocated at 10R behind the rotor blade with a computational domain cylinder diameter of 7D. These dimensions of the computational domain wereexpected to allow wake expansion and avoid unfavorable boundary dependency effects [33,34]. The rotational frame wasconsidered in the computational-domain to account for the rotor speed, with the blade taken as a non-slip wall boundary condition. Inflation layers of the prismatic type wereused forthe rotor face together with unstructured meshes for the computational-domain to get an improved resolution of the flow in boundary layers to get the wall distance y+ to be less than one [34]. The wall distance y+ is a non-dimensional distance from the wall to the first mesh node used to describe the height of the first grid element next to a wall in a CFD simulation. An iterative optimization process has been undertaken in two steps using Matlab code. The first step was the optimization process in the BEM. Here the tangential and axial induction factors have been the constraints required to converge to an ideal value within an acceptable tolerance level. The flow angle or the angle of attack, force coefficients, and element forces were variables to be determined and the tip speed ratio was considered as constant when designing for selected values. Taking the results of the first step as input, the second step considered the design power as the constraint. In the second step, the number of blades, wind velocity, density, and viscosity were the constants while the radius, angular velocity, and chord length were considered as variables.

CFD Modeling and Simulation
Ansys-Fluent, a widely used CFD modeling software, is used for the CFD simulation of the rotor. The CFD modeling considers one-third of a sector of the full cylindrical domain using one blade only with periodic boundary conditions to reduce computation time considerably. Figure 1 shows the boundary conditions and computational-domain of the simulation. The cylindrical computational-domain at the upstream inlet waslocated in front of the blade at three times the blade radius (3R), and its diameter wasfive times the diameter of the blade (5D). On the downstream side, it waslocated at 10R behind the rotor blade with a computational domain cylinder diameter of 7D. These dimensions of the computational domain wereexpected to allow wake expansion and avoid unfavorable boundary dependency effects [33,34]. The rotational frame wasconsidered in the computational-domain to account for the rotor speed, with the blade taken as a non-slip wall boundary condition. Inflation layers of the prismatic type wereused forthe rotor face together with unstructured meshes for the computational-domain to get an improved resolution of the flow in boundary layers to get the wall distance y+ to be less than one [34]. The wall distance y+ is a non-dimensional distance from the wall to the first mesh node used to describe the height of the first grid element next to a wall in a CFD simulation. In CFD simulations of wind turbine blades, the fluid-flow was analyzed based on the Navier-Stokes equations. RANS equation is investigated usually with either two-equation turbulence models of the k-ε and k-ω series or four-equation shear stress transport k-ω SST turbulence model. The k-ε turbulence model could not depict flows at the boundary layers accurately. It is, therefore, suitable for far-field flow simulations. On the other hand, the k-ω turbulence model could not represent far field flows precisely. Hence, it is appropriate for modeling  In CFD simulations of wind turbine blades, the fluid-flow was analyzed based on the Navier-Stokes equations. RANS equation is investigated usually with either two-equation turbulence models of the k-ε and k-ω series or four-equation shear stress transport k-ω SST turbulence model. The k-ε turbulence model could not depict flows at the boundary layers accurately. It is, therefore, suitable for far-field flow simulations. On the other hand, the k-ω turbulence model could not represent far field flows precisely. Hence, it is appropriate for modeling the boundary layers. The transition k-ω SST turbulence uses the k-ε model and k-ω model alternatively depending on the flow condition. It has the capability of flow separation forecasting on the rotor and reflects the effect of free-stream pressure gradients and turbulence [33,35,36]. Therefore, a four-equation shear stress transport, k-ω SST turbulence model, was employed in the CFD analysis. As per the assessment and investigations of different researchers, the k-ω SST turbulence model was reliable for flow analysis with good agreement with experimental results. However, at large angles of attack, it was confirmed that the CFD simulation results tend to slightly underpredict the power and overpredict the thrust [37][38][39][40][41]. An incompressible analysis was employed, as most of the wind turbine operation is in the subsonic region, with a constant local air density and viscosity. A coupled algorithm for pressure, which considerably increases the convergence rate, was considered in solving the incompressible Reynold-averaged Navier-Stokes (RANS) model [33,42].
Residuals and net mass imbalance were engaged to limit the running time of the solver and ensure convergence. Residuals of the specific dissipation ω, the kinetic energy turbulence k, velocities, and continuity were limited to less than 10 −4 during the solution process. In addition to these, the net imbalance was also checked to be below 0.1% [33,37,43], and an integral static pressure surface monitor to remain constant for a considerable amount of iterations after convergence.
Furthermore, ASHES (simis.io) was employed as an off-design performance evaluation tool. It is a time-domain wind-turbine aerodynamics module, which uses AeroDyne of FAST with a rich user interface. Parameters like a chord, twist angle, pitch axis, airfoil type, and shapefiles of the airfoils at their corresponding radial positions from Matlabwere exported to ASHES as input in the table format. It internally implements xfoil to determine the lift and drag constants for the shapefiles to get the necessary parameters in the evaluation process.

Wind Distribution Characteristics and Design Parameters
The wind characteristics of the four sites in the Geba catchment area as bins of power density distribution and wind speed distribution are shown in Figure 2a-d. The figures show also graphs of the Weibull probability wind distribution and power density distribution. Figure 2e shows weighted average distribution over the catchment based on the data from the four sites. Figure 2f shows the corresponding cumulative distribution graphs of the sites the weighted average over the catchment. The results in Figure 2e,f show that the weighted average is close to the wind conditions of the Mayderhu site.
Energies 2020, 13, x FOR PEER REVIEW 8 of 21 the boundary layers. The transition k-ω SST turbulence uses the k-ε model and k-ω model alternatively depending on the flow condition. It has the capability of flow separation forecasting on the rotor and reflects the effect of free-stream pressure gradients and turbulence [33,35,36]. Therefore, a four-equation shear stress transport, k-ω SST turbulence model, was employed in the CFD analysis. As per the assessment and investigations of different researchers, the k-ω SST turbulence model was reliable for flow analysis with good agreement with experimental results. However, at large angles of attack, it was confirmed that the CFD simulation results tend to slightly underpredict the power and overpredict the thrust [37][38][39][40][41]. An incompressible analysis was employed, as most of the wind turbine operation is in the subsonic region, with a constant local air density and viscosity. A coupled algorithm for pressure, which considerably increases the convergence rate, was considered in solving the incompressible Reynold-averaged Navier-Stokes (RANS) model [33,42]. Residuals and net mass imbalance were engaged to limit the running time of the solver and ensure convergence. Residuals of the specific dissipation ω, the kinetic energy turbulence k, velocities, and continuity were limited to less than 10 −4 during the solution process. In addition to these, the net imbalance was also checked to be below 0.1% [33,37,43], and an integral static pressure surface monitor to remain constant for a considerable amount of iterations after convergence.
Furthermore, ASHES (simis.io) was employed as an off-design performance evaluation tool. It is a time-domain wind-turbine aerodynamics module, which uses AeroDyne of FAST with a rich user interface. Parameters like a chord, twist angle, pitch axis, airfoil type, and shapefiles of the airfoils at their corresponding radial positions from Matlabwere exported to ASHES as input in the table format. It internally implements xfoil to determine the lift and drag constants for the shapefiles to get the necessary parameters in the evaluation process.

Wind Distribution Characteristics and Design Parameters
The wind characteristics of the four sites in the Geba catchment area as bins of power density distribution and wind speed distribution are shown in Figure 2a-d. The figures show also graphs of the Weibull probability wind distribution and power density distribution. Figure 2e shows weighted average distribution over the catchment based on the data from the four sites. Figure 2f shows the corresponding cumulative distribution graphs of the sites the weighted average over the catchment. The results in Figure 2e,f show that the weighted average is close to the wind conditions of the Mayderhu site. Based on the wind distribution and equations in Sections 2.1 and 2.2, the parameters for the design of the representative blade were amean wind speed of 6.66 m/s, a standard deviation of 2.94 m/s, a shape factor , of 2.43, and scale factor , of 7.52 m/s. As can be seen from the plots in Figure  2e, the maximum power density wasat a wind speed of around 9.5 m/s, while the most frequent wind wasat 6 m/s. Hence, wind speed of 10 m/s wasselected as a design wind speed. The range between cut-in wind speed of 2 m/s and cut-out wind speed of 18 m/s was considered as the usable operating range based on the power density distribution. A summary of the design parameters considered isshown in Table 2.  Based on the wind distribution and equations in Sections 2.1 and 2.2, the parameters for the design of the representative blade were amean wind speed of 6.66 m/s, a standard deviation of 2.94 m/s, a shape factor k, of 2.43, and scale factor c, of 7.52 m/s. As can be seen from the plots in Figure 2e, the maximum power density was at a wind speed of around 9.5 m/s, while the most frequent wind was at 6 m/s. Hence, wind speed of 10 m/s was selected as a design wind speed. The range between cut-in wind speed of 2 m/s and cut-out wind speed of 18 m/s was considered as the usable operating range based on the power density distribution. A summary of the design parameters considered is shown in Table 2. The output of the blade design analyzed using Matlab code shown in Figure 3a in a 3D plot. Figure 3b shows the same rotor blade after exporting to Inventor as a step file and processed to get the blade surface for further use in the CFD analysis. The output of the blade design analyzed using Matlab code shown in Figure 3a in a 3D plot. Figure 3b shows the same rotor blade after exporting to Inventor as a step file and processed to get the blade surface for further use in the CFD analysis.

CFD Analysis and Validation
The results of the CFD analysis of the blade sections are shown in Figure 4 at different radial positions from root to the tip.

CFD Analysis and Validation
The results of the CFD analysis of the blade sections are shown in Figure 4 at different radial positions from root to the tip. Figure 4 presents velocity streamlines in the left and pressure contours in the right for each section. The flow was fully attached at most of the radial positions except at around the root section. The root part started at r/R = 0.042 and the hub region assumed was up to around 0.115. Along the trailing edge, flow separated at around these radial positions. It was due to the modification of the twist angle out of the optimized values to create a smooth connection to the hub plate. The figure also clearly shows that the angle of attack was close to zero at the tip region and increases towards the root region being maximum at around r/R = 0.115 as a result of the optimum twist angle employed. An almost perpendicular angle of attack was used at the root because of positioning the blade root for proper connection with the hub plate.
Furthermore, the surface pressure contour of the turbine blade shows that pressure distribution was found highest on the pressure side from the tip to around two-thirds of the span along the leading edge. Besides, a lowest pressure was observed on the suction side along the leading edge at similar radial positions. Hence, these differences in surface pressure between the pressure and suction surfaces of the blade created the driving torque with the highest values at around the tip region and decreased while moving towards the root part resulting in smaller pressure differences. It also shows that flow velocity was high on the low-pressure sides and low on the high-pressure sides.
Energies 2020, 13, x FOR PEER REVIEW 10 of 21 The output of the blade design analyzed using Matlab code shown in Figure 3a in a 3D plot. Figure 3b shows the same rotor blade after exporting to Inventor as a step file and processed to get the blade surface for further use in the CFD analysis.

CFD Analysis and Validation
The results of the CFD analysis of the blade sections are shown in Figure 4 at different radial positions from root to the tip.   Furthermore, the surface pressure contour of the turbine blade shows that pressure distribution was found highest on the pressure side from the tip to around two-thirds of the span along the leading edge. Besides, a lowest pressure was observed on the suction side along the leading edge at similar radial positions. Hence, these differences in surface pressure between the pressure and suction surfaces of the blade created the driving torque with the highest values at around the tip region and decreased while moving towards the root part resulting in smaller pressure differences. It also shows that flow velocity washigh on the low-pressure sides and low on the high-pressure sides. Figure 5 illustrates the limiting streamlines on the pressure and suction sides of the blade at the design conditions. On the pressure side, almost all flows weresmooth, and on the suction side, most of the streamlines remained attached to the blade surface except at around the root part. Separation started at the root end and reattached back around r/R = 0.115, which wasthe point where the optimized twist modification starts creating a connection to the hub plate. Then, the flow will be smooth again along the reattachment line near the trailing edge, as shown in the suction side of the blade. It indicates that the advantage of using twisted small wind turbines to ensure smooth flow.   Figure 5 illustrates the limiting streamlines on the pressure and suction sides of the blade at the design conditions. On the pressure side, almost all flows were smooth, and on the suction side, most of the streamlines remained attached to the blade surface except at around the root part. Separation started at the root end and reattached back around r/R = 0.115, which was the point where the optimized twist modification starts creating a connection to the hub plate. Then, the flow will be smooth again along the reattachment line near the trailing edge, as shown in the suction side of the blade. It indicates that the advantage of using twisted small wind turbines to ensure smooth flow.
As a validation process, Figure 6 shows the power coefficient for the variable speed characteristics of the wind turbine applied between the cut-in and cut-out wind speeds of 2-18 m/s, respectively. The corresponding rotational speeds werearound 6.5 rad/s and 32.5 rad/s, respectively. The figure presents the C P plotof CFD and BEM results at the rated power of 5 kW and rated wind speed of 10 m/s. Accordingly, the CFD results werein good agreement to the BEM results at most of the wind speeds with a slight decrease at high wind speeds. This small deviation of the power coefficient in CFD was due to an increased angle of attack at higher wind speeds, and is an expected kind of variation as discussed in Section 2.7, which was verified by different researchers using experiments. Thus a maximum power coefficient of around 51.8% was achieved at the design/rated wind speed as shown in the figure. As a validation process, Figure 6 shows the power coefficient for the variable speed characteristics of the wind turbine applied between the cut-in and cut-out wind speeds of 2-18 m/s, respectively. The corresponding rotational speeds werearound 6.5 rad/s and 32.5 rad/s, respectively. The figure presents the CP plotof CFD and BEM results at the rated power of 5 kW and rated wind speed of 10 m/s. Accordingly, the CFD results werein good agreement to the BEM results at most of the wind speeds with a slight decrease at high wind speeds. This small deviation of the power coefficient in CFDwasdue to an increased angle of attack at higher wind speeds, and is an expected kind of variation as discussed in Section 2.7, which was verified by different researchers using experiments. Thus a maximum power coefficient of around 51.8% was achieved at the design/rated wind speed as shown in the figure.
Pressure side Suction side  As a validation process, Figure 6 shows the power coefficient for the variable speed characteristics of the wind turbine applied between the cut-in and cut-out wind speeds of 2-18 m/s, respectively. The corresponding rotational speeds werearound 6.5 rad/s and 32.5 rad/s, respectively. The figure presents the CP plotof CFD and BEM results at the rated power of 5 kW and rated wind speed of 10 m/s. Accordingly, the CFD results werein good agreement to the BEM results at most of the wind speeds with a slight decrease at high wind speeds. This small deviation of the power coefficient in CFDwasdue to an increased angle of attack at higher wind speeds, and is an expected kind of variation as discussed in Section 2.7, which was verified by different researchers using experiments. Thus a maximum power coefficient of around 51.8% was achieved at the design/rated wind speed as shown in the figure.

Performance Evaluation
It is essential to look into different indicator plots to assess the main features of the performance of the designed blade for a range of rotor speeds at different wind speed values. Therefore, Figures 7 and 8 show computed power and torque versus rotor speeds and power coefficient for varying wind speeds. power curve following the maximum power coefficient Cp-max points until it reached a rated value at the rated wind speed of 10 m/s and the corresponding rotor speed of 307.2 rpm. Beyond the rated wind speed, the power remained at a constant value, which wasthe rated power, moving horizontally to the left, which lowered the rotor speed considering the variable speed operated wind turbine throughout its range. It also shows the maximum power-tracking graph applicable for variable speed generators.     Figure 7a represents power versus rotational speed for a range of wind speeds. It includes a power curve following the maximum power coefficient C p-max points until it reached a rated value at the rated wind speed of 10 m/s and the corresponding rotor speed of 307.2 rpm. Beyond the rated wind speed, the power remained at a constant value, which wasthe rated power, moving horizontally to the left, which lowered the rotor speed considering the variable speed operated wind turbine throughout its range. It also shows the maximum power-tracking graph applicable for variable speed generators.
Similarly, Figure 7b shows torque versus rotor speed for the range of wind speeds and operating points of power. Likewise, the operation coincided with the C P-max locus until it grasped the rated wind speed. All along this curve, rotor speed roseto the rated value. Above the rated wind speed, the operating point moved together with the rated power locus to the stall region. Based on these, it waspossible to generate the power, power coefficient, and rotor speed curves of the turbine to evaluate the performance in terms of AEP and CF of the designed blade.
As a result, Figure 8 shows the variable-speed characteristics of the rotor using the power, power coefficient, and rotor speed against the assumed usable range of wind speeds. In this figure, two different regimes are shown, a maximum-C p and power regulated regions. The maximum-C p section spanned from the cut-in wind speed to the rated wind speed of 10 m/s, with maximum turbine C p at each wind speed. Hence, the C p was 41.4% at the wind speed of 2 m/s and reached around 51.8% at the rated wind speed. Likewise, the power controlled section ranged from the rated to the cutout wind speed with a constant rated generator power of 5 kW realized by decreasing the rotational speed to a regulated constant value. In this region, the C p decreased in a parabolic curve to about 8.9% at a wind speed of 18 m/s. The performance parameters weredetermined by implementing the equations discussed in Section 2.2 (Equations (5), (6), and (10)) to the power curve of the designed blade in Figure 8 and the weighted average wind distribution Figure 2e. The calculated parameters werefound to be the average machine power, P w of 1956.1 W, annual energy production, AEP of 17.14 MWh, and power factor, CF of 38.17% considering mechanical efficiency of the blade only, which is C p . As direct-drive generators wereused in most cases for small wind turbines, assuming a generator efficiency of 97%, for instance, the AEP and CF became 16.63 MWh and 37.03%, respectively.
Consequently, the selected design conditions and processes had led to increased performance of the rotor. Such high performance of small wind rotors, in terms of power coefficient and capacity factor, could be achieved by implementing maximum power tracking control mechanisms together with variable speed generators applying variable speed fixed-pitch operation with assisted stall regulation control strategy.

Selection of Rated Power for Low Cost of Energy
Now consider the use of the designed blade for the four sites at the designed rated power of 5 kW, and rated wind speed of 10 m/s. The AEP and CF for the sites weredetermined based on the power curve of Figure 8 and the respective wind characteristics of the sites in Figure 2a-d. The results of calculation show AEP of around 8. 43, 9.16, 22.55, and 17.01 in MWh with corresponding CF of 18.8%, 20.4%, 50.2%, and 37.9% respectively, for Mekelle, Dera, Hagereselam, and Mayderhu. Then, examining the relationship between AEP, CF, and CoE against rated power for each wind-site, parameters for the optimal use of the blade were assessed. Thus, plots of AEP and CF against rated-power werefound as shown in Figure 9. Generally, the Hagereselam site performed the highest in both parameters; the Mayderhu site delivered moderate performance while Dera and Mekelle delivered low performance. As a general pattern, the AEP increased with increased rated power while the CF decreased. It is crucial to select appropriate rated-power that can give an advantage to the sites in minimizing the CoE. The optimal power for the respective sites wasdetermined based on the AEP for various rated powers and the relative-CoEconsidering the suggested cost scaling of the power electronics and blade discussed in Section 2.4. As a result, Figure 10 shows relative-CoE against rated power. Using the turbine at the rated power of 5 kW and rated wind speed of 10 m/s at each site, the relative-CoE increase in comparison to the nominal site value was3%, 7%, and 1% for Mekelle, Dera, and Hagereselam sites, respectively, and with no change for the Mayderhu site. Therefore, to find the rated power that would give minimum CoE, relative-CoE very close to one was identified for the sites. As a general pattern, the AEP increased with increased rated power while the CF decreased. It is crucial to select appropriate rated-power that can give an advantage to the sites in minimizing the CoE. The optimal power for the respective sites was determined based on the AEP for various rated powers and the relative-CoE considering the suggested cost scaling of the power electronics and blade discussed in Section 2.4. As a result, Figure 10 shows relative-CoE against rated power. Using the turbine at the rated power of 5 kW and rated wind speed of 10 m/s at each site, the relative-CoE increase in comparison to the nominal site value was 3%, 7%, and 1% for Mekelle, Dera, and Hagereselam sites, respectively, and with no change for the Mayderhu site. Therefore, to find the rated power that would give minimum CoE, relative-CoE very close to one was identified for the sites. Accordingly, using Figure 10, the recommended range of rated power with corresponding rated wind speeds, AEP, and CF of the four sites are shown in Table 3. Still the Hagereselam site performed at the top in terms of all the performance parameters, and the Mekelle site wasthe lowest in terms of annual energy production and capacity factor. The Mayderhu site remained at around the design values. On top of the advantage of selecting the rated power that can give low CoE, the use of a single blade design concept had an advantage in minimizing the cost of design and production. If specific blades were considered, for the four sites, four different designs and molds would be leading to at least four times the single blade design and mold cost. Hence, the application of the single blade concept together with selection of appropriate rated power for the respective site minimizes the cost of energy. In future work, the structural analysis based on the wind conditions considered should be studied, especially for the low wind sites, focusing on determining efficient structure to minimize material and production costs, besides the safety requirements.

Conclusions
Catchment based aerodynamically optimized small wind turbine was designed using representative wind conditions. The proposed weighting average was based on the wind power density of the wind sites within the catchment area. The blade design was applying the single blade Accordingly, using Figure 10, the recommended range of rated power with corresponding rated wind speeds, AEP, and CF of the four sites are shown in Table 3. Still the Hagereselam site performed at the top in terms of all the performance parameters, and the Mekelle site wasthe lowest in terms of annual energy production and capacity factor. The Mayderhu site remained at around the design values. On top of the advantage of selecting the rated power that can give low CoE, the use of a single blade design concept had an advantage in minimizing the cost of design and production. If specific blades were considered, for the four sites, four different designs and molds would be leading to at least four times the single blade design and mold cost. Hence, the application of the single blade concept together with selection of appropriate rated power for the respective site minimizes the cost of energy. In future work, the structural analysis based on the wind conditions considered should be studied, especially for the low wind sites, focusing on determining efficient structure to minimize material and production costs, besides the safety requirements.

Conclusions
Catchment based aerodynamically optimized small wind turbine was designed using representative wind conditions. The proposed weighting average was based on the wind power density of the wind sites within the catchment area. The blade design was applying the single blade concept and low cost of energy. In designing, BEM was implemented using Matlab code to optimize the local chords, twist angles, radius, and rotational speed of the turbine and exported to Ansys for CFD analysis. Ansys-Fluent was applied to validate the results of BEM, and good agreement found for the power coefficient for a different combination of wind speeds.
Accordingly, the maximum power coefficient achieved was around 51.8% at a design wind speed of 10 m/s. Considering the single blade concept for all the sites within the catchment area, the performance of the blade at a rated power of 5 kW resulted in capacity factors ranging from 18.8% to 50.2%. Therefore, analyzing the relationship between rated power against annual energy production, capacity factor, and relative cost of energy, a range of appropriate rated power that would give a low cost of energy at higher performances recommended for each site. It can be concluded that for Mekelle and Dera sites, which are low wind sites, rated-powers of 3.5 kW and 3 kW are appropriate at rated wind speeds of 8.8 m/s and 8.4 m/s, performing at capacity factors around 25.1% and 31.8% and corresponding annual energy production around 7.3 and 8.4 in MWh, respectively. While optimum rated-power of 5 kW would be proper for the Mayderhu site at 10 m/s rated wind speed, executing with a capacity factor around 38.6% and annual energy production around 16.9 MWh. On the other hand, as the Hagereselam site is a relatively high wind site, rated power of 6 kW would be proper at 10.5 m/s rated wind speed that would achieve a capacity factor of 45.9% and annual energy production of 24.1 MWh.
Even though the performance of the blade at the Mekelle site was low relative to other wind sites in the catchment, it was moderate to high compared to the most commonly reported wind turbines. On the other hand, the performance at the Hagereselam site was the highest in the catchment area considered and high compared to the existing wind rotors. As the Mayderhu wind condition is very close to the representative weighted average of all sites used for designing the blade, the performance was found very close to the results at design conditions. In conclusion, it is possible to achieve high performance at a low cost of energy using a single blade concept at appropriate rated power, implementing proper design conditions and procedures for different sites with different wind conditions.