Next Article in Journal
A Volunteered Geographic Information-Based Environmental Decision Support System for Waste Management and Decision Making
Next Article in Special Issue
Investigations of HAWT Airfoil Shape Characteristics and 3D Rotational Augmentation Sensitivity Toward the Aerodynamic Performance Improvement
Previous Article in Journal
Eco-Friendly Cellulose Nanofiber Extraction from Sugarcane Bagasse and Film Fabrication
Previous Article in Special Issue
A CFD Based Application of Support Vector Regression to Determine the Optimum Smooth Twist for Wind Turbine Blades
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Nonlinear Optimization Models for Wind Power Plants Using Box-Behnken Design of Experiment: A Case Study for Turkey

1
Department of Industrial Engineering, Yozgat Bozok University, Yozgat 66200, Turkey
2
Department of Environmental Engineering, Marmara University, Istanbul 34722, Turkey
3
Department of Industrial Engineering, Gaziantep Islam Science and Technology University, Gaziantep 27010, Turkey
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(15), 6017; https://doi.org/10.3390/su12156017
Submission received: 8 June 2020 / Revised: 16 July 2020 / Accepted: 22 July 2020 / Published: 27 July 2020
(This article belongs to the Special Issue Consistent Computational Approaches for Wind Energy Applications)

Abstract

:
This study aims to develop an optimization model for obtaining the maximum benefit from wind power plants (WPPs) to help with reducing external dependence in terms of energy. In this sense, design of experiment and optimization methods are comprehensively combined in the wind energy field for the first time. Existing data from installed WPPs operating in Turkey for the years of 2017 and 2018 are analyzed. Both the individual and interactive effects of controllable factors, namely turbine power (MW), hub height (m) and rotor diameter (m), and uncontrollable factor as wind speed (m/s) on WPPs are investigated with the help of Box-Behnken design. Nonlinear optimization models are utilized to estimate optimum values for each decision variable in order to maximize the amount of energy to be produced for the future. Based on the developed nonlinear optimization models, the optimum results with high desirability value (0.9587) for the inputs of turbine power, hub height, rotor diameter and wind speed are calculated as 3.0670 MW, 108.8424 m, 106.7597 m, and 6.1684 m/s, respectively. The maximum energy output with these input values is computed as 9.952 million kWh per unit turbine, annually. The results of this study can be used as a guideline in the design of new WPPs to produce the maximum amount of energy contributing to supply escalating energy needs by more sustainable and clean ways for the future.

Graphical Abstract

1. Introduction

Energy in different forms, such as conventional fossil fuels, nuclear, and renewable energy, is utilized throughout the world for human needs [1]. Considering the adverse effects of traditional sources, the significance of renewable energy about decreasing polluting emissions for diminishing the impact on human health, environment and climate change is undeniable [2]. The World Energy Forum estimated that the time remaining for the depletion of fossil-based oil, coal and gas reserves is no more than 100 years. As alternative energy resources can reproduce themselves, exhaustion is no longer a worry for policymakers in response to ascending energy needs [3]. Other reasons behind the inclination for renewable energy resources include fluctuations in oil prices, desire to provide energy security without interruption, and environmental impacts resulting from emissions of carbon dioxide [4]. The U.S. Energy Information Administration (EIA) proposes that renewables including wind, solar and hydroelectric power will be the fastest growing resources of energy by 2050, compared to the conventional resources such as coal, natural gas, and nuclear energy [5].
The importance of renewable energy sources for electricity production is increasing worldwide [6]. On a global scale, the portion of alternative energy resources in the total energy usage was 15% in 2018 and anticipated to rise by 28% in 2050. Among the renewables, wind energy is predicted to be the main contributor to this increase [7]. Wind turbines convert the linear kinetic energy of wind into rotational kinetic energy, which is then converted into electrical energy [8,9].
Wind is an accessible and readily obtainable energy source all over the world. The International Panel on Climate Change (IPCC) Special Report emphasizes that it is necessary to reduce emissions by at least 49%, as compared to 2017’s value throughout the world by 2030 and reach net-zero carbon emissions by 2050 to prevent the impacts of climate change [10]. Contribution of wind energy is essential for controlling emissions on a global scale at the targeted time, considering that the energy sector is responsible for about 40% of CO2 emissions, globally. It was estimated in the 4th IPCC Assessment Report that the capacity of wind energy will escalate to 1000 GW by 2020, which corresponds to 2600 TWh of electricity that avoids the emission of 1500 million tons of carbon dioxide annually [7].
With the extensive use of wind power, its share in electricity production is also increasing [11]. Approximately 14% of the EU’s consumed energy per annum was supplied by wind energy in 2018. Wind energy has become the most suitable choice for a wide range of places all over the world, owing to its quickly decreasing cost, both onshore and offshore [12]. Wind itself is the fuel used to produce electricity, differently from most power plants requiring a source [13]. As an eco-friendly energy source, wind energy does not cause any emissions of nitrous oxides (NOx) and sulfur dioxide (SO2) which are emitted from conventional energy sources [14]. Additionally, unlike fossil fuels, wind energy does not release any CO2 during its operation [15]. Air pollution, caused by those pollutants, is not only a concern for our climate but also a severe threat for human health. Indeed, pollution of air was reported to be responsible for one out of every nine deaths in 2012 [16]. In this context, wind energy is significant for the positive development of air quality, too [7]. Generating electricity with this technology may contribute to the prohibition of millions of tons of carbon and other emissions while saving several billion barrels of oil. Wind energy, which is both easy and fast to set up and operate, even in areas that are rugged and difficult to reach, offers a feasible source of electricity for developing countries [17]. The power from the wind is neither exhausted nor does the price rise [18], and finally, wind is present both at night and day with high conversion performance to electrical energy [19].
Turkey has an important geopolitical position by establishing a bridge between Asia and Europe. The growth of the industry and the economy with the increase in the population in Turkey causes the demand for energy also to rise. Wind energy offers an excellent alternative source of energy regarding energy security, considering the high energy dependency of Turkey on external resources [20]. By the end of 2019, the energy sources of Turkey based on installed capacity, comprised of hydraulic energy with 31.4%, natural gas with 28.6%, coal with 22.4%, wind with 8.1%, solar with 6.2%, geothermal with 1.6% and other resources with 1.7% [21]. Undoubtedly, the utilization of alternative energy resources plays a significant role in being energy independent. The predicted potential of wind energy of Turkey is 48,000 MW that makes it a very favorable sustainable option for the future as compared to other renewables [20]. In this respect, the primary aim of this study is to obtain maximum benefit from the wind energy potential to help minimize the dependence on foreign fossil fuel sources besides dealing with the other environmental concerns explained above. This work represents the first application of the design of experiment (DOE) method and optimization technique combined for the field of wind energy, comprehensively. Firstly, the approach of design of experiment is used to determine and analyze the factors that are effective for wind power generation. Secondly, the optimization technique is employed to calculate the optimum values of the factors in order to maximize the amount of electrical energy produced.
In the literature, different optimization models were developed to get maximum energy or minimum cost subject to only one factor or one constraint such as wind speed to maximize produced power, material types to minimize cost, etc. Negm and Maalawi compared five separate optimization models for wind turbine tower design to achieve optimum solutions for the productivity of the system [22]. Christie and Bradley analyzed the optimization of land utilization for wind farms. The goal was to maximize power density, the power per unit area, of needed land for which the wake effect was taken into account [23]. Gao et al. investigated wind turbine layout optimization to maximize energy output and minimize the capital cost, for which three different conditions related to the wind speed and wind direction were implemented for a 2 km × 2 km wind farm [24]. Chan et al. studied blade shape optimization of Savonius wind turbine to enhance its power coefficient, in which three models were developed to optimize the blade shape of this specific turbine [25]. Differently from the examples mentioned above, in this research, we applied an optimization technique considering four factors, namely turbine hub height, rotor diameter, turbine power and wind speed, for the maximum output of the wind power. In addition, we determined how these variables individually and interactively have influenced the amount of energy produced using the design of experiment method. In this work, the existing data from the wind power plants operating in Turkey were analyzed to reach optimum values of the decision variables for the maximized energy output.
The rest of the study is structured as follows. Section 2 delivers the methodology of the study, which consists of the proposed three main stages that combine statistical analysis, Box-Behnken design of experiment and nonlinear optimization models. Section 3 discusses the details of the results of the methodology using WPPs’ data in Turkey. Section 4 provides the conclusions of the study.

2. Materials and Methods

The methodology of the study was developed to optimize the wind power plants (WPPs) in order to produce maximum energy output using the historical data of operational WPPs in Turkey. There are 155 WPPs located in Turkey from which the data belonging to 2017 and 2018 were collected [21]. The wind turbines of Turkey are installed in six of the seven physiographic regions of Turkey. However, they are more concentrated in the Aegean and Marmara regions as quantitively displayed in Figure 1. The main reasons for the establishment of a large number of WPPs in these areas are the large population, heavy industry, high-altitude lands having the wind speed higher than the average of Turkey, etc. The total number of operational WPPs is 63 in Aegean Region, 51 in Marmara region, 21 in Mediterranean region, 12 in Central Anatolian region, 6 in Black Sea region, and 2 in Southeast Anatolian region. There have been no WPPs installed in Eastern Anatolian region so far. Located in the Aegean region, İzmir is the city with the highest number of wind farms with 38 WPPs. Balıkesir, in the Marmara Region, ranks the second with 21 WPPs.
The data of 43 WPPs out of 155 were ignored due to lack or contradictory data and the fact that some wind farms have not started production yet. The contradictory data, especially for the hub height and the amount of annual produced energy, were encountered in different sources, so these WPPs were left out of the scope of this study. Additionally, some WPPs have different types of turbines (i.e., distinct turbine power, hub height and rotor diameter). Nevertheless, the data of the energy generation amount is available for the entire facility in total. Thus, these were not taken into account in order to avoid confusion in the calculations. In general, the areas with high wind speed are preferred to gain energy from wind effectively, but this is not the only factor to consider. We used hub height, rotor diameter, turbine power, and wind speed at 50 m as factors or input variables and we obtained the annual energy production amount for the years of 2017 and 2018 of the installed 112 facilities as a response or output variable (see Figure 2). Due to the fact that the number of turbines varies between the years in some WPPs, the average annual energy production per unit turbine was taken into consideration in each phase of the study.
The wind speed has an important effect on the power produced by the wind turbine. According to the study of wind speed variability and wind power potential over Turkey, Aslan et al. showed that the wind speed fitted the Weibull distribution based on hourly data of 335 stations covering 34 years [2]. Secondly, the hub height of the turbine should be taken into account concerning energy generation because of its impact on wind speed. The wind speed increases as hub height increases. The wind power-law helps to extrapolate wind speed data at different hub heights as a practical tool utilized in wind power applications. The wind profile power-law relationship is expressed as [26] (Equation (1)):
v v 0 = ( h h 0 ) α
where v is the wind speed at hub height h , v 0 is the wind speed at the initial hub height, h 0 , and α is an empirically derived wind speed power-law coefficient depending on the stability of the air. The equation that shows the relationship between rotor diameter and wind speed is expressed as follows [27] (Equation (2)):
v = v 0 [ 1 2 α ( l 2 l 2 + α x ) 2 ]
where l / 2 is the radius of the rotor, and x is axial induction factor [28]. With the aid of the combination of the equations explained above, the extractable wind power by the turbine rotor is given by (Equation (3)):
P m a x E = 1 2 ρ A v 3 c p
where A donates the swept area of the turbine, c p symbolizes the power coefficient. The power coefficient is defined as the percent of the extraction of the energy by a wind turbine which is 59.3 % at the most, theoretically. The value of cp is given as [29] (Equation (4)):
c p = 0.593 ,   referred   to   as   the   Betz   limit
The swept area of the turbine (A) in Equation (4) can be computed with the help of the equation for the area of a circle considering the length of the turbine blades as the radius of the circle during rotation [30] (Equation (5)):
A = π l 2
where l represents the value of the blade length as the radius of the rotor.
The flow chart in Figure 3 illustrates the phases of the study in detail. There are three main phases to build the methodology of the study, which are statistical analysis, Box-Behnken design of experiment (DOE) and nonlinear optimization model. Each phase constitutes the infrastructure of the next phase.

2.1. Statistical Analysis

The production of energy amount by the WPPs operated in the provinces of Turkey was organized for the analysis. Statistical analysis is the basis for constructing the next steps of the methodology of this study. The first phase before applying the DOE with the optimization method is to arrange the statistical data descriptively. Quartiles (denoted as q1, q2 and q3) were used in order to separate the whole data set, belonging to 112 WPPs operating in Turkey, into three parts that help to have an idea about the distribution of the data. Table 1 shows the descriptive statistics of the factors and responses.
The Box plots of the four factors, namely turbine power, hub height, rotor diameter and wind speed, were generated to visually show the spread of the data set by making use of minimum value, q1 (first quartile), q2 (second quartile or median), q3 (third quartile) and maximum value of each variable. As seen from Figure 4, the data of turbine power are mostly centered in the 1st quartile (2.1 MW–2.75 MW). Similarly, the data of rotor diameter and wind speed are almost concentrated in the q1. However, hub height values are distributed more intensely in the 3rd quartile ranging from 92 m to 117 m. The intervals of quartiles for each factor are tabulated in the design part in detail.

2.2. Box-Behnken Design of Experiment

Designed experiments are useful to ascertain cause-and-effect relationships for the process of interest, which in turn help on the optimization of the response variable. DOE assists in maximizing the correctness of the result or response of a process conversely with minimizing the quantity of experiments, to escalate the system performance. DOE not only helps to specify the most important factors that influence the output but also points out the interactions among them [31]. The components of the design of experiment approach are composed of factors or variables, levels and responses [32,33]. Factors, as controllable and uncontrollable, are the input variables of the process [34]. Levels are the settings related to each factor of the system. Outputs or responses are the constituents that determine the process performance.
Models are frequently utilized in the engineering field to formulate and solve any kind of problem. The model employed for DOE to show the relationship between the response (y) (i.e., produced energy annually) and the input variables (x’s) (i.e., hub height, turbine power, rotor diameter and wind speed) is expressed as (Equation (6)):
y = β 0 + i = 1 p β i x i + i = 1 p j = 1 p β i j x i x j + i = 1 p j = 1 p k = 1 p β i j k x i x j x k                 (   i j k ) .
where β 0 is the overall mean response, β i is the main effects for each factor ( i = 1 , 2 , 3 , . , p ) , β i j is the two-way interaction between i th and j th factors, β i j k is the three-way interaction between i th , j th and k th factors.
Box-Behnken design (BBD), one of the subsets of DOE methodology, helps to develop higher-order models that need fewer runs as compared to other factorial designs. Thus, application of BBD not only saves time but also reduces cost with the help of decreasing the number of experiments. BBD aims to optimize the output of a system with fewest runs by combining the mathematical and statistical methods necessary for the optimization modeling and analysis of experiments where output is affected by multiple factors [33].
In this study, the BBD tool was applied to improve wind energy production process and point out the interactions among considered parameters. This design type was chosen because of the fact that it requires less design studies. More specifically, in the full-factorial experimental design with three factors having three levels of each factor, the number of designs is 162 (81 × 2) with two replicates, while the number of required designs is only 54 (27 × 2) in BBD tool. Based on BBD method, the process yield (y) is given as (Equation (7)):
y = f   ( x 1   ,     x 2   ,     x 3   ,     ,   x n ) + ϵ
where x i   denotes the numbers of factors maximizing the yield, y, of a process and 𝜖 symbolizes the noise or error observed in the response y.
The limits for the levels of variables required for the BBD design were individually determined using Box plots of the factors. The number of levels for each parameter was chosen as 3. The levels and limit values of factors are depicted in Table 2.
In this research, four different factors, each having three levels, were considered with regard to their effects on wind energy production. While the number of basic designs is 27 in Box-Behnken design, the total number of data used is 54 by replicating the study twice. Table 3 represents all possible combinations of factors according to the Box-Behnken design of experiment both with coded levels and actual levels of each variable. Factor 1, factor 2, factor 3 and factor 4 terms are the basic names used in Box-Behnken design of experiment, and their values are denoted as −1.0, 0 and 1.0 for level 1, level 2 and level 3, respectively. The factors used for this study were replaced by actual values corresponding to these notations.

2.3. Nonlinear Optimization Models

Optimization models are developed so as to determine the optimum decision variables that meet the objective function(s) (by maximizing and/or minimizing) of any problem to be examined [35]. In many interesting maximization and minimization problems, the objective(s) or constraint(s) may not follow a linear trend, making them hard to solve. This type of optimization model is called to be nonlinear [36]. We observed that the optimization model obtained in this study fits into a nonlinear mathematical modeling type. The reason is that the factors were taken into account both singularly and interactively in the optimization model. The objective function of the nonlinear optimization model is formulated as the following [37] (Equation (8)):
max o r min   Z = f   ( x 1     ,   x 2     ,   x 3     ,   ,   x n ) .
The set of all points ( x 1   ,   x 2   ,   x 3   ,   ,   x n ) such that x i is a real number is R n . The following subsets of R n (called intervals) will be of particular interest (Equation (9)):
[ l i ,   u i ] = all   x   satisfying   l i x i u i
where l i ,   u i represent lower limit and upper limit for the constraints, respectively.
The aim of this experiment is to get the maximum amount of energy produced in wind power plants. The objective of the proposed nonlinear programming model is to maximize the estimated objective function subject to the boundary constraints over the experimental region. The nonlinear model that contains all factors neglecting significance levels of the model is given as below (see Appendix A for detailed equation) (Equation (10)):
y k W h = β 0 + β 1 A + β 2 B + β 3 C + β 4 D + β 11 A 2 + β 22 B 2 + β 33 C 2 + β 44 D 2 + β 12 A B + β 13 A C + β 14 A D + β 23 B C + β 24 B D + β 34 C D .
This regression equation is defined as an objective function in this study, where y k W h refers to the output of the WPPs (as million kilowatt-hours). Then, the full optimization model is expressed as (Equation (11)):
t o   m a x i m i z e   y k W h = { [ T P ( A ) ,   h ( B ) ,   l ( C ) ,   v ( D ) ] 1 ,   [ T P ( A ) ,   h ( B ) ,   l ( C ) ,   v ( D ) ] 2 ,   ,   [ T P ( A ) ,   h ( B ) ,   l ( C ) ,   v ( D ) ] j } s u b j e c t   t o 0.8 MW x A j 3.6 MW   the   value   of   factor   A   for   the   j t h   WPP 64 m x B j 117 m   the   value   of   factor   B   for   the   j t h   WPP 52.9 m x C j 126 m       the   value   of   factor   C   for   the   j t h   WPP 5.4 m / s x D j 7.75 m / s       the   value   of   factor   D   for   the   j t h   WPP x i j ,   n o n n e g a t i v e x i j R n
where i = { A , B , C , D } , and j = { t h e   a s s i g n e d   o f   W P P } . According to the BBD of experiment, a total of 27 × 2 different WPPs have been designed with different factor values for optimum WPP design, generating maximum annual energy. The limits of the constraints have been determined by BBD.
Furthermore, the development of additional optimization models for the mean and standard deviation of the response variable are required simultaneously to improve the optimization process. The idea behind this concept is to maximize the mean value and minimize the standard deviation of outputs symbolized by y [38]. The use of BBD method necessitates to minimize the standard deviation to reach the target value for industrial problems. Thus, we have extended the developed optimization models with this approach since wind power generation is a real industrial issue worldwide. The polynomial models for the mean ( y μ ) and standard deviation ( y σ ) of the response were formed by Vining and Myers as following [39] (Equations (12) and (13)):
y μ = β 0 + i = 1 k β i x i + i = 1 k β i i x i 2 + i = 1 k j = 1 k β i j x i x j + ε μ               i < j ,      
y σ = γ 0 + i = 1 k γ i x i + i = 1 k γ i i x i 2 + i = 1 k j = 1 k γ i j x i x j + ε σ               i < j      
where β 0 ,   β 1 ,   β 2 ,   ,   β k , and γ 0 ,   γ 1 ,   γ 2 ,   ,   γ k are regression coefficients,   x i   is the number of the factors, and ε μ and ε σ are the random errors of the mean and standard deviation of the response variable, respectively. Determining the values which minimize performance deviation from the target value is preferable. The additional optimization models were introduced as below (see Equations (14) and (15)):
M i n i m i z e         y ^ σ = f   ( x 1   ,     x 2   ,     x 3   ,     ,   x n ) σ S u b j e c t   t o       f   ( x 1   ,     x 2   ,     x 3   ,     ,   x n ) μ = τ μ τ σ   ,   x i ,   n o n n e g a t i v e x i R n
M a x i m i z e         y ^ μ = f   ( x 1   ,     x 2   ,     x 3   ,     ,   x n ) μ S u b j e c t   t o           f   ( x 1   ,     x 2   ,     x 3   ,     ,   x n ) σ = τ σ τ σ ,   x i ,   n o n n e g a t i v e x i R n
where τ μ and τ σ are the target values of the mean and standard deviation of outputs. y ^ μ and y ^ σ represent the functions of regressions. Target values represent the points corresponding to the constraints of the optimization models developed for the study.

3. Results and Discussion

This study involves the application of Box-Behnken design that is part of the design of experiment with a total of 54 runs (27 base runs with two replications), while considering four variables, each with three levels, for the wind power plants to provide maximum energy produced. The factors, which are turbine power, hub height, rotor diameter and wind speed denoted as A, B, C and D, respectively, were analyzed statistically. Furthermore, R2, the value of the coefficient of multiple terms for the regression showing the statistical measure of how close the data are to the fitted regression line, was calculated as 0.9179 (R2 adj is 0.8885). With regard to the statistical analysis, the individual and interactive effects of factors were determined by calculating F-value, t-Ratio and p-value. According to Table 4, it is seen that the overall model, the factors (linear) and the interactions with themselves (square) are important for the production of wind energy (p-value ≤ 0.05). However, when the 2-way interactions are considered, it is understood that the interactions of factor A (turbine power) with factors B (hub height) and C (rotor diameter) are insignificant in terms of the amount of energy generated (p-value ≥ 0.05).
The main effects plot was made to show how the mean of response changed concerning the levels of each factor. The impact of factors (A, B, C and D) on the response variable (Y) is separately demonstrated in Figure 5. It has been observed that the response variable increases up to a certain point with the increase of factors A and C. Increasing the level of factor B had a very low positive effect on the outcome. As seen from the figure, the increase in the level of wind speed (factor D) escalated the result variable up to a certain point and then started decreasing. The similarity between these separate curves that show the effect of each factor on energy production in WPPs was that they all followed a nonlinear trend.
After investigating individual and interactive impacts of each factor on response with their significance degrees by statistical analysis, a nonlinear optimization model (defined in the previous section) was established to find out the optimum values of input and output variables. The objective function of the initial model is to maximize the annual energy output per unit turbine, which is defined as Target0.
The desirability values check the validity of the results obtained in the optimization models. Desirability value, ranging from 0.0 to 1.0, evaluates how well the combination of variables meets the objectives defined for the answers. The closer this value to 1.0, the better the desired objective function value is. The optimum results with high desirability value (0.9587) for factors A, B, C, D were calculated as 3.0670 MW, 108.8424 m, 106.7597 m, 6.1684 m/s, respectively, to produce the maximum amount of annual energy which is estimated to be 9.952 million kWh per unit turbine. The best ten feasible solutions, according to their desirability values, are represented in descending order in Table 5. This table includes the top ten feasible results for turbine power (A), hub height (B), rotor diameter (C) and wind speed (D) and gives the annual average amount of energy per unit turbine with these configurations.
The optimization process was improved with additional optimization models developed for the mean and standard deviation of the response variable. The optimization results when the mean of output is maximized while keeping the standard deviation constant is stated as Target1 (max y k W h of output subject to constant y ^ σ ). The optimization results of the mean of output for the factors A, B, C, D and the response were calculated as 2.9395 MW, 101.1193 m, 94.7760 m, 6.4038 m/s, and 9.390 million kWh per unit turbine, respectively, with the desirability value of 0.9246.
The second optimization model presents the situation in which standard deviation is minimized while holding mean of response at a constant level that is expressed as Target2 (min y ^ σ subject to constant y k W h ). The optimum results, with high desirability value of 0.9144, for the factors A, B, C, D and the response were computed as 1.7665 MW, 106.2929 m, 95.4418 m, 6.8004 m/s, and 7.143 million kWh per unit turbine, respectively. The minimized standard deviation value of response was computed as 0.24680 million kWh.
The multi-objective optimization (max y k W h and min y ^ σ ) of output results, in which mean of response is maximized (Target1) while standard deviation is minimized (Target2), simultaneously, are shown in Table 6. It can be evaluated as the combinational representation of the results of the optimization models concerning both mean and standard deviation of output. The optimum results, having desirability value of 0.9156, for the factors A, B, C, D and the output were calculated as 2.7100 MW, 94.5152 m, 112.0857 m, 6.6581 m/s, and 9.684 million kWh per unit turbine, respectively. The minimized standard deviation value of output was 0.28630 million kWh.
The results of all optimization models with their desirability values were displayed together in Table 6. The averages of all optimum results were computed for the factors A, B, C, D as 2.62075 MW, 102.6925 m, 102.2658 m, and 6.50767 m/s, respectively, with a high average desirability value of 0.92833, to get maximum energy production which was 9.952 million kWh in Target0, 9.537 million kWh in Target1 per unit turbine. The average of the standard deviation of the response variable was calculated as 0.26655 million kWh in Target2.
In terms of the validation of the method used, we have optimized each decision variable separately by fixing the other three decision variables according to the mean value of parameters which was determined as the reference case. Table 7 shows the comparative results of decision variables of multi-objective functions’ optimization models and the findings after validation process. The percent difference between the optimum values of each variable is way below 1%, demonstrating that the applied method and the results are more than 99% confident.
The European Wind Energy Union classifies wind speeds as nearly good for 6.5 m/s, good for 7.5 m/s and very good for 8.5 m/s which is suitable for energy generation from wind [40]. All of the optimization models developed propose that even in areas in nearly good wind speed group, the higher amount of wind energy can be produced with optimizing other variables. In terms of comparison, Lee et al. investigated the optimization of hub height to get maximum annual net profit. They limited the hub height value between 30 m and 300 m in the optimization model developed for WPPs. The optimum hub height value they find is 125.8 m with V = 7 m/s and α = 0.25 for case 1 and 219.2 m with V = 5 m/s and α = 0.3 for case 2 [41]. One of the most important factors affecting the initial capital cost in WPPs is the hub height. Therefore, we can assume that the hub height value we calculated in this study, which is 105 m, is more reasonable in terms of cost. Lee et al. calculated the optimum hub height value as 75 m, which is calculated by considering the cost of establishment. However, the amount of energy produced according to the calculated hub height is low in their study [41]. Cheng et al. studied the optimization of rotors for higher power capture from the wind. They tried to optimize Darrieus-type vertical axis wind turbine rotor having a constant turbine power of 5 MW. The optimum value was computed as 120.96 m for rotor diameter and 143 m for hub height. Moreover, the reference height for wind speed was considered as 79.78 m [42].
In short, it was observed that the desirability values, which indicated how much the combination of the variables have achieved the target defined for the response variable(s), for all of the developed optimization models were close to 1 and very similar to each other. This closeness supports the feasibility of our study, which might assist not only operators but also researchers, entrepreneurs and governments for the future.

4. Conclusions

This study focuses on gaining maximum benefit from wind energy that has the potential to meet increasing energy demands in a cost-effective and environmentally friendly way without depending on fossil fuel sources. In this respect, the data from operating wind power plants in Turkey were collected to apply the design of experiment and optimization techniques together in the wind energy field for the first time in such a comprehensive manner. Turbine power, hub height, rotor diameter and wind speed (at 50 m) were selected as decision variables whose limits required for the design were determined using Box plots. Both the individual and interactive effects of controllable and uncontrollable factors on wind power generation were investigated with the help of Box-Behnken design. The results showed that all the input factors were significant for wind energy production based on p-value, which is the sign of the importance of the results. The 2-way interactions of the design parameters, that is, the interaction of factors with each other, were found to be significant, too, except for the interactions of turbine power with hub height and rotor diameter.
After analyzing individual and interactive impacts of each factor on the response with significance degrees, a nonlinear optimization model was established in order to find out the optimum values of input and output variables. The optimum results with high desirability value (0.9587) for the factors A, B, C, D were calculated as 3.0670 MW, 108.8424 m, 106.7597 m, 6.1684 m/s, respectively, to produce the maximum amount of energy which was 9.952 million kWh per unit turbine, annually. The optimization process was improved with additional optimization models developed for the mean and standard deviation of the response variable in the way of maximizing mean of output, minimizing standard deviation of response and combining them both. It was observed that the desirability values for all of the developed optimization models were close to 1 and very similar to each other. The similarity between these results advocates the feasibility of our study to gain the maximum amount of energy from wind power plants. The results may serve as guidelines to contribute to satisfying the escalating energy needs of the future by more sustainable and clean ways.

Author Contributions

Conceptualization, Y.A.A., M.T. and A.A.; methodology, Y.A.A. and A.A.; software, Y.A.A.; validation, M.T., K.E. and A.A.; investigation, Y.A.A.; data curation, Y.A.A.; writing—original draft preparation, Y.A.A., M.T. and K.E.; writing—review and editing, Y.A.A. and M.T.; visualization, A.A.; supervision, M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the Ministry of Energy and Natural Resources for sharing WPPs’ annual energy production amount and the features of WPPs’ structures, and the Turkish State Meteorological Service for providing the wind data for all stations.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The nonlinear model that contains all factors neglecting significance levels of the model is given as below (see Equation (A1)):
y k W h = 51897234 848278 A + 134692 B + 177710 C + 14333926 D 682673 A 2 596 B 2 363 C 2 935391 D 2 + 8582 A B + 8227 A C + 752165 A D + 801 B C 19659 B D 26156 C D
This regression equation is defined as an objective function in this study, where y k W h refers to the output of the wind power plants (as kilowatt-hour). Then, the full optimization model is expressed as (Equation (A2)):
M a x i m i z e   y k W h = 51897234 848278 A + 134692 B + 177710 C + 14333926 D 682673 A 2 596 B 2 363 C 2 935391 D 2 + 8582 A B + 8227 A C + 752165 A D + 801 B C 19659 B D 26156 C D
The optimization model for the minimum standard deviation of the response was given as (Equation (A3)):
M i n i m i z e     y ^ σ = 17411963 143235 A 66017 B + 10071 C 4195339 D 3560 A 2 + 281 B 2 172 C 2 + 304941 D 2 + 2765 A B + 2036 A C 26994 A D + 44 B C 94   B D + 1642 C D S u b j e c t   t o y ^ μ = 51897234 1848278 A + 134692 B + 177710 C + 14333926 D 682673 A 2 596 B 2 363 C 2 935391 D 2 + 8582 A B + 8227 A C + 752165 A D + 801 B C 19659 B D 26156 C D = 6885648
and the optimization model with the maximum mean of the response variable was given as (Equation (A4)):
M a x i m i z e   y ^ μ = 51897234 1848278 A + 134692 B + 177710 C + 14333926 D 682673 A 2 596 B 2 363 C 2 935391 D 2 + 8582 A B + 8227 A C + 752165 A D + 801 B C 19659 B D 26156 C D S u b j e c t   t o y ^ σ = 17411963 143235 A 66017 B + 10071 C 4195339 D 3560 A 2 + 281 B 2 172 C 2 + 304941 D 2 + 2765 A B + 2036 A C 26994 A D + 44 B C 94 B D + 1642 C D = 461089

References

  1. Asif, M.; Muneer, T. Energy supply, its demand and security issues for developed and emerging economies. Renew. Sustain. Energy Rev. 2007, 11, 1388–1413. [Google Scholar] [CrossRef]
  2. Arslan, H.; Baltaci, H.; Akkoyunlu, B.O.; Karanfil, S.; Tayanc, M. Wind speed variability and wind power potential over Turkey: Case studies for Çanakkale and İstanbul. Renew. Energy 2020, 145, 1020–1032. [Google Scholar] [CrossRef]
  3. Toklu, E. Overview of potential and utilization of renewable energy sources in Turkey. Renew. Energy 2013, 50, 456–463. [Google Scholar] [CrossRef]
  4. Coetzee, R.A.M.; Mwesigye, A.; Huan, Z. A numerical analysis and optimization of the dynamic performance of a multipurpose solar thermal system for residential applications. Sci. Technol. Built Environ. 2018, 24, 1156–1173. [Google Scholar] [CrossRef]
  5. U.S. Energy Information Administration. Annual Energy Outlook 2019 with Projections to 2050; Office of Energy Analysis, U.S. Department of Energy: Washington, DC, USA, 2019.
  6. Matayoshi, H.; Howlader, A.M.; Datta, M.; Senjyu, T. Control strategy of PMSG based wind energy conversion system under strong wind conditions. Energy Sustain. Dev. 2018, 45, 211–218. [Google Scholar] [CrossRef]
  7. Kaygusuz, K. Wind energy status in renewable electrical energy production in Turkey. Renew. Sustain. Energy Rev. 2010, 14, 2104–2112. [Google Scholar] [CrossRef]
  8. Çapik, M.; Yilmaz, A.O.; Çavuşoĝlu, I. Present situation and potential role of renewable energy in Turkey. Renew. Energy. 2012, 46, 1–13. [Google Scholar] [CrossRef]
  9. Talebi, J.; Ganjefar, S. Fractional order sliding mode controller design for large scale variable speed wind turbine for power optimization. Environ. Prog. Sustain. Energy. 2018, 37, 2124–2131. [Google Scholar] [CrossRef]
  10. Chen, Z.; Li, W.; Guo, J.; Bao, Z.; Pan, Z.; Hou, B. Projection of Wind Energy Potential over Northern China Using a Regional Climate Model. Sustainability 2020, 12, 3979. [Google Scholar] [CrossRef]
  11. Chen, L.; Ponta, F.L.; Lago, L.I. Perspectives on innovative concepts in wind-power generation. Energy Sustain. Dev. 2011, 15, 398–410. [Google Scholar] [CrossRef]
  12. Hales, D. REN21, Renewables 2018 Global Status Report; Renewable Energy Policy Network: Paris, France, 2018. [Google Scholar]
  13. Masukume, P.M.; Makaka, G.; Mukumba, P. Optimization of the Power Output of a Bare Wind Turbine by the Use of a Plain Conical Diffuser. Sustainability 2018, 10, 2647. [Google Scholar] [CrossRef] [Green Version]
  14. Cullen, J. Measuring the environmental benefits of wind-generated electricity. Am. Econ. J. Econ. Policy 2013, 5, 107–133. [Google Scholar] [CrossRef] [Green Version]
  15. Hossain, F.M.; Hasanuzzaman, M.; Rahim, N.A.; Ping, H.W. Impact of renewable energy on rural electrification in Malaysia: A review. Clean Technol. Environ. Policy 2015, 17, 859–871. [Google Scholar] [CrossRef]
  16. World Health Organization. Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease; World Health Organization: Geneva, Switzerland, 2016. [Google Scholar]
  17. Panwar, N.L.; Kaushik, S.C.; Kothari, S. Role of renewable energy sources in environmental protection: A review. Renew. Sustain. Energy Rev. 2011, 15, 1513–1524. [Google Scholar] [CrossRef]
  18. Kaplan, Y.A. Overview of wind energy in the world and assessment of current wind energy policies in Turkey. Renew. Sustain. Energy Rev. 2015, 43, 562–568. [Google Scholar] [CrossRef]
  19. Santhakumar, S.; Palanivel, I.; Venkatasubramanian, K. Building a low cost wind turbine in highways for rural house electricity demand. Environ. Prog. Sustain. Energy 2019, 38, 278–285. [Google Scholar] [CrossRef] [Green Version]
  20. Emeksiz, C.; Demirci, B. The determination of offshore wind energy potential of Turkey by using novelty hybrid site selection method. Sustain. Energy Technol. Assess. 2019, 36, 100562. [Google Scholar] [CrossRef]
  21. Ministry of Energy and Natural Resources. Wind Energy Potential Report; Republic of Turkey Ministry of Energy and Natural Resources: Ankara, Turkey, 2020. Available online: https://www.enerji.gov.tr/en-US/Pages/Wind (accessed on 6 June 2020).
  22. Negm, H.M.; Maalawi, K.Y. Structural design optimization of wind turbine towers. Comput. Struct. 2000, 74, 649–666. [Google Scholar] [CrossRef]
  23. Christie, D.; Bradley, M. Optimising land use for wind farms. Energy Sustain. Dev. 2012, 16, 471–475. [Google Scholar] [CrossRef]
  24. Gao, X.; Yang, H.; Lin, L.; Koo, P. Wind turbine layout optimization using multi-population genetic algorithm and a case study in Hong Kong offshore. J. Wind Eng. Ind. Aerodyn. 2015, 139, 89–99. [Google Scholar] [CrossRef]
  25. Chan, C.M.; Bai, H.L.; He, D.Q. Blade shape optimization of the Savonius wind turbine using a genetic algorithm. Appl. Energy 2018, 213, 148–157. [Google Scholar] [CrossRef]
  26. Gul, M.; Tai, N.; Huang, W.; Nadeem, M.H.; Yu, M. Evaluation of Wind Energy Potential Using an Optimum Approach based on Maximum Distance Metric. Sustainability 2020, 12, 1999. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, L.; Zuo, M.J.; Xu, J.; Zhou, Y.; Tan, A.C. Optimizing wind farm layout by addressing energy-variance trade-off: A single-objective optimization approach. Energy 2019, 189, 116149. [Google Scholar] [CrossRef]
  28. Miao, H.; Dong, D.; Huang, G.; Hu, K.; Tian, Q.; Gong, Y. Evaluation of Northern Hemisphere surface wind speed and wind power density in multiple reanalysis datasets. Energy 2020, 200, 117382. [Google Scholar] [CrossRef]
  29. Schubel, P.J.; Crossley, R.J. Wind turbine blade design. Energies 2012, 5, 3425–3449. [Google Scholar] [CrossRef] [Green Version]
  30. Yang, Y.; Li, H.; Yao, J.; Gao, W. Research on the characteristic parameters and rotor layout principle of dual-rotor horizontal axis wind turbine. Energy 2019, 189, 116270. [Google Scholar] [CrossRef]
  31. Wahdame, B.; Candusso, D.; Franc, X.; Harel, F. Design of experiment techniques for fuel cell characterisation and development. Int. J. Hydrog. Energy 2009, 34, 967–980. [Google Scholar] [CrossRef]
  32. Montgomery, D.C. Design and Analysis of Experiments; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef] [Green Version]
  33. Panda, D.; Satapathy, A.K.; Sarangi, S.K. Thermo-hydrodynamic analysis and optimal design of a GM cycle cryorefrigerator using response surface methodology and particle swarm optimization. Sci. Technol. Built Environ. 2019, 25, 1467–1481. [Google Scholar] [CrossRef]
  34. Antony, J. Design of Experiments for Engineers and Scientists; Elsevier: London, UK, 2014. [Google Scholar] [CrossRef]
  35. Bayraktar, M.; Fabrizio, E.; Perino, M. The extended building energy hub: A new method for the simultaneous optimization of energy demand and energy supply in buildings. HVAC R Res. 2012, 18, 67–87. [Google Scholar] [CrossRef]
  36. Winston, W. Operations Research-Applcations and Algorithms 4e; Brooks/Cole-Thomson: Belmont, CA, USA, 2004. [Google Scholar] [CrossRef]
  37. Yurkiewicz, J. Operations Research: Applications and Algorithms; Winston, W.L., Ed.; Duxbury Press: Boston, MA, USA, 1987; p. 1025. [Google Scholar]
  38. Ding, R.; Lin, D.K.J.; Wei, D. Dual-response surface optimization: A weighted MSE approach. Qual. Eng. 2004, 16, 377–385. [Google Scholar] [CrossRef]
  39. Vining, G.G.; Myers, R.H. Combining Taguchi and Response Surface Philosophies: A Dual Response Approach. J. Qual. Technol. 1990, 22, 38–45. [Google Scholar] [CrossRef]
  40. Fazelpour, F.; Soltani, N.; Rosen, M.A. Wind resource assessment and wind power potential for the city of Ardabil, Iran. Int. J. Energy Environ. Eng. 2015, 6, 431–438. [Google Scholar] [CrossRef] [Green Version]
  41. Lee, J.; Kim, D.R.; Lee, K.S. Optimum hub height of a wind turbine for maximizing annual net profit. Energy Convers. Manag. 2015, 100, 90–96. [Google Scholar] [CrossRef]
  42. Cheng, Z.; Wang, K.; Ong, M.C. Assessment of performance enhancement of a semi-submersible vertical axis wind turbine with an optimized Darrieus rotor. Eng. Struct. 2018, 167, 227–240. [Google Scholar] [CrossRef]
Figure 1. The map of the wind power plants’ locations in Turkey’s 6 out of 7 physiographic regions and 30 provinces.
Figure 1. The map of the wind power plants’ locations in Turkey’s 6 out of 7 physiographic regions and 30 provinces.
Sustainability 12 06017 g001
Figure 2. Parameters used for wind turbine structure.
Figure 2. Parameters used for wind turbine structure.
Sustainability 12 06017 g002
Figure 3. Flow chart of the proposed methodology.
Figure 3. Flow chart of the proposed methodology.
Sustainability 12 06017 g003
Figure 4. Box plots showing data distribution of factors: (a) Turbine Power; (b) Hub Height; (c) Rotor Diameter; (d) Wind Speed.
Figure 4. Box plots showing data distribution of factors: (a) Turbine Power; (b) Hub Height; (c) Rotor Diameter; (d) Wind Speed.
Sustainability 12 06017 g004aSustainability 12 06017 g004b
Figure 5. Main effects of input factors on the response variable.
Figure 5. Main effects of input factors on the response variable.
Sustainability 12 06017 g005
Table 1. Descriptive statistics of inputs and responses.
Table 1. Descriptive statistics of inputs and responses.
VariablesUnitNMeanSDMinq1q2 (Median)q3Max
FactorsTurbine PowerMW1122.62770.55250.80002.10002.75003.00003.6000
Hub Heightm11286.2949.52064.00080.00085.00091.000117.000
Rotor Diameterm112101.9713.9352.90090.00103.00113.00126.00
Wind speedm/s1126.80000.50755.40006.50006.80007.00007.7500
ResponseAnnual energy production
in 2018
million kWh1126.82842.42340.19285.11726.84628.719511.9324
Annual energy production
in 2017
million kWh1125.82193.09180.32883.45485.76958.236612.5206
Note 1: Unit: the measurement unit; MW: megawatts; m: meters; m/s: meter per second; kWh: kilowatt-hour. Note 2: N: the number of WPPs (wind power plants) operated; Mean: the average value, Min: the minimum value, Max: the maximum value: SD (standard deviation) is a criterion to calculate the dispersion or variation quantity belonging to a data set. In statistics, data are divided into three parts to form quartiles in order to understand the distribution of the data set; q1: 25% of the data set is less than or equal to this value; q2 (median): 50% of the data set is less than or equal to this value; q3: 75% of the data set is less than or equal to this value.
Table 2. Coded factors levels and limits.
Table 2. Coded factors levels and limits.
FactorsNotationsLevel 1Level 2Level 3
Coded Level −1Coded Level 0Coded Level 1
ValueLimitValueLimitValueLimit
Turbine PowerA0.80.8–2.12.22.11–3.03.63.1–3.6
Hub HeightB6464–8090.581–9111792–117
Rotor DiameterC52.9052.90–9089.4591–11312692–126
Wind SpeedD5.405.4–6.56.5756.6–77.757.1–7.75
Table 3. Box-Behnken design of experiment according to coded factor levels.
Table 3. Box-Behnken design of experiment according to coded factor levels.
Proposed
Types
of WPP
The Features of WPPs
Coded LevelsActual Levels
Factor 1
Unit 1
Factor 2
Unit 2
Factor 3
Unit 3
Factor 4
Unit 4
A
(MW)
B
(m)
C
(m)
D
(m/s)
WPP 1.0−1.00−1.000.000.000.864.089.456.575
WPP 2.01.00−1.000.000.003.664.089.456.575
WPP 3.0−1.001.000.000.000.8117.089.456.575
WPP 4.01.001.000.000.003.6117.089.456.575
WPP 5.00.000.00−1.00−1.002.290.552.905.400
WPP 6.00.000.001.00−1.002.290.5126.05.400
WPP 7.00.000.00−1.001.002.290.552.907.750
WPP 8.00.000.001.001.002.290.5126.07.750
WPP 9.0−1.000.000.00−1.000.890.589.455.400
WPP 10.01.000.000.00−1.003.690.589.455.400
WPP 11.0−1.000.000.001.000.890.589.457.750
WPP 12.01.000.000.001.003.690.589.457.750
WPP 13.00.00−1.00−1.000.002.264.052.906.575
WPP 14.00.001.00−1.000.002.2117.052.906.575
WPP 15.00.00−1.001.000.002.264.0126.06.575
WPP 16.00.001.001.000.002.2117.0126.06.575
WPP 17.0−1.000.00−1.000.000.890.552.906.575
WPP 18.01.000.00−1.000.003.690.552.906.575
WPP 19.0−1.000.001.000.000.890.5126.06.575
WPP 20.01.000.001.000.003.690.5126.06.575
WPP 21.00.00−1.000.00−1.002.264.089.455.400
WPP 22.00.001.000.00−1.002.2117.089.455.400
WPP 23.00.00−1.000.001.002.264.089.457.750
WPP 24.00.001.000.001.002.2117.089.457.750
WPP 25.00.000.000.000.002.290.589.456.575
WPP 26.00.000.000.000.002.290.589.456.575
WPP 27.00.000.000.000.002.290.589.456.575
Table 4. Analysis of variance for factors of wind power.
Table 4. Analysis of variance for factors of wind power.
SourceModel TypesF-Value (1.5017)t-Ratio (2.02)Prob>|t|
Prob>|F|
Status
Model31.16029.00<0.001Significant
ALinear237.8715.42<0.001Significant
B4.6400−2.150.037Significant
C62.0407.88<0.001Significant
D12.040−3.470.001Significant
AASquare37.440−6.12< 0.002Significant
BB3.6600−1.1910.063 *Significant *
CC4.9100−2.22<0.033Significant
DD34.880−5.91<0.001Significant
AB2-way Interaction1.59001.260.215Insignificant
AC2.78001.670.103Insignificant
AD24.0103.07<0.001Significant
BC9.44004.90<0.004Significant
BD5.8800−2.42<0.020Significant
CD19.790−4.45<0.001Significant
Note: p-value < 0.05, if the individual and interactive effects of the factors that could influence the response and are statistically significant; * at the margin of statistical provisionally significance (p-value < 0.1); Prob>|t| the p-value associated with a two-tailed test.
Table 5. Best ten solutions for Target0 ( y k W h ) based on desirability value.
Table 5. Best ten solutions for Target0 ( y k W h ) based on desirability value.
Feasible
Solutions
A
(MW)
B
(m)
C
(m)
D
(m/s)
Target0
(Million kWh)
Desirability
Values
1.03.06700108.842106.7596.168409.95520.95870
2.03.60000103.627125.3197.286359.83430.94530
3.03.2516177.4070116.9045.921589.52930.91058
4.01.68100116.695126.0005.400009.27830.88203
5.03.6000068.0000118.8677.332829.16640.86929
6.03.6000064.000056.09507.731378.89580.83849
7.02.1308874.5810126.0005.496828.69170.81527
8.02.2000090.500089.45006.575008.45550.78839
9.03.16695112.38381.28106.186598.34590.77593
10.02.2712364.383079.53406.478678.26120.76629
Table 6. The results of all optimization models.
Table 6. The results of all optimization models.
Input FactorsTarget0
y k W h
Target1
y ^ μ
Target2
y ^ σ
Optimization
Models
A
(MW)
B
(m)
C
(m)
D
(m/se)
Maximize
(Million kWh)
Maximize
(Million kWh)
Minimize
(Million kWh)
Composite
Desirability
(0.0 < CD> 1.0)
Target03.0670108.8424106.75976.16849.952ConstantConstant0.95872
Target12.9395101.119394.77606.4038-9.390Constant0.92460
Target21.7665106.292995.44186.8004-Constant0.246800.91440
Multi-Objective2.710094.5152112.08576.6581-9.6840.286300.91560
Average2.62075102.6925102.26586.507679.9529.5370.266550.92833
Table 7. The comparison of the decision variables’ results.
Table 7. The comparison of the decision variables’ results.
A
(MW)
B
(m)
C
(m)
D
(m/s)
Optimum 1Optimum 2Differences (%)
A (MW)0.8–3.6Fixed (86.294)Fixed (101.97)Fixed (6.80)2.7102.7230.4790%
B (m)Fixed (2.627)64–117Fixed (101.97)Fixed (6.80)94.515295.050.5658%
C (m)Fixed (2.627)Fixed (86.294)52.9–126Fixed (6.80)112.0857112.7090.5560%
D (m/s)Fixed (2.627)Fixed (86.294)Fixed (101.97)5.4–7.756.65816.61060.7100%
Optimum 1: The optimum results of the decision variables based on multi-objective functions’ optimization models. Optimum 2: The optimum results of the decision variables based on a reference case in which one constraint is optimized at a time and three other constrains are fixed.

Share and Cite

MDPI and ACS Style

Ayaz Atalan, Y.; Tayanç, M.; Erkan, K.; Atalan, A. Development of Nonlinear Optimization Models for Wind Power Plants Using Box-Behnken Design of Experiment: A Case Study for Turkey. Sustainability 2020, 12, 6017. https://doi.org/10.3390/su12156017

AMA Style

Ayaz Atalan Y, Tayanç M, Erkan K, Atalan A. Development of Nonlinear Optimization Models for Wind Power Plants Using Box-Behnken Design of Experiment: A Case Study for Turkey. Sustainability. 2020; 12(15):6017. https://doi.org/10.3390/su12156017

Chicago/Turabian Style

Ayaz Atalan, Yasemin, Mete Tayanç, Kamil Erkan, and Abdulkadir Atalan. 2020. "Development of Nonlinear Optimization Models for Wind Power Plants Using Box-Behnken Design of Experiment: A Case Study for Turkey" Sustainability 12, no. 15: 6017. https://doi.org/10.3390/su12156017

APA Style

Ayaz Atalan, Y., Tayanç, M., Erkan, K., & Atalan, A. (2020). Development of Nonlinear Optimization Models for Wind Power Plants Using Box-Behnken Design of Experiment: A Case Study for Turkey. Sustainability, 12(15), 6017. https://doi.org/10.3390/su12156017

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop