A hybrid SA-MFO algorithm for function optimization and engineering design problems

This paper presents a hybrid algorithm based on using moth-flame optimization (MFO) algorithm with simulated annealing (SA), namely (SA-MFO). The proposed SA-MFO algorithm takes the advantages of both algorithms. It takes the ability to escape from local optima mechanism of SA and fast searching and learning mechanism for guiding the generation of candidate solutions of MFO. The proposed SA-MFO algorithm is applied on 23 unconstrained benchmark functions and four well-known constrained engineering problems. The experimental results show the superiority of the proposed algorithm. Moreover, the performance of SA-MFO is compared with well-known and recent meta-heuristic algorithms. The results show competitive results of SA-MFO concerning MFO and other meta-heuristic algorithms.


Introduction
Optimization is defined as the process of finding global or near global optimum solutions for a given problem. Many problems in the real world can be observed as optimization problems. Over the past several decades, several optimization algorithms are proposed to solve many optimization problems. Simulated annealing (SA) algorithm is one of the most popular iterative meta-heuristic optimization algorithm. It used to solve continuous and discrete problems [16]. The key privilege of SA is that its ability to escape from local optima called metropolis process. However, the major shortage of SA is that its efficiency for finding the global optimum is unsatisfactory. This is due to, when SA generates a new candidate solution, it does not learn intelligently from its searching history [37]. Several studies have been introduced for enhancing SA performance. Some of these are based on changing the acceptance or generation mechanisms. How-B Gehad Ismail Sayed GehadIsmail_FCI@yahoo.com Aboul Ella Hassanien aboitcairo@gmail.com 1 Faculty of Computers and Information, Cairo University, Cairo, Egypt 2 Scientific Research Group in Egypt (SRGE), Cairo, Egypt ever, this kind of modifications in SA schema does not usually inherit the ability of SA for escaping from local optima [11]. On the other hand, population-based meta-heuristics algorithms equipped with the ability to guide their search through using intelligent learning mechanism. Thus, the populationbased meta-heuristics algorithms showed better efficiency. Some of these algorithms are particle swarm optimisation (PSO) [15], ant colony optimization (ACO) [4], ant lion optimizer (ALO) [22], differential evolution (DE) [35], Social group optimization (SGO) [30], and moth-flame optimization (MFO) [24].
In this paper, MFO is selected to be studied and analyzed. MFO algorithm is one of recent meta-heuristic optimization algorithms proposed in 2015. The main inspiration of MFO came from the navigation method of moths in nature called transverse orientation. The inventor of this algorithm, Mirjalili, showed that MFO obtained very competitive results compared with other meta-heuristic optimization algorithms. Also, MFO shows competitive results for automatic mitosis detection in breast cancer histology images [32]. However, MFO like other meta-heuristic algorithms suffers from entrapping at local optima and low convergence rate. Thus, many studies are proposed to solve this problem. Several studies have been proposed to enhance the performance of MFO. Li et al. [17] proposed Lévy-flight moth-flame optimization (LMFO) algorithm to improve the performance of MFO. The proposed LMFO used Lévy-flight strategy in the searching mechanism of MFO. The experimental results on 19 unconstrained benchmark functions and two other constrained problems showed the superior performance of LMFO. Emery et al. [5] employed chaos parameter in the spiral equation of updating the position of moths. They applied their approach in feature selection application. Recently, many studies hybridize two or more algorithms to obtain optimal solutions for optimization problems [31]. Some of these studies hybridize differential evolution with biogeographybased optimization to solve global optimization problem [8]. Some hybridize particle swarm optimization with differential evolution for solving constrained numerical and engineering optimization problems [19]. Moreover, some of them hybridize chaotic with meta-heuristic algorithms for solving feature selection problem [33]. Also, some hybridize genetic algorithm and particle swarm optimization for multimodal problems [13]. Furthermore, some researchers considering introducing new methods to improve meta-heuristic algorithms. Author in [25] introduces new mutation rule to be combined with the basic mutation strategy of differential evolution algorithm. He named his algorithm as enhanced adaptive differential evolution (EADE). The experimental results show that EADE is very competitive algorithm for solving large-scale global optimization problems. Authors in [28] combined two binary bat algorithm (BBA) and with local search scheme (LSS)to improve the diversity of BBA and enhance the convergence rate. The results show that the proposed hybrid algorithm obtained promising results for solving large-scale 01 knapsack problems. Authors in [36] introduced a mutation operator guided by preferred regions to enhance the performance of an existing set-based evolutionary many-objective optimization algorithm. The results obtained from 21 instances of seven benchmark show that their proposed algorithm is superior. The main drivers for this kind of hybridization are to avoid local optima and improve search results from global optimization. This becomes an essential task in multidimensional and multiobjective problems.
Therefore, in this paper, we tried another kind of hybridization to solve the main problems of MFO. To best of our knowledge, this kind of hybridization of MFO with SA is not used before. The salient features of SA and MFO are hybrid to create a new approach, namely simulated annealing moth-flame algorithm (SA-MFO). In this paper, simulated annealing method is applied during the updating positions of flames of the standard version of MFO to further enhance the performance of MFO. The proposed SA-MFO is applied to solve both constrained and unconstrained benchmark problems. The proposed mimetic algorithm is mainly based on MFO algorithm, whereas, metropolis mechanism of SA is used to slow down the convergence rate of MFO and increases the probability of moths to reach the global optima. Also, metropolis mechanism is used to increase the diversity of the population. The proposed SA-MFO has the advantages of both SA and MFO algorithms: the ability to escape from local optima mechanism of SA and fast searching ability and learning mechanism for guiding the generation of candidate solutions of MFO.
The proposed mimetic algorithm applied on 23 benchmark functions and four well-known engineering problems. These problems are applied previously on the original MFO. Thus, we almost make a fair comparison between the proposed SA-MFO and MFO in terms of evaluation criteria, number of iterations, number of independent runs, population size, and number of dimensionality. The simulation results demonstrate that the proposed SA-MFO can significantly improve the performance of MFO. Also, the experimental results show that SA-MFO is superior to the other algorithms.
The rest of this paper is organized as follows: Section 2 provides a brief overview of the standard moth flame optimization algorithm, followed by a brief introduction to simulated annealing. Section 3 describes in detail the proposed hybrid algorithm, and then Sect. 4 discusses the main results and evaluate the performance of the proposed algorithm. Finally, Sect. 5 concludes the paper.

Inspiration
In this section, one of the recent stochastic population-based algorithms is employed namely moth flame optimization (MFO) [24]. The main inspiration of MFO came from the navigation method of moths in nature. Moths are fancy insects that are very similar to the butterfly family. In nature, there are greater than 160,000 various species of this insect. Larvae and adult are the two main milestones in their lifetime. The larva is converted to moth by cocoons. Special navigation methods in the night are the most interesting fact about moths. They used a mechanism called transverse orientation for their navigation. Moths fly using a fixed angle with respect to the moon, which is a very efficient mechanism for long traveling distances in a straight line.

Mathematical model of MFO
Let the candidate's solutions are moths, and the problem's variables are the position of moths in the space. P is the spiral function where moths move around the search space. Each moth updates his position with respect to flame using the following equations: where M i indicates the ith moth and F j is j t h flame. There are other types of spiral functions can be utilized respect to the following rules: 1. The initial point of spiral should start from the moth. 2. The final point of spiral should be the position of the flame. 3. The Fluctuation range of spiral shouldn't exceed the search space.
where D i is the distance of the ith moth for the jth flame, t is a random number in [−1, 1] and b is a constant for defining the shape of the P. D is calculated using the following equation.
where M i is the ith moth, F j indicates the jth flame, and D i indicates the distance of the ith moth to the jth flame. Another concern, the moths update their position with respect to n different locations in the search space which can degrade the best promising solutions exploitation. Therefore, the number of flames adaptively decreases over the course of iterations using the following formula: where I is the current number of iterations, I N is the maximum number of iterations and FN is the maximum number of flames. MFO utilizes Quicksort method and the computational complexity of this sort method is O(nlogn) and O(n 2 ) in the best and worst case, where n denotes the number of moths.

Simulated annealing (SA)
Simulated annealing (SA) was proposed by Metropolis et al. [21] in 1953. The main inspiration of SA came from the simulation modeling of a molecular movement in the materials during annealing (physical annealing principles). The process of cooling and heating a material to recrystallized is called annealing. At high temperature, the particles move in disorder and when the temperature decreases gradually, the particles converge to the lowest state of energy. SA uses probability acceptance mechanisms for considering the current optimal solution. SA was successfully applied to solve the optimization problems by Kirkpatrick et al. [15] in 1983. Through using the probability of acceptance of metropolis process, it can be simulated the annealing mechanism in finding the lowest energy state which is the optimum solution. The main parameters of SA are metropolis acceptance mechanism, initial condition, and cooling scheduling method. They are described as the following: Initial condition at the beginning of the searching process of SA algorithm starts with initializing the initial temperature T 0 , which have a significant impact on the performance. The higher value of T 0 will increase the computational time of the optimization process. On the contrary, if the T 0 is too low, it will cause that SA will not effectively explore the search space. Thus it is needed to be selected carefully to obtain the optimal solution.
Metropolis acceptance mechanism It is one of key affecting on obtaining a near optimal solution quickly. The metropolis rules used to indicate whether the nearby solution with low fitness value is accepted or not as the current solution. This mechanism affects the capability of SA to escape from the local optimum. Let F C represents the fitness value of the current solution, and F N denotes the fitness value of the neighboring solution. When F(Y 0 ) is better than F(Y ), then SA will use the acceptance probability mechanism to indicating either to consider the neighbor solution as the current solution or not. The probability of acceptance mechanism is defined as follows: where Pr is defined as the probability of acceptance and T k is temperature value at time k.
Cooling schedule method During the physical annealing process, as the temperature decreases, the cooling rate decreases too to reach the stable ground state. This it's needed that the system at the beginning to be cooled faster and slower as gradually decreased of the temperature. The cooling schedule parameter is defined as follows: where α denotes the temperature coefficient, T k is the initial temperature value and T k+1 is temperature at time k

The proposed SA-MFO algorithm
According to the MFO algorithm, moths' positions are randomly initialized. Then moth swarm updates their position movement using Eq. (2). Through this, the next moth position is determined. The standard version of MFO automatically accepts the current moth position as the new moth position. However, SA-MFO uses the metropolis acceptance mechanism of SA. This mechanism indicates whether to accept the new moth position or not. In case of the new moth, the position is not accepted, another candidate position is recalculated. The selection is based on its fitness value. Using this mechanism, a moth solution able to escape from the local optima. Also, the quality of the solution is enhanced with fastest convergence rate [34]. Based on temperature, cooling schedule parameters and the fitness value difference of the metropolis mechanism, the optimization process will explore the search space effectively to find an optimal solution. This process will be repeated until a new position is accepted using metropolis acceptance rule, or the termination criterion is reached.

Parameters initialization
In the beginning, SA-MFO starts with setting MFO parameters and randomly initialized moths' positions within the search space. Each position represents a solution in the search space. The initial parameter settings are presented in Table 1.

Fitness function
At each iteration, each moth position is evaluated using a predefined fitness function f (z). Then based on this value, the metropolis rule is applied to determine whether to accept or not.

Updating positions
If the current moth accepted as the new moth, then the new moths change its position using Eq. (2).

Termination criteria
The optimization process terminates when it reaches the maximum number of iterations or when the best solution is found [7]. In our case, we used the maximum number of iterations is 1000 iterations (see Table 1). The overall SA-MFO algorithm is proposed in Algorithm (1). Moreover, the summarize flowchart of the proposed simulated annealing moth-flame algorithm is given in Fig. 1. Calculate Calculate D using Equ. (3). 12: Update M(i, j) using Equ. (1) and (2) with respect to moth. 13: else 14: Set

Results and discussion
In this paper, 23 optimization benchmark problems with different characteristics and four other engineering problems are employed to evaluate the performance of the proposed SA-MFO algorithm. The employed 23 benchmark problems divided into two families: unimodal and multimodal. Tables 11 and 12 in Appendix 1 and 2 describes the mathematical formulation and the properties of the benchmark functions used in this paper. In the table, lb and ub represent the lower and upper bound of the search space respectively, dim indicates the number of dimensions, and opt indicates the optimum solution of the function. More information about the adopted benchmark functions can be found in [18]. The main characteristics of the first family, namely unimodal test functions have only one global optimum and no local optima. These functions are used to evaluate the exploitation and convergence rate of the adopted algorithms. However, the second family, namely multimodal test functions have multiple global optimum and multiple local optima. The characteristic of this family is used to evaluate the capability of the algorithm in avoiding the local optima and the explorative ability of an algorithm.
In this section, three main experiments were conducted to evaluate the performance of the SA-MFO algorithm. The first experiment aims to determine the optimal parameter value for both T 0 and α. The second experiment aims to evalu- ate and compare the performance of SA-MFO with MFO and five well-known optimization algorithms: particle swarm optimization (PSO) [15], ant bee colony (ABC) [14], moth flame optimization (MFO) [24], multi verse optimization (MVO) [23], and ant lion optimizer (ALO) [22] in solving 23 numerical optimization problems using different statistical measurements. The parameters settings for all meta-heuristic optimization algorithms are shown in Table 2. The rest of parameters including, the maximum number of iterations for all functions is set to 1000, the number of search agents is set to 50 individuals for all comparisons and the initialization of the search agents is same for all adopted algorithms. The reason behind this that we want to make almost a fair comparison for these algorithms. The third and last experiment aims to evaluate the performance of SA-MFO on solving engineering problems, also to compare the best-obtained results with other meta-heuristic algorithms. All the obtained results for both benchmark and engineering problems are calculated on average for 30 independent runs. Moreover, all these experiments are performed on the same PC with Core i3 and RAM 2 GB on OS windows 7.

Parameters optimization experiment
Due to the stochastic behavior of the optimization algorithms, there is a need to set certain parameters carefully. These parameters can have a significant influence on the performance of the algorithm. The main parameters of SA are T 0 and α which affecting the searching efficiency to find the optimal solution. These two parameters need to be optimized effectively. In this paper, the optimal parameters values for T 0 and α are identified based on the highest success rate for 50 The significance bold means indicates the optimal values for T 0 , which are determined based on highest success rate for 50 runs of F1; F10 and F18 The significance bold means indicates the optimal values for α, which are determined based on highest success rate for 50 runs of F1; F10 and F18 runs of F 1 , F 10 and F 18 . This paper uses the same mechanism proposed in [34] which depends on optimizing one parameter per simulation. When an optimal value of the parameter is determined, the optimization process is repeated till to find an optimal value of the next parameter. All the obtained results reach 1 × E −3 are observed as a success. Table 3 shows the simulation results obtained for selecting the optimal value of T 0 . Changing the value of T 0 can highly affect on the initial acceptance rate of SA and SA-MFO. As, if the T 0 is too high that will lead the algorithm to spend much time in exploiting solution space and if the α is too low that will cause increasing the search time. In this paper, T 0 set with 10% and increases by 10% each time till to reach 90% as shown in Table 3. As it can be seen, T 0 with value 90% improves the success rate (Table 3). α can highly affect by the acceptance of the poor solution which may lead to exploring the search space without finding an optimal solution. If α is too low, it may lead to stuck in the local optima. In this paper, α set with 10% and increases by 10% each time till to reach 90% as shown in (Table 4). As it can be seen, α with value 10% improves the success rate.

Comparison using numerical benchmark functions experiment
In this subsection, the performance of SA-MFO is analyzed and compared with other meta-heuristic algorithms in terms of mean, standard deviation (std.) for 30 independent runs.  Moreover, the p value from Wilcoxon rank sum test is used as well. Table 5 shows the statistical results obtained for 23 benchmark functions using SA-MFO, PSO, ALO, ABC, MFO, and MVO. As it can be observed from Table 5, SA-MFO outperforms the other algorithms for both unimodal and multimodal test function in most of the functions. The second best results belong the MFO and PSO algorithms. Considering the characteristic of both multi-modal and uni-modal test functions and the obtained results, it may be concluded that SA-MFO algorithm has a high exploration capability for exploring the promising area in the search space and the capability for avoiding the local optima. For further analyzing the performance of the proposed SA-MFO algorithm, convergence curves for all the used benchmark functions are shown in Fig. 2. In this figure, the performance of the proposed mimetic algorithm is compared with MFO, ABC, MVO, PSO, and ALO. As it can be seen, the proposed SA-MFO algorithm overtakes the other algorithms in most of the cases in terms of both stability and high performance. Moreover, it can be observed that ABC in the most case provides the lowest score. PSO and MVO are in second place. Figure 3 shows the mathematical representation of some benchmark functions, average fitness values calculated for whole moths in each iteration, trajectory changes for the first moth in the first dimension during optimization. As it can be observed, the fluctuations of the first moth are gradually decreased over the course of the iteration. Through this behavior, it can be guaranteed the transition between exploration and exploitation of the search space. Also, Fig. 3 compares the convergence curves of the proposed algorithm with original MFO. From all these figures, it can be observed that using metropolis mechanism can significantly enhance the performance of the original MFO.
From all the obtained results, it can be observed that the proposed SA-MFO in some cases obtained similar results as the original MFO. The explanation for these results is that Metropolis mechanism is only applied in case of the current solution is worse than the best solution obtained so far. Thus, in some cases, based on the probability of acceptance, the current solution is not accepted. Therefore, SA-MFO obtains similar results as original MFO.
A non-parametric statistical test namely Wilcoxon rank sum test is conducted in this paper. This test is used to verify if a pair of two solutions is statistically different or not [40]. In this paper, Wilcoxon rank sum test is used to compare and analyze the difference of the solution (population) of the proposed algorithm with other algorithms. The best values of p when p value < 0.05 which can be considered sufficient evidence against the null hypothesis. The obtained p values produced by Wilcoxon's test of 23 benchmark functions for the pair-wise comparison of the best score obtained for all 1000 iterations with a 5% significance level from the employed statistical test are reported in Table 6. Such pair-wise comparison is constituted by SA-MFO vs. MFO, SA-MFO vs. ALO, SA-MFO vs. PSO, SA-MFO vs. ABC, and SA-MFO vs. MVO. In this table, N/A denotes "not applicable," which mean that the proposed algorithm could not be compared with the other algorithm in the rank-sum test. As it can be observed from this table, SA-MFO outperforms the other algorithms for both unimodal and multimodal test function. As the SA-MFO provides p value less than 0.05 for most of the functions. Also, it can be seen that SA-MFO reject the null hypothesis for F 16 and F 17 for most of the competitive algorithms. From all the obtained results, it can be concluded that the SA-MFO algorithm has high exploitation and exploration compared with the original version MFO and other meta-heuristic algorithms. Regarding the p value, it can be concluded that SA-MFO is statistically significant.

Comparison using engineering design problems experiment
In this section, the performance of the proposed SA-MFO is evaluated using four well-known engineering problems, namely three-bar truss, welded beam, pressure vessel and tension/compression spring design problems. The detailed description of the adopted engineering problem is presented in Appendix 2. Table 7 presents the optimal solution for three-bar truss design problem obtained by particle swarm optimization with differential evolution (PSODE) [20], dynamic stochastic selection (DEDS) [42], mine blast algorithm (MBA) [29], water cycle algorithm (WCA) [6], MFO and SA-MFO. As it can be seen from this table, SA-MFO obtained better results compared with the standard MFO. Table 8 compares the optimal solution obtained for pressure vessel problem by SA-MFO with hybrid particle swarm optimization (HPSO) [10], co-evolutionary particle swarm optimization using Gaussian distribution (CPSOGD) [27], GA based on using of dominance-based tour tournament selection (GA-TTS) [3], co-evolutionary differential evolution (CDE) [12], PSO [15], hybrid Nelder-mead simplex search and particle swarm optimization (NMPSO) [41], Gaussian quantum-behaved particle swarm optimization (GQPSO) [1], and MFO. As it can be observed, SA-PSO and ACO obtained the optimal solutions compared with the other algorithms. Table 9 compares the best obtained results for tension/compression spring design problem of SA-MFO and CPSO, DEDS, NM-PSO, GA-TTS, hybrid evolutionary algorithm and adaptive constraint handling technique (HEAA) [39], mine blast algorithm (MBA) [30], WCA, differential evolution with level comparison (DELC) [38], and MFO. In this table, NA indicates "not available" value. As it can be seen, NM-PSO obtained the optimal solution. However, it can be observed that SA-MFO obtained better results compared with original MFO. The best solutions for welded beam prob-

Computational complexity of the proposed SA-MFO algorithm
In this section, the computational complexity of the proposed hybrid algorithm is computed. The computational complexity of a meta-heuristic algorithm is mainly depended on the number of variables, the number of moths, the flames' sorting mechanism in each iteration and the maximum number of iterations [24]. Since the original MFO uses Quicksort mechanism, the computational complexity is O(n 2 ) in the worst case. In addition, the temperature in simulated annealing goes through O(log(m)) [9]. The overall computational where T is the maximum number of iterations and v is the number of variables.

Conclusions
This paper introduced a new hybrid algorithm based on using MFO and SA. The metropolis acceptance mechanism of SA is employed into the standard MFO to enhance the performance of MFO. The simulation results using 23 benchmark functions show that the proposed SA-MFO is superior compared with the original version MFO and other Table 8 Comparison of best obtained results using various optimization algorithms for pressure vessel problem z GA-TTS

Appendix 1: List of 23 benchmark optimization functions
See Tables 11 and 12.