An Arithmetic-Trigonometric Optimization Algorithm with Application for Control of Real-Time Pressure Process Plant

This paper proposes a novel hybrid arithmetic–trigonometric optimization algorithm (ATOA) using different trigonometric functions for complex and continuously evolving real-time problems. The proposed algorithm adopts different trigonometric functions, namely sin, cos, and tan, with the conventional sine cosine algorithm (SCA) and arithmetic optimization algorithm (AOA) to improve the convergence rate and optimal search area in the exploration and exploitation phases. The proposed algorithm is simulated with 33 distinct optimization test problems consisting of multiple dimensions to showcase the effectiveness of ATOA. Furthermore, the different variants of the ATOA optimization technique are used to obtain the controller parameters for the real-time pressure process plant to investigate its performance. The obtained results have shown a remarkable performance improvement compared with the existing algorithms.


Introduction
In the era of advanced technological development, the complexity of implementing real-time processes and applications requires effective metaheuristic optimization tuning techniques to meet the global demands with more reliability [1]. The traditional optimization techniques determine the local minima based on the analytical calculations using the relevant processes models, which generally produce a single optimal solution at each run [2,3]. On the other hand, the metaheuristic techniques overcome the conventional methods due to the effective gradient-free mechanisms and the excellent handling of local optima of the respective functions [4,5]. In most metaheuristic algorithms, the optimal solution is found by increasing or decreasing the desired objective function [6,7]. Every optimization algorithm has its challenges in facing the real-time applications that are diverse in nature, multimodal system structure, complexity in implementation, bounded space limitations, longer processing time, flexibility in parameter modifications, and application-specific requirements [8][9][10].
In general, there are different categories of metaheuristic optimization algorithms available. The categorization is based on inspiration: nature-based, swarm and evolutionary intelligence, physics-based, and human activities-based techniques [4,11]. Most metaheuristic optimization techniques identify the optimum solution area by incorporating the two critical search phases: exploration and exploitation [12]. The exploration stage attain the convergence rate and achieve better optimal zone finding in the exploration and exploitation phases.
Some of the significant contributions in this research article can be listed as follows: 1.
The proposed ATOA technique avoids premature convergence and accelerates the search mechanism using the trigonometric functions (i.e., sin, cos, and tan).

2.
The different combinations of the proposed trigonometric function quickly relocate its position from one local minimum to another without getting stuck and reducing the computational complexity. 3.
The proposed optimization is simulated and validated with 33 different benchmark functions to determine its rate of convergence and optimal solution zone identification performance. 4.
The ATOA optimized controller is implemented on the real-time pressure process plant to validate the proposed optimization technique.
The paper's organization is given as follows: A comprehensive design procedure of the proposed ATOA technique is discussed in Section 2. The performance analysis of the proposed ATOA and conventional algorithms on different benchmark functions is given in Section 3. The experimental results of ATOA optimized controller parameters on the real-time pressure process plant are compared to the conventional techniques in Section 4. Summary, concluding remarks, and future directions of the proposed work are provided in Section 5.

The Arithmetic-Trigonometric Optimization Algorithm
This section gives the development of the proposed ATOA based on the inspirations from the SCA and AOA algorithms. Selecting these two specific algorithms is due to their flexibility in modifying the design structure and parameters based on the system requirements. In addition, their simple and easy implementable structure allows it to be used for different engineering problems.

Sine Cosine Algorithm
The SCA algorithm is also a population-based metaheuristic optimization technique. The optimization will be initiated using a group of populations with random solutions. This random set needs a continuous evaluation to reach the optimal solution based on the desired objective function. The probability of finding the best optimal solution in the single run of this metaheuristic algorithm is very low. Still, the chances will be increased if there are a sufficient number of iterations and random set solutions available. The SCA also consists of the exploration and exploitation phases to reach the global minima by locating the optimal search regions. The position update equation during this exploration phase is given as [17] X t+1 Once the area is located by exploration, the fluctuations that occur in the random solutions will be reduced in the exploitation phase. The position update equation during this exploitation phase is given as Equation (2): In Equations (1) and (2), • X t i and X t+1 i are the positions of i th solution at t th and (t + 1) th iterations. • P t i is the point of destination of i th solution at t th iteration. • r 1 , r 2 , r 3 , and r 4 are the random variables.
The parameter r 1 guides the direction of movement. This direction is the next position towards the region between the solution and the destination. The guiding principle of r 1 is defined as where t and T are the number of current and maximum iterations, respectively. The parameter r 2 determines how far the movement of the particles has to go either towards or away from the destination. Similarly, the parameter r 3 specifies the goal to define the distance movement using the random weights to emphasize or de-emphasize them. The equations for generating these r 2 and r 3 values are given as Furthermore, the choice of emphasizing or de-emphasizing effect based on the value of r 3 is as follows: Finally, the parameter r 4 helps switch between the sin and cos functions during the exploration and exploitation. The range of r 4 lies in [0, 1], and this switching will be performed as follows: The pictorial representation of switching conditions of the SCA is given in Figure 1. Here, the random variable r 1 movement towards the destination is based on the current solution and position.

Destination (P)
Solution (X) Next position r 1 1 Next position r 1 >1 Figure 1. Position update mechanism in SCA algorithm.

Arithmetic Optimization Algorithm
The design of the AOA starts with the usage of four arithmetic operators as per the order given as follows: division (÷), multiplication ( × ), subtraction (−), and addition (+), as shown in Figure 2. The conditions for selecting the exploration and exploitation phases of the algorithm are depicted in Figure 2a. The hierarchical order of these four mathematical operators is illustrated in Figure 2b. In contrast, the mechanism of searching for the optimal solution area is shown in Figure 2c [1]. In AOA, three phases are involved in finding the best optimal solution, and they are given as follows: initialization, exploration, and exploitation. Comparing Figure 1 with Figure 2 clearly shows that the AOA follows the footsteps of the SCA in finding the optimal solution.

Initialization
The optimization commences with the randomly generated initial candidate solutions X, which are obtained using Equation (7) with the order of n × d. Here, n is the number of solutions, and d is the number of variables. The best variable obtained from each iteration of the matrix X is considered the best optimal solution, which is given as The math optimizer accelerated (MOA) function determines the next search phase (i.e., exploration or exploitation) of the algorithm with the help of division and multiplication operators. The equation to find the MOA coefficient is given as where "min" is the minimum value of MOA, "max" is the maximum value of MOA, t is the current iteration, and T is the maximum number of iterations.
The next operation in the AOA is categorized into two phases based on the MOA function value. The phase search operation is used to commence the next phase of AOA, which is decided based on the MOA value using the equation as follows:

Exploration
In the exploration operation, division and multiplication operators produce practical search operation values or results, which lead to the identification of near-optimal values for the next phase (exploitation). Here, the exploration is started based on the condition given by MOA in (9). The first operator, ÷, will initiate the searching performance task if the value of r 2 ≤ 0.5. The adjoining operator × will be ignored during this process until the current job is completed. Otherwise, the second operator × will start the first task if r 2 > 0.5 instead of the ÷ operator. The position update of the exploration phase is obtained using the equation, as given as follows.
In the exploration operation, division and multiplication operators produce effective search operation values or results, which lead to the identification of near-optimal values for the next phase (exploitation). Here, the exploration is started based on the condition given by MOA (r 1 ≤ MOA). The first operator, ÷, will initiate the searching performance task at the beginning if r 2 ≤ 0.5. In the meantime, the adjoining operator, × , will be ignored during this process until the current job is completed. Otherwise, the second operator, × , will start the first task (r 2 > 0.5) instead of the ÷ operator. The position update of the exploration phase is obtained using the following equation: where MOP is the math optimizer probability coefficient, which is calculated using In Equation (11), • α is a sensitive parameter that defines the exploitation accuracy. • t and T are the current and maximum number of iterations, respectively.

Exploitation
The operators used in exploitation are subtraction and addition, which have a low dispersion rate, leading to a faster and more accessible approach to reach the target value, which yields a lower number of iterations in the optimization process of the exploration phase. Here, these exploitation operators have enhanced communication which allows better coordination between them. The exploitation search phase is conditioned with the help of the MOA value, as given in Equation (9).
In this stage, the subtraction operator will start performing the searching task if r 3 ≤ 0.5. The addition operator will be ignored during this searching period until the current job is completed. Otherwise, the addition operator will start the task instead if r 3 > 0.5 (see Figure 2c). Furthermore, exploitation search operators will search deeper in the denser regions to find a better optimal solution. It is worth mentioning that, due to the low dispersion rate of the addition and subtraction operators, the probability of getting stuck in the same search region is very minimal. If it happens, they will recover faster than other operators from that local search region. The position update in the exploitation phase is obtained using the following equation:

Arithmetic-Trigonometric Optimization Algorithm
The above two metaheuristic techniques inspired us to hybridize them to create this new ATOA technique. Unlike SCA, the proposed ATOA uses four arithmetic operators instead of the random variables. The ATOA also maintains the same structure and position update mechanism as the existing AOA. On top of that, the proposed method uses trigonometric functions (i.e., sin, cos, and tan) to reposition the solution around the desired solution. This condition ensures the exploration and exploitation search spaces remain around the desired destination area. By combining all these techniques, the design of the proposed ATOA algorithm is explained underneath.
The cyclic and infinite vertical patterns of sin, cos, and tan functions allow the discovered solution to readjust around the nearby solution. This approach guarantees the desired space needs during the exploitation phase in between the two adjacent solutions. During the exploration phase, by adjusting the range of the trigonometric functions, the solutions yield the capability to search outside the space between their corresponding optimum solution zone. A position update is achieved by changing ranges of the sin, cos, and tan functions. Here, the solution repositions inside or outside their search space between itself and the adjacent key. The position update, either inside or outside, is determined by the random number r 2 used in the ATOA. Hence, a balance between switching exploration and exploitation is essential, and it is regulated based on MOA. In addition, the range of the trigonometric function is controlled by using MOP. The diverse patterns of the trigonometric functions in the range of [0, 3π] are illustrated in Figure 3.  The random numbers r 1 and r 2 determine the algorithm phase and operators switching based on its value. During the standalone implementation of sin, cos, and tan functions in the proposed ATOA technique, only one trigonometric function is used in all the phases (i.e., exploration and exploitation) in between the bounded values. Thus, the position update of the standalone sin function in the exploration and exploitation phases of ATOA is obtained using the following Equations (13) and (14), respectively. Similarly, for the following trigonometric function, i.e., cos, the existing sin is replaced with cos during both phases, and the same procedure is repeated for tan.
The different combinations of the trigonometric functions used during the exploration and exploitation phase are determined based on the MOA coefficient. The motivation to select the various trigonometric functions combinations is to address the lead-lag characteristics sin and cos. In addition, it is worth mentioning that the non-synchronous nature of the tan function will produce asymmetric values when it combines with the sin or cos. Thus, this article has experimented with the standalone usage of the tan in the exploration and exploitation phases. The range variations in the trigonometric functions are carried out using the MOP coefficient.
The term ATOAsc represents that the sin function is used in the exploration phase and the cos function during the exploitation phase. Thus, the position update of the ATOAsc algorithm during exploration and exploitation phases is given as follows: Similarly, in ATOAcs, the cos function will be automatically assigned as the initial function during the exploration phase and the sin operation in the exploitation phase. Thus, the position update of this ATOAcs algorithm is given as follows: In Equations (13) to (18), • x i,j (t + 1) represents the i th solution in the (t + 1) th iteration at the j th position. • best(x j ) is best solution obtained at j th position. • UB j and LB j are the upper and lower boundaries at j th position. • is the constant integer. • µ is the search control parameter. Table 1 gives the list of developed algorithms and the functions used during the exploration and exploitation phases. The implementation of all the developed algorithms is the same, and only the functions in the corresponding position update equation will change. For an example case, the pseudocode for the implementation of the proposed ATOAcs algorithm is described in Algorithm 1. In this case, the exploration phase uses the cos function, and the exploitation phase uses the sin function. On the other hand, during the standalone implementation of ATOAs, the sin function will be used in the exploration and exploitation phases. Calculate the best(x j ), MOA(t) and MOP(t). 4: for i = 1 : size(X, 1) do 5: for j = 1 : size(X, 2) do 6: Generate if r 1 ≤ MOA then 8: Generate if r 2 ≤ 0.5 then 10: Compute the x i,j (t + 1) using the first case in (17).

Performance Analysis on Benchmark Functions
This section discusses the simulation and comparative analysis of the proposed ATOA algorithm with conventional AOA. Simulations over 33 different benchmark functions are carried out to validate the proposed ATOA efficiency. The corresponding benchmark functions and their equations and the desired range limits are given in Table 2. Here, the abbreviated terms cat. and func. represents the category and benchmark function, respectively. In all the cases, the different variants of the proposed ATOA algorithm and the existing AOA algorithm are compared to prove the effectiveness of the proposed algorithm. Furthermore, all the simulations are performed in MATLAB/Simulink software (2021a) using the 3.10 GHz Intel(R) Xeon PC with 16.00 GB of RAM.

Selection of Benchmark Functions
For effective comparison analysis, the benchmark functions are selected based on the independence of one another in multiple dimensions, different local minima, and diverse boundary values. The number of iterations and the population size of the search solutions are also kept constant, respectively, at 300 and 30. During the simulation analysis, the optimization parameters used in all the compared algorithms are min = 0.2, max = 1.0, α = 5, µ = 0.499, and = 2.2204 × 10 −16 . A single global minima value function is known as the unimodal function. Therefore, they are employed to test the exploitation precision capabilities of the optimization algorithm. Such unimodal benchmark functions used for the simulation are F1-F6, as given in Table 2.

Cat. Func. Description Range
F32 The proposed method is also validated using the multimodal and hybrid composition functions in single and multidimension. The multimodal functions have multiple local minima points used to test the algorithms' exploration ability, i.e., to check the ability to switch from the local minima to the global minima without getting caught in the same position. F7-F10 are the multimodal benchmark functions, and their equations are given in Table 2. In addition, the functions F11-F18 shown in Table 2 have multiple local minima but with fixed dimensions. These functions will help to test the stability of the algorithm. It is worth mentioning that the first 13 benchmark functions use higher dimension values of (30/100/500/1000) to evaluate the proposed ATOA. Other considered benchmark functions F19-F33 are hybrid, as given in Table 2. In the unimodal functions, the probability of getting caught in the same search space is minimal, as the optimal search solutions have a much clearer position update to move from their current position to the desired position. Additionally, in the case of hybrid functions from F19-F33, the search space identification and the next position movement are more challenging due to multiple local minima values.
The surface plots of all the considered functions are shown in Figure 4. These plots will help to visualize required search space requirements for the single and multi global minima functions. In these functions, many categories were present based on their surface search area and shape. The classifications are bowl-shaped, plate-shaped, valley-shaped, single, and multiple-local minima functions. These functions have multiple local minima values with an ample search space, wide range, numerous layers, and multidimensional characters.

Numerical Analysis on Benchmark Functions
In this section, the different variants of the proposed ATOA algorithm are numerically compared with the conventional AOA technique. The comparison results of the various benchmark functions (F1-F33) are given in Table 3. Here, the abbreviated terms refer to the different combinations of the trigonometric functions as mentioned in the methodology section. They are ATOAs, ATOAc, ATOAt, ATOAsc, and ATOAcs, which refer to ATOAsine, ATOAcos, ATOAtan, ATOAsinecos, and ATOAcossine, respectively. The performance of all the algorithms is evaluated in terms of mean value, the best solution, worst value, and the standard deviation (std. dev.). The obtained results have shown a tough competition between the proposed ATOAcs and the AOA in the unimodal benchmark functions (F1-F6) in the mean, best, and worst value. During this performance, the proposed ATOAcs obtained significantly closer convergence to the global minima. While noticing the standard deviation results of the unimodal function, both the ATOAcs and AOA had the balanced solution-attaining ability. Among the different proposed variants, the ATOAcs had better performance than AOA due to the leading characteristics of the cos function. The other proposed combinations had a performance setback while compared with the AOA in the unimodal functions. In the multimodal and the fixed-dimension benchmark functions (F7-F18), all the proposed ATOA variants outperformed the conventional AOA in all the parameters, and among these, ATOAcs took the first position. This performance shows the efficiency of the proposed ATOA in the exploration and exploitation capabilities, which resulted in converging very closely or strictly at global minima value. It is worth mentioning that in the performance in the hybrid benchmark functions, the proposed ATOA variants produced the same results as the previous stage with a higher and faster convergence rate.  Additionally, in order to compare the performance of the different techniques quantitatively, the Friedman ranking test is used [45]. This method ranks the functions based on the nonparametric analysis of variance produced in one complete run of the algorithm. The Friedman ranking test for the different compared algorithms for all the benchmark functions is given in Table 4. The ranking is based on the best mean value obtained by each algorithm in the respective function. For instance, rank 1 denotes the closest value to the global minima, and rank 6 denotes the farthest convergence value.
In the initial benchmark functions, the existing AOA had dominance, but in the latter, the proposed ATOA variants showed outstanding search space identification, position update, and convergence performance. Regarding the final mean values, the ATOAcs had the least mean value of 2.454 and secured the first rank. The conventional AOA had a mean value of 3.696 that was 1.242 times higher than ATOAcs. ATOAs took second place by obtaining the value of 2.636. The respective values of ATOAsc and ATOAc were 3.242 and 2.878. The proposed ATOAt had a value of 3.484, which is almost close to the conventional AOA. Thus, the proposed technique validates their abilities to attain the best mean value for the numerical analysis, and the same results are reflected in the convergence performance.

Convergence Analysis
The convergence speed plot of the different algorithms in the benchmark functions are shown in Figure 5. In the initial benchmark functions (F1-F6), the ATOA had a slow search rate due to the distribution of the obtained solutions to the different search regions instead of aggregating them in the current search area. However, the proposed ATOA overcomes those issues in the benchmark functions F7-F10 and shows faster convergence at the global minima with fewer iterations. In addition, in functions F11-F18, the proposed ATOA variants can converge all over the iterations continuously. This indicates that the proposed combinations of the trigonometric functions in this research enhance the exploration and exploitation search space of the original AOA.
The proposed trigonometric combinations assisted the existing AOA to tackle all issues by making the ATOA search mechanism faster, using every iteration's current best optimal value. These made the proposed ATOA algorithms find those global minima values in lower iterations than the existing AOA algorithm in most benchmark functions. It is worth mentioning that, in the hybrid benchmark functions F19-F29, ATOA found the optimal positions quickly even though these functions contain multiple local minima. The benchmark functions results of F20, F21, and F30-F33 show that the proposed ATOA experienced multiple local minima positions. However, the ATOA recovered quickly from those local minima positions and found the desired global minima with few iterations using the proposed trigonometric functions. The proposed ATOA accomplished better accuracy in most of the benchmark functions as opposed to the compared algorithms. Finally, sometimes the AOA is more effective than the ATOA in some of the initial functions.
The overall ranking and the performance of the proposed ATOA combinations produced a better optimal solution.

Performance Analysis on Control of Real-time Pressure Process Plant
In the first part of this section, detailed information about the real-time pressure process plant, including its schematics and P&I diagram, is given, along with its mathematical modeling design procedure. Then, the selection of the controller and its implementation along with their optimized controller parameters are given. Finally, the performance analysis on the process plant is discussed. Figure 6a shows the experimental setup of the pressure process plant with its peripherals. The tank VL 202 is the main process tank supplied with a maximum pressure of 10 bar from the centralized air compressor. The hand valve HV 202 acts as the primary input pressure control for the process plant. Once this hand valve is opened, the input pressure to the buffer tank is controlled using the process control valve PCV 202. The pressure inside the VL 202 is measured using the digital pressure transmitter PT 202 attached to the process tank. The PT 202 measures the pressure values and converts them to digital voltage signals in the range of 0 to 5 V. The pressure-indicating controller PIC 202 receives these voltage signals and feeds them to the host PC via I/O interface boards to generate the control signal for the PCV 202. Additionally, an analog pressure gauge is available to indicate the real-time pressure changes inside the tank.

Industrial-Scale Setup of Real-Time Pressure Process Plant
The process tank contains a pressure-release hand valve attached at the bottom to release the pressure during an emergency. This hand valve provides a safer way to discharge the pressurized air from the VL 202 in case of a pneumatic process control valve (PCV 202) failure. The pressure inside the buffer tank is regulated by releasing the excess air via a gas outlet available on top of the process tank. This outlet gas pipe is connected with another process control valve, PCV 203. The PCV 303 is maintained at 50% opening during the experimentation to avoid excessive pressure build-up inside the VL 202. The PCV 202 and PCV 203 receive their control signals from the host PC via the I/O interface boards.
The piping and instrumentation diagram of the pressure process plant is illustrated in Figure 6b. The pressure process plant is controlled using "Remote Desktop Connection" mode from the central control room for safety. The host PC in the process plant uses the MATLAB/Simulink software to generate control signals based on the input set-point and the actual process variable. Data transmission and communication between the host PC and the process field devices (i.e., control valve actuator and pressure transmitter) happens via the peripheral component interconnect (PCI) cards. The PCI card contains three separate modules, namely 1713U, 1720U, and 1751. PCI-1720U performs the function of sending the control signals to the host PC using its 12-bit, 4-channel analog output port. In addition, all the interfacing cards are provided with isolation protection of 2500 V DC between the outputs of the PCI bus. The analog input from the process plant is received via a 12-bit, 32-channel analog input card PCI-1713U with a sampling rate of 100 k samples. The PCI-1751 card carries out the digital signals' data transmission from the PT 202 to the host PC and vice versa. It contains 48 bits of parallel digital input/output, enabling the process plant to be controlled remotely.
The mathematical modeling of the process plant is obtained using the open-loop step response characteristics. This response contains all the behavioral data and process dynamics such as process gain K, process dead-time L p , and process time constant T. Thus an open-loop step response of the pressure process is carried out, and the plant dynamics are obtained as follows: K = 0.866, L p = 1, and T = 1.365. Based on the first-order plus dead-time system characteristic equation, the process model of the pressure process plant is constructed using the obtained plant dynamics. Thus, the transfer function of the pressure process plant is given as follows:

ATOA-Based Fractional-Order Predictive PI Control of Pressure Process Plant
PI controllers are the most well-known and widely adopted controllers in most industries. However, they perform poorly in an environment having delay, noise, networkinduced delay, and high-frequency noise [47,48]. On the other hand, advanced control strategies, such as model predictive controller and generalized predictive controller, are computationally and structurally complex for implementation. Thus, the controllers need a simple design, effective dead-time compensation, and high disturbance-rejection characteristics to overcome the above issues. Such a type of controller is proposed by Arun et al. [46], which has the dead-time compensating capability of the Smith predictor and robustness nature of the fractional-order controllers. Therefore, in this article, such an efficient fractionalorder predictive PI (FOPPI) controller designed by them is utilized. The selection of pressure process application for this research is due to its nonlinearity and sensitivity [49]. The FOPPI controller generates an effective control signal which is free from load variations and plant uncertainties. However, their controller parameters are obtained analytically, and those values are insufficient to mitigate the real-time process disturbances [50]. Thus, the proposed ATOA and conventional AOA techniques are utilized in this research to obtain effective controller parameters. The adequately tuned parameters are used for experimentation on the real-time pressure process plant.
The control signal u(s) of the adopted FOPPI controller is given as where K p = 1 K , and T i = T. In addition, u(s) and e(s) are the control and error signals, K p is the proportional gain, T i is the integral time, λ is the fractional-order integrator, and L p is the process dead-time.
For the optimization algorithms, the objective is to obtain the effective controller parameters (i.e., K p , K i , and λ), the K i is the integral gain of the controller, and it is obtained In the analytical design technique, the FOPPI controller parameters are obtained directly from its design. The closed-loop system containing the FOPPI controller and the obtained pressure process model G p (s) in the presence of the ATOA is shown in Figure 7. The proposed ATOA is used to tune the FOPPI controller parameters using the integral time absolute error (ITAE) value. The ITAE value is assigned as the objective function for the proposed ATOA and AOA. The calculation of the ITAE value is obtained as

Performance Analysis
The performance analysis of the FOPPI under various types of ATOA and AOA optimizations are compared in terms of process rise time (t r ), settling time before (t s1 ) and (t s2 ), overshoot (%OS), and disturbance rejection. A 30% disturbance is injected at 150 s in the process feedback loop to investigate the set-point tracking capability of the FOPPI controller. In the analytical design, the controller parameters are obtained based on the FOPPI design structure itself. The different optimized FOPPI controllers using the various optimization algorithms and their numerical performance analysis are given in Table 5.
While observing the numerical analysis results given in Table 5, the proposed ATOA algorithm variants performed better than the AOA and the analytical design. Among all the compared algorithms, the existing analytical method had the faster rise time of 0.7665, followed by AOA, ATOASC, ATOAt, ATOAs, ATOAc, and ATOAcs, respectively. In this, the first-ranked ATOAcs had the slowest rise time, of 2.4432 s, compared to other methods. The proposed ATOAcs settled faster than all at 61.3744 s before the disturbance injection, even with the slowest rise time. Second place is secured by ATOAs with the settling time at 65.8043 s, followed by ATOAsc at 69.5402 s, ATOAc at 72.1039 s, AOA at 75.8397 s, ATOAt at 77.8114s, and analytical design last with 83.0137s. In this performance, the Friedman ranking order of the algorithms is also reflected, except for ATOAt. After the disturbance injection, the AOA and analytical design managed to track the set-point, but both of them settled slower than the proposed ATOA at 276.7360 s and 280.4273 s, respectively. Here, the ATOAt managed to settle faster than the AOA and analytical design with a difference of 2.7189 s and 6.4102 s, respectively. In this case, the ATOAcs recovered effectively from the disturbance and settled 22.7199 s faster than the analytical design. The controller parameters obtained by the proposed algorithms significantly reduced the peak overshoot values, which indicates the excellent objective function-finding capability of the ATOA techniques. In the peak overshoot performance, ATOAc produced the most negligible value of 2.8020%, which is 19.427% less than the analytical design. ATOAt had the second-least value of 3.1546%, followed by ATOAs, ATOAsc, ATOACS, AOA, and analytical technique.
The disturbance rejection, peak overshoot, and set-point tracking performance of the compared algorithms are shown in Figure 8. The obtained results indicate the excellent tracking and faster recovering ability of the proposed ATOA optimized controller parameters. Zoomed regions of Figure 8 during the before and after disturbance are given in Figure 9. Regarding the results in region A, the AOA and analytical method produced more oscillations, which will significantly deteriorate the lifetime and performance of the control valve actuator. The same trend is repeated during the performance after the disturbance injection (see region B). It is worth mentioning that the control signals generated from the ATOA optimized values are less aggressive, and smooth compared with AOA and analytical design (see regions C and D). Once again, a similar pattern of Friedman ranking results is obtained in the ITAE performance, and the ATOAcs had the minimal value of 1.4224. Then, ATOAs had the second-least value of 1.7391, followed by ATOAsc with 2.0154. Here, the conventional AOA almost had a comparable value to the proposed ATOAt. Lastly, the analytical design had a value of 3.7159, 2.2935 higher than the proposed ATOAcs.

Summary and Conclusions
In the initial part of this section, the contributions of the proposed research are summarized. Then, the concluding remarks and the future directions of the research are provided.

Summary
The different combinations of the proposed ATOA, AOA, and the conventional analytical design algorithms helped us to investigate their performance on various benchmark functions and the real-time pressure process. The usage of the basic arithmetic and trigonometric functions of the ATOA provides a straightforward design and easy implementation. However, the advanced functions such as modulus, cubic, and sigmoid may catch up with the proposed technique's performance. At the same time, even after the addition of the new functions, the proposed ATOA maintains the simple structure of the conventional SCA and AOA. All the above claims are validated based on the obtained results, which shows the faster convergence ability of ATOA in lesser iterations. Furthermore, the Friedman ranking of the compared algorithms given in Table 4 validates the effectiveness of using the trigonometric functions in the existing algorithms. In addition, the ATOA technique dramatically improves the search mechanism, which leads to the identification of the efficient, optimized FOPPI controller parameters shown in Table 5. The following points highlight some of the essential contributions of this article:

1.
The proposed ATOA outperformed all the compared algorithms in most of the benchmark functions in terms of mean, best, and standard deviation.

2.
The proposed ATOA variants produced the best global minima in fewer number of iterations. Among them, ATOAcs achieved phenomenal results in all the comparative research analysis. 3.
The proposed ATOAcs and ATOAs had an efficient global optima search mechanism and yielded better performance, and it is proven by the Friedman ranking test given in Table 4.

4.
The ATOA-optimized FOPPI controller parameters performed effectively by reducing the peak overshoot, actively tracking the set-point, and efficiently minimizing the disturbance impacts on the process.

5.
The control signals of the optimized FOPPI controller greatly smooth the control actions by filtering out the undesired stochastic disturbances.

Conclusions
This paper develops an enhanced arithmetic-trigonometric optimization algorithm by incorporating different trigonometric functions, namely sin, cos, and tan. The performance of the proposed metaheuristic ATOA algorithm was validated on thirty-three different optimization benchmark functions. Furthermore, the ATOAs were compared for the mean, global best, worst, and standard deviation performance to validate their effectiveness. The convergence performance results showed the proposed ATOA to achieve faster global minima in fewer iterations. Friedman ranking of the comparison of the different algorithms showed the improvements obtained by the proposed ATOA algorithm. Finally, experimentation on the real-time pressure process was carried out to prove the dynamic abilities of the proposed algorithm in obtaining the best optimal solution. In future, newer evolutionary and additional arithmetic operators will be considered. Furthermore, an attempt to hybridize the proposed ATOA with the other population-based metaheuristic optimization algorithms will be made.