Multi ‐ Objective Optimization and Optimal Airfoil Blade Selection for a Small Horizontal ‐ Axis Wind Turbine (HAWT) for Application in Regions with Various Wind Potential

: The type of airfoil with small wind turbine blades should be selected based on the wind potential of the area in which the turbine is used. In this study, 10 low Reynolds number airfoils, namely, BW ‐ 3, E387, FX 63 ‐ 137, S822, S834, SD7062, SG6040, SG6043, SG6051, and USNPS4, were selected and their performance was evaluated in a 1 kW wind turbine in terms of the power coeffi ‐ cient and also the startup time, by performing a multi ‐ objective optimization study. The blade ele ‐ ment momentum technique was utilized to perform the calculations of the power coefficient and startup time and the differential evolution algorithm was employed to carry out the optimization. The results reveal that the type of airfoil used in the turbine blade, aside from the aerodynamic performance, completely affects the turbine startup performance. The SG6043 airfoil has the highest power coefficient and the BW ‐ 3 airfoil presents the shortest startup time. The high lift ‐ to ‐ drag ratio of the SG6043 airfoil and the low inertia of the turbine blades fitted with the BW ‐ 3 airfoil make them suitable for operation in windy regions and areas with low wind speeds, respectively.


Introduction
In today's developed world, energy is a cardinal factor for the sustainable development and economic growth of a country [1], and, in recent years, continuous access to inexpensive and reliable energy has become an essential right for humanity [2].However, the rise in energy demand and dependence on fossil fuels has led to energy crises, declining resources, and environmental pollution [3].Although optimizing energy systems [4,5] and employing energy storage devices [6] are considered by scholars, the use of renewable energy sources [7] is one of the most effective solutions to cope with these difficulties.Among renewable energy sources, wind energy is the most common and sustainable type of energy and, currently, the application of this free energy source is being considered more than ever [8].Harnessing wind energy and converting it to electrical power is carried out with the aid of wind turbines.
The shape of a wind turbine blade has a substantial role in its performance and should be determined according to the design goal or goals.Compared to the experimental and simulation studies conducted on wind turbine blades, optimization algorithms present more feasible and cost-effective solutions for acquiring blade sections with suitable performance in different scenarios.The enhancements obtained with these algorithms result from their searching strategy.Contrary to the other techniques, these algorithms do not depend on searching and testing the entire possibilities.These algorithms make an effort to acquire the optimal solutions by trying the model performance at each trial by the means of an objective function.
Recently, nature-inspired metaheuristic strategies for solving different optimization problems, especially in computer science and engineering, have rapidly increased.Soni et al. [9] compared the performance of 12 nature-inspired optimization algorithms, including Bat, Lion, Particle swarm, Water wave, Elephant herding, Optics inspired, Cuckoo search, Flower, Genetic, Differential evolution (DE), Harmony, and Simulated annealing algorithms.They discussed the speed, accuracy, performance, convergence, efficiency, and complexity of these algorithms.Their findings revealed that compared to the other selected algorithms, the differential evolution algorithm presents fast speed and performance, as well as the best accuracy.
Tušar and Filipič [10] compared the performance of genetic algorithms and differential evolution in solving multi-objective optimization problems.Their results showed that the DE algorithm provides better and more stable solutions than genetic algorithms.Similar results were observed in the research of Lilla et al. [11].
According to the International Electrotechnical Commission (IEC) 61400-2 [12], a small horizontal-axis wind turbine has a swept area of less than 200 m 2 , which corresponds to a blade radius and a power output of less than 8 m and around 50 kW, respectively.In a further subdivision, small wind turbines (SWTs) are divided into three groups, namely micro, mid, and mini, with output power less than 1, 5, and 50 kW, respectively [13].Although large wind turbines have gained significant attention from researchers [14], SWTs have not received the research interest they deserve.These turbines can be employed to supply electricity for domestic and agricultural purposes, and their application is economically viable [15].Unlike large wind turbines which are used in windy areas, these turbines can also be installed in places where there is no potential for continuous high-speed wind [16] .
An important difference between large and SWTs is the lack of pitch controllers in SWTs for reducing costs [17].Thus, in these turbines, optimally positioning the blade sections in the wind direction is not possible.This point is significantly important when the SWTs are beginning to rotate from a stationary state, since throughout the startup process, the high angles of attack along the blade, result in the reduction in the lift force and aerodynamic torque, and, consequently, the blade rotates slowly and power generation is delayed [18].
The time required for the turbine to accelerate from a stationary state to a tip speed ratio of around one is called the startup time (Ts) [19].In addition to achieving the maximum power coefficient (Cp), which is the main goal in designing wind turbine blades, reducing this time has also been one of the interesting topics for researchers in recent years.In this regard, the optimization of SWT blades has been performed to maximize the Cp and minimize Ts [19].The results showed that most of the turbine power is obtained from the blade tip part while the root section can facilitate the startup process.By properly distributing the chord length and twist angle in these two parts of the blade, a significant improvement was observed in the Cp and Ts of these turbines.Pourrajabian et al. [20] compared the performance of solid and hollow blades in terms of Cp and Ts in an SWT.It was observed that compared to the solid blades, the hollow blades exhibit a higher Cp and lower Ts.Rahgozar et al. [21] studied the influence of linear and non-linear distributions of the twist angle and chord length on the Cp and Ts of a small horizontal-axis wind turbine.The results showed that, while the use of non-linear distribution can increase the Cp in windy areas, with linear distribution the blades have a better Ts in areas with low wind speeds.Minimizing the aerodynamic noise [22], mass, and cost [23] are other goals of designing and optimizing SWT blades that have been studied by researchers.
The shape of the airfoil determines the drag and lift coefficients (Cd and Cl) and plays a crucial role in the aerodynamic torque generated by the wind turbine blade and its output power.In this regard, researchers have always attempted to use airfoils with high aerodynamic efficiency in blade design [24].
In general, in large turbines, several airfoils are used in different parts of the blades, each of which has its own task.Thin airfoils with better aerodynamic performance are used at the tip of the blade to provide more power, while thick airfoils are employed at the root of large blades to resist structural stresses [25].Unlike large blades, to reduce construction costs as much as possible, the airfoil type of the SWT blades does not change [26].The operating Reynolds number (Re) range of SWTs is less than 500,000 [27].Therefore, there are not a variety of airfoils to use in these turbines.In this range of Re, the flow on the upper surface of the airfoil is often laminar, so the formation of a laminar separation bubble is possible.To solve this problem and to reduce the adverse pressure gradient on the upper surface of the airfoil, airfoils thinner than the traditional ones are recommended for operation in SWT blades [28].Giguere and Selig [29] designed a group of SG series airfoils (SG6040-SG6043) specifically for use in SWT blades.Studies have shown that increasing the leading edge nose radius, as well as cusping of the trailing edge, improves the aerodynamic performance of thin airfoils, and, in this regard, some airfoils have been proposed specifically for use in SWTs [30].
Given the importance of the airfoil type used in the blade for the aerodynamic performance of a wind turbine, its proper selection, in terms of the startup process, needs further study.The previous works that were performed in this field are mostly focused on the distribution of chord length and twist angle along the blade, and the effect of the airfoil type used in the blade has not been investigated from the startup viewpoint, which is discussed in the present study.For this purpose, the performance of 10 different airfoils that are specifically designed for SWT blades is examined in terms of the Cp and Ts.
The remainder of this study is organized as follows.The selected airfoils are introduced in Section 2. Section 3 presents the numerical techniques that are employed, as well as their validations.The results of the study which include the performance evaluation of the airfoils in windy areas and areas with low wind speed are presented in Section 4. Finally, the conclusions are discussed in Section 5.

The Selected Airfoils
Despite the variety of airfoils available for use in large wind turbine blades, there are not many airfoils designed for SWTs.By examining previous research studies and also considering the operating Re of the turbine blade investigated in the present study, 10 airfoils that are used in the industry of SWTs, including, BW-3 [31], E387 [30], FX 63-137 [30], S822 [31], S834 [30], SD7062 [31], SG6040 [31], SG6043 [31], SG6051 [32], and USNPS4 [31] were selected for examination.Their aerodynamic coefficients, which are presented in Figure A1 of Appendix A, were tabulated and coded for various Re and angle of attack values, and they were employed in the numerical method.The geometry of the selected airfoils has been shown in Figure 1.In this figure, x/c indicates the distance along the chord line, while y/c shows the airfoil thickness normalized with the chord length.

Multi-Objective Optimization
A multi-objective optimization was used to determine how to optimally distribute the twist angle and chord length along the blade for each of the selected airfoils to maximize the Cp and minimize Ts.For this purpose, in the present study, the differential evolution (DE) algorithm which was first introduced by Storn et al. [33] in 1997 was employed.DE is a population-based method that, among the recently presented metaheuristic search algorithms (MSAs), has become one of the most popular ones for solving optimization problems [34].Presenting better global convergence and robustness, simpler and more straightforward implementation, as well as having good performance in solving multi-objective optimization problems, are the main factors that make this algorithm superior to other MSAs [35].In this algorithm, the initial population (blades) is generated completely randomly, based on the upper and lower limits of the design variables.Then, using mutation and crossover operators, children (new blades) are generated.Subsequently, with the help of the selection operator, each child is compared with the parent which generated it.If the child has a better objective evaluation, it replaces the parent, and the next generation is generated.This process is repeated in the same way until the termination criteria are met.In the present study, the termination criterion was the number of generations.Although raising this number (500 generations in this study) increases the x/c computation time, it guarantees to reach the global solution .It is worth mentioning that to ensure that the obtained answer is a global one, the values of the objective function were monitored.Since its values remained unchanged for 100 consecutive generations, it was made sure that the answer is global and the solution can be terminated.Employing the weighted coefficient method [36], to maximize Cp and minimize Ts, the following objective function was considered [20]: using which, the objective function tries to simultaneously maximize the Cp (the first goal) and minimize the Ts (the second goal).The weighting factor, n (0  1) determines the influence of these two goals on the objective function.Values of max  and min  represent the blades with the highest Cp and also the fastest Ts, respectively, which are updated in each generation.In other words, in each generation, the blades are compared and those with the best Cp and Ts values become candidates and form max  and min  for that generation.When the solution proceeds to the next generation, a new set of max  and min  are specified.It is noteworthy that these values remain constant after a sufficiently large number of generations.This way, it is made sure that the DE algorithm approaches the optimal solution.

Calculating Design Goals
To calculate the  , the first objective of the design, the blade element momentum (BEM) technique was employed.This theory, which is the most widely used engineering model for simulating wind turbine blades and propellers [37], was proposed in 1935 by Glauret [38].BEM is the combination of the blade element method and momentum theory in which the blade is divided into several independent elements, each element having a specific twist angle and chord length.Figure 2 shows the velocities on a blade element, where  is the wind velocity, and   is the total velocity, which is obtained according to the value of the velocity at the blade section ( 1 ) and also the angular velocity of the turbine ().The angles shown in this figure are the angle of attack (α), twist angle (θ), and flow angle ().Based on the figure,  can be calculated from the following equation: where a and  represent the axial and rotational induction factors, respectively, whose values are calculated in an iterative process [19].Figure 3 shows the axial ( ) and tangential ( ) forces on a blade element.Using the continuity equation along with the momentum equations for each element, the local values of axial force and aerodynamic torque  are calculated from: In the above equations,  is the air density,  is the radial distance from the rotor hub,  is the chord length, and   and   are the lift and drag coefficients for the blade element, respectively.By summing the torque generated by each element, the total aerodynamic torque Q is calculated.Finally, the   of the turbine is calculated from [39]: in which, S is the swept area of the blades.
To take the tip losses into account, Prandtl's model () is employed.The purpose of this work is to determine the effect of blade tip losses along the blade.This parameter, which is a function of the number of blades (), flow angle (ϕ), and also the radial coordinate along the blade (), is defined as follows [40]: By using the BEM theory and using analytical correlations for aerodynamic coefficients at high angles of attack, it is possible to determine the   from the stationary state.Wood [41] showed that at high α values, from an aerodynamic viewpoint, airfoils behave like a two-dimensional flat plate, and to calculate their lift and drag coefficients, the following equations can be used with acceptable accuracy and regardless of the airfoil type: By implementing Equations ( 8) and ( 9) in Equation ( 4) and normalizing c and r with the blade tip radius (), and also all speeds (,   , and  1 ) with the wind velocity () the startup torque   is calculated as follows: where,  is the hub radius, and    ⁄ represents the local tip speed ratio.From the assumption that no power is generated throughout the startup time and the aerodynamic torque is merely used to rotate the blades from the stationary state, the variation in the tip speed ratio can be computed using the following ordinary differential equation [20]: in which,    ⁄ is the tip speed ratio,  shows the generator resistive torque, and  indicates the total rotational inertia, which is the sum of the blade inertia ( ) and the generator inertia ( ).By considering a specific value for the tip speed ratio (λ), this equation is solved and the time to reach  from the stationary state ( 0) to the specific value (which is considered in advance) is calculated.The value of  1 is considered the ending point of the startup process and the above equation is solved using this method.It is necessary to explain that the parameters such as R,  , and also U, are constant values in the present study.
By normalizing c and r by R, the moment of inertia of the blade is determined from [19]: (12) in which,  is the density of the blade, and A is the surface area of the airfoil used in the blade, assuming that the unit is chord length.
It should be noted that the startup performance of a 5 kW wind turbine measured by Sessrego and Wood [42] has confirmed the accuracy of Equation (11).This reference indicates that the startup stage involves two important steps: (1) the long idling step, in which the blade has a slow rotational acceleration and continues until the angles of attack along the blade fall to less than 20°, and (2) the rapid acceleration step, in which the blade reaches an operational angular velocity.Therefore, to minimize the startup time, the initial phase should be minimized.Integrating Equation (11) up to λ = 1 properly captures this step.To calculate Equation ( 11), the Adams-Moulton method was used and the trapezoidal integration method was employed to calculate the blade inertia of Equation ( 12).

Adjusting the Input Parameters
The flowchart of the used optimization algorithm is illustrated in Figure 4.A code that was developed in MATLAB software was employed for the optimization process and calculation of the design objectives.Based on the values suggested in the literature [19,43], the blade was divided into 15 elements.Table 1 shows the minimum and maximum values of the design variables according to the limitations and construction problems for the twist angle/chord length distribution.Additionally, Table 2 lists the input parameters considered in the DE algorithm.
Table 1.The specified limits for the design parameters [20].

Parameter Minimum Maximum
Twist, (°) −5 25 Chord, c/R 0.01 0.2 The selected small horizontal-axis wind turbine was introduced and numerically evaluated by Pourrajabian et al. [44].The technical specifications of this turbine, which include the output power (P), number of blades (N), blade radius (R), hub radius ( ), angular velocity (), rated tip speed ratio ( ), rated wind speed ( ), generator inertia (( ), generator resistive torque ( ), blade density ( ), startup wind speed ( ), and the type of airfoil used in the blade are given in Table 3.

Validation of the BEM Code
To evaluate the accuracy of the developed BEM code, the results of the present study have been validated with the experimental results from the wind tunnel test of a small wind turbine with a radius of 1.5 m.This two-bladed turbine was examined by Anderson et al. [45], and then further investigated by Wood [19] and Spera [37].The geometry of this turbine, which uses NACA4412 airfoil as the blade cross-section, is shown in Figure 5.The aerodynamic coefficients of the NACA4412 airfoil, including the lift and drag coefficients at various angles of attack and different Reynolds numbers, are also depicted in Figure 6.As mentioned in Section 2, all aerodynamic coefficients of the airfoils used in the present study are coded, and, based on the Re number and the angle of attack of each blade element, the corresponding values of Cl and Cd are computed and called during the BEM calculations and optimization steps.
Figure 7 shows the changes of Cp at different λ values.It can be observed that the results obtained from the developed numerical code are well-compatible with the experimental results with a maximum error of 6.47%.It is noteworthy that the mean squared error (MSE) of the obtained numerical values was found to be 0.000228562.More detailed data for this comparison are given in Table A1 in Appendix B.

Validation of the Optimization Code
The ideal equations proposed by Burton et al. [46] were employed to evaluate the performance of the optimization code.By ignoring the tip losses and drag, the following correlations determine the twist angle and the chord length along the blade to achieve the maximum power coefficient: Considering the BW-3 airfoil for the 1 kW turbine blade of the present study, Figure 8 shows the distribution of the twist angle (θ) and chord length (c) obtained from the optimization algorithm (n = 1) and the ideal Equations of (13 ) and (14 ).As can be observed, there is a good agreement between the distributions along the blade.Therefore, it is ensured that the final answer of the optimization code is optimal.Of course, there are some discrepancies in the distributions considered by the optimization algorithm in the root section, the reason for which is revealed in detail in the next section.It should be noted that, to use the ideal Equations ( 13) and ( 14), it is necessary to specify the α and Cl values of the airfoil at the maximum lift-to-drag ratio, which, according to Figure 9, are 5.39° and 1.034, respectively, for the BW-3 airfoil at the Re of 200,000.

Discussion of Results
Similar to the specified ranges for the design variables (mentioned in Table 1) and also the value of Tables 2 and 3, n is an input parameter for running the optimization.In multi-objective optimization problems, there is no single answer.In this study, the answers are determined by changing the value of n in the objective function (Equation ( 1)), which leads to the Pareto front.Here, Pareto front members are the blades whose at least one component of the objective function has a larger Cp or 1/Ts value than the other blades.Generally, larger values are considered for n, because the output power is a more important goal than the startup time [19].Imposing low weighting factors (n) causes poor aerodynamic performance of the blade and makes it practically unusable.Selecting the final blade from the Pareto front depends on the wind potential.Although n = 1 is used in windy areas to achieve the maximum power coefficient, n = 0.8 is recommended in areas with low wind speeds where the startup performance of turbines is more important [20,21].
By considering the input parameters mentioned in Section 3.3, as well as the explanations provided above, the optimization process was used to evaluate the performance of the selected airfoils, the results, which include Cp and Ts, are summarized in Table 4.It is necessary to explain that the details of the optimization results were also obtained for n = 0.6, and the tabulated results can be found and compared in Table A1 of Appendix C.

Investigating the Performance of Airfoils in Windy Areas
In areas with high wind speeds (n = 1) where Ts is removed from the objective function, the SG6043 airfoil yields the highest Cp while the S822 and SG6040 airfoils present the lowest Cp (see Table 4).The lift-to-drag ratio of airfoils is an important parameter in the efficiency of wind turbine blades.Figure 10 shows the maximum lift-to-drag ratio of the selected airfoils over a wide range of Re numbers.
It can be observed that, compared to the other airfoils, this value is higher for the SG6043 airfoil.It should be noted that the type of airfoil completely influences the distribution of twist and chord along the blade.Indeed, the optimization algorithm considers the best geometry for the blade based on the aerodynamic coefficients and geometric characteristics of each airfoil.Figure 11 shows the distribution of twist and chord of the selected airfoils for n = 1.By focusing on this figure, it can be seen that the values of the chord length for the blade that was fitted with the FX 63-137 airfoil are smaller than the values considered for the other selected airfoils.Due to this, despite the high lift-to-drag ratio of this airfoil (Figure 10), the Cp is not so high.Because the smaller chord length leads to a lower Re and, hence, reduces the lift-to-drag ratio.The opposite is true for the USNPS4 airfoil.According to Figure 10, it can be observed that this airfoil does not have a high lift-to-drag ratio, but the larger chord length values of the blade fitted with this airfoil increase the Re and, consequently, raise the Cp.
It should be noted that for all selected airfoils, the twist angles and chord lengths are not smooth at the root elements (as seen in Figure 8).This is due to the low contribution of these elements in the total aerodynamic torque and, hence, the Cp of the turbine, which makes the optimization algorithm focus on the evolution of the design variables in the middle and tip elements and creates a smoother distribution in these elements.To take a closer look at this point, Figure 12 illustrates the torque generation distribution along the blade throughout the power generation for the SG6043 and S822 airfoils.As can be seen, regardless of the airfoil type, a significant portion of the aerodynamic torque is generated by the middle and tip elements of the blade, and the share of the root elements is considerably smaller.It should be noted that the slight decrease in the aerodynamic torque at the tip of the blade is because of considering the blade tip losses in the BEM calculations.

Investigating the Performance of Airfoils in Areas with Low Wind Speed
In areas with low wind speeds (n = 0.8), small wind turbines need to react quickly to the wind and start generating power in the shortest possible time.In this regard, Worasinchai et al. [18] have shown that reducing startup time increases energy capture.
According to Table 4, the startup performance of the BW-3 airfoil is better than the other selected airfoils.Thus, the turbine that uses this airfoil has the shortest Ts.The reason is the low inertia of the blades fitted with this airfoil (see Figure 13).

Airfoil
The S822, S834, and SG6040 airfoils have the highest Ts and their application in areas with low wind speeds is not recommended at all.The common aspect of these three airfoils is their high surface area compared to the other selected airfoils, which is an effective factor in the blade moment of inertia, aside from the chord length (Figure 14).But it is also important to note that a thinner airfoil does not always present a better performance at low wind speeds.For example, the E387 airfoil is thinner than the SG6043 airfoil and has a smaller surface area (see Figure 14), but it has a longer Ts (Table 4 for n=0.8).The reason is the higher moment of inertia and the lower startup torque (the startup torque is examined at the end of the current section) of the blades fitted with this airfoil compared to the blades fitted with the SG6043 airfoil.
Figure 15 shows the optimal distribution of twist and chord for the selected airfoils for n = 0.8.By comparing this figure with Figure 11, it can be observed that regardless of the airfoil type, raising the twist and chord values in the root part of the blade decreases the Ts.It should be noted that during the startup process, the angles of attack along the blade are high, which decrease by raising the twist angle, and this increases the startup torque.However, since the blade root elements have a smaller share in power generation, the optimization algorithm applies these increments (twist and chord) only in this area so that the Cp does not decline much.Figure 16 shows the variations of startup torque during the startup time of the selected airfoils for n = 0.8.Regardless of the airfoil type, as the blade begins to rotate, the startup torque initially decreases, and then it starts to increase.When the blade is stationary, the drag force is perpendicular to the direction of rotation and does not affect the startup torque.As the blade begins to rotate, the drag force reduces the startup torque, but then the startup torque starts to increase.Figure 16 shows that the blade that uses the USNPS4 airfoil produces the highest startup torque, but having a high moment of inertia raises its Ts.The lowest startup torque is produced in the blade with the BW-3 airfoil, but, as mentioned, its low moment of inertia is the advantageous aspect of this airfoil in areas with low wind speed.It is noteworthy that the startup torque is a function of two parameters, including the twist angle and the chord length (see Equation (10)).The difference observed in the values of startup torque for the selected airfoils is due to the different twist angle and chord length values that the optimization algorithm has considered for each airfoil.Finally, Figure 17 summarizes the performance of the selected airfoils in terms of Cp and Ts.The arrangement of these airfoils is based on exhibiting the highest Cp in windy areas and having the shortest Ts in areas with low wind speeds.

Conclusions
In the present study, the effect of the airfoil type was investigated on the power coefficient, as well as the startup time of a 1 kW small wind turbine by developing a numerical code.For this purpose, 10 airfoils, namely, BW-3, E387, FX 63-137, S822, S834, SD7062, SG6040, SG6043, SG6051, and USNPS4, were selected and the differential evolution (DE) optimization technique was employed to determine the best geometric shape of the blades.The twist angle and chord length were considered as design variables while maximizing the power coefficient and minimizing the startup time were defined as design goals.The blade element moment (BEM) method was used to calculate the power coefficient and the startup time.The BEM code and the optimization tools were both validated.The main results can be summarized as:  Regardless of the type of airfoil, using ideal equations to determine the twist angle and chord length to maximize the power coefficient, raises the turbine startup time;  The SG6043 airfoil has the highest power coefficient while the S822 and SG6040 airfoils have the lowest power coefficients.The reason for the superiority of SG6043 is its high lift-to-drag ratio.It is highly recommended to use this airfoil in windy areas where the purpose of designing small wind turbines is to achieve the maximum power coefficient;  Among the optimal blades for achieving the maximum power coefficient, the blades with the FX 63-137 and USNPS4 airfoils have the shortest and longest chord lengths, respectively.This has caused the power coefficient not to be as high as expected, despite the high lift-to-drag ratio of the FX 63-137 airfoil;  Regardless of the airfoil type, raising the twist angle and chord length in the root section reduces the turbine startup time;  From the startup viewpoint, the BW-3 airfoil has the best performance among the selected airfoils.This is due to the low inertia of the blades fitted with this airfoil.Therefore, in areas with low wind speeds where having a low startup time is greatly important, the use of this airfoil is highly recommended;  The S822, S834, and SG6040 airfoils have the highest startup time.The common aspect of these three airfoils is their high surface area;  Although the thinness of the airfoil is an advantage for reducing the blade moment of inertia and hence obtaining a better performance of the turbine at low wind speeds, the blade fitted with a thinner airfoil does not necessarily have a lower startup time than the blade with a thicker airfoil.This is because the airfoil type completely affects the distribution of twist angle and chord length.This is accompanied by fundamental changes in the startup torque and moment of inertia, both of which play an influential role in the startup process of the turbine;  Regardless of the airfoil type, when the blade begins to rotate, the startup torque first decreases slightly and then starts to increase;  The highest startup torque is produced by the blade fitted with the USNPS4 airfoil and the lowest startup torque is produced by the blade fitted with the BW-3 airfoil.

Appendix A
Figure A1 illustrates the changes in Cl and Cd values for different airfoil blades which were investigated throughout this research.These values were fed into the developed numerical code as a database for the current numerical computations.It is necessary to explain that, in addition to the relevant references, these coefficients are also available on the website maintained by Professor Michael Selig's group at the University of Illinois at www.ae.illinoise.edu/m-seliq/(accessed on 5 August 2022).Numerous researchers in this field have employed the data acquired from this website for their study [48][49][50][51], and the use of this website as a valid source for obtaining lift and drag coefficient values was recommended in a reference book [19] [31], (b) E387 [30], (c) FX 63-137 [30], (d) S822 [31], (e) S834 [30], (f) SD7062 [31], (g) SG6040 [31], (h) SG6043 [31], (i) SG6051 [31], and (j) USNPS4 [31].

Appendix B
Table A1 presents the Cp values from the current study and the experimental research of Anderson et al. [45].This table also presents the deviation between these two sets of results.

Appendix C
Table A2 shows the optimal values of Cp, Ts, Jb, and QS0 for n = 1, 0.8, and 0.6 weighting factors.The values of n = 1 and n = 0.8 are two accepted values in the literature for the design of small wind turbine blades to operate in windy areas and areas with low wind speed, respectively.In this regard, the work of Pourrajabian et al. [20] and Rahgozar et al. [21] can be referred to.Considering small values for n results in impractical blades whose aerodynamic performance is poor.The table below summarizes the optimization results of the current research with n = 0.6, along with n = 1 and n = 0.8.As can be noticed from the table, regardless of the airfoil type, considering n = 0.6 drastically reduces the power coefficient, which is indeed the main goal of wind turbine blade design.

Figure 1 .
Figure 1.The shape of the selected airfoils.

Figure 2 .
Figure 2. Velocities on an element of the blade.

Figure 3 .
Figure 3. Axial and tangential forces on a blade element.

Figure 5 .
Figure 5.The twist angle and chord length distributions of the turbine in the work of Anderson et al. [45].

Figure 7 .
Figure 7. Changes in the Cp at different λ values for the present study and the research of Anderson et al. [45].

Figure 8 .
Figure 8.The twist angle (top) and chord length (bottom) distributions of the blade fitted with the BW-3 airfoil.

Figure 12 .
Figure 12.The torque generated by the blade elements during power generation for n = 1.

Figure 14 .
Figure 14.The surface area (A) of the selected airfoils.

Figure 16 .
Figure 16.Variations of the startup torque for the selected airfoils (n = 0.8).

Figure 17 .
Figure 17.Comparison of the performance of the selected airfoils.

Table 2 .
The considered input parameters in the DE algorithm.

Table 4 .
Optimal values of Cp and Ts for the selected airfoils. .

Table A1 .
Detailed changes of Cp at different λ values for the current research and the study of Anderson et al. [45].

Table A2 .
The optimal values of Cp, Ts, Jb, and QS0 for selected airfoils for n = 1, 0.8, and 0.6.