1. Introduction
Network Control Systems (NCS) have become increasingly popular in mining, process, and manufacturing industries that produce various goods such as food and beverages, different metals, multiple chemicals, pulp and paper, automobiles, textiles, crude oil refineries, and power generation plants [
1]. NCS monitors and controls field instruments effectively, ensuring efficient and accurate operations [
2]. In recent times, WMN has gained widespread acceptance in the industry, particularly with the wide acceptance of different wireless networking protocols, namely WIA-PA, Zigbee, WirelessHART, Bluetooth, and ISA 100 wireless. One of the primary advantages of these wireless standards is eliminating clunky cabling—a costly and time-consuming process to repair and maintain [
3,
4]. Additionally, WMN enables the expansion of network capability to the regions where it is difficult or impossible to install the wired cables, thus making monitoring and controlling operations in remote locations possible. Moreover, using WMNs increases data reliability, reduces data loss or interference, and offers improved security measures, protecting sensitive data from unauthorized access or theft [
5]. It is essential to have a thorough knowledge of several factors such as transmission power, network topology, and wireless node density. Ignoring these factors while deploying WMN routers can lead to inadequate client coverage and connectivity, substandard network range, control-loop failures, and transmission loss [
6]. It is crucial to consider the technical constraints of the deployment location and the basic topology before deploying wireless mesh routers. These constraints include factors such as radio signal strength, frequency range, and interference from other wireless networks [
7,
8].
WMNs offer many benefits for all the process control and automation industries [
9]. However, some issues must be addressed to enhance the network’s performance, including network coverage, router connectivity, network traffic, compatibility, data security, etc. While comparing the outdoor regions, the indoor area is challenging while implementing the WMN for process monitoring and control because of the stochastic interferences [
10]. In some outdoor areas, installations require fewer clients than in indoor regions. Furthermore, WMN can be implemented in multiple places with the availability of various control services in indoor and outdoor conditions. Determining the appropriate location for installing WMN routers for real-world applications is crucial. Therefore, optimization is essential for ensuring reliable and efficient WMN performance [
11]. Optimization involves fine-tuning various parameters, such as router placement, transmission power, network topology, and routing protocols, to achieve the best possible network performance [
12]. Optimizing these parameters allows the network to ensure that the WMN provides seamless coverage, reliable connectivity, and low latency, even in challenging industrial environments. Additionally, optimizing industrial WMNs can help improve network scalability and reduce maintenance costs [
13].
At the same time, a feedback controller is essential in WMNs to ensure the system’s stability. In addition, the feedback controllers will help the system to maintain stable and accurate control over the process variables [
14]. However, to achieve optimal performance, the controller must be optimized based on the system characteristics, process model, and external disturbances of the system it controls. Optimization techniques such as trial and error, Ziegler–Nichols, and the Cohen–Coon methods are commonly used to optimize feedback controllers. These conventional techniques involve adjusting the controller’s parameters to find the optimal values that result in suboptimal, unstable, and ineffective control. Thus, many researchers proposed using metaheuristic-based techniques to optimize the controller parameters for different applications [
15,
16,
17,
18,
19,
20].
The population-based evolutionary algorithms use mutations and crossovers, which produce new solutions using nature-inspired concepts from evolution. Swarm-based optimization techniques also employ these procedures [
21]. However, they utilize two core processes for organizing initial swarms: scattering and intensifying them to specific high-potential regions [
22]. Simultaneously, it should be noted that the No Free Lunch theorem affirms that no such existing metaheuristic optimization technique can be employed to solve every category of the existing issues in simulation or real-time conditions. In addition, this theorem highlights the need for a thorough understanding of the issue and the appropriate optimization technique selection to achieve optimal results [
23]. Thus, a brief review of the Harris Hawk Optimization (HHO) and Arithmetic Optimization Algorithm (AOA) for different types of applications was carried out and is presented in the following sections for identifying the research gap for the proposed methodology.
1.1. Harris Hawk Optimization
Initial development of the Harris Hawk Optimization, a swarm-based optimization technique, was proposed by Heidari et al. [
24]. This method is motivated by the cooperation and hunting strategies of Harris hawks in nature, which involve team action and reaction, as well as prey evasion. The main objective of HHO is to discover solutions to single-objective problems. The HHO method uses hawks’ chasing actions as search agents and the prey as the best position. This approach has the potential to solve a variety of real-world optimization problems, such as engineering design, pattern and speech recognition, manufacturing optimization, bio-mechanical engineering, power production quality, feature selection, and medical image segmentation [
25,
26,
27,
28,
29]. As with other metaheuristic techniques, HHO uses exploration and exploitation phases to identify and catch the prey. However, the HHO differs from others in the exploitation phase, where it catches the prey based on soft and hard besiege and rapid dives on soft and hard besiege movements. Detailed information on the HHO algorithm can be found in these articles [
24,
25,
28]. As a result of the benefits provided by HHO, there has been a significant increase in HHO research, with the number of studies continuing to grow rapidly to numerous applications and creating new HHO variations.
The authors of [
30] aimed to accelerate the convergence of the HHO by proposing an improved version of the algorithm that employs two strategies. The first strategy involves improving the exploration phase of HHO by incorporating opposition-based learning and logarithmic spiral techniques. The second strategy combines the Modify Rosenbrock method to enhance the HHO’s local searching capability and convergence rate. The authors tested their algorithm by conducting experiments on 23 standard benchmark functions. They also compared their results with various state-of-the-art metaheuristic algorithms. In [
31], two different variants of the HHO, namely binary HHO and quadratic binary HHO, were used to improve the feature selection problem at various data classification conditions. Binary HHO has a transfer function that can be either S-shaped or V-shaped, effectively converting a continuous variable into a binary one. The quadratic binary HHO is designed to improve the binary HHO by incorporating a binary quadratic model.
Dokeroglu et al. [
32] developed a binary variant of the multi-objective HHO approach to address a classification task. Their novel discrete operators for exploitation (besieging) and exploration (perching) were introduced. Their study employed four machine learning methods applied to a COVID-19 dataset for binary classification problems to increase prediction accuracy. Reference [
33] proposed a novel optimization algorithm, called the modified HHO, to reconfigure photovoltaic modules and disperse shadow regularly to increase the power generated. The algorithm was tested on different shade patterns with 9 × 9, 6 × 4, and 6 × 20 photovoltaic arrays based on the available simulation data. The outcomes were analyzed and compared with varying techniques of reconfiguration based on metrics such as mismatch power, fill factor, power loss, and power enhancement. Hussain et al. designed an improved variant of the HHO named long-term memory HHO that enhances the diversity of search agents by maintaining a broad search region for the exploration phase resulting in better convergence results. The proposed technique uses different engineering optimization test problems, including power flow optimization for the power generation systems [
34].
The improved HHO algorithm, as introduced by Houssein et al. [
35], addresses the challenge of optimizing large-scale wireless sensor networks. The approach involves utilizing Prim’s shortest path algorithm to establish minimum transmission paths from the sink node to all other sensor nodes in order to reconstruct the network. The researchers have proposed energy-efficient wireless sensor networks using a hybrid Harris hawk-salp swarm optimization algorithm in [
36]. The approach uses a hybrid optimization algorithm combining the HHO and salp swarm techniques to optimize the network’s energy consumption. Additionally, a mobile sink strategy reduces energy consumption by allowing the destination node to move toward the sensor nodes instead of the other way around. For PID controller optimization, researchers used HHO to control the speed of the DC motor [
37]. Using the HHO algorithm for tuning the controllers improves performance in terms of steady-state error, rise time, overshoot, and settling time, resulting in better efficiency, accuracy, and stability of the motor control system. An overall detailed review of the HHO literature is given in
Table 1.
1.2. Arithmetic Optimization Algorithm
The AOA is a recently designed mathematics-inspired metaheuristic optimization technique that draws inspiration from the distribution behavioral nature of arithmetic functions from the mathematics proposed by Abualigah et al. [
40]. It uses basic arithmetic operations, such as addition, subtraction, multiplication, and division, as the primary operator to explore the desired objective function and the optimal solution. AOA is also a population-based optimization technique, a set of solutions, called individuals, are generated randomly and evaluated based on their fitness concerning the considered problem. The best-fit individuals are chosen and subjected to various mathematical operations to generate new individuals based on the desired objective values in both phases. These new individuals are then evaluated for their fitness and the process is repeated until the desired optimization level is reached [
41].
A key advantage of AOA is that it can prevent getting stuck in local minima, a common issue many optimization algorithms face. AOA achieves this by employing a mechanism called elitism, which ensures that the best solutions in the population are always preserved and passed over to the next generation, allowing the algorithm to explore different areas of the solution and converge to the global minimum efficiently. Premkumar et al. [
42] proposed a newly developed approach called Multi-Objective AOA (MOAOA) to solve real-world constrained multi-objective optimization problems. The latest version of the AOA algorithm has been improved by adding two new functions—elitist non-dominance sorting and a crowding distance-based mechanism. This new algorithm has been tested on 35 optimization problems with constraints and five test problems without constraints. Its effectiveness has been compared with four other advanced multi-objective optimization techniques. Additional performance metrics are also considered to evaluate the MOAOA for higher accuracy and effective convergence, in which MOAOA achieved superior performance.
In [
43,
44,
45], the researchers developed an enhanced form of artificial neural network which uses a two-step process to identify, locate, and measure the extent of damage in plate structures made of functionally graded materials, petroleum products, and other electrical systems. In the first stage, a damage indicator based on the frequency response function is utilized to forecast which components of the material have been affected. In the second stage, the networks are employed to quantify the damage and it eliminates healthy elements from the numerical model and utilizes information from the defective components to estimate the degree of damage. Abualigah et al. [
46] proposed a new method for various multilevel thresholding in data analysis and image segmentation using the AOA and differential evolution technique. The DAOA algorithm was assessed using traditional examination images from nature and CT COVID-19 images. The accuracy of the segmented images was measured using the peak signal-to-noise ratio and structural similarity index test.
Zheng et al. [
47] developed a new hybrid algorithm called DESMAOA, combining two meta-heuristic algorithms—the Slime Mold Algorithm (SMA) and the AOA. The DESMAOA algorithm was designed to improve the optimization capability of the existing techniques. The SMAOA algorithm is first used to improve the SMA algorithm, and then two strategies from SMA and AOA are combined to create the DESMAOA. In [
48], a new method is developed to control an automatic voltage regulator using a robust Model Predictive Controller (MPC). The suggested approach was crafted to address the challenge of uncertain automatic voltage regulation parameters. Frequency domain conditions are obtained through the Hermite–Biehler theorem to ensure stability in the face of perturbations. The tuning of MPC parameters is performed using the AOA technique while accommodating stability constraints. A time-domain objective is established to optimize the voltage regulator’s performance by minimizing voltage peak overshoot and fastening the process settling.
In [
49], the authors introduced a forced switching mechanism for AOA termed IAOA that helps search agents to switch between various local optima. The effectiveness of the IAOA was tested on multiple benchmark functions and real-world test problems and the results showed that it outperforms other optimization algorithms in most cases. Different types of metaheuristic optimization were compared for the various engineering problems in [
50]. The researchers proposed an enhanced version of the AOA called nAOA. The nAOA technique uses logarithmic and exponential mathematical operators to improve algorithm performance. The overall summary of the AOA technique is given in
Table 2.
Some of the main contributions of this research article are listed as follows:
The primary motivation of HHAOA is to solve the problem of placing the routers optimally to achieve adequate coverage and connectivity, along with finding the best parameter for the Fractional-Order Predictive PI (FOPPI) controller.
The proposed HHAOA technique makes use of all the arithmetic operator combinations, enabling smooth and efficient relocations between local minima without getting stuck. This approach effectively minimizes computational complexity, resulting in a more streamlined process outcome.
To achieve a better convergence rate and reach the desired optimal solution location, the proposed optimization was simulated and validated using 33 different benchmark functions.
The proposed technique was implemented to enhance nodes’ placement and mitigate congestion in wireless mesh networks (WMNs). The results reveal that the HHAOA technique is highly effective in reducing the number of node placements, leading to significant cost savings and improved network congestion compared to alternative algorithms.
The HHAOA-optimized FOPPI controller was implemented on a real-time pressure process plant to validate our proposed optimization algorithm, and the results show a better performance than conventional controllers.
2. Proposed Hybrid Harris Hawks-Arithmetic Optimization Algorithm
The primary focus of this study was developing a hybrid HHAOA technique to enhance the coverage and connectivity of a WMN, minimize network congestion, and optimize parameters for the FOPPI controller [
54]. The proposed technique hybridizes the HHO and AOA techniques. Despite the impressive capabilities of HHO and AOA, some drawbacks need to be addressed, including the risk of premature convergence, becoming trapped in multiple local optima, and its phase-switching mechanism. The proposed HHAOA approach improves convergence behavior, the position-switching mechanism, and solution quality. When implementing the hybrid method for searching, the process becomes notably more thorough and effective. This is due to the ability to navigate throughout the desired search area and avoid becoming trapped in local optima. As a result, a diverse range of potential solutions can be generated, increasing the likelihood of finding the optimal outcome. Furthermore, the HHAOA will be used to optimize the placement of WMN routers via simulation. Additionally, real-time experiment experimentation on the pressure process plant using the HHAOA-optimized FOPPI controller was compared with traditional techniques.
2.1. Proposed Hierarchical Structure
The proposed HHAOA hierarchical structure is depicted in
Figure 1. The system comprises a top layer (primary layer) with
M HHO search agents and a bottom layer (secondary layer) with groups containing an
N AOA population. The AOA execution in the bottom layer initiates the process of updating the search agents’ positions. In order to find the most optimal solution, it is crucial to ensure that the positions of all search agents in the upper layer are updated with the best solution discovered by the corresponding group in the lower layer. New equations can be formulated for both the exploitation and exploration phases, leading to an even more effective solution. This approach allows for a comprehensive problem analysis and significantly improves the overall outcome.
2.2. Exploration
This section proposes an exploration mechanism for hybridizing the HHO and AOA algorithms. The algorithm draws inspiration from the hunting behavior of Harris hawks, capitalizing on their exceptional eyesight to effectively track and recognize prey. However, the prey may not be immediately visible in some cases, prompting the hawks to patiently wait, observe, and monitor the objective location for several hours until detecting a potential target.
The algorithm expertly utilizes Harris hawks to represent potential solutions, with the ultimate goal of comparing the best solution to the ideal outcome. The algorithm flawlessly imitates the hawks’ actions by strategically placing them in different locations and using various approaches until the desired solution is found. Each strategy is meticulously employed in equal measure to ensure maximum efficiency. The hawks rest closer to the prey (rabbit) based on the positions of other family members during
. When the condition is met with
, the hawks will randomly perch on tall trees within the range of their family members. Instead of selecting the random tall trees, the best position nearest to the rabbit will be identified using the AOA algorithm, which will further enhance the chances of reaching and identifying the target much faster. Therefore, the exploration strategy of the proposed HHAOA is obtained using the equation given below.
where
describes the location of
t-th solution in the top layer (HHO) corresponding to the
j-th search solution in the AOA layer.
j represents the current iteration of the algorithm.
,
,
,
,
c, and
q are the random numbers that lies in the span of [0,1].
,
,
,
,
, MOP are the random hawk selected, the best location obtained so far in the present iteration, average mean of the hawk’s position, upper and lower bound ranges of the variables, and the math-optimizer probability coefficient. The MOP and average mean of the hawk’s location values are calculated based on the equation given below.
where
is a crucial parameter factor that determines the precision of the exploitation phase,
t, and
T represent the current and the highest possible number of iterations, respectively. The position of every hawk in iteration j is represented by
, and the total number of hawks is denoted by
N.
2.3. Transition Stage
The HHO algorithm can alternate between exploration and exploitation modes depending on the amount of energy the prey has left while attempting to escape. The prey’s energy decreases notably during its escape. This energy is determined by using the equation below:
The escaping energy of the prey is represented as E and the maximum number of iterations is represented by T with the current iteration being represented as t. Additionally, indicates the initial energy state of the prey, which randomly changes within the interval of (−1,1) at each iteration in HHO. When decreases from 0 to −1, the prey is physically exhausted, while an increase in from 0 to 1 demonstrates that the prey is becoming more powerful. During the iterations, the escaping energy (E) gradually decreases. When the escaping energy is |E| , the hawks explore various locations to search for the prey, thereby performing the exploration phase. Conversely, when |E| < 1, the algorithm focuses on exploiting the neighbourhood of the solutions.
2.4. Exploitation
During the attacking phase, Harris hawks employ a surprise pounce maneuver to take down their intended prey, which is identified in the previous stage. However, the prey may attempt to evade danger, resulting in various chasing styles observed in actual real-life scenarios.
Based on the prey’s fleeing behavior and the Harris hawks’ pursuit strategies, the HHO proposes four potential approaches to simulate the attacking phase. As preys always try to escape threatening situations, they have a chance (r) of successfully escaping () or not () before the surprise pounce. Hawks employ a hunting technique known as hard or soft encirclement to capture their prey which involves surrounding the prey from various angles, considering its energy level. The hawks gradually move closer to their prey to increase their chances of a successful surprise attack in practical settings. As the fleeing prey loses energy over time, the hawks intensify the besiege process to capture the exhausted prey effortlessly. Parameter E models this strategy and enables the HHO to alternate between soft and hard besiege processes. Similar to the exploration, the AOA will reach the solution (prey) closer in this phase.
2.4.1. Soft Besiege
The rabbit retains sufficient energy to flee through random and deceptive jumps during
and |E|
. However, its efforts prove to be unsuccessful, and the Harris hawks gradually encircle it, causing the rabbit to become increasingly tired before ultimately launching a surprise attack. The following Equation (
5) represents this surprise attack of HHAOA:
where
.
2.4.2. Soft Besiege with Progressive Rapid Dives
When the rabbit’s energy is sufficient (|E|
) for successful escape and the hawk’s movement is rapid (r < 0.5), a gentle ambush is carried out before the surprise attack. This approach is more intelligent than the previous method. To mathematically model the prey’s escape patterns and the predator’s leapfrog movements, the HHAOA algorithm utilizes the idea of incorporating the Levy flight (LF) pattern. LF simulates the erratic zigzag movements of prey (especially rabbits) in the fleeing phase and hawk’s irregular, sudden, and rapid dives as they approach the targeted prey. The LF mechanism used here is given below.
where
and
are random numbers that lie in the range of inside [0,1],
is a constant value that is set at 1.5.
During the soft besiege, hawks perform multiple rapid dives around the rabbit, constantly adjusting their position and direction to match the deceptive movements of the prey. This mechanism is also observed in other competitive situations in nature. The movement of the hawks in this phase is obtained using the following position update rule given in Equation (
7).
where
S = Random vector in the dimension D .
2.4.3. Hard Besiege
When the conditions meet
and |E|
, the rabbit will become exhausted and have less energy to flee. As a result, the Harris hawks will attempt to capture the rabbit by surrounding it aggressively and launching surprise attacks. The following equation can describe this intense attack phase.
2.4.4. Hard Besiege with Progressive Rapid Dives
If |E|
and
r , the rabbit lacks the necessary energy to flee and must resort to a sudden, aggressive dive technique to catch its prey off guard. This maneuver places the prey in a position similar to a soft dive. Still, the hawks work to minimize the distance between themselves and the prey’s average location while attempting to escape. The movement of this phase is obtained as follows:
where,
The HHAOA proposal consists of two distinct phases: exploration and exploitation. The initial phase is responsible for finding a new position to enhance the optimization process. At the same time, the latter stage utilizes four different strategies to maximize the position update to converge effectively. Detailed information about when to switch between these two phases is presented in
Table 3.
2.5. Pseudocode of Proposed Algorithm
The application of HHAOA for the industrial wireless mesh networks (WMN) and FOPPI controller parameter optimization is comprehensively explained in the pseudocode for the proposed technique presented in Algorithm 1. The respective implementation of the HHAOA using the flowchart is illustrated in
Figure 2.
2.6. Algorithm Complexity
It is crucial to note that the complexity of the proposed HHAOA algorithm is heavily influenced by three critical processes: initialization of the transition stage, fitness function evaluation, and position updating of Hawks. The initialization process has a complexity of C(n) when there are n hawks. Updating the location vector of all hawks involves searching for the optimal location, which can be complex. This involves C(T × n) + C(T × n × D) operations, where T represents the maximum number of iterations and D is the dimension of the current particular problem in hand. Therefore, it is imperative to understand that the computational complexity of HHAOA can be expressed as C(n × (T + TD + 1)).
Algorithm 1 Pseudocode of HHAOA technique. |
- Input:
Random WMN routers positions and manually calculated FOPPI controller parameters. - Output:
Optimal WMN connection and FOPPI parameters - 1:
Initialize the search agents positions (j = 1, 2, 3, ..., N) - 2:
Check search space boundary and initiate the case - 3:
while ( T) do - 4:
Attain the initial solution - 5:
for (All hawks ()) do - 6:
Amend the position update - 7:
- 8:
Phase: Exploration - 9:
if & then - 10:
if then - 11:
Update the solution using Equation ( 1) condition 1 - 12:
else if then - 13:
Update the solution using Equation ( 1) condition 2 - 14:
end if - 15:
else if & then - 16:
Update the solution using Equation ( 1) condition 3 - 17:
end if - 18:
Phase: Exploitation - 19:
if < 1 then - 20:
if ( 0.5 and r ≥ 0.5) then - 21:
if then - 22:
Solution update using Equation ( 5) condition 1 - 23:
else if then - 24:
Solution update using Equation ( 5) condition 2 - 25:
end if - 26:
else if ( 0.5 and r < 0.5) then - 27:
if or then - 28:
Solution update using Equation ( 7) condition 1 - 29:
else if or then - 30:
Solution update using Equation ( 7) condition 2 - 31:
end if - 32:
else if ( & r ≥ 0.5) then - 33:
if then - 34:
Solution update using Equation ( 8) condition 1 - 35:
else if then - 36:
Solution update using Equation ( 8) condition 2 - 37:
end if - 38:
else if ( 0.5 & r < 0.5) then - 39:
if or then - 40:
Solution update using Equation ( 9) condition 1 - 41:
else if or then - 42:
Solution update using Equation ( 9) condition 2 - 43:
end if - 44:
end if - 45:
end if - 46:
end for - 47:
t = t + 1 - 48:
end while - 49:
Return for optimal WMN connectivity with all the routers and optimal FOPPI controller parameters for real-time pressure process plant.
|
4. Results and Discussion
This section presents the benchmark functions and wireless mesh network simulation analysis in the first and second subsections. Lastly, the HHAOA-optimized FOPPI controller was experimented on with the pressure process plant. During the simulation of the benchmark functions, the Hawks population was kept at 100, and ran for 300 iterations, while in the WMN the simulation ran for 500 iterations, keeping the Hawks population constant with dimensions of 1000 m by 1000 m. The experimentation and simulation were conducted in the Intel(R) Xeon PC 3.10 GHz and 16.00 GB RAM using the MATLAB/Simulink software (2021a). Additionally, extra parameters for the algorithms are listed as min = 0.2, max = 1.0, convergence constant , initial and escaping energies = [0 1] and E = [−1 1], = 5, = 1.5, and .
4.1. Performance Analysis on Benchmark Functions
A simulation analysis for 30 benchmark functions (refer to
Table 4 and
Table 5) using AOA, MFO, SCA, GWO, WOA, HHO, and the proposed HHAOA was conducted. The selected benchmark functions were independent of each other in multiple dimensions and had various local minima, global minima (
), different ranges, and diverse boundary values to ensure a reliable comparative analysis. Furthermore, the algorithms were comprehensively tested using multimodal and hybrid composition functions in single and multiple dimensions.
4.1.1. Benchmark Functions
The algorithm’s accuracy was tested using unimodal functions (F1 to F7). As for multimodal functions, they contain multiple local minimum points, which serve as a measure of the algorithm’s exploration capacity. This capacity pertains to its ability to move seamlessly from local to global minima without getting stuck in a single position.
Table 4 shows the F8 to F13 multimodal benchmark functions that have many local minima. The F14 to F21 functions also have multiple local minima with fixed dimensions. These functions help to evaluate the stability of optimization algorithms. It is necessary to mention that the first 13 benchmark functions have higher dimension values of 30, 100, 500, and 1000 to evaluate the proposed HHAOA.
The remaining functions, from F22 to F30, are hybrid functions as detailed in
Table 5. While working with unimodal functions, finding the optimal solutions was comparatively effortless as they were conveniently located and easily reachable from the current location. The likelihood of being trapped in a particular search area is minimal. Finding the optimal position and movement between hybrid functions F19 and F33 is very difficult because there are many local minimum values to consider in the search space.
Figure 6 displays the surface plots for all functions, providing a clear visualization of the search space needed to identify the global minimum for both single and multi-minima functions. Functions can be classified into different categories based on their surface area and shape. These categories include bowl-shaped, plate-shaped, and valley-shaped functions. There are also functions with single or multiple local minima. The latter type of function has a broad search space, wide range, numerous layers, and multi-dimensional features due to its multiple local minima values.
4.1.2. Convergence Analysis
The benchmark procedures were evaluated through 300 iterations, with 100 for the loop condition, to determine their convergence properties. Convergence, in this context, indicates the point at which an algorithm locates the smallest possible fitness value within the specified number of maximum iterations.
Figure 7 illustrates the convergence analysis plot for all the benchmark functions in the presence of all the optimization algorithms.
In the unimodal functions (F1–F4), the conventional HHO had the fastest convergence in reaching the global minima. However, in the remaining functions, the proposed HHAOA relocated the positions faster, reaching the objective function in the first place. In this performance, the MFO had the major setback of exploring desired search regions, leading to inadequate performance and securing the last position. In the multimodal functions (F8–F13), using the hawk’s population effectively, the proposed HHAOA had a better transition phase from exploration to exploitation while noticing the convergence speed, which has a more significant speed of reaching the global minima with fewer iterations. In addition, the amplitude reduction rate is directly proportional to its convergence speed in the HHAOA performance in these functions. These results, in turn, illustrate the ability of HHAOA to focus more on the desired search locations during the iterations.
It is important to note that, in F9, many algorithms were unable to reach the desired value because of the existence of multiple minimum values. (see
Figure 6). Still, the HHAOA explicitly performed well and had the best and faster convergence near the global minima. SCA and MFO had the least convergence speed in these functions, followed by GWO, WOA, and AOA, respectively. Notably, the conventional HHO had performance comparable to that of the proposed HHAOA but could not outperform it because of the premature convergence towards reaching the objective function. In the fixed-dimension multimodal functions (F14–F21), the proposed HHAOA continued the first position by effectively finding the best minimum value. In these functions, AOA has the slowest convergence rate. Surprisingly, the HHO had the least performance and lost the ability to relocate the position to find the desired search regions within the desired number of iterations.
As the number of dimensions increases, the quality of outcomes and the efficiency of alternative approaches deteriorate noticeably. This indicates that HHO can effectively uphold a favorable equilibrium between exploration and exploitation behaviors in scenarios involving multiple variables. The performances in the hybrid functions (F22–F30) have a similar results trend. Notably, in F24, MFO has the best convergence performance, followed by the HHAOA. The proposed HHAOA performs best in non-zero global minima functions by converging nearer to the global minima with fewer iterations, even at less than 50 iterations in most of the functions. Lastly, the proposed HHAOA can deliver exceptional outcomes across all dimensions and consistently outperforms other methods when dealing with problems involving many variables.
4.1.3. Quantitative Analysis
In this section, the efficiency of each algorithm was evaluated by measuring how closely the statistical data match the global minimum of the benchmark function being tested. HHAOA performs better than other algorithms, as it has significant performance outcomes. However, this also means that HHAOA can identify a better optimal solution with fewer iterations and is closer to the objective function, as evidenced by its lower mean value. In contrast, a lower standard deviation (Stdv) signifies more excellent convergence stability, effectiveness, and reliability. Consequently, HHAOA is capable of avoiding local optima with great success.
Table 6 show the numerical comparison for all the optimization techniques. Here, the abbreviated terms are
—Benchmark Functions,
—Global minima, and
—Standard deviation.
The unimodal function results show that the conventional HHO had a lead in reaching the global minima value more than the other algorithms. Surprisingly, GWO and WOA have an average performance compared to the newly developed AOA and SCA algorithms. The proposed HHAOA has notable performances in F4–F7 due to the arithmetic operators’ and hawks’ ability to narrow down the desired global minima and avoid further searching for the exploration locations. Meanwhile, the MFO performs poorly in these functions, which shows its inability to reach even for the singular objective functions.
The HHAOA method is able to reach the global minimum in both multimodal and fixed-dimension functions better than other techniques. HHO reached the desired value in most functions, but HHAOA converged at the exact minimum values in many cases. It’s important to note that almost all algorithms converged at the desired global minimum values in F9 and F16. The MFO competed with the proposed HHAOA in mean and best value in the fixed dimension functions. At the same time, the well-performed GWO and WOA in the unimodal function have a significant setback in this multi-objective function. They struggled to reach the best value in most of the and had an almost massive difference in the final convergence values than the desired global minima.
Lastly, the same results’ trend has been repeated in the hybrid functions. Notably, the HHO and SCA had poor performance and these algorithms quickly progressed in the first few iterations, but then the rate of improvement slowed down. The GWO and MFO had a closer convergence value in most functions, showing the same approach to finding the best optimum values (see
Figure 7), whereas the proposed HHAOA had a stable and smooth value closer to the global minima, it suggests that the algorithm converges stably. However, in F23, the HHAOA has a lower standard deviation, indicating that the algorithm consistently finds a better solution for a stable convergence. Overall, these statistical measurements prove that the convergence behavior of the HHAOA helps to identify effective search locations for improving the exploration and exploitation phases of the algorithm.
The Friedman ranking test, which compares algorithms based on their best mean values in numerical analysis, is displayed in
Table 7. In this comparative analysis, the Friedman ranking test is used to evaluate the performance of various metaheuristic approaches. The algorithms were ranked according to their mean value for the corresponding functions and the resulting list clearly compares each approach’s effectiveness. The top-performing algorithm is awarded a rank of 1 as it is the closest to the global minimum. On the other hand, the algorithm with the lowest mean value is assigned a rank of 5, indicating its significant deviation from the desired global minimum value. In this analysis, the proposed HHAOA has the smallest average value of 1.967, securing the first rank, which signifies the technique’s effectiveness in identifying the best optimum value with fewer iterations in most benchmark functions. In addition, the HHAOA has a massive 181.342% faster and best mean value compared with the SCA, which has a final average value of 5.534 with the last rank. The HHAOA has 35.587%, 64.412%, and 79.664% of increased performance than the HHO, GWO, and respectively. Lastly, the MFO has an average of 4.2 and secures fifth rank, followed by AOA with the value of 5.067, securing the second-to-last rank.
4.2. Simulation Analysis of Industrial WMNs
Understanding the convergence behavior of optimization algorithms is crucial for obtaining optimal solutions. In the performance evaluation of the proposed HHAOA algorithm against other algorithms, the convergence behavior was analyzed using 500 iterations with 100 runs for 100 search agent sets, shown in
Figure 8.
Table 8 shows the statistical analysis of the convergence. By observing the convergence ability of these algorithms, their performance can be determined clearly. The objective function was utilized to find the minimum fitness value, which demonstrates improved connectivity and coverage.
In the numerical analysis, the HHAOA method yielded the best outcomes with the least mean value of 0.499. Comparatively, the AOA came in second place with a mean value of 0.531. While the HHAOA method is 6.412% faster than AOA, it is worth noting that the WOA method is lagging with a minimal difference of 0.007. This resulted in WOA reaching a third place with a mean value of 0.456, which is 8.617% higher than the HHAOA method. The GWO, HHO, SCA, and MFO methods were followed in the remaining ranking order, which makes it clear that the proposed HHAOA method produced better results in the numerical analysis.
Figure 9 shows the network connectivity and coverage region for the WMN for all the algorithms and the initial WMN topology.
Figure 9a illustrates that the initial iteration resulted in insufficient mesh client coverage and high network congestion due to multiple overlapping mesh routers, which increased deployment costs. Even with redundant mesh routers, their connectivity did not reach the maximum capacity of 100%.
Figure 9b–h shows the optimized WMN connectivity and coverage of the AOA, MFO, SCA, GWO, WOA, HHO, and the proposed HHAOA, respectively. In these, SCA has the maximum number of disconnected clients and more number of overlapped router placements. It has a 150% less efficient connectivity than the HHO, which secured second place, with four unconnected clients in its optimized WMN. Additionally, the AOA has eight disconnected clients, followed by MFO, WOA, and GWO, with seven, five, and four disconnected clients in the WMN network. However,
Figure 9h shows that the proposed HHAOA has dramatically improved the client coverage and optimally deployed the mesh routers to achieve full network connectivity with the least number of routers and had only one disconnected client. Likewise, the proposed HHAOA significantly reduced network congestion by reducing the number of mesh routers by 31.7% while maintaining high coverage and connectivity. Furthermore, other comparing algorithms produced network topologies where mesh routers overlapped, resulting in more significant interference. In contrast, the topology formed by the HHAOA approach was more widely dispersed, resulting in improved client coverage.
4.3. Performance Evaluation on Pressure Process Control
In the numerical analysis, the comparison was carried out for the process rise time (
), settling time before (
) disturbance, overshoot (
), and settling time after (
) disturbance injection. In order to evaluate the ability of the FOPPI to decrease the stochastic disturbance and to track the set-point effectively, an external disturbance of 35% was injected at 100 s in the process feedback loop.
Table 9 shows the quantitative analysis of the optimized FOPPI controller using different optimization techniques, and the comparative analysis is shown in
Figure 10. Here, controller parameters
,
and
were obtained based on the ITAE value and integral time was obtained by
, where
is the integral gain.
The numerical results from the table show that the MFO has the fastest rise time of 0.7552 s, followed by the AOA, HHO, WOA, and SCA with values of 0.8404, 0.8411, and 0.9381 s, respectively. Here the GWO has the slowest rise time of 5.1078 s which is 576.351% slower performance than the fastest MFO. The proposed HHAOA optimized FOPPI had the second-last rise time of 1.4504 s which is almost two times slower than the MFO.
During the settling time before disturbance (), even with the slowest rise time, the HHAOA-optimized FOPPI managed to settle faster at 37.2265 s which is 29.7129 s faster than the slowest settled SCA with the settling at 66.9394 s. The result shows that the proposed HHAOA has 79.816% faster performance, which is twice as fast as SCA. The HHO comes second fastest, settling at 49.2073 s, followed by GWO, WOA, MFO, and AOA with the respective values of 52.9052, 54.2561, 57.5289, and 59.6059 s. While observing the settling time after disturbance injection performance, once again, the HHAOA settled faster at 119.2158, which is 34.7982 s ahead of SCA (151.0140s). The same settling time trend is also observed in this, with the HHO settling at 130.0259 s followed by GWO settling at 135.1132 s. The second last rank is for AOA settling at 142.2151 s. In the peak overshoot performance, the HHAOA has 4.7601%, corresponding to a 396.037% reduction of the value, while comparing the SCA at 23.3340%. Amazingly, the GWO has the least peak overshoot value of 2.5131%, followed by WOA, AOA, HHO, and HHAOA, respectively. In this numerical performance analysis, in most of the cases, the Friedman ranking order is observed.
Based on the comparison results, it is evident that the FOPPI controller optimized by HHAOA is more efficient in producing a stable control signal and rejecting disturbances. However, it is worth noting that the HHAOA may have a slower rise time when tracking the initial set-point, which is noticeable. In the same period, the HHAOA has the fastest settling ahead of others with minimal peak overshoot, are shown in
Figure 11, section A. Later, in the external disturbance injection at 100 s, the HHAOA optimized FOPPI had the fastest disturbance recovery and better set-point tracking ability (see,
Figure 11, section B). The control actions of the FOPPI for different algorithms are shown in
Figure 11, sections C and D. Here, the HHAOA has the robust and effectively generated control signal at 4.3 itself, while others are starting at above 6.0. This is also true in the case of the after-disturbance injection scenario. It is clear from the results that finding the optimized controller parameters for the real-time process plant is essential in order to improve its performance.
5. Summary and Conclusions
This section presents the overall summary, which includes the proposed techniques and their advancements given in the first part. Later, the future scope of the current research and the concluding remarks are given.
5.1. Summary
This study analyzed different optimization techniques: AOA, MFO, SCA, GWO, WOA, HHO, and the proposed HHAOA. The tests were conducted on various benchmark functions and optimal router placement and connectivity for WMN benchmark functions. Additionally, the optimized FOPPI controller for a real-time pressure process plant was examined. The HHAOA technique outperformed all other algorithms in every comparison analysis, including the rigorous Friedman ranking test. This statistical evidence unequivocally proves the superiority of the HHAOA method over existing techniques. Moreover, it is crucial to highlight the notable contributions that further emphasize the effectiveness of the proposed technique.
Utilizing a range of multi-hopping methodologies, the proposed HHAOA can effectively detect the global minima with fewer attempts throughout most benchmark functions, resulting in a more efficient and accurate optimization process. These can be evidenced by its consistently higher mean, best, and standard deviation scores in benchmark functions testing.
Through various multi-hopping techniques, the proposed HHAOA is able to identify the global minima in significantly fewer attempts in most of the benchmark functions.
In the WMN, the HHAOA demonstrated a highly competent desired objective searching technique, producing the most optimal path for routing network traffic. The algorithm significantly reduces congestion by minimizing the data transmitted across the network, resulting in improved network performance.
HHAOA carefully selects the appropriate access points to connect to the network, ensuring that clients are always connected to nearby routers. This makes it an essential optimization technique for successfully deploying wireless mesh networks for the best performance.
Notable improvements in producing the robust and smoother control signal of the pressure process were achieved by optimizing the FOPPI controller parameters using HHAOA. These include smaller peak overshoot, dynamic set-point tracking, and practical disturbance rejection ability.
5.2. Conclusions
This research article introduced a novel optimization technique called hybrid HHAOA, which combines two existing algorithms, HHO and AOA, to achieve a better performance. In order to assess the efficacy of the HHAOA algorithm, tests were conducted on a total of 33 optimization benchmark functions. The analysis involved a comparison of performance results, which were based on various measures, including mean, global best, worst, and standard deviation. The convergence performance of the HHAOA algorithm is faster at achieving global minima with fewer iterations than other algorithms. The comparison results were evaluated using Friedman ranking, which showed that the proposed HHAOA algorithm significantly outperforms various algorithms with a 181.342% increased performance ranking in terms of the final mean value. In addition, the best connectivity, network overlapping minimization, and optimal router placement for WMN using the proposed HHAOA were simulated for 500 iterations with 100 search agents. The HHAOA produced the most satisfactory desired performance by creating the best client router connectivity with only one client disconnected in the network. Additionally, the network overlapping was significantly reduced, with a 31.7% reduction in the WMN routers substantially minimizing the operational cost. Experimentation was conducted on a real-time pressure process to further demonstrate the HHAOA algorithm’s effectiveness. The findings showed that the proposed algorithm performed best in the dynamic processes. Furthermore, using HHAOA-optimized FOPPI resulted in a more reliable, smooth, and robust control signal, leading to quicker settling and reduced peak overshoot which, in turn, significantly minimizes the wear and tear on the control valve. As part of future research, newer evolutionary algorithms with different mathematical operators will be investigated to enhance the algorithm’s performance. Additionally, attempts will be made to hybridize the HHAOA algorithm with other metaheuristic optimization algorithms to widen its applicability to more complex, real-time industrial and engineering problems.