A new bio-inspired metaheuristic algorithm for solving optimization problems based on walruses behavior

This paper introduces a new bio-inspired metaheuristic algorithm called Walrus Optimization Algorithm (WaOA), which mimics walrus behaviors in nature. The fundamental inspirations employed in WaOA design are the process of feeding, migrating, escaping, and fighting predators. The WaOA implementation steps are mathematically modeled in three phases exploration, migration, and exploitation. Sixty-eight standard benchmark functions consisting of unimodal, high-dimensional multimodal, fixed-dimensional multimodal, CEC 2015 test suite, and CEC 2017 test suite are employed to evaluate WaOA performance in optimization applications. The optimization results of unimodal functions indicate the exploitation ability of WaOA, the optimization results of multimodal functions indicate the exploration ability of WaOA, and the optimization results of CEC 2015 and CEC 2017 test suites indicate the high ability of WaOA in balancing exploration and exploitation during the search process. The performance of WaOA is compared with the results of ten well-known metaheuristic algorithms. The results of the simulations demonstrate that WaOA, due to its excellent ability to balance exploration and exploitation, and its capacity to deliver superior results for most of the benchmark functions, has exhibited a remarkably competitive and superior performance in contrast to other comparable algorithms. In addition, the use of WaOA in addressing four design engineering issues and twenty-two real-world optimization problems from the CEC 2011 test suite demonstrates the apparent effectiveness of WaOA in real-world applications. The MATLAB codes of WaOA are available in https://uk.mathworks.com/matlabcentral/profile/authors/13903104.

these initial solutions are improved. Finally, the best solution found during the implementation of the algorithm is introduced as the solution to the problem 4 . However, none of the metaheuristic algorithms guarantee that they will be able to provide the optimal global solution. This insufficiency is due to the nature of random search in these types of optimization approaches. Hence, the solutions derived from metaheuristic algorithms are known as quasi-optimal solutions 5 .
Exploration and exploitation capabilities enable metaheuristic algorithms to provide better quasi-optimal solutions. Exploration refers to the ability to search globally in different areas of the problem-solving space to discover the best optimal area. In contrast, exploitation refers to the ability to search locally around the available solutions and the promising areas to converge to the global optimal. Balancing exploration and exploitation is the key to the success of metaheuristic algorithms in achieving effective solutions 6 . Achieving better quasioptimal solutions has been the main challenge and reason for researchers' development of various metaheuristic algorithms 7,8 .
The main research question is that despite the numerous metaheuristic algorithms introduced so far, is there still a need to develop new algorithms? The No Free Lunch (NFL) theorem 9 answers the question that the optimal performance of an algorithm in solving a set of optimization problems gives no guarantee for the similar performance of that algorithm in solving other optimization problems. The NFL theorem concept rejects the hypothesis that a particular metaheuristic algorithm is the best optimizer for all optimization applications over all different algorithms. Instead, the NFL theorem encourages researchers to continue to design newer metaheuristic algorithms to achieve better quasi-optimal solutions for optimization problems. This theorem has also motivated the authors of this paper to develop a new metaheuristic algorithm to address optimization challenges. This paper's novelty and contribution are in designing a new metaheuristic algorithm called the Walrus Optimization Algorithm (WaOA), which is based on the simulation of walrus behaviors in nature. The main contributions of this article are as follows: • The natural behaviors of walruses inspire WaOA's design in feeding when migrating, fleeing, and fighting predators. • WaOA is mathematically modeled in three phases: exploration, exploitation, and migration.
• The efficiency of WaOA in handling optimization problems is tested on sixty-eight standard objective functions of various types of unimodal, multimodal, the CEC 2015 test suite, and the CEC 2017 test suite. • WaOA performance is compared with the performance of ten well-known metaheuristic algorithms.
• The success of WaOA in real-world applications is challenged in addressing four engineering design issues and twenty-two real-world optimization problems from the CEC 2011 test suite.
The rest of the paper is as follows. The literature review is presented in the "Literature review" section. The proposed WaOA approach is introduced and modeled in the "Walrus Optimization Algorithm" section. Simulation studies are presented in the "Simulation studies and results" section. The efficiency of WaOA in solving engineering design problems is evaluated in the "WaOA for real world-application" section. Conclusions and future research directions are included in the "Conclusions and future works" section.

Literature review
Metaheuristic algorithms are based on the inspiration and simulation of various natural phenomena, animal strategies and behaviors, concepts of biological sciences, genetics, physics sciences, human activities, rules of games, and any evolution-based process. Accordingly, from the point of view of the main inspiration used in the design, metaheuristic algorithms fall into five groups: evolutionary-based, swarm-based, physics-based, human-based, and game-based.
Evolutionary-based metaheuristic algorithms have been developed using the concepts of biology, natural selection theory, and random operators such as selection, crossover, and mutation. Genetic Algorithm (GA) is one of the most famous metaheuristic algorithms, which is inspired by the process of reproduction, Darwin's theory of evolution, natural selection, and biological concepts 10 . Differential Evolution (DE) is another evolutionary computation that, in addition to using the concepts of biology, random operators, and natural selection, uses a differential operator to generate new solutions 11 .
Swarm-based metaheuristic algorithms have been developed based on modeling natural phenomena, swarming phenomena, and behaviors of animals, birds, insects, and other living things. Particle Swarm Optimization (PSO) is one of the first introduced metaheuristics methods and was widely used in optimization fields. The main inspiration in designing PSO is the search behaviors of birds and fish to discover food sources 12,13 . Ant Colony Optimization (ACO) is a swarm-based method inspired by the ability and strategy of an ant colony to identify the shortest path between the colony to food sources 14 . Grey Wolf Optimization (GWO) is a metaheuristic algorithm inspired by grey wolves' hierarchical structure and social behavior while hunting 15 . Marine Predator Algorithm (MPA) has been developed inspired by the ocean and sea predator strategies and their Levy flight movements to trap prey 16 . The strategy of the tunicates and their search mechanism in the process of finding food sources and foraging have been the main inspirations in the design of the Tunicate Swarm Algorithm (TSA) 17 . Some other swarm-based methods are White Shark Optimizer (WSO) 18 , Reptile Search Algorithm (RSA) 19 , Raccoon Optimization Algorithm (ROA) 20 , African Vultures Optimization Algorithm (AVOA) 21 , Farmland Fertility Algorithm (FFA) 22 , Slime Mould algorithm (SMA) 23 , Mountain Gazelle Optimizer (MGO) 24 , Sparrow Search Algorithm (SSA) 25 , Whale Optimization Algorithm (WOA) 26 , Artificial Gorilla Troops Optimizer (GTO) 27 , and Pelican Optimization Algorithm (POA) 28 .
Physics-based metaheuristic algorithms have been inspired by physics' theories, concepts, laws, forces, and phenomena. Simulated Annealing (SA) is one of the most famous physics-based methods, the main inspiration

Walrus Optimization Algorithm
In this section, employed fundamental inspiration and the theory of the proposed Walrus Optimization Algorithm (WaOA) is stated, then its various steps are modeled mathematically.

Inspiration of WaOA.
Walrus is a big flippered marine mammal with a discontinuous distribution in the Arctic Ocean and subarctic waters of the Northern Hemisphere around the North Pole 52 . Adult walruses are easily identifiable with their large whiskers and tusks. Walruses are social animals who spend most of their time on the sea ice, seeking benthic bivalve mollusks to eat. The most prominent feature of walruses is the long tusks of this animal. These are elongated canines seen in both male and female species that may weigh up to 5.4 kg and measure up to 1 m in length. Males' tusks are slightly thicker and longer and are used for dominance, fighting, and display. The most muscular male with the longest tusks dominates the other group members and leads them 53 . An image of walrus is presented in Fig. 1. As the weather warms and the ice melts in late summer, walruses prefer to migrate to outcrops or rocky beaches. These migrations are very dramatic and involve massive aggregations of walruses 54 . The walrus has just two natural predators due to its large size and tusks: the polar bear and the killer whale (orca). Observations show that the battle between a walrus and a polar bear is very long and exhausting, and usually, polar bears withdraw from the fight after injuring the walrus. However, walruses harm www.nature.com/scientificreports/ polar bears with their tusks during this battle. In the fight against walruses, killer whales can hunt them successfully, with minimal and even no injuries 55 .
The social life and natural behaviors of walruses represent an intelligent process. Of these intelligent behaviors, three are the most obvious: (i) Guiding individuals to feed under the guidance of a member with the longest tusks.
Tracking the best population member in the search process directs the algorithm toward promising areas. In the social life of walruses, the most potent walrus, which can be recognized as having the longest tusk, is responsible for guiding the other walruses. Moving walruses in this process leads to significant changes in their position. Simulating these large displacements increases the algorithm's ability in global search and exploration ability.
(ii) Migration of walruses to rocky beaches. One of the natural behaviors of walruses is their migration due to warming weather in summer. In this process, walruses make big changes in their position by moving towards outcrops or rocky beaches. In the WaOA simulation for a walrus, the position of other walruses are assumed as migration destinations, one of these positions is randomly selected, and the walrus moves towards it. In the design of WaOA, imitating this strategy, global search and discovery capabilities are improved. The difference between the migration strategy and the foraging process under the guidance of the strongest walrus is that in this process, the population update process is prevented from relying on a particular member, such as the best member of the population. This updating process prevents early convergence and the algorithm from getting stuck in local optima.
(iii) Fight or escape from predators. The fighting strategy of walruses in the face of their predators, such as the polar bear and the killer whale, is a long chase process. This chasing process takes place in a small area around the walrus position and causes small changes in the walrus position. Therefore, simulating the small displacements of the walrus by aiming at better positions during the fight leads to an increase in WaOA's ability to search locally and exploit to converge to better solutions.
Mathematical modeling of these behaviors is the primary inspiration for developing the proposed WaOA approach.

Algorithm initialization.
WaOA is a population-based metaheuristic algorithm in which the searcher members of this population are walruses. In WaOA, each walrus represents a candidate solution to the optimization problem. Thus, the position of each walrus in the search space determines the candidate values for the problem variables. Therefore, each walrus is a vector, and the population of walruses can be mathematically modeled using so-called the population matrix. At the beginning of WaOA implementation, populations of walruses are randomly initialized. This WaOA population matrix is determined using (1).
where X is the walruses' population, X i is the i th walrus (candidate solution), x i,j is the value of the j th decision variable suggested by the i th walrus, N is the number of walruses, and m is the number of decision variables.
As mentioned, each walrus is a candidate solution to the problem, and based on its suggested values for the decision variables, the objective function of the problem can be evaluated. The estimated values for the objective function obtained from walruses are specified in (2).
where F is the objective function vector and F i is the value of the objective function evaluated based on the i th walrus.
Objective function values are the best measure of the quality of candidate solutions. The candidate solution that results in the evaluation of the best value for the objective function is known as the best member. On the other hand, the candidate solution that results in the worst value for the objective function is called the worst member. According to the update of the values of the objective function in each iteration, the best and worst members are also updated.
Mathematical modelling of WaOA. The process of updating the position of walruses in the WaOA is modeled in three different phases based on the natural behaviors of this animal.
Phase 1: feeding strategy (exploration). Walruses have a varied diet, feeding on more than sixty species of marine organisms, such as sea cucumbers, tunicates, soft corals, tube worms, shrimp, and various mollusks 57 . However, walrus prefers benthic bivalve mollusks, particularly clams, for which it forages by grazing around the www.nature.com/scientificreports/ sea floor, seeking and detecting food with its energetic flipper motions and susceptible vibrissae 58 . In this search process, the strongest walrus with the tallest tusks guides the other walrus in the group to find food. The length of the tusks in the walruses is similar to the quality of the objective function values of the candidate solutions. Therefore, the best candidate solution with the best value for the objective function is considered the strongest walrus in the group. This search behavior of the walruses leads to different scanning areas of the search space, which improves the exploration power of the WaOA in the global search. The process of updating the position of walruses is mathematically modeled based on the feeding mechanism under the guidance of the most vital member of the group, using (3) and (4). In this process, a new position for walrus is first generated according to (3). This new position replaces the previous position if it improves the objective function's value; this concept is modeled in (4).
where X P 1 i is the new generated position for the i th walrus based on the 1st phase, x P 1 i,j is its j th dimension, F P 1 i is its objective function value, rand i,j are random numbers from the interval [0, 1] , SW is the best candidate solution which is considered as the strongest walrus, and I i,j are integers selected randomly between 1 or 2. I i,j is used to increase the algorithm's exploration ability so that if it is chosen equal to 2, it creates more significant and broader changes in the position of walruses compared to the value of 1, which is the normal state of this displacement. These conditions help improve the algorithm's global search in escaping from the local optima and discovering the original optimal area in the problem-solving space.
Phase 2: migration. One of the natural behaviors of walruses is their migration to outcrops or rocky beaches due to the warming of the air in late summer. This migration process is employed in the WaOA to guide the walruses in the search space to discover suitable areas in the search space. This behavioral mechanism is mathematically modeled using (5) and (6). This modeling assumes that each walrus migrates to another walrus (randomly selected) position in another area of the search space. Therefore, the proposed new position is first generated based on (5). Then according to (6), if this new position improves the value of the objective function, it replaces the previous position of walrus.
where X P 2 i is the new generated position for the i th walrus based on the 2nd phase, x P 2 i,j is its j th dimension, F P 2 i is its objective function value, X k , k ∈ {1, 2, . . . , N} and k � = i, is the location of selected walrus to migrate the i th walrus towards it, x k,j is its j th dimension, and F k is its objective function value.
Phase 3: escaping and fighting against predators (exploitation). Walruses are always exposed to attacks by the polar bear and the killer whale. The strategy of escaping and fighting these predators leads to a change in the position of the walruses in the vicinity of the position in which they are located. Simulating this natural behavior of walruses improves the WaOA exploitation power in the local search in problem-solving space around candidate solutions. Since this process occurs near the position of each walrus, it is assumed in the WaOA design that this range of walrus position change occurs in a corresponding walrus-centered neighborhood with a certain radius. Considering that in the initial iterations of the algorithm, priority is given to global search in order to discover the optimal area in the search space, the radius of this neighborhood is considered variable so that it is first set at the highest value and then becomes smaller during the iterations of the algorithm. For this reason, local lower/upper bounds have been used in this phase of WaOA to create a variable radius with algorithm repetitions. For simulation of this phenomenon in WaOA, a neighborhood is assumed around each walrus, which first is generated a new position randomly in this neighborhood using (7) and (8). then if the value of the objective function is improved, this new position replaces the previous position according to (9).  www.nature.com/scientificreports/ where X P 3 i is the new generated position for the i th walrus based on the 3rd phase, x P 3 i,j is its j th dimension, F P 3 i is its objective function value, t is the iteration contour, lb j and ub j are the lower and upper bounds of the j th variable, respectively, lb t local,j and ub t local,j are local lower and local upper bounds allowable for the j th variable, respectively, to simulate local search in the neighborhood of the candidate solutions.
Repetition process, pseudocode, and flowchart of WaOA. After updating the walruses' position based on the implementation of the first, second, and third phases, the first WaOA iteration is completed, and new values are calculated for the position of the walruses and the objective functions. Update and improve candidate solutions is repeated based on the WaOA steps according to Eqs. (3)-(9) until the final iteration. Upon completion of the algorithm execution, WaOA introduces the best candidate solution found during execution as the solution to the given problem. The WaOA implementation flowchart is presented in Fig. 2, and its pseudocode is specified in Algorithm 1.

Input all information of optimization problem.
Create initial population. Set = 1 and = 1. Phase 2: Select randomly as the migration destination for using Eq. (5).

<
Save the best candidate solution so far.
Output the best quasi-optimal solution of the objective function found by WaOA.  , and TLBO has a computational complexity equal to to O(Nm(1 + 2T)) . Therefore, it is clear that the proposed WaOA approach has higher computational complexity than all algorithms used for comparison. However, to make a fair comparison, we used the population size of each metaheuristic algorithm in the simulation analysis so that the total number of function evaluations is the same for all employed algorithms.

Simulation studies and results
In this section, WaOA simulation studies on optimization applications are presented. The efficiency of WaOA in providing the optimal solution has been tested on sixty-eight standard objective functions, including unimodal, high-dimensional multimodal, fixed-dimensional multimodal, the CEC 2015 test suite, and the CEC 2017 test suite. The information on these test functions is specified in the Appendix and Tables A1 to A5.
The reasons for choosing these benchmark functions are as follows. Unimodal functions F1 to F7 are suitable for evaluating the exploitation ability of metaheuristic algorithms in convergence towards the global optimal as they do not have a local optimum. Multimodal functions F8 to F23 are suitable options for evaluating the exploration ability of metaheuristic algorithms due to having multiple local optimal. The CEC 2015 and the CEC 2017 test suites have complex benchmark functions that are suitable for evaluating the ability of metaheuristic algorithms to balance exploration and exploitation during the search process. WaOA performance is compared with ten well-known GA, PSO, GSA, TLBO, GWO, MVO, MPA, TSA, RSA, and WSO algorithms to determine the quality of WaOA results. The values set for the control parameters of the employed algorithms are specified in Table 1. The WaOA and mentioned competitor algorithms had been implemented on F1 to F23, each in twenty independent runs containing a thousand iterations (i.e., T = 1000 ). In this study, parameter N is considered equal to 20 for WaOA, 30 for TLBO, and 60 for other competitor algorithms to equalize the number of function evaluations. In this case, considering the computational complexity of each algorithm, the number of function evaluations for each metaheuristic algorithm is equal to 60,000.
Optimization results are reported using four statistical indicators: mean, best, standard deviation, and median. In addition, each algorithm's rank in handling each objective function is determined based on the average criterion.
Evaluation unimodal objective function. Unimodal  Evaluation high-dimensional multimodal objective functions. High dimensional multimodal functions with several local and globally optimal solutions have been selected to evaluate WaOA exploration capability in global search. The optimization results of F8 to F13 functions using WaOA and competitor algorithms are reported in Table 3. What can be deduced from the results of this table is that WaOA has converged to the global optimal in optimizing F9 and F11. WaOA is also the best optimizer for optimizing F10, F12, and F13. TSA is the best optimizer for the F8 objective function, while WaOA is the second-best optimizer for this objective function. Analysis of the simulation results shows that WaOA has an acceptable performance in optimizing high-dimensional multimodal objective functions and has provided a superior outcome compared to ten competitor algorithms.
Evaluation fixed-dimensional multimodal objective function. The fixed-dimensional multimodal functions, which have fewer local solutions than functions F8 to F13, have been selected to evaluate WaOA's ability to balance exploration and exploitation. The optimization results of F14 to F23 functions are reported in Table 4. The results show that WaOA ranks first as the best optimizer in handling all F14 to F23 functions. Furthermore, analysis of the simulation results shows the superiority of WaOA over ten compared algorithms due to the high power of WaOA in balancing exploration and exploitation. The performances of WaOA and competitor algorithms in solving F1 to F23 functions are presented as boxplot diagrams in Fig. 3. Intuitive analysis of these boxplots shows that the proposed WaOA approach has provided superior and more effective performance than competitor algorithms by providing better results in statistical indicators in most of the benchmark functions.
Statistical analysis. In this subsection, the superiority of WaOA over competitor algorithms is statistically analyzed to determine whether this superiority is significant or not. To perform statistical analysis on the obtained results, Wilcoxon signed-rank test 59 is utilized. Wilcoxon signed-rank test is a non-parametric test that is used to detect significant differences between two data samples. The results of statistical analysis using this test are presented in Table 5. What can be seen from the study of the simulation results is that WaOA has a significant statistical superiority over the competitor algorithm in cases where the p-value is less than 0.05. www.nature.com/scientificreports/ Sensitivity analysis. WaOA is a population-based optimizer that performs the optimization process in a repetitive-based calculation. Accordingly, the parameters N (the number of members of the population) and T (the total number of iterations of the algorithm) are expected to affect the WaOA optimization performance. Therefore, WaOA's sensitivity analysis to parameters T and N is presented in this subsection. For analyzing the sensitivity of WaOA to the parameter N , the proposed algorithm for different values of the parameter N equal to 20, 30, 50, and 100 is used to optimize the functions of F1 to F23. Optimization results are given in Table 6, and WaOA's convergence curves under this analysis are presented in Fig. 4. What is evident from the analysis of WaOA's sensitivity to the parameter N is that increasing the searcher agents improves WaOA's search capability in scanning the search space, which enhances the performance of the proposed algorithm and reduces the values of the objective function. www.nature.com/scientificreports/ For analyzing the sensitivity of the proposed algorithm to the parameter T , WaOA for different values of the parameter T equal to 200, 500, 800, and 1000 is used to optimize the functions of F1 to F23. Optimization results are in Table 7, and the WaOA's convergence curves under this analysis are presented in Fig. 5. Based on the obtained results, it is found that increasing values of T gives the algorithm more opportunity to converge to better solutions based on exploitation ability. Therefore, it can be seen that with increasing values of T , the optimization process has become more efficient, and as a result, the values of the objective function have decreased.

Evaluation of the CEC 2015 test suite. The optimization results of the CEC 2015 test suite, including
C15-F1 to C15-F15 using WaOA and competitor algorithms, are released in Table 8. The simulation results show that WaOA is the best optimizer for C15-F1 to C15-F8, C15-F10, C15-F13, and C15-F14 functions. In addition, in solving C15-F9 after MVO, in C15-F11 after WSO, C15-F12, and C15-F15 after GSA, the proposed WaOA is the second-best optimizer. Analysis of simulation results shows that WaOA provides better results in most functions of the CEC 2015 test suite, and in total, with the first rank of the best optimizer in handling the CEC 2015 test suite, has provided superior performance compared to competitor algorithms.  Table 9. What can be seen from the analysis of the simulation results is that WaOA is the first best optimizer for C17-F1 to C17-F6, C17-F8 to C17-F30 functions. in solving C17-F7, proposed WaOA after GSA is the second-best optimizer. Comparison of simulation results shows that WaOA has provided better results in most functions of CEC 2017 test suite, and has provided superior performance in solving this test suite compared to competing algorithms.

Informed consent.
Informed consent was not required as no human or animals were involved.
Ethical approval. This article does not contain any studies with human participants or animals performed by any of the authors.

WaOA's application to real-world problems
Metaheuristic algorithms are one of the most widely used techniques in dealing with real-world applications. This section tests WaOA performance in optimizing four engineering design challenges and twenty-two constrained optimization problems from the CEC 2011 test suite. It should be noted that to model the constraints of optimization problems, the penalty function has been used. Thus, if a solution does not meet any of the constraints of the problem, a penalty coefficient is added to the value of its objective function corresponding to each noncompliance of the constraint, and as a result, it is known as an inappropriate solution.
tension/compression spring design optimization problem. Tension/compression spring design is a challenge in real-world applications with the aim of minimizing the weight of tension/compression spring. A schematic of this design is shown in Fig. 6 59 . The tension/compression spring problem formulation is as follows: The results of using WaOA and competing algorithms in optimizing the Tension/compression spring design variables are presented in Table 10. The simulation results show that WaOA has provided the optimal solution to this problem with the values of the variables equal to (0.0519693, 0.363467, 10.9084) and the corresponding objective function value equal to 0.012672. The statistical results obtained from the performance of WaOA and competitor algorithms are reported in Table 11, which shows the superiority of WaOA in providing better values for statistical indicators. The WaOA convergence curve while achieving the solution for tension/compression spring is shown in Fig. 7.

Welded beam design.
Welded beam design is a real global challenge in engineering sciences whose main goal in design is to reduce the fabrication cost of the welded beam. A schematic of this design is shown in Fig. 8 60 .

With
WaOA and competing algorithms are implemented on the welded beam design problem, and the results are presented in Table 12. Based on these results, WaOA has provided the optimal solution to this problem with the values of the variables equal to (0.20573, 3.470489, 9.036624, 0.20573) and the corresponding objective   Table 13. This table shows that WaOA performs better in terms of statistical indicators. The convergence curve from the WaOA implementation on the welded beam design is shown in Fig. 9.

Pressure vessel design.
Pressure vessel design is a real-world optimization challenge that aims to reduce design costs. A schematic of this design is shown in Fig. 12 63 . The formulation of pressure vessel design problem is as follows:    Table 17. Statistical results indicate that WaOA has effectively optimized the pressure vessel design challenge by providing more favorable values for statistical indicators. The WaOA convergence curve in achieving the optimal solution is shown in Fig. 13.

Evaluation of twenty-two real-world optimization problems from the CEC 2011 test suite.
In this subsection, the performance of WaOA in handling real-world applications is challenged on twenty-two constrained optimization problems from the CEC 2011 test suite. This test suite has twenty-two optimization problems, namely: parameter estimation for frequency-modulated (FM) sound waves, the Lennard-Jones potential problem, the bifunctional catalyst blend optimal control problem, optimal control of a nonlinear stirred tank reactor, the Tersoff potential for model Si (B), the Tersoff potential for model Si (C), spread spectrum radar polyphase code design, transmission network expansion planning (TNEP) problem, large-scale transmission pricing problem, circular antenna array design problem, and the ELD problems (which consist of DED instance 1, DED instance 2, ELD Instance 1, ELD Instance 2, ELD Instance 3, ELD Instance 4, ELD Instance 5, hydrothermal scheduling instance 1, hydrothermal scheduling instance 2, and hydrothermal scheduling instance 3), the Messenger spacecraft trajectory optimization problem, and the Cassini 2 spacecraft trajectory optimization problem. Full details and description of the CEC 2011 test suite are available at 64 . The results of employing WaOA and competitor algorithms on these real-world optimization problems are presented in Table 18. The boxplot diagrams obtained from the performance of metaheuristic algorithms in handling CEC 2011 test suite optimization problems are drawn in Fig. 14. Based on the simulation results, WaOA is the first best optimizer to solve all C11-F1 to C11-F22 optimization problems. Based on the simulation results, the proposed WaOA approach has provided better results in most of the optimization problems and has provided superior performance in handling the CEC 2011 test suite in competition with competing algorithms. Also, the results obtained from the statistical analysis for p-value show that WaOA has a significant statistical superiority compared to competitor algorithms.

Conclusions and future works
In this study, a new bio-inspired metaheuristic algorithm called the Walrus Optimization Algorithm (WaOA) was developed based on the natural behaviors of walruses. Feeding, escaping, fighting predators, and migrating are the primary sources of inspiration used in the design of WaOA. Therefore, the WaOA theory was explained, and its mathematical modeling was presented in three phases: (i) feeding strategy, (ii) migration, and (iii) escaping and fighting against predators. Sixty-eight standard benchmark functions of various types of unimodal, multimodal, the CEC 2015 test suite, and the CEC 2017 test suite, were employed to analyze WaOA performance in providing solutions. The information on these test functions is specified in the Appendix and Tables A1 to A5. The optimization results of unimodal functions showed the high ability of WaOA exploitation in local search to converge towards global optimal. The optimization results of multimodal functions indicated the high ability of WaOA exploration in global search and not to be caught in locally optimal solutions. WaOA's performance results were compared with the ten well-known metaheuristic algorithms. The simulation and comparison results showed that the proposed WaOA approach has a high ability to balance exploration and exploitation and is much 0 ≤ x 1 , x 2 ≤ 100, and 10 ≤ x 3 , x 4 ≤ 200.       www.nature.com/scientificreports/ superior and more competitive against ten competitor metaheuristic algorithms. In addition, the results of the WaOA implementation in addressing the four design issues and twenty-two real-world optimization problems from the CEC 2011 test suite demonstrates the effectiveness of the proposed approach in real-world applications.
Although it was observed that WaOA had provided superior results in most of the benchmark functions, the proposed approach has some limitations. The first limitation facing all metaheuristic algorithms is that it is always possible to design newer algorithms that can provide better results than existing algorithms. The second limitation of WaOA is that the proposed method may fail in some optimization applications. The third limitation of WaOA is that the nature of random search in this algorithm leads to the fact that there is no guarantee of achieving the global optimum. Moreover, the authors do not claim that the proposed WaOA approach is the

Data availability
All data generated or analyzed during this study are included directly in the text of this submitted manuscript. There are no additional external files with datasets.