Next Article in Journal
Development Trends and Research Frontiers of Preferential Flow in Soil Based on CiteSpace
Next Article in Special Issue
Monitoring of Urban Black-Odor Water Using UAV Multispectral Data Based on Extreme Gradient Boosting
Previous Article in Journal
An Eco-Study for a Feasible Project: “Torun and Its Vistula Stretch—An Important Green Navigation Spot on a Blue Inland Waterway”
Previous Article in Special Issue
Research on Water-Level Recognition Method Based on Image Processing and Convolutional Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameters Tuning of Fractional-Order Proportional Integral Derivative in Water Turbine Governing System Using an Effective SDO with Enhanced Fitness-Distance Balance and Adaptive Local Search

1
School of Water Conservancy and Hydropower, Hebei University of Engineering, Handan 056038, China
2
Illinois State Water Survey, Prairie Research Institute, University of Illinois at Urbana-Champaign, Champaign, IL 61820, USA
*
Author to whom correspondence should be addressed.
Water 2022, 14(19), 3035; https://doi.org/10.3390/w14193035
Submission received: 2 August 2022 / Revised: 22 September 2022 / Accepted: 23 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Application of AI and UAV Techniques in Urban Water Science)

Abstract

:
Supply-demand-based optimization (SDO) is a swarm-based optimizer. However, it suffers from several drawbacks, such as lack of solution diversity and low convergence accuracy and search efficiency. To overcome them, an effective supply-demand-based optimization (ESDO) is proposed in this study. First, an enhanced fitness-distance balance (EFDB) and the Levy flight are introduced into the original version to avoid premature convergence and improve solution diversity; second, a mutation mechanism is integrated into the algorithm to improve search efficiency; finally, an adaptive local search strategy (ALS) is incorporated into the algorithm to enhance the convergence accuracy. The effect of the proposed method is verified based on the comparison of ESDO with several well-regarded algorithms using 23 benchmark functions. In addition, the ESDO algorithm is applied to tune the parameters of the fractional-order proportional integral derivative (FOPID) controller of the water turbine governor system. The comparative results reveal that ESDO is competitive and superior for solving real-world problems.

1. Introduction

With the development of social economy and technology, many complex optimization problems have appeared in the fields of communication, transportation, machinery, e-commerce, automation, materials and economics [1,2,3,4,5,6,7,8,9,10]. Meta-heuristic algorithms, an effective tool, simulate one or some natural processes in nature and have unique advantages in solving complex optimization problems. The No-Free−Lunch (NFL) theorem [11] proves that any optimizers cannot provide the best solutions for all different optimization problems. Therefore, plenty of meta-heuristics have sprung up in recent years according to different inspirations. Meta-heuristic algorithms can be divided into three categories: evolutionary-based (EB), physis-based (PB), and swarm-based (SB).
The most representative EB is the genetic algorithm (GA) [12], which simulates the evolution process of biological groups in nature. With the development of GA, many improved versions and variants have emerged, but most of them obtain high-quality solutions through mutation, crossover and selection steps. Other EBs are genetic programming (GP) [13], evolution strategies (ES) [14] and evolutionary programming (EP) [15].
PBs simulate the physical laws in nature. The annealing algorithm (SA) is a classical PB. The theoretical idea behind SA is the motion of molecules in a solid material as it cools gradually from high temperature. In recent years, numerous new physics-inspired algorithms have been proposed, including gravitational search algorithm (GSA) [16], atom search algorithm (ASO) [17], electromagnetic mechanism-like algorithm (EM) [18], central force optimization (CFO) [19], hysteretic optimization (HO) [20], charged system search (CSA) [21], Lichtenberg algorithm (LA) [22], henry gas solubility optimizer (HGSO) [23], optics inspired optimization (OIO) [24], and so on.
SBs simulate the collective behaviors of natural or artificially formed organizations. Particle swarm algorithm (PSO) [25], a well-known optimizer, has received a high degree of attention and widespread concern from scholars of many various fields. Since its emergence, other new SBs have been developed, including bacteria foraging optimization (BFO) [26], manta ray foraging optimization (MRFO) [27], virus colony search (VCS) [28], cuckoo search algorithm (CS) [29], artificial ecosystem-based optimization (AEO) [30], trees social relations optimization (TSR) [31], tuna swarm optimization (TSO) [32], black widow optimization algorithm (BWOA) [33], etc. Among three types of meta-heuristic algorithms, since SBs have an inherent advantage in exploration and exploitation, they have a larger family and are more competitive in solving complex optimization problems.
Inspired by supply-demand relation of commodities in economics, a swarm-based optimizer called supply-demand-based optimization (SDO) is proposed [34]. It simulates the cobweb model in market economy mathematically. Due to its simplicity and feasibility, it can achieve a balance between exploration and exploitation. Therefore, since the emergence of SDO, it has attracted a wide attention and has been applied in different fields. Ali et al. [35] used SDO for dielectric strength test and obtained the optimal value of SIR/TiO2 filler parameter percentage under different conditions. Ginidi et al. [36] extracted the electrical parameters of different photovoltaic models through SDO, and the comparisons with other algorithms showing that SDO is highly competitive. Alturki et al. [37] applied SDO to the selection problem of hybrid energy system, it can obtain lower loss of power supply probability (LPSP) and lower annualized system cost (ASC). Jing et al. [38] proposed an improved SDO algorithm for the prediction of foundation pit deformation, and the results demonstrated its feasibility and effectiveness. Ibrahim et al. [39] used SDO to estimate the parameters of induction motor and obtained a small deviation.
The operating state of the water turbine governor system is closely related to the stability of hydropower station operation and the safety of power consumption. Therefore, the parameter selection of the water turbine governor system becomes an important subject worthy of long-term study. At the beginning of the study, the parameters are obtained mainly through a large number of empirical formulas derived from field tests. However, the application value obtained by this method is not proportional to the time and financial consumption. The emergence of PID regulation law improves this situation, and is widely used because of its simple operation and good robustness. With the in-depth study of computer technology, it is found that traditional integer-order PID control law cannot be well modeled and analyzed for some complex non-linear engineering problems. The water turbine governor system has strong memory dependence, and its dynamic process can be better described by fractional-order calculus. The most intuitive difference between fractional-order PID (FOPID) controller and traditional integer-order PID controller is that there are two more adjustable parameters in FOPID, which can achieve better control effect. Therefore, the FOPID controller has strong application value.
The rise of intelligent optimization techniques provides a more effective method for parameter adjustment of FOPID controller. Vanchinathan et al. [40] applied the artificial bee colony algorithm to tuning for brushless DC motor. Oguzhan Karahan [41] used FOPID controller in core power control in molten salt reactors and adjusted the parameters by using the cuckoo algorithm. Munagala [42] performed robustness analysis on FOPID controller in an AVR system by using chaotic black widow algorithm. These studies demonstrated the superiority of FOPID in effectiveness of different control systems.
Although SDO as a swarm-based optimizer has better optimization ability than some other optimizers when handling optimization problems, it still has some big promotion space. SDO has the limitations of solution diversity, search efficiency and solution accuracy when tracking complex problems. Therefore, to overcome the drawbacks, this study designs an effective SDO termed ESDO; it incorporates the EFDB and Levy flight, mutation mechanism and adaptive local search strategy into SDO to perform global search. The effect of ESDO is verified based on the comparative results of ESDO with several well-regarded algorithms on a set of benchmark functions and tuning FOPID controller of water turbine governor. The results verify the efficiency and superiority of ESDO. The major difference of our study and its competitors is that ESDO makes three different strategies to improve its overall optimization performance in terms of solution diversity, convergence accuracy and search efficiency.
The rest of this study is as follows. Section 2 gives the main structure of SDO and provides the proposed ESDO by combining several strategies. In Section 3, the experimental results on some functions are investigated to assess the effectivity of the proposed ESDO. Section 4 provides an application of ESDO in tuning the FOPID controller of a water turbine governor. Section 5 gives some conclusions of the study.

2. Effective Supply-Demand-Based Optimization (ESDO)

2.1. Supply-Demand-Based Optimization (SDO)

In SDO [34], there are d commodity prices and quantities as candidate solutions and possible candidate solutions, respectively. After evaluation, the one with better fitness is selected as the current candidate solution. The mathematical expressions for commodity price and quantity are:
X = x 1 x 2 x n = x 1 1 x 1 2 x 1 d x 2 1 x 2 2 x 2 d x n 1 x n 2 x n d ,
Y = y 1 y 2 y n = y 1 1 y 1 2 y 1 d y 2 1 y 2 2 y 2 d y n 1 y n 2 y n d ,
where x i and y i     i = 1 n represent the price vectors and quantity vectors in each market, and x i j and   y i j i = 1 n ; j = 1 d represent the price and quantity of the j th commodity in the i th market, respectively. There are n markets in total.
Then, the fitness values of all prices and quantities are evaluated by the following functions:
F x = F x 1 F x 2 F x 3 F x n T ,
F y = F y 1 F y 2 F y 3 F y n T ,
where T represents transposition of the matrix.
The vector of the equilibrium price and quantity are represented by:
M i = F x i 1 n i = 1 n F x i P = M / i = 1 n M i x 0 = r 1 1 n i = 1 n x i r < 0.5 x k , k = R ( P ) r 0.5 ,
N i = F y i 1 n i = 1 n F y i Q = N / i = 1 n N i y 0 = y k , k = R ( Q ) ,
where r and r 1 are the random numbers in 0 , 1 , R P and R Q are the Roulette Wheel Selection.
The expressions of demand function and supply function are given by, respectively:
x i t + 1 = x 0 β y i t + 1 y 0 ,
y i t + 1 = y 0 + α x i t x 0 ,
where x i t and y i t respectively represent the i th price and quantity of commodity at time t , and β and α are respectively the demand weight and supply weight.
Substituting (7) into (8), the demand equation can be rewritten into the following form:
x i t + 1 = x 0 α β x i t y 0 ,
Therefore, the commodity price is updated by adjusting the values of α and β , and it is updated according to the equilibrium price relative to the current price. α and β can be expressed as:
α = 2 T t + 1 T s i n 2 π r ,
β = 2 c o s 2 π r ,
where, α β < 1 is equivalent to the stability mode in the supply and demand mechanism, emphasizing the exploitation ability, α β > 1 is equivalent to the instability mode and emphasizing the overall exploration ability. Figure 1 shows the two modes of SDO.

2.2. Proposed Method

Since SDO easily suffers from low search efficiency and misses some better solutions, ESDO is an enhanced version of SDO in order to overcome them. For SDO, the fitness value of each solution and the distance between each solution and the best solution are very two important factors affecting the search efficiency of the algorithm for the optimal solution. In this study, inspired by a fitness-distance balance (FDB) [43], an enhanced fitness-distance balance (EFDB) is developed to replace the selection for the equilibrium quantity and price of commodities. Meanwhile, to strengthen exploration ability of the algorithm in the search space, a Levy flight strategy is introduced to the weight to improve the convergence ability. A mutation mechanism is employed to enhance the search efficiency of SDO. Meanwhile, an adaptive local search strategy is used to improve the convergence accuracy of the algorithm. ESDO is a newly presented optimizer and not yet applied in any real-world application.

2.2.1. Combining Enhanced Fitness-Distance Balance (EFDB) and Levy Flight

In [44], Kati et al. proposed an improved version of SDO, in which the equilibrium quantity is replaced with the commodity quantity selected by the FDB method to provide diversity. In this study, according to [44], we propose the EFDB method, in which the equilibrium quantity and price of commodities are replaced with the commodity quantity and price by using the FDB method, respectively. This method can further strengthen the solution diversity. The following is a specific description to the EFDB method. The fitness value of commodity price is calculated by [43]:
i = 1 N F x i = if   goal   is   minimization   1 - n o r m G i if   goal   is   maximization   n o r m G i
F x = [ f x 1 , , f x N ] 1 × N
i = 1 n , i b e s t , D i = ( x i 1 x b e s t 1 ) 2 + + ( x i d x b e s t d ) 2
D x = [ d x 1 , , d x n ] 1 × n
i = 1 n , M i = 0.5 ( 1 + I t / M a x I t ) × n o r m F x i + 0.5 ( 1 I t / M a x I t ) × n o r m D x i
The equilibrium price vector x0 is redefined by:
P = M i i = 1 n M i
As an analogy, the equilibrium quantity vector y0 is redefined by:
Q = N i i = 1 n N i
where
i = 1 n , M i = 0.5 ( 1 + I t / M a x I t ) × n o r m F y i + 0.5 ( 1 I t / M a x I t ) × n o r m D y i
where normFyi is the normalized fitness value of Fyi, normDyi is the normalized distance value of the ith commodity, It is the current iteration and MaxIt is the maximum number of iterations. This selection method combines the fitness and the distance to calculate the score for each individual in the population. Therefore, this strategy can effectively improve solution diversity and avoid local solutions. In addition, inspecting Equation (16), in the early iterations, the selection for the equilibrium quantity and price of commodities takes into account the largest fitness value and the furthest distance from the current optimal individual so far, it will contribute to exploration; in the later iterations, the selection focuses more on the fitness value rather than the distance, it will be dedicated to exploitation.
The Levy flight [45,46] whose step-width obeys non-uniform Levy distribution, is a random walk; thus, it has the superior ability to enhance exploring space search [47]. The step-width of the Levy flight is produced by [48,49]:
L(s) ~ tλ (1≤ λ ≤ 3)
s = u σ u | v | 1 b
where λ is a stability/tail index, s is the step-width and u and v obey the normal distribution, respectively:
u ~ N(0,1), v ~ N(0,1),
σ u = Γ ( 1 + b ) sin ( π b 2 ) Γ ( 1 + b 2 ) β 2 b 1 2 1 b
where Γ denotes the standard Gamma function and b = 1.5.
The weight α in Equation (10) is reformulated by:
α = 5 u σ u 2 | v | 1 β ( T t + 1 ) T sin ( 2 π r )

2.2.2. The Mutation Mechanism

To improve the search efficiency of the algorithm, the mutation mechanism is employed in this study. Although some different mutation strategies are introduced in the literature [50,51,52], the Gaussian mutation is one of the most frequently used mutation methods since it is more effective and simpler to implement [53]. So, in ESDO, the supply function and demand function are modified by, respectively:
y i ( t + 1 ) = y 0 + α ( x i ( t ) x 0 ) + r o u n d ( 0.5 ( 0.01 + r a n d ) ) r n
x i ( t + 1 ) = x 0 β ( y i ( t + 1 ) y 0 ) + r o u n d ( 0.5 ( 0.01 + r a n d ) ) r n
where
rn ~ N(0,1)

2.2.3. Adaptive Local Search (ALS) Strategy

Local search strategy is an important way to improve the current best solution. The chaotic local search is a classic local search (CLS) method [54], which used the chaotic map to improve the solution quality by searching the neighborhood around the best solution so far. However, the step-size for this local search cannot decrease as iterations go on, which will affect the solution accuracy and search efficiency. So, to dynamically adjust the step-size of the local search, an adaptive local search (ALS) strategy is given by:
x b e s t ' ( t ) = x best ( t ) + s t e p ( 2 r a n d 1 ) ( U b L b )
s t e p = e 20 t T
x b e s t ' ( t ) is a new current solution generated at time t.
The update of the current best solution is given by:
x b e s t ( t + 1 ) = x b e s t ' ( t )   if   f i t ( x b e s t ' ( t ) ) < f i t ( x b e s t ( t ) ) x b e s t ( t + 1 ) = x best ( t )   else ,
If the fitness value of x b e s t ' ( t ) is better than that of x b e s t ( t ) , the current best solution is replaced with the new one, or it remains unchanged. It can be observed from Equation (28) that a bigger step-size contributes to exploration in the early iterations, with the increase of iterations and, a small step-size is greatly dedicated to exploitation. In addition to improving the convergence accuracy, this adaptive search strategy also strengthens the balance between exploration and exploitation to some extent.

2.2.4. The Proposed ESDO Algorithm

By introducing the EFDB method and Levy flight, the mutation mechanism, and the ALS strategy to strengthen the optimization performance, the ESDO algorithm is proposed. The pseudocode of ESDO algorithm is given in Figure 2.

3. Experimental Results and Analysis

3.1. Test Functions and Parameter Setting

To assess the performance of the ESDO algorithm, a classical suit of benchmark set, containing 23 test functions (see Table A1 in Appendix A for details), are employed. The 23 benchmark functions include 7 unimodal functions (UFs) (F1–F7), 6 multimodal functions (MFs) (F8–F13) and 10 low-dimensional multimodal functions (LMFs) (F14–F23). Meanwhile, several competitive optimizers, including whale optimization algorithm (WOA) [55], gray wolf optimizer (GWO) [56], and GSA, are used and their results are provided for a comparison. For all the considered optimizers, the population size and the maximum number of iterations are set to 50 and 500, respectively. The ESDO algorithm is firstly analyzed qualitatively based on exploration and exploitation. Then, the Wilcoxon signed-rank test and Friedman test are statistically analyzed, respectively. The experimental results are based on 30 independent runs. The other parameter settings of all the considered optimizers are described in Table 1.

3.2. Exploitation Analysis

The functions F1–F7 having only one extremum are used to assess exploitation of algorithms. The comparison results of the algorithms on these UNs are listed in Table 2, in which the underline indicates the best value among all the algorithms. From this table, ESDO provides better solutions on functions F1–F5 and F7 in terms of mean and Std. ESDO performs as well as SDO and GWO on function F6. Therefore, ESDO obtains better results on most of UNs. Figure 3 shows the convergence curves of the algorithms on F1–F7. ESDO exhibits superior convergence rate over other optimizers in exploiting the optimal solution.

3.3. Exploration Analysis

The functions F8–F23 having multiple extrema are used to evaluate exploration of algorithms. The comparison results of the algorithms on these functions are shown in Table 3 and Table 4. Inspecting Table 3, ESDO obviously outperforms other optimizers on function F12 and F13. ESDO, SDO and WOA provide the same results in terms of mean on functions F9 and F10. In addition, ESDO is third best algorithm on function F7. The performance of ESDO ranks only second to WOA on function F8. Observing Table 4, for F15 and F20, ESDO obtains the second-best results, which are only inferior to those of SDO. For functions F16–F19 and F22–F23, ESDO, SDO and one of other algorithms offer the same best results. For F21, ESDO provides the best results. The convergence curves of the algorithms on F8–F23 are depicted in Figure 4 and Figure 5, which manifests that ESDO has better convergence performance with comparison to other algorithms when tackling these test functions.

3.4. Statistical Analysis

To evaluate the overall performance of ESDO and rank it statistically, the Wilcoxon signed-rank test (WSRT) [57] and Friedman test (FT) [58] are employed, respectively. Table 5 and Table 6 provide the results of WSRT. In the table, “=” indicates there is no significant difference between ESDO and its competitor for a considered problem, “+” indicates that the performance of ESDO is better than that of a competitor for a considered problem and “−” is the opposite. Table 7 summarizes the results of WSRT in Table 5 and Table 6. In Table 7, The WSRT results reveal that ESDO is superior to SDO, GWO, WOA and GSA on functions 7, 18, 22 and 17 out of 23, respectively, indicating the superior performance of ESDO to its competitors statistically. Figure 6 depicts the rank of each function for the comparative optimizers, Figure 7 gives the mean of these ranks. From Figure 7, ESDO ranks the first among the considered optimizers, demonstrating that it exhibits the best optimization ability compared to its counterparts.

4. ESDO for FOPID Controller of Water Turbine Governor System

4.1. System Description

The traditional integer-order PID controller has only three parameters ( K p , K i , K d ) [59], while the FOPID controller used in this experiment has five parameters ( K p , K i , K d , λ , μ ), so as to achieve more accurate control. Theoretically, λ and μ can take any number. In the case of λ = μ = 1 , FOPID is converted to integer-order PID.
The transfer function of FOPID is described by:
G s ( s ) = K p + K i s λ + K d s μ T v s + 1 ,
where K p , K i and K d represent the proportional coefficient, integral coefficient and differential coefficient, respectively, s represents the Laplace operator, and λ and μ represent the exponents of the integral operator and the differential operator, respectively. T v represents the differential time constant.
The variation of the unit load produces a deviation e. The FOPID controller will convert this difference into an adjustment signal. After receiving the signal, the mechanical hydraulic system will operate to adjust the opening of the guide vane, and then adjust the flow and restore the speed. Figure 8 is the FOPID model of the water turbine governor system.

4.2. Experimental Results and Analysis

This study is simulated by MATLAB/Simulink. The fitness function uses the time multiplication error integration criterion and is defined by:
f ITAE = 0 t s t e ( t ) d t ,
where e ( t ) is the deviation of the actual output from the expected output.
The specific experimental settings and parameter values are shown in Table 8.
In this experiment, in order to be more practical, four different working conditions under 0–20% load are adopted, and the controller performance from ESDO is compared with that from SDO, WOA and GSA under each working condition. Each algorithm runs 20 times, and the average optimal fitness value, overshoot and adjustment time based on 20 experiments are compared to verify the parameter setting effect of ESDO on the FOPID controller. Under different load conditions, the overshoot is represented by the maximum value. The simulation results are shown in Table 9. The convergence curve of average fitness, convergence curves of average FOPID parameters and the speed response curves obtained of the average FOPID parameters under different load conditions are shown in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20, respectively.
Under 4% load, the fitness value of the FOPID controller offered by ESDO is the smallest among the five algorithms, which shows that ESDO can jump out of local optimum well and find the global optimum; meanwhile, it effectively improves the optimization accuracy. From the convergence curves under 4% load in Figure 9, the convergence rate of ESDO is also faster than other algorithms, which indicates that ESDO has the strongest ability to find the best solution among these algorithms. The overshoot and adjustment time are two important indicators to measure whether a control system is stable. It can be seen from Table 9 that under the condition of 4% load, the overshoot and adjustment time obtained are the minimum after the governor of the water turbine is tuned by ESDO. The trend of the response curves in Figure 11 shows that after the peak value, the oscillation amplitude of the curve obtained by ESDO is the smallest, and it can reach a stable value quickly. On the contrary, the oscillation amplitude of the response curves obtained by SDO, GWO, WOA and GSA are very large, and the state is very unstable. It proves that the dynamic regulation ability of the FOPID controller offered by ESDO is stronger than those offered by the other four algorithms, and the turning parameters of the turbine governor with ESDO can make the control system more stable and efficient.
Under 8% load, the fitness value of the FOPID controller offered by ESDO is also the smallest, showing that the global exploration ability and local exploitation ability of ESDO are significantly enhanced. The convergence curves of average fitness in Figure 12 demonstrate that, under the 8% load condition, ESDO can quickly find the optimal value in the optimization process, greatly improving the efficiency of the algorithm. After the governor of the water turbine is tuned by different algorithms, the overshoot of the FOPID controller offered by ESDO is only second to that offered by GSA and very close to that offered by GSA, but the adjustment time of the FOPID controller offered by ESDO is far better than that offered by GSA. It can be found from Figure 14 that the FOPID controller from ESDO also has the smallest oscillation amplitude and the smallest duration, the FOPID controller from GSA has the largest oscillation amplitude and the longest duration, and the FOPID controllers from SDO, GWO and WOA are in between. Considering comprehensively, under 8% load condition, ESDO also has the strongest optimization ability for the FOPID controller among the five algorithms.
When the load is 12%, the FOPID controller offered by ESDO also gets the minimum fitness value. From Figure 15, although the convergence speed of ESDO is slower than that of SDO at the beginning of the iteration, in the middle of the iteration, ESDO outperforms SDO. At this time, the convergence speed of ESDO is the fastest, and it tends to jump out of the local extreme value to find a smaller value. Moreover, in the regulation system, the overshoot of the FOPID controller offered by ESDO is relatively good, and the adjustment time is also the smallest. From the speed response curves in Figure 17, after the first peak, the ESDO curve is very stable with little fluctuation, while the response curves of the FOPID controller offered by other algorithms fluctuate greatly and last a long time.
When the load is 16%, the fitness value of the FOPID controller offered by ESDO is still the smallest. From Figure 18, there is a fast convergence speed at the beginning of iteration, and its ability to jump out of local extreme value is also the strongest among all algorithms. According to the speed response curves in Figure 20, the overshoot of of the FOPID controller offered by ESDO is very close to that offered by GSA, which is better than those offered by other algorithms. The adjustment time of the FOPID controller offered by ESDO is the shortest among all algorithms, after falling from the peak, the fluctuation weakens rapidly and finally tends to be stable.
Therefore, through tuning the FOPID parameters of the water turbine governor under different load conditions, these convincing results reveal that ESDO is excellent in tackling real-world engineering applications.

5. Conclusions

To better handle optimization problems, an effective supply-demand-based optimization (ESDO) is proposed, it combines three strategies, including the enhanced fitness-distance balance (EFDB) with Levy flight, mutation mechanism and adaptive local search, to improve solution diversity, and convergence accuracy and search efficiency. The experimental results of ESDO with several well-regarded algorithms on 23 benchmark functions discover that ESDO has superior optimization ability. In addition, the practicability of ESDO is also verified by tuning the parameters of the FOPID controller of the water turbine governor system. The experimental results show that the turbine governor system tuned by ESDO is better than those tuned by other algorithms in terms of response time and overshoot.
In this study, some secondary and high-level terms are ignored when establishing the hydraulic turbine simulation model, and higher-order models can be established in the future. The implementation method of FOPID can also be improved.
In future studies, ESDO could be specifically intended to handle a variety of problems in hydraulic engineering, such as optimal allocation of power station [60], intelligent fault diagnosis of turbine [61] and optimal design of gate [62,63,64].

Author Contributions

Conceptualization, W.Z. and L.W.; methodology, H.Z. and W.Z.; software, Z.Z.; validation, Z.Z., K.Z. and H.Z.; formal analysis, H.Z.; investigation, writing—original draft preparation, W.Z and H.Z.; writing—review and editing, H.Z.; visualization, L.W.; supervision, L.W.; project administration L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Unimodal test functions.
Table A1. Unimodal test functions.
NameFunctionDRange f opt
Sphere f 1 ( x ) = i = 1 n x i 2 30 100 , 100 D 0
Schwefel2.22 f 2 ( x ) = i = 1 n x i + i = 1 n x i 30 10 , 10 D 0
Schwefel1.2 f 3 ( x ) = i = 1 n ( j = 1 i x j ) 2 30 100 , 100 D 0
Schwefel2.21 f 4 ( x ) = max i x i , 1 i n 30 100 , 100 D 0
Rosenbrock f 5 ( x ) = i = 1 n 1 ( 100 ( x i + 1 x i ) 2 + ( x i 1 ) 2 ) 30 30 , 30 D 0
Step f 6 ( x ) = i = 1 n ( x i + 0.5 ) 2 30 100 , 100 D 0
Quartic f 7 ( x ) = i = 1 n i x 4 + r a n d o m [ 0 , 1 ) 30 1.28 , 1.28 D 0
Table A2. Multimodal test functions.
Table A2. Multimodal test functions.
NameFunctionDRange f opt
Schwefel f 8 ( x ) = i = 1 n ( x i sin ( x i ) ) 30 500 , 500 D 12 , 569.5
Rastrigin f 9 ( x ) = i = 1 n ( x i 2 10 cos ( 2 π x i ) + 10 ) 2 30 5.12 , 5.12 D 0
Ackley f 10 ( x ) = 20 exp ( 0.2 1 n i = 1 n x i 2 ) exp ( 1 n i = 1 n cos 2 π x i ) + 20 + e 30 32 , 32 D 0
Griewank f 11 ( x ) = 1 4000 i = 1 n x i 100 2 i = 1 n cos ( x i 100 i ) + 1 30 600 , 600 D 0
Penalized f 12 ( x ) = π n { 10 sin 2 ( π y 1 ) + i = 1 n = 1 ( y i 1 )   2 [ 1 + 10 sin 2 ( π y i + 1 ) ] + ( y n 1 )   2 } + i = 1 30 u ( x i , 10 , 100 , 4 ) 30 50 , 50 D 0
Penalized2 f 13 ( x ) = 0.1 { sin 2 ( 3 π x 1 ) + i = 1 29 ( x i 1 ) 2 p [ 1 + sin 2 ( 3 π x i + 1 ) ] + ( x n 1 )   2 [ 1 + sin 2 ( 2 π x 30 ) ] } + i = 1 30 u ( x i , 5 , 10 , 4 ) 30 50 , 50 D 0
Table A3. Low-dimension Multimodal test functions.
Table A3. Low-dimension Multimodal test functions.
NameFunctionDRange f opt
Foxholes f 14 ( x ) = [ 1 500 + j = 1 25 1 j + j = 1 2 ( x i a ij ) 6 ] 1 2 65.536 , 65.536 D 0.998
Kowalik f 15 ( x ) = i = 1 11 a i x 1 ( b i 2 + b i x 2 ) b i 2 + b i x 3 + x 4 2 4 5 , 5 D 3.075 × 10 4
Six Hump Camel f 16 ( x ) = 4 x 1 2 2.1 x 1 4 + 1 3 x 1 6 + x 1 x 2 4 x 2 2 + 4 x 2 4 2 5 , 5 D 1.036
Branin f 17 ( x ) = ( x 2 5.1 4 π 2 x 1 2 + 5 π x 1 6 ) 2 + 10 ( 1 1 8 π ) cos x 1 + 10 2 5 , 10 × 0 , 15 0.398
GoldStein-Price f 18 ( x ) = [ 1 + ( x 1 + x 2 + 1 )   2 ( 19 14 x 1 + 3 x 1 2 14 x 2 + 6 x 1 x 2 + 3 x 2 2 ) ] × [ 30 + ( 2 x 1 + 1 3 x 2 )   2 ( 18 32 x 1 + 12 x 1 2 + 48 x 2 36 x 1 x 2 + 27 x 2 2 ) ] 2 2 , 2 D 3
Hartman 3 f 19 ( x ) = i = 1 4 exp [ j = 1 3 a ij ( x j p ij ) 2 ] 3 0 , 1 D 3.86
Hartman 6 f 20 ( x ) = i = 1 4 exp [ j = 1 6 a ij ( x j p ij ) 2 ] 6 0 , 1 D 3.322
Shekel 5 f 21 ( x ) = i = 1 5 ( x i a i ) ( x i a i ) T + c i 1 4 0 , 10 D 10.1532
Shekel 7 f 22 ( x ) = i = 1 7 ( x i a i ) ( x i a i ) T + c i 1 4 0 , 10 D 10.4028
Shekel 10 f 23 ( x ) = i = 1 10 ( x i a i ) ( x i a i ) T + c i 1 4 0 , 10 D 10.5364

References

  1. Lu, R.R.; Wang, J.Y.; Fu, X.T.; Lin, S.H.; Wang, Q.; Zhang, B. Performance analysis and optimization for UAV-based FSO communication systems. Phys. Commun. 2022, 51, 101594. [Google Scholar] [CrossRef]
  2. Beus, M.; Krpan, M.; Kuzle, I.; Pandžić, H.; Parisio, A. A model predictive control approach to operation optimization of an ultracapacitor bank for frequency control. IEEE Trans. Energy Convers. 2021, 36, 1743–1755. [Google Scholar] [CrossRef]
  3. Chang, C.; Han, M. Production scheduling optimization of prefabricated building components based on dde algorithm. Math. Probl. Eng. 2021, 2021, 6672753. [Google Scholar] [CrossRef]
  4. Hu, G.; Du, B.; Wang, X.; Wei, G. An enhanced black widow optimization algorithm for feature selection. Knowl.-Based Syst. 2022, 235, 107638. [Google Scholar] [CrossRef]
  5. Yoo, S.; Oh, S. Flow analysis and optimization of a vertical axis wind turbine blade with a dimple. Eng. Appl. Comp. Fluid Mech. 2021, 15, 1666–1681. [Google Scholar] [CrossRef]
  6. Wang, L.; Cao, Q.; Zhang, Z.; Mirjalili, S.; Zhao, W. Artificial rabbits optimization: A new bio-inspired meta-heuristic algorithm for solving engineering optimization problems. Eng. Appl. Artif. Intell. 2022, 114, 105082. [Google Scholar] [CrossRef]
  7. Yildiz, A.R.; Mehta, P. Manta ray foraging optimization algorithm and hybrid Taguchi salp swarm-Nelder–Mead algorithm for the structural design of engineering components. Mater. Test. 2020, 64, 706–713. [Google Scholar] [CrossRef]
  8. Zeinalzadeh, A.; Pakatchian, M.R. Evaluation of novel-objective functions in the design optimization of a transonic rotor by using deep learning. Eng. Appl. Comp. Fluid Mech. 2021, 15, 561–583. [Google Scholar] [CrossRef]
  9. Pei, H.; Cui, Y.; Kong, B.; Jiang, Y.; Shi, H. Structural parameters optimization of submerged inlet using least squares support vector machines and improved genetic algorithm-particle swarm optimization approach. Eng. Appl. Comp. Fluid Mech. 2021, 15, 503–511. [Google Scholar] [CrossRef]
  10. Hu, G.; Zhong, J.; Du, B.; Wei, G. An enhanced hybrid arithmetic optimization algorithm for engineering applications. Comput. Meth. Appl. Mech. Eng. 2022, 394, 114901. [Google Scholar] [CrossRef]
  11. Wolpert, D.H.; Macready, W.G. No free lunch theorems for optimization. IEEE Trans. Evol. Comput. 1997, 1, 67–82. [Google Scholar] [CrossRef]
  12. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  13. Koza, J.R. Genetic programming as a means for programming computers by natural selection. Stat. Comput. 1994, 4, 87–112. [Google Scholar] [CrossRef]
  14. Huning, A. Evolutionsstrategie. optimierung technischer systeme nach prinzipien der biologischen evolution. Arch. Philos. Law Soc. Philos. 1976, 62, 298–300. [Google Scholar]
  15. David, B.F. Artificial Intelligence through Simulated Evolution. In Evolutionary Computation: The Fossil Record; Wiley-IEEE Press: Hoboken, NJ, USA, 1998; pp. 227–296. [Google Scholar] [CrossRef]
  16. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A gravitational search algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  17. Zhao, W.; Wang, L.; Zhang, Z. Atom search optimization and its application to solve a hydrogeologic parameter estimation problem. Knowl.-Based Syst. 2019, 163, 283–304. [Google Scholar] [CrossRef]
  18. Xing, B.; Gao, W.J. Electromagnetism-like Mechanism Algorithm. In Innovative Computational Intelligence: A Rough Guide to 134 Clever Algorithms; Intelligent Systems Reference Library; Springer: Cham, Switzerland, 2014; Volume 62, pp. 347–354. [Google Scholar] [CrossRef]
  19. Formato, R.A. Central force optimization: A new metaheuristic with applications in applied electromagnetics. Prog. Electromagn. Res. 2007, 77, 425–491. [Google Scholar] [CrossRef]
  20. Pál, K.F. Hysteretic optimization for the Sherrington-Kirkpatrick spin glass. Phys. A Stat. Mech. Its Appl. 2006, 367, 261–268. [Google Scholar] [CrossRef]
  21. Kaveh, A.; Talatahari, S. A novel heuristic optimization method: Charged system search. Acta Mech. 2010, 213, 267–289. [Google Scholar] [CrossRef]
  22. Pereira, J.L.J.; Francisco, M.B.; Diniz, C.A.; Oliver, G.A.; Cunha, S.S., Jr.; Gomes, G.F. Lichtenberg algorithm: A novel hybrid physics-based meta-heuristic for global optimization. Expert Syst. Appl. 2021, 170, 114522. [Google Scholar] [CrossRef]
  23. Hashim, F.A.; Houssein, E.H.; Mabrouk, M.S.; Al-Atabany, W.; Mirjalili, S. Henry gas solubility optimization: A novel physics-based algorithm. Futur. Gener. Comp. Syst. 2019, 101, 646–667. [Google Scholar] [CrossRef]
  24. Kashan, A.H. A new metaheuristic for optimization: Optics inspired optimization (OIO). Comput. Oper. Res. 2015, 55, 99–125. [Google Scholar] [CrossRef]
  25. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science, IEEE, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar] [CrossRef]
  26. Passino, K.M. Biomimicry of bacterial foraging for distributed optimization and control. IEEE Control Syst. Mag. 2002, 22, 52–67. [Google Scholar] [CrossRef]
  27. Zhao, W.; Zhang, Z.; Wang, L. Manta ray foraging optimization: An effective bio-inspired optimizer for engineering applications. Eng. Appl. Artif. Intell. 2020, 87, 103300. [Google Scholar] [CrossRef]
  28. Li, M.D.; Zhao, H.; Weng, X.W.; Han, T. A novel nature-inspired algorithm for optimization: Virus colony search. Adv. Eng. Softw. 2016, 92, 65–88. [Google Scholar] [CrossRef]
  29. Yang, X.S.; Deb, S. Cuckoo search via Levy flights. In Proceedings of the World Congress on Nature & Biologically Inspired Computing (NaBIC 2009), Coimbatore, India, 9–11 December 2009; IEEE Publications: Piscataway, NJ, USA; pp. 210–214. [Google Scholar] [CrossRef]
  30. Zhao, W.; Wang, L.; Zhang, Z. Artificial ecosystem-based optimization: A novel nature-inspired meta-heuristic algorithm. Neural Comput. Appl. 2020, 32, 9383–9425. [Google Scholar] [CrossRef]
  31. Alimoradi, M.; Azgomi, H.; Asghari, A. Trees social relations optimization algorithm: A new Swarm-Based metaheuristic technique to solve continuous and discrete optimization problems. Math. Comput. Simul. 2022, 194, 629–664. [Google Scholar] [CrossRef]
  32. Xie, L.; Han, T.; Zhou, H.; Zhang, Z.R.; Han, B.; Tang, A. Tuna swarm optimization: A novel swarm-based metaheuristic algorithm for global optimization. Comput. Intell. Neurosci. 2021, 2021, 9210050. [Google Scholar] [CrossRef]
  33. Peña-Delgado, A.F.; Peraza-Vázquez, H.; Almazán-Covarrubias, J.H.; Torres Cruz, N.; García-Vite, P.M.; Morales-Cepeda, A.B.; Ramirez-Arredondo, J.M. A novel bio-inspired algorithm applied to selective harmonic elimination in a three-phase eleven-level inverter. Math. Probl. Eng. 2020, 2020, 8856040. [Google Scholar] [CrossRef]
  34. Zhao, W.; Wang, L.; Zhang, Z. Supply-Demand-Based Optimization: A Novel Economics-Inspired Algorithm for Global Optimization. IEEE Access 2019, 7, 73182–73206. [Google Scholar] [CrossRef]
  35. Ali, A.M.; Nasrat, L.; Hassan, M.H.; Elsayed, S.K.; Kamel, S. Supply Demand-Based Optimization Algorithm for Estimating Break Down Voltage of Silicon Rubber Insulators. In Proceedings of the 2021 IEEE CHILEAN Conference on Electrical, Electronics Engineering, Information and Communication Technologies (CHILECON), Valparaíso, Chile, 6–9 December 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1–9. [Google Scholar] [CrossRef]
  36. Ginidi, A.R.; Shaheen, A.M.; El-Sehiemy, R.A.; Elattar, E. Supply demand optimization algorithm for parameter extraction of various solar cell models. Energy Rep. 2021, 7, 5772–5794. [Google Scholar] [CrossRef]
  37. Alturki, F.A.; Al-Shamma’a, A.A.; Farh, H.M.; AlSharabi, K. Optimal sizing of autonomous hybrid energy system using supply-demand-based optimization algorithm. Int. J. Energy Res. 2021, 45, 605–625. [Google Scholar] [CrossRef]
  38. Jing, C.; Wang, H.; Li, H. Deformation Prediction of Foundation Pit Based on Exponential Power Product Model of Improved Algorithm. Geofluids 2021, 2021, 7055693. [Google Scholar] [CrossRef]
  39. Ibrahim, S.A.; Kamel, S.; Hassan, M.H.; Elsayed, S.K.; Nasrat, L. Developed Algorithm Based on Supply-Demand-Based Optimizer for Parameters Estimation of Induction Motor. In Proceedings of the 2021 IEEE International Conference on Automation/XXIV Congress of the Chilean Association of Automatic Control (ICA-ACCA), Valparaíso, Chile, 22–26 March 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1–6. [Google Scholar]
  40. Vanchinathan, K.; Selvaganesan, N. Adaptive fractional order PID controller tuning for brushless DC motor using artificial bee colony algorithm. Results Control Optim. 2021, 4, 100032. [Google Scholar] [CrossRef]
  41. Karahan, O. Design of optimal fractional order fuzzy PID controller based on cuckoo search algorithm for core power control in molten salt reactors. Prog. Nucl. Energy 2021, 139, 103868. [Google Scholar] [CrossRef]
  42. Munagala, V.K.; Jatoth, R.K. Improved fractional PIλDμ controller for AVR system using Chaotic Black Widow algorithm. Comput. Electr. Eng. 2022, 97, 107600. [Google Scholar] [CrossRef]
  43. Kahraman, H.T.; Aras, S.; Gedikli, E. Fitness-distance balance (FDB): A new selection method for meta-heuristic search algorithms. Knowl.-Based Syst. 2020, 190, 105169. [Google Scholar] [CrossRef]
  44. Kati, M.; Kahraman, H.T. Improving supply-demand-based optimization algorithm with fdb method: A comprehensive research on engineering design problems. Mühendislik Bilimleri Ve Tasarım Derg. 2020, 8, 156–172. [Google Scholar] [CrossRef]
  45. Hu, G.; Chen, L.; Wang, X.; Wei, G. Differential Evolution-Boosted Sine Cosine Golden Eagle Optimizer with Lévy Flight. J. Bionic Eng. 2022, 2022, 1–36. [Google Scholar] [CrossRef]
  46. Lai, V.; Huang, Y.F.; Koo, C.H.; Ahmed, A.N.; El-Shafie, A. Optimization of reservoir operation at Klang Gate Dam utilizing a whale optimization algorithm and a Lévy flight and distribution enhancement technique. Eng. Appl. Comp. Fluid Mech. 2021, 15, 1682–1702. [Google Scholar] [CrossRef]
  47. Liu, J.; Shi, J.; Hao, F.; Dai, M. A novel enhanced global exploration whale optimization algorithm based on Lévy flights and judgment mechanism for global continuous optimization problems. Eng. Comput. 2022, 2022, 1–29. [Google Scholar] [CrossRef]
  48. Mantegna, R.N. Fast, accurate algorithm for numerical simulation of Levy stable stochastic processes. Phys. Rev. E 1994, 49, 4677–4683. [Google Scholar] [CrossRef] [PubMed]
  49. Viswanathan, G.M.; Afanasyev, V.; Buldyrev, S.V.; Murphy, E.J.; Prince, P.A.; Stanley, H.E. Lévy flight search patterns of wandering albatrosses. Nature 1996, 381, 413–415. [Google Scholar] [CrossRef]
  50. Higashi, N.; Iba, H. Particle swarm optimization with Gaussian mutation. In Proceedings of the 2003 IEEE Swarm Intelligence Symposium. SIS’03 (Cat. No. 03EX706), Indianapolis, IN, USA, 26 April 2003; pp. 72–79. [Google Scholar] [CrossRef]
  51. Gharehchopogh, F.S. An Improved Tunicate Swarm Algorithm with Best-random Mutation Strategy for Global Optimization Problems. J. Bionic Eng. 2022, 19, 1177–1202. [Google Scholar] [CrossRef]
  52. Alireza, A.L.F.I. PSO with adaptive mutation and inertia weight and its application in parameter estimation of dynamic systems. Acta Autom. Sin. 2011, 37, 541–549. [Google Scholar] [CrossRef]
  53. Lu, H.; Sriyanyong, P.; Song, Y.H.; Dillon, T. Experimental study of a new hybrid PSO with mutation for economic dispatch with non-smooth cost function. Int. J. Electr. Power Energy Syst. 2010, 32, 921–935. [Google Scholar] [CrossRef]
  54. Tang, Z.; Tao, S.; Wang, K.; Lu, B.; Todo, Y.; Gao, S. Chaotic Wind Driven Optimization with Fitness Distance Balance Strategy. Int. J. Comput. Intell. Syst. 2022, 15, 46. [Google Scholar] [CrossRef]
  55. Mirjalili, S.; Lewis, A. The whale optimization algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  56. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  57. Derrac, J.; García, S.; Molina, D.; Herrera, F. A practical tutorial on the use of nonparametric statistical tests as a methodology for comparing evolutionary and swarm intelligence algorithms. Swarm Evol. Comput. 2011, 1, 3–18. [Google Scholar] [CrossRef]
  58. Friedman, M. The use of ranks to avoid the assumption of normality implicit in the analysis of variance. J. Am. Stat. Assoc. 2012, 32, 674–701. [Google Scholar] [CrossRef]
  59. Zhao, W.; Shi, T.; Wang, L.; Cao, Q.; Zhang, H. An adaptive hybrid atom search optimization with particle swarm optimization and its application to optimal no-load PID design of hydro-turbine governor. J. Comput. Des. Eng. 2021, 8, 1204–1233. [Google Scholar] [CrossRef]
  60. Hemeida, M.G.; Ibrahim, A.A.; Mohamed, A.A.A.; Alkhalaf, S.; El-Dine, A.M.B. Optimal allocation of distributed generators DG based Manta Ray Foraging Optimization algorithm (MRFO). Ain Shams Eng. J. 2021, 12, 609–619. [Google Scholar] [CrossRef]
  61. Chen, C.; Mo, C. A method for intelligent fault diagnosis of rotating machinery. Digit. Signal Prog. 2004, 14, 203–217. [Google Scholar] [CrossRef]
  62. Liu, J.; Li, D.; Wu, Y.; Liu, D. Lion swarm optimization algorithm for comparative study with application to optimal dispatch of cascade hydropower stations. Appl. Soft. Comput. 2020, 87, 105974. [Google Scholar] [CrossRef]
  63. Lu, P.; Zhou, J.; Wang, C.; Qiao, Q.; Mo, L. Short-term hydro generation scheduling of Xiluodu and Xiangjiaba cascade hydropower stations using improved binary-real coded bee colony optimization algorithm. Energy Conv. Manag. 2015, 91, 19–31. [Google Scholar] [CrossRef]
  64. Liu, J.; Liu, X.; Wu, Y.; Yang, Z.; Xu, J. Dynamic multi-swarm differential learning Harris Hawks Optimizer and its application to optimal dispatch problem of cascade hydropower stations. Knowl.-Based Syst. 2022, 242, 108281. [Google Scholar] [CrossRef]
Figure 1. (A) stability mode, (B) instability mode [34].
Figure 1. (A) stability mode, (B) instability mode [34].
Water 14 03035 g001
Figure 2. Pseudocode of ESDO algorithm.
Figure 2. Pseudocode of ESDO algorithm.
Water 14 03035 g002
Figure 3. Convergence comparisons of UNs.
Figure 3. Convergence comparisons of UNs.
Water 14 03035 g003
Figure 4. Convergence comparisons of MNs (F8–F13).
Figure 4. Convergence comparisons of MNs (F8–F13).
Water 14 03035 g004
Figure 5. Convergence comparisons of LMNs (F14–F23).
Figure 5. Convergence comparisons of LMNs (F14–F23).
Water 14 03035 g005
Figure 6. Ranks of different algorithms for each function.
Figure 6. Ranks of different algorithms for each function.
Water 14 03035 g006
Figure 7. Mean of ranks of different algorithms.
Figure 7. Mean of ranks of different algorithms.
Water 14 03035 g007
Figure 8. FOPID block diagram of the water turbine governor system.
Figure 8. FOPID block diagram of the water turbine governor system.
Water 14 03035 g008
Figure 9. Convergence curves of average fitness with 4% load.
Figure 9. Convergence curves of average fitness with 4% load.
Water 14 03035 g009
Figure 10. Convergence curves of FOPID parameters with 4% load.
Figure 10. Convergence curves of FOPID parameters with 4% load.
Water 14 03035 g010
Figure 11. Speed response curves with 4% load.
Figure 11. Speed response curves with 4% load.
Water 14 03035 g011
Figure 12. Convergence curves of average fitness with 8% load.
Figure 12. Convergence curves of average fitness with 8% load.
Water 14 03035 g012
Figure 13. Convergence curves of FOPID parameters with 8% load.
Figure 13. Convergence curves of FOPID parameters with 8% load.
Water 14 03035 g013
Figure 14. Speed response curves with 8% load.
Figure 14. Speed response curves with 8% load.
Water 14 03035 g014
Figure 15. Convergence curves of average fitness with 12% load.
Figure 15. Convergence curves of average fitness with 12% load.
Water 14 03035 g015
Figure 16. Convergence curves of FOPID parameters with 12% load.
Figure 16. Convergence curves of FOPID parameters with 12% load.
Water 14 03035 g016
Figure 17. Speed response curves with 12% load.
Figure 17. Speed response curves with 12% load.
Water 14 03035 g017
Figure 18. Convergence curves of average fitness with 16% load.
Figure 18. Convergence curves of average fitness with 16% load.
Water 14 03035 g018
Figure 19. Convergence curves of FOPID parameters with 16% load.
Figure 19. Convergence curves of FOPID parameters with 16% load.
Water 14 03035 g019
Figure 20. Speed response curves with 16% load.
Figure 20. Speed response curves with 16% load.
Water 14 03035 g020
Table 1. Parameter setting of each algorithm.
Table 1. Parameter setting of each algorithm.
AlgorithmParametersValues
GSAGravitational constant; decreasing coefficient100; 20
WOAControl parameter[0, 2]
SDOConvergence factorlinearly decreases from 2 to 0
GWOControl parameter[0, 2]
ESDOConvergence factorlinearly decreases from 2 to 0
Table 2. Comparison results of UNs (F1–F7).
Table 2. Comparison results of UNs (F1–F7).
NO.IndexESDOSDOWOAGWOGSA
F1Mean9.09 × 10−2311.49 × 10−1652.62 × 10−831.87 × 10−333.66 × 10−17
Std009.62 × 10−832.21 × 10−331.15 × 10−17
F2Mean4.93 × 10−1161.30 × 10−712.40 × 10−538.85 × 10−203.21 × 10−8
Std2.70 × 10−1157.04 × 10−718.38 × 10−535.55 × 10−206.07 × 10−9
F3Mean3.86 × 10−1943.66 × 10−12830167.46591.01 × 10−75.38 × 102
Std0.00 × 1001.99 × 10−1278562.6735624.27 × 10−72.49 × 102
F4Mean3.28 × 10−1133.22 × 10−7536.644947911.83 × 10−83.31 × 100
Std1.79 × 10−1131.76 × 10−7429.450251311.45 × 10−81.56 × 100
F5Mean2.40 × 1002.59 × 1012.75 × 1012.64 × 1013.85 × 101
Std1.11 × 1007.19 × 10−14.09 × 10−16.35 × 10−13.18 × 101
F6Mean003.33 × 10−206.33 × 10−1
Std001.83 × 10−108.09 × 10−1
F7Mean7.49 × 10−59.85 × 10−52.14 × 10−31.14 × 10−33.09 × 10−2
Std7.58 × 10−57.61 × 10−52.91 × 10−35.78 × 10−41.73 × 10−2
Table 3. Comparison results of MNs (F8–F13).
Table 3. Comparison results of MNs (F8–F13).
NO.IndexESDOSDOWOAGWOGSA
F8Mean−8525.069636−8498.031439−11,263.28556−5998.877219−2735.739471
Std5.93 × 1027.29 × 1021.55 × 1031.19 × 1033.67 × 102
F9Mean0001.51785667818.27407681
Std0002.13 × 1004.47 × 100
F10Mean8.88 × 10−168.88 × 10−165.27 × 10−154.21 × 10−145.05 × 10−9
Std002.41 × 10−154.53 × 10−158.67 × 10−10
F11Mean0005.48 × 10−31.76 × 101
Std0009.90 × 10−34.99 × 100
F12Mean5.77 × 10−52.62 × 10−45.68 × 10−32.94 × 10−26.31 × 10−1
Std8.05 × 10−54.09 × 10−44.10 × 10−31.75 × 10−25.01 × 10−1
F13Mean2.32 × 10−31.93 × 10−22.21 × 10−13.40 × 10−12.59 × 100
Std4.53 × 10−33.96 × 10−21.23 × 10−11.97 × 10−13.77 × 100
Table 4. Comparison results of LMNs (F14–F23).
Table 4. Comparison results of LMNs (F14–F23).
NO.IndexESDOSDOWOAGWOGSA
F14Mean0.9980038380.9980038382.4323008372.6078010745.100884714
Std003.36 × 1002.44 × 1002.80 × 100
F15Mean3.07 × 10−43.07 × 10−47.54 × 10−42.48 × 10−33.28 × 10−3
Std1.03 × 10−141.30 × 10−176.96 × 10−46.07 × 10−31.83 × 10−3
F16Mean−1.031628453−1.031628453−1.031628453−1.031628439−1.031628453
Std6.78 × 10−166.78 × 10−161.05 × 10−101.46 × 10−85.05 × 10−16
F17Mean0.3978873580.3978873580.3978893410.3979133460.397887358
Std004.92 × 10−61.34 × 10−40
F18Mean333.0000197863.0000124093
Std1.12 × 10−151.18 × 10−155.81 × 10−51.44 × 10−53.24 × 10−15
F19Mean−3.862782148−3.862782148−3.859813445−3.861222768−3.862782148
Std2.70 × 10−152.71 × 10−153.69 × 10−32.62 × 10−32.31 × 10−15
F20Mean−3.302179651−3.310105859−3.229898305−3.259834789−3.321995172
Std4.51 × 10−23.63 × 10−29.89 × 10−27.99 × 10−21.37 × 10−15
F21Mean−10.15319968−9.983266281−9.389732298−9.060734145−7.353288325
Std7.01 × 10−159.31 × 10−12.01 × 1002.26 × 1003.54 × 100
F22Mean−10.40294057−10.40294057−8.736752324−10.40158908−10.40294057
Std1.04 × 10−169.90 × 10−162.85 × 1006.32 × 10−49.33 × 10−16
F23Mean−10.53640982−10.53640982−7.471755875−10.53516103−10.53640982
Std1.21 × 10−151.98 × 10−153.59 × 1005.47 × 10−42.29 × 10−15
Table 5. WSRT results of ESDO for SDO and WOA.
Table 5. WSRT results of ESDO for SDO and WOA.
Fun.SDO vs. ESDOWOA vs. ESDO
p-ValueT+T−Winnerp-ValueT+T−Winner
F11.73 × 10−60465+1.73 × 10−60465+
F21.73 × 10−60465+1.73 × 10−60465+
F31.73 × 10−60465+1.73 × 10−60465+
F41.73 × 10−60465+1.73 × 10−60465+
F51.73 × 10−60465+1.73 × 10−60465+
F610465=10465=
F73.68 × 10−23341311.73 × 10−60465+
F89.92 × 10−1233232=6.89 × 10−542639
F910465=10465=
F1010465=4.34 × 10−60465+
F1110465=0.50465=
F124.49 × 10−2135330+1.73 × 10−60465+
F132.43 × 10−2123342+1.73 × 10−60465+
F1410465=1.73 × 10−60465+
F155.73 × 10−1166299=1.73 × 10−60465+
F1610465=1.73 × 10−60465+
F1710465=1.73 × 10−60465+
F18112453=1.73 × 10−60465+
F1910465=1.73 × 10−60465+
F204.53 × 10−119446=4.86 × 10−535430+
F2110465=1.73 × 10−60465+
F2210465=1.73 × 10−60465+
F2310465=1.73 × 10−60465+
Table 6. WSRT results of ESDO for GWO and GSA.
Table 6. WSRT results of ESDO for GWO and GSA.
Fun.GWO vs. ESDOGSA vs. ESDO
p-ValueT+T−Winnerp-ValueT+T−Winner
F11.73 × 10−60465+1.73 × 10−60465+
F21.73 × 10−60465+1.73 × 10−60465+
F31.73 × 10−60465+1.73 × 10−60465+
F41.73 × 10−60465+1.73 × 10−60465+
F51.73 × 10−60465+1.73 × 10−60465+
F610465=1.56 × 10−20465+
F71.73 × 10−60465+1.73 × 10−60465+
F81.73 × 10−60465+1.73 × 10−60465+
F94.63 × 10−60465+1.73 × 10−60465+
F101.31 × 10−60465+1.73 × 10−60465+
F113.91× 10−30465+1.73 × 10−60465+
F122.35 × 10−63462+1.73 × 10−60465+
F131.92 × 10−61464+1.24 × 10−520445+
F141.73 × 10−60465+1.73 × 10−60465+
F151.73 × 10−60465+1.73 × 10−60465+
F161.73 × 10−60465+10465=
F171.73 × 10−60465+10465=
F181.73 × 10−60465+5.34 × 10−715450+
F191.73 × 10−60465+10465=
F202.41 × 10−454411+6.25 × 10−215450=
F211.73 × 10−60465+1.75 × 10−40465+
F221.73 × 10−60465+10465=
F231.73 × 10−60465+10465=
Table 7. Statistics results of WSRT for ESDO.
Table 7. Statistics results of WSRT for ESDO.
Function TypesESDO vs. SDO
(+/=/−)
ESDO vs. WOA
(+/=/−)
ESDO vs. WOA
(+/=/−)
ESDO vs. GSA
(+/=/−)
Unimodal5/1/16/1/06/1/07/0/0
Multimodal2/4/02/3/16/0/06/0/0
Low-dimensional0/10/010/0/010/0/04/6/0
Total7/15/118/4/122/1/017/6/0
Table 8. The specific experimental settings and parameter values.
Table 8. The specific experimental settings and parameter values.
Population SizeMaximum Numbers of Iterations K p K i K d λ μ
30100 [ 0 , 20 ] [ 0 , 20 ] [ 0 , 20 ] [ 0.5 , 1.9 ] [ 0.5 , 1.9 ]
T v T y T y 1 K y T w T a e n
0.20.30.310.5101
Table 9. Experimental results at 4%, 8%, 12%, and 16% loads.
Table 9. Experimental results at 4%, 8%, 12%, and 16% loads.
LoadAverageAlgorithms
ESDOSDOGWOWOAGSA
4%ITAE1.50281.70091.92374.07465.5928
Kp12.648712.045711.739210.319910.678
Ki4.58884.25214.51355.02794.0414
Kd3.34033.02823.37053.69213.3575
λ1.02281.02301.01530.99871.0830
μ1.50791.45191.42801.34161.3857
Overshoot0.49920.51000.50700.51130.5072
Ts(s)10.1610.8115.2021.7828.5100
8%ITAE3.18913.43784.39428.249019.0396
Kp11.429510.735210.29928.599610.2509
Ki4.29844.15324.12794.63324.2505
Kd3.392813.29583.18013.74713.7841
λ1.02341.02111.01560.96791.0563
μ1.43411.35631.36111.15681.3633
Overshoot1.00731.02741.02721.06011.0064
Ts(s)10.0210.7710.7727.950028.7300
12%ITAE5.11165.25966.252012.297433.1630
Kp10.653310.264910.29007.78788.9427
Ki4.03613.86553.76924.46303.3708
Kd3.77963.19633.12253.99084.0885
λ1.02331.02291.02360.96531.0339
μ1.40561.37351.42441.14031.3098
Overshoot1.50731.53381.51351.58061.5024
Ts(s)10.3210.4510.4222.5828.09
16%ITAE8.43689.09789.945518.871040.3035
Kp10.15489.62699.23277.42389.2411
Ki3.46103.31663.37393.48103.1995
Kd3.77432.74253.26793.78483.9186
λ1.02431.01971.01700.98321.0734
μ1.39361.38591.33151.12471.3524
Overshoot2.06532.10872.10052.15382.0577
Ts(s)11.0911.1513.4027.2927.90
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, W.; Zhang, H.; Zhang, Z.; Zhang, K.; Wang, L. Parameters Tuning of Fractional-Order Proportional Integral Derivative in Water Turbine Governing System Using an Effective SDO with Enhanced Fitness-Distance Balance and Adaptive Local Search. Water 2022, 14, 3035. https://doi.org/10.3390/w14193035

AMA Style

Zhao W, Zhang H, Zhang Z, Zhang K, Wang L. Parameters Tuning of Fractional-Order Proportional Integral Derivative in Water Turbine Governing System Using an Effective SDO with Enhanced Fitness-Distance Balance and Adaptive Local Search. Water. 2022; 14(19):3035. https://doi.org/10.3390/w14193035

Chicago/Turabian Style

Zhao, Weiguo, Hongfei Zhang, Zhenxing Zhang, Kaidi Zhang, and Liying Wang. 2022. "Parameters Tuning of Fractional-Order Proportional Integral Derivative in Water Turbine Governing System Using an Effective SDO with Enhanced Fitness-Distance Balance and Adaptive Local Search" Water 14, no. 19: 3035. https://doi.org/10.3390/w14193035

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