*Article* **Load Frequency Control and Automatic Voltage Regulation in a Multi-Area Interconnected Power System Using Nature-Inspired Computation-Based Control Methodology**

**Tayyab Ali <sup>1</sup> , Suheel Abdullah Malik <sup>1</sup> , Ibrahim A. Hameed 2,\* , Amil Daraz <sup>3</sup> , Hana Mujlid <sup>4</sup> and Ahmad Taher Azar 5,6,7,\***


**Citation:** Ali, T.; Malik, S.A.; Hameed, I.A.; Daraz, A.; Mujlid, H.; Azar, A.T. Load Frequency Control and Automatic Voltage Regulation in a Multi-Area Interconnected Power System Using Nature-Inspired Computation-Based Control Methodology. *Sustainability* **2022**, *14*, 12162. https://doi.org/10.3390/ su141912162

Academic Editors: Herodotos Herodotou, Sheraz Aslam and Nouman Ashraf

Received: 19 August 2022 Accepted: 20 September 2022 Published: 26 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

**Abstract:** The stability control of nominal frequency and terminal voltage in an interconnected power system (IPS) is always a challenging task for researchers. The load variation or any disturbance affects the active and reactive power demands, which badly influence the normal working of IPS. In order to maintain frequency and terminal voltage at rated values, controllers are installed at generating stations to keep these parameters within the prescribed limits by varying the active and reactive power demands. This is accomplished by load frequency control (LFC) and automatic voltage regulator (AVR) loops, which are coupled to each other. Due to the complexity of the combined AVR-LFC model, the simultaneous control of frequency and terminal voltage in an IPS requires an intelligent control strategy. The performance of IPS solely depends upon the working of the controllers. This work presents the exploration of control methodology based on a proportional integral–proportional derivative (PI-PD) controller with combined LFC-AVR in a multi-area IPS. The PI-PD controller was tuned with recently developed nature-inspired computation algorithms including the Archimedes optimization algorithm (AOA), learner performance-based behavior optimization (LPBO), and modified particle swarm optimization (MPSO). In the earlier part of this work, the proposed methodology was applied to a two-area IPS, and the output responses of LPBO-PI-PD, AOA-PI-PD, and MPSO-PI-PD control schemes were compared with an existing nonlinear threshold-accepting algorithm-based PID (NLTA-PID) controller. After achieving satisfactory results in the two-area IPS, the proposed scheme was examined in a three-area IPS with combined AVR and LFC. Finally, the reliability and efficacy of the proposed methodology was investigated on a three-area IPS with LFC-AVR with variations in the system parameters over a range of Â ± 50%. The simulation results and a comprehensive comparison between the controllers clearly demonstrates that the proposed control schemes including LPBO-PI-PD, AOA-PI-PD, and MPSO-PI-PD are very reliable, and they can effectively stabilize the frequency and terminal voltage in a multi-area IPS with combined LFC and AVR.

**Keywords:** PI-PD controller; load frequency control; automatic voltage regulator; nature-inspired optimization; multi-area interconnected power system

#### **1. Introduction**

Research efforts and specializations in power systems are increasing day by day to acquire reliable power with nominal voltage and frequency. In a power system, the main

goal is to provide nominal voltage and frequency to all consumers without any interruption. The simultaneous control of load frequency and terminal voltage in an interconnected electrical power system is the fundamental area of research for all practicing engineers. The mutilation of frequency or voltage can spoil the performance and life expectancy of equipment associated with IPS [1]. The active and reactive powers can change with load demands in IPS. The active power can be adjusted by a speed governor in an LFC loop, whereas reactive power can be controlled by an exciter in an AVR loop. In order to fulfill the active power demand, a turbine input is continuously regulated in LFC, or else the changing frequency will vary the machine's speed. In AVR, terminal voltage remains within the prescribed limit if the excitation of generators is regulated properly to match the reactive power demand. A lot of literature is available on individual AVR or LFC systems; however, relatively less research work has been carried out on combined LFC-AVR due to its complex design. The PID controller was extensively used in multi-area IPS due to its simple design and easier installation. For instance, the artificial electric field algorithmbased hybridized approach to tune the fuzzy PID controller was suggested for combined LFC and AVR with the incorporation of different energy storage devices [1]. A particle swarm-optimized Ziegler–Nicholas (PSO-ZN)-based PID controller was examined for AVR-LFC control in PV integration and a conventional power system [2]. PI and PID with filter (PIDF) controllers based on the sine cosine algorithm were also inquired for a two-area, two-source IPS. The redox flow batteries were assimilated for further improvements in the system dynamics [3]. The doctor and patient optimization (DPO)-based accelerating PID controller (PIDA) was proposed for the LFC-AVR problem in a multi-area IPS with renewable energy sources [4]. The PID controller was employed for collective AVR-LFC in a two-area IPS. A nonlinear threshold-accepting algorithm was explored to find the optimum parameters of the PID controller [5]. PI and I controllers for AVR and an LFC loop were also investigated for a single-area IPS [6]. In [7], due to the inclusion of deregulated environments in IPS, a fuzzy logic controller (FLC) was recommended for a two-area LFC-AVR problem. A fractional order controller (PIDµF) based on the lightning search algorithm (LSA) was also proposed for LFC-AVR with wind and a reheat thermal plant as the generating companies (GENCOs) of area-1,and with diesel and a nonlinear reheat thermal plant as the GENCOs of area-2 under deregulated environments [8].The PID controller was optimized with the hybridization of the artificial electric field algorithm and differential evolution for a two-area IPS with a joint LFC-AVR [9,10]. In [11], PID with the firefly algorithm was employed for a two-area IPS with AVR-LFC. The moth flame optimization (MFO)-based fractional order PID controller was proposed for both LFC and AVR loops [12]. For a single-area synchronous generator, the combined LFC-AVR was explored using a hardware environment [13]. In [14,15], the authors inspected the firefly algorithm, particle swarm optimization, and the genetic algorithm-based PID controller for AVR-LFC loops. The novel state-observer (SO)-based integral double-derivative controller based on magneto-tactic-bacteria optimization (MBO) was presented for voltage–frequency control in a hybrid IPS [16]. The model predictive controller (MPC) was also used to improve AVR-LFC responses [17]. In [18], the heuristic computation-based two degrees of freedom state-feedback PI controller was exploited for the AVR loop in synchronous generators. A combination of the bacterial foraging optimization algorithm and particle swarm optimization was utilized to tune the PI controller for the AVR system with a static synchronous compensator [19]. In [20], a sliding mode controller with the addition of a gene ralized extended state observer was successfully explored to optimize the LFC loop in a multi-area IPS. The PID controller tuned with the many optimizing liaisons (MOL) algorithm was applied to a two-area IPS with non-reheat thermal sources in the presence of GDB [21]. Moreover, a comprehensive research work was presented for individual LFC loops as presented in [22–37]. A brief literature summary of AVR-LFC is provided in Table 1. It can be seen that much less attention has been given to the combined LFC-AVR problem in multi-area IPS due to its complex structure. The literature survey also depicts that modified forms of the PID controller were explored very rarely for combined AVR-LFC. Different

modified forms such as PI, PIDF, PIDµF, and FO-IDF have been explored due to their excellent time response characteristics with fast convergence, but the PI-PD controller has not been employed for combined LFC-AVR multi-area IPS. Due to its modified structure having a control branch in the feedback path, complex systems can be well optimized with PI-PD as compared to classical control schemes such PI and PID, etc. To obtain optimal controller parameters, an intelligent tuning algorithm is needed, which can optimize the controller with minimum error/fitness. In the past, nature-inspired optimization algorithms have received a lot of attention from researchers because of their strengths and abilities to tackle a variety of complex optimization issues in engineering. These strategies have also been used successfully to obtain optimal controller parameters. The classical nature-inspired computing techniques have shown very satisfactory performances for both individual and combined LFC-AVR. Moreover, researchers have also presented some novel nature-computing algorithms such as dandelion optimizer [38], modified particle swarm optimization (MPSO), bald eagle search (BES) [39], the transient search algorithm (TSO) [40], learner performance-based behavior optimization (LPBO) [41], the Archimedes optimization algorithm (AOA) [42], etc. These recently introduced techniques such as MPSO, LPBO, and AOA have not been considered for the optimal tuning of the PI-PD control scheme. It will be worth choosing these nature-inspired techniques for the optimization of multi-area IPS with combined LFC-AVR. Keeping in mind the existing research gap, the nature-inspired computation-based PI-PD control scheme is proposed in this research for multi-area IPS with combined AVR-LFC. The main contributions of this work are:



**Table 1.** Literature on ALR-LFC.


**Table 1.** *Cont.*

Table 2 demonstrates the nomenclature used in this study. This research paper is organized in following way: The power system model is described in Section 2. The proposed control methodology is presented in Section 3. Section 4 contains a description and flow charts of nature-inspired computation algorithms including LPBO, AOA, and MPSO. The implementation and results of the proposed techniques are summarized in Section 5. Lastly, conclusions and future guidelines are given in Section 6.



#### **2. Power System Model**

The multi-area IPS model under study is shown in Figure 1. The terminal voltage was maintained at nominal value by stabilizing the generator fields, while the load frequency was regulated by controlling real power. Figure 1a represents the AVR-LFC model of a power system for a single area, where *i* and *j* represent area-1 and area-2, respectively [5].

**2. Power System Model**

1( ) ( ) *K s <sup>i</sup>* ), amplifier ( ( )

( ) 1 *a i a i*

*K sT*

[5].

( ( )

nor ( ( )

*K sT*

( ) 1 *s i s i*

> *K sT*

( ) 1 *G i G i*

area-*i* and area-*j* is represented by T*ij*.

The multi-area IPS model under study is shown in Figure 1. The terminal voltage was maintained at nominal value by stabilizing the generator fields, while the load frequency was regulated by controlling real power. Figure 1a represents the AVR-LFC model of a power system for a single area, where *i* and *j* represent area-1 and area-2, respectively

*Vt i* ( ) , *Ve i* ( ) , *Vref i* ( ) ,and *Vs i* ( ) refer to the terminal output, error, reference, and sensor voltage in area-1, respectively. The AVR system of area-*i* consists of a controller (

> *K sT*

( ) 1 *g i g i*

), exciter ( ( )

( ) 1 *p i p i*

*K sT*

( ) 1 *e i e i*

> *K sT*

( ) 1 *t i t i*

). ( )*<sup>i</sup> f* denotes fre-

), and sensor

), gover-

*K sT*

), generator ( ( )

), speed regulation (*Ri*), and generator/load ( ( )

).Area-1′s LFC system has <sup>a</sup> controller 2 ( ) ( ) *K s <sup>i</sup>* , turbine ( ( )

quency deviation (Hz), *X G i* ( ) shows the valve position of the governor (p.u.MW), *PG i*( ) represents the deviation in the output of the generator (p.u.MW), *PD i*( ) (p.u.MW) denotes the deviation in load, speed regulation is represented by *R*( )*<sup>i</sup>* (Hz p.u.MW−1), and *Ptie i*( ) is the tie–line power. The purpose of tie–line is to interconnect multiple areas in IPS. Figure 1b shows the tie–line connections. The synchronization coefficient between

**3. Proposed Control Methodology** The proportional integral derivative (PID) controller is commonly utilized in industrial applications owing to its easier implementation and simpler structure. The PID controller provides a satisfactory performance in most of the systems; however, the modified forms of the PID control structure have shown improved performance in many control systems, such as the AVR-LFC interconnected power system. The proportional integral– proportional derivative controller (PI-PD) is a modified version of PID, which is designed in such a way to eliminate system errors with optimum transient and steady state response [43]. The PI part of PI-PD exists in a feed forward path and directly responds to the error signal coming from the summing junction. The PD part is located in the feedback path, and it is unaffected by sudden changes in the set point specification. The closed-loop response can be improved significantly with the addition of a controller part in the feed-*Vt*(*i*) , *Ve*(*i*) , *Vre f*(*i*) , and *Vs*(*i*) refer to the terminal output, error, reference, and sensor voltage in area-1, respectively. The AVR system of area-*i* consists of a controller (*K*1(*i*) (*s*)), amplifier ( *<sup>K</sup>a*(*i*) *sTa*(*i*)+1 ), generator ( *<sup>K</sup>g*(*i*) *sTg*(*i*)+1 ), exciter ( *<sup>K</sup>e*(*i*) *sTe*(*i*)+1 ), and sensor ( *<sup>K</sup>s*(*i*) *sTs*(*i*)+1 ). Area-10 s LFC system has a controller *K*2(*i*) (*s*), turbine ( *<sup>K</sup>t*(*i*) *sTt*(*i*)+1 ), governor ( *<sup>K</sup>G*(*i*) *sTG*(*i*)+1 ), speed regulation (*R<sup>i</sup>* ), and generator/load ( *<sup>K</sup>p*(*i*) *sTp*(*i*)+1 ). ∆*f* (*i*) denotes frequency deviation (Hz), ∆*XG*(*i*) shows the valve position of the governor (p.u.MW),∆*PG*(*i*) represents the deviation in the output of the generator (p.u.MW), ∆*PD*(*i*) (p.u.MW) denotes the deviation in load, speed regulation is represented by *R*(*i*) (Hz p.u.MW−1), and ∆*Ptie*(*i*) is the tie–line power. The purpose of tie–line is to interconnect multiple areas in IPS. Figure 1b shows the tie–line connections. The synchronization coefficient between area-*i* and area-*j* is represented by T*ij*.

#### back path. The PI-PD controller has been successfully employed in the recent past in dif-**3. Proposed Control Methodology**

and reference signal.

ferent applications [44–50]. The proposed control methodology with the combined LFC-AVR system is given in Figure 2. The transfer function of PI-PD controllers is represented as: 1 2 ( ) ( ) ( ) ( ) ( ) *<sup>i</sup> p p d <sup>K</sup> U s K E s K K s Y s s* (1) The proportional integral derivative (PID) controller is commonly utilized in industrial applications owing to its easier implementation and simpler structure. The PID controller provides a satisfactory performance in most of the systems; however, the modified forms of the PID control structure have shown improved performance in many control systems, such as the AVR-LFC interconnected power system. The proportional integral–proportional derivative controller (PI-PD) is a modified version of PID, which is designed in such a way

where *U*(*s*), *Y*(*s*), *R*(*s*), and *E*(*s*) denote the control, output, reference, and error signals,

*E s Y s R s* ( ) ( ) ( ) (2)

to eliminate system errors with optimum transient and steady state response [43]. The PI part of PI-PD exists in a feed forward path and directly responds to the error signal coming from the summing junction. The PD part is located in the feedback path, and it is unaffected by sudden changes in the set point specification. The closed-loop response can be improved significantly with the addition of a controller part in the feedback path. The PI-PD controller has been successfully employed in the recent past in different applications [44–50]. The proposed control methodology with the combined LFC-AVR system is given in Figure 2. The transfer function of PI-PD controllers is represented as:

$$\mathcal{U}L(\mathbf{s}) = (\mathcal{K}\_{p1} + \frac{\mathcal{K}\_i}{\mathcal{s}})E(\mathbf{s}) - (\mathcal{K}\_{p2} + \mathcal{K}\_d \mathbf{s})\mathcal{Y}(\mathbf{s})\tag{1}$$

$$E(\mathbf{s}) = \mathbf{Y}(\mathbf{s}) - \mathbf{R}(\mathbf{s}) \tag{2}$$

where *U*(*s*), *Y*(*s*), *R*(*s*), and *E*(*s*) denote the control, output, reference, and error signals, respectively. The cost function (J) is minimized to obtain the best possible parameters of the controllers. J depends upon *E*(*s*), which is basically the difference between the output and reference signal. *Sustainability* **2022**, *14*, x FOR PEER REVIEW 7 of 32

**Figure 2.** Proposed control methodology with combined LFC−AVR system. **Figure 2.** Proposed control methodology with combined LFC−AVR system.

In order to minimize the error signal, different types of performance indices can be used such as the integral of the squared value of the error signal (ISE), the integral of the time multiplied with the absolute value of the error signal (ITAE),the integral of the time multiplied with the squared value of the error signal (ITSE), and the integral of the absolute value of error (IAE) represented by the following equations: In order to minimize the error signal, different types of performance indices can be used such as the integral of the squared value of the error signal (ISE), the integral of the time multiplied with the absolute value of the error signal (ITAE), the integral of the time multiplied with the squared value of the error signal (ITSE), and the integral of the absolute value of error (IAE) represented by the following equations:

JITSE,two-area = 2 2 2 2 2

*t t t tie tie tie f f f V V V P P P dt* (7)

*t t t tie tie tie t f f f V V V P P P dt* (9)

*t t t ptie ptie ptie t f f f V V V P P P dt* (10)

$$\mathcal{J}\_{\text{ISE,two-area}} = \int\_0^T [\Delta f\_1^2 + \Delta f\_2^2 + \Delta V\_{l1}^2 + \Delta V\_{l2}^2 + \Delta P\_{\text{tie12}}^2]dt\tag{3}$$

1 2 1 2 12 <sup>0</sup> [ ]

*t t tie t f f V V P dt* (5)

$$\mathcal{J}\_{\text{ITAE}, \text{two}-\text{area}} = \int\_0^T t \left[ |\Delta f\_1| + |\Delta f\_2| + |\Delta V\_{t1}| + |\Delta V\_{t2}| + |\Delta P\_{ptic12}| \right] dt \tag{4}$$

JISE,three-area = 2 2 2 2 2 2 2 2 2

JIAE,three-area =

1 2 3 1 2 3 1 2 3 <sup>0</sup> [| | | | | | | | | | | <sup>|</sup> | | | | | |] *<sup>T</sup>*

JITAE,three-area <sup>=</sup>1 2 3 1 2 3 1 2 3 <sup>0</sup> [| | | | | | | | | | | | | | | | | |] *<sup>T</sup>*

JITSE,three-area = 2 2 2 2 2 2 2 2 2

For three-area IPS, we can write:

*T*

*T*

where,

*T*

1 2 3 1 2 3 1 2 3 <sup>0</sup> [ ]

*t t t ptie ptie ptie f f f V V V P P P dt* (8)

1 2 3 1 2 3 1 2 3 <sup>0</sup> [ ]

$$\mathcal{J}\_{\text{TISE,two-area}} = \int\_0^T t \left[ \Delta f\_1^2 + \Delta f\_2^2 + \Delta V\_{t1}^2 + \Delta V\_{t2}^2 + \Delta P\_{tiel12}^2 \right] dt \tag{5}$$

$$\mathbf{J}\_{\text{IAE,two-area}} = \int\_0^T [|\Delta f\_1| + |\Delta f\_2| + |\Delta V\_{t1}| + |\Delta V\_{t2}| + |\Delta P\_{ptie12}|]dt\tag{6}$$

For three-area IPS, we can write:

$$\mathbf{J}\_{\text{ISE, three-area}} = \int\_0^T [\Delta f\_1^2 + \Delta f\_2^2 + \Delta f\_3^2 + \Delta V\_{t1}^2 + \Delta V\_{t2}^2 + \Delta V\_{t3}^2 + \Delta P\_{\text{tie1}}^2 + \Delta P\_{\text{tie2}}^2 + \Delta P\_{\text{tie3}}^2]dt\tag{7}$$

$$\begin{aligned} \mathbf{J}\_{\text{IAE,three}} &= \\ \int\_0^T [|\Delta f\_1| + |\Delta f\_2| + |\Delta f\_3| + |\Delta V\_{I1}| + |\Delta V\_{I2}| + |\Delta V\_{I3}| + |\Delta P\_{\text{pitch}1}| + |\Delta P\_{\text{pitch}2}| + |\Delta P\_{\text{pitch}3}|] dt \end{aligned} \tag{8}$$

$$\mathcal{J}\_{\text{TISE, three-area}} = \int\_0^T t \left[ \Delta f\_1^2 + \Delta f\_2^2 + \Delta f\_3^2 + \Delta V\_{t1}^2 + \Delta V\_{t2}^2 + \Delta V\_{t3}^2 + \Delta P\_{\text{tie1}}^2 + \Delta P\_{\text{tie2}}^2 + \Delta P\_{\text{tie3}}^2 \right] dt \tag{9}$$

$$\mathbf{J}\_{\text{ITAE},\text{fittree}-\text{area}} = \int\_0^T t \left[ |\Delta f\_1| + |\Delta f\_2| + |\Delta f\_3| + |\Delta V\_{\text{II}}| + |\Delta V\_{\text{II}}| + |\Delta V\_{\text{B}}| + |\Delta P\_{\text{pitch}}| + |\Delta P\_{\text{pitch}}| + |\Delta P\_{\text{pitch}}| \right] dt \tag{10}$$

where,

$$
\Delta V\_{t1} = V\_{ref} - V\_{t1} \\
\Delta V\_{t2} = V\_{ref} - V\_{t2} \\
\Delta V\_{t3} = V\_{ref} - V\_{t3} \tag{11}
$$

$$
\Delta P\_{\text{ptie1}} = \Delta P\_{\text{ptie1}2} + \Delta P\_{\text{ptie13}} \\
\Delta P\_{\text{ptie2}} = \Delta P\_{\text{ptie21}} + \Delta P\_{\text{ptie23}} \\
\Delta P\_{\text{ptie3}} = \Delta P\_{\text{ptie31}} + \Delta P\_{\text{ptie32}} \tag{12}
$$

When the cost function is minimized, the algorithm returns the best optimum parameters of the controller. To optimize the cost function (J), nature-inspired computation algorithms including LPBO, AOA, and MPSO were adapted.

#### **4. Nature-Inspired Computation Algorithms**

Due to their ability to solve complex valued problems, nature-inspired computation algorithms have gained brilliant attention in IPS. Keeping in view their remarkable contribution, an effort was made in this research to optimize the combined LFC and AVR-based IPS using nature-inspired computation techniques.

#### *4.1. Learner Performance-Based Behavior Optimization*

Rashid and Rahman presented a novel nature-inspired learner performance-based behavior optimization (LPBO) technique in 2020. The basic concept behind this algorithm is based on the fact that how students are admitted to different departments of a university is based on their high school performance. After admission, students must be able to improve their intellectual level to improve their skills. In this way, both exploitation and exploration phases are preserved. In this algorithm, a random population is generated with various ranges of grade point average (GPA). The applications of some of these learners will be rejected or accepted based on their fitness. After that, the population is divided in to subpopulation. Fitness is calculated and is then sorted into separate groups. The new population's structure is changed using crossover and mutation operators. A specified number of learners is acquired by different departments based on the minimum GPA criteria. This rejection and acceptance process is continued until all departments have their vacancies filled. Population fitness is improved in each iteration based on group learning, intellectual level, and teaching level [41]. Figure 3 presents the flow chart of the LPBO algorithm. Note that the LPBO population represents the PI-PD controller's parameters in this case.

**Figure 3.** Flow chart of LPBO algorithm. **Figure 3.** Flow chart of LPBO algorithm.

#### *4.2. Archimedes Optimization Algorithm 4.2. Archimedes Optimization Algorithm*

The Archimedes optimization algorithm (AOA) is a new state-of-the-art algorithm based on the Archimedes principle. It deals with both convex and non-convex problems. It was invented in 2021 by Fatmaand Houssein. It defines the relationship between a buoyant force and an object submerged in water. The object will sink if the displaced fluid weight is less than the weight of the object. Similarly, if the displaced fluid and object weight are equal, the object floats on the fluid. An object has volume, acceleration, and density that results in the buoyancy force, as a result fluid's net force is always zero. AOA is a very effective nature-inspired algorithm in a way that it analyzes a problem with a global optimum solution. AOA fences in both exploitation and exploration phases since it is a global optimization algorithm. A comprehensive area must be examined to identify the global optimum solution of a given problem. Firstly, the fluid's random position is initialized, and then AOA evaluates the initial population fitness to discover the best possible solution until the selection criteria are met. The density and volume of each object changes at each AOA iteration. The new density, volume, and acceleration are obtained using the object's fitness. The AOA population represents the PI-PD controller's parameters [42]. Figure 4 presents a flow chart diagram of AOA. The Archimedes optimization algorithm (AOA) is a new state-of-the-art algorithm based on the Archimedes principle. It deals with both convex and non-convex problems. It was invented in 2021 by Fatmaand Houssein. It defines the relationship between a buoyant force and an object submerged in water. The object will sink if the displaced fluid weight is less than the weight of the object. Similarly, if the displaced fluid and object weight are equal, the object floats on the fluid. An object has volume, acceleration, and density that results in the buoyancy force, as a result fluid's net force is always zero. AOA is a very effective nature-inspired algorithm in a way that it analyzes a problem with a global optimum solution. AOA fences in both exploitation and exploration phases since it is a global optimization algorithm. A comprehensive area must be examined to identify the global optimum solution of a given problem. Firstly, the fluid's random position is initialized, and then AOA evaluates the initial population fitness to discover the best possible solution until the selection criteria are met. The density and volume of each object changes at each AOA iteration. The new density, volume, and acceleration are obtained using the object's fitness. The AOA population represents the PI-PD controller's parameters [42]. Figure 4 presents a flow chart diagram of AOA.

**Figure 4. Figure 4.** Flow Flow chart of AOA. chart of AOA.

#### *4.3. Modified Particle Swarm Optimization (PSO) 4.3. Modified Particle Swarm Optimization (PSO)*

Drawing inspiration from swarm intelligence, Eberhart and Kennedy proposed the particle swarm optimization (PSO) algorithm in 1995. In PSO, the movement of particles (candidate solutions) over a defined search space depends upon their velocity and position. The movement of particles is incited by the best possible positions known as local bests. These local bests lead particles toward the best possible position [51]. In modified particle swarm optimization (MPSO), the global learning coefficient is updated using a combination of existing local and global learning coefficients. The modification in the PSO algorithm is being made to improve the convergence characteristics of the controller. Figure 5 depicts the flow chart of the MPSO algorithm. Remember that in this research work, the particles represent the PI-PD controller's parameters. Drawing inspiration from swarm intelligence, Eberhart and Kennedy proposed the particle swarm optimization (PSO) algorithm in 1995. In PSO, the movement of particles (candidate solutions) over a defined search space depends upon their velocity and position. The movement of particles is incited by the best possible positions known as local bests. These local bests lead particles toward the best possible position [51]. In modified particle swarm optimization (MPSO), the global learning coefficient is updated using a combination of existing local and global learning coefficients. The modification in the PSO algorithm is being made to improve the convergence characteristics of the controller. Figure 5 depicts the flow chart of the MPSO algorithm. Remember that in this research work, the particles represent the PI-PD controller's parameters.

**Figure 5.** Flow chart of MPSO algorithm. **Figure 5.** Flow chart of MPSO algorithm.

#### **5. Implementation and Results Discussion 5. Implementation and Results Discussion**

Multiple simulations were carried out in MATLAB/Simulink to express the validation of the proposed control methodology. Firstly, a two-area, two-source IPS with combined LFC and AVR was optimized using LPBO, AOA, and MPSO. ITSE was chosen as the error criterion, due to efficient error convergence characteristics. After achieving Multiple simulations were carried out in MATLAB/Simulink to express the validation of the proposed control methodology. Firstly, a two-area, two-source IPS with combined LFC and AVR was optimized using LPBO, AOA, and MPSO. ITSE was chosen as the error criterion, due to efficient error convergence characteristics. After achieving successful

results, the proposed methodology was applied to a three-area, three-source IPS with LFC-AVR loops. IPS with LFC-AVR loops. *5.1. Optimization of Two-Area Interconnected Power System*

successful results, the proposed methodology was applied to a three-area, three-source

#### *5.1. Optimization of Two-Area Interconnected Power System* The two-area IPS model under investigation with a collective LFC-AVR system is

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 12 of 32

The two-area IPS model under investigation with a collective LFC-AVR system is shown in Figure 6. The system parameters of the two-area IPS are specified in Appendix A. The system parameters of area-1 and area-2 were chosen from [5] for a direct comparison of the proposed methodology with the NLTA-PID controller. The parameters of optimization algorithms such as MPSO, LPBO, and AOA used in simulations are given in Table 3. The tie–line connection between area-1 and are-2 can be established using Figure 7. The optimal parameters of MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control schemes are given in Table 4. For the sake of the assessment of the proposed control schemes, the evaluation of the time response of each schemes was carried out and comparisons were made with the results of NLTA-PID [5]. Further, a comparison between the proposed control schemes such as MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD is also presented in detail in this section. shown in Figure 6. The system parameters of the two-area IPS are specified in Appendix A. The system parameters of area-1 and area-2 were chosen from [5] for a direct comparison of the proposed methodology with the NLTA-PID controller. The parameters of optimization algorithms such as MPSO, LPBO, and AOA used in simulations are given in Table 3. The tie–line connection between area-1 and are-2 can be established using Figure 7. The optimal parameters of MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control schemes are given in Table 4. For the sake of the assessment of the proposed control schemes, the evaluation of the time response of each schemes was carried out and comparisons were made with the results of NLTA-PID [5]. Further, a comparison between the proposed control schemes such as MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD is also presented in detail in this section.

**Figure 6.** *Cont*.

**Figure 6.** Two-area IPS with combined LFC−AVR. **Figure 6.** Two-area IPS with combined LFC−AVR.

**Table 3.** Parameters of optimization techniques.


centage

centage

tants

centage

centage

tants

**MPSO LPBO AOA** Parameter Value Parameter Value Parameter Value Population size 20 Population size 20 Population size 20 Iterations 10 Iterations 10 Iterations 10

efficient 2.88 Mutation Rate 0.03 C3 (constant) <sup>2</sup>

0.7 C1 (constant) 2

0.3 C2 (constant) 6

6 C4 (constant) 0.5

0.7 C1 (constant) 2

0.3 C2 (constant) 6

6 C4 (constant) 0.5

ization (u,l) 0.9, 0.1

ization (u,l) 0.9, 0.1

spring <sup>14</sup> Range of Normal-

spring <sup>14</sup> Range of Normal-

$$\begin{aligned} \xleftarrow{\Delta f\_1} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star} \xleftarrow{\star}} \end{aligned}$$

Global Learning Coefficient 2.88 Mutation Rate 0.03 C3 (constant) <sup>2</sup> Max. Velocity Limit 0.2 Number of Mu-

**Figure 7.** Tie–line connection. **Figure 7.** Tie–line connection.

Inertia Weight

Personal Learning

Global Learning Co-

Coefficient 2.74 Mutation Per-

Max. Velocity Limit 0.2 Number of Mu-

Min. Velocity Limit −0.2 Number of Off-

Min. Velocity Limit −0.2 Number of Off-

**Table 3.** Parameters of optimization techniques.


**Table 4.** Optimal values of controller parameters (area-2).

Figure 8 shows the frequency deviation curves of area-1 and area-2 using NLTA-PID [5], MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It can be seen that the proposed control schemes provided a very satisfactory frequency deviation response. For the area-1 LFC, the settling time of NLTA-PID [5] was lower than the proposed schemes but at the cost of a high undershoot. NLTA-PID provided an undershoot of −0.285, whereas the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD provided−0.130, −0.135, and −0.115, respectively. It can be noticed that the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD provided 54%, 52.6%, and 60%, respectively, better undershoot responses as compared to the NLTA-PID controller in area-1. For area-2, NLTA-PD provided a quick settling, but it provided an undershoot of −0.275, whereas the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD provided−0.135, −0.170, and −0.120, respectively. It was verified that the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD provided 51%, 38%, and 56%, respectively, better undershoot responses as compared to the NLTA-PID controller. The percentages of overshoots and steady state (s-s) errors were almost zero with each proposed technique. *Sustainability* **2022**, *14*, x FOR PEER REVIEW 15 of 32

**Figure 8.** LFC response with PI−PD control scheme. (**a**) ∆*f*1; (**b**) ∆*f*2. **Figure 8.** LFC response with PI−PD control scheme. (**a**) ∆*f* <sup>1</sup> ; (**b**) ∆*f* <sup>2</sup> .

Figure 9 shows the terminal voltage of area-1 and area-2 using the NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respec-

overshoot in area-1 and area-2, respectively, but the proposed technique provided a negligible overshoot percentage at the cost of the settling time with all tuning techniques. It can be observed that the proposed LPBO-PI-PD and AOA-PI-PD control schemes pro-

Figure 10 shows the tie–line power response using NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It can be observed from the results that LPBO-PI-PD and AOA-PI-PD provided tie–line power responses with no undershoot; however, this was at the cost of a slightly small overshoot. In addition, the tie–line power responses yielded by MPSO-PI-PD, LPBO-PI-PD, and

duced settling times approximately the same as those achieved by NLTA-PID.

(**a**) (**b**) **Figure 9.** AVR response with PI-PD control scheme. (**a**) Vt1; (**b**) Vt2.

AOA-PI-PD were satisfactory.

Figure 9 shows the terminal voltage of area-1 and area-2 using the NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It is clear that the proposed control schemes provided a very satisfactory transient response in both area-1 and area-2. It is identified that NLTA-PID provided 18% and 17% overshoot in area-1 and area-2, respectively, but the proposed technique provided a negligible overshoot percentage at the cost of the settling time with all tuning techniques. It can be observed that the proposed LPBO-PI-PD and AOA-PI-PD control schemes produced settling times approximately the same as those achieved by NLTA-PID. Figure 9 shows the terminal voltage of area-1 and area-2 using the NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It is clear that the proposed control schemes provided a very satisfactory transient response in both area-1 and area-2. It is identified that NLTA-PID provided 18% and 17% overshoot in area-1 and area-2, respectively, but the proposed technique provided a negligible overshoot percentage at the cost of the settling time with all tuning techniques. It can be observed that the proposed LPBO-PI-PD and AOA-PI-PD control schemes produced settling times approximately the same as those achieved by NLTA-PID.

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 15 of 32

(**a**) (**b**) **Figure 8.** LFC response with PI−PD control scheme. (**a**) ∆*f*1; (**b**) ∆*f*2.

**Figure 9.** AVR response with PI-PD control scheme. (**a**) Vt1; (**b**) Vt2. **Figure 9.** AVR response with PI-PD control scheme. (**a**) Vt1; (**b**) Vt2.

Figure 10 shows the tie–line power response using NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It can be observed from the results that LPBO-PI-PD and AOA-PI-PD provided tie–line power responses with no undershoot; however, this was at the cost of a slightly small overshoot. In addition, the tie–line power responses yielded by MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD were satisfactory. Figure 10 shows the tie–line power response using NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It can be observed from the results that LPBO-PI-PD and AOA-PI-PD provided tie–line power responses with no undershoot; however, this was at the cost of a slightly small overshoot. In addition, the tie–line power responses yielded by MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD were satisfactory. *Sustainability* **2022**, *14*, x FOR PEER REVIEW 16 of 32

**Figure 10.** Tie–line power response in two−area IPS with combined LFC−AVR. **Figure 10.** Tie–line power response in two−area IPS with combined LFC−AVR.

**Error**

Tables 5 and 6 show the summary of LFC and AVR responses using NLTA-PID, MPSO-PI-PD, LPBO-PI-PD and AOA-PI-PD control schemes in a two-area IPS, respectively. Tables 5 and 6 show the summary of LFC and AVR responses using NLTA-PID, MPSO-PI-PD, LPBO-PI-PD and AOA-PI-PD control schemes in a two-area IPS, respectively.

**Settling** 

**Time % Overshoot Undershoot s-s**

**Time % Overshoot s-s**

**Error**

**Error**

**Area-1 Area-2**

**Area-1 Area-2**

**Rise Time**

Figure 11 shows the graphical comparison of the performance parameters of NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It is very clear that the proposed PI-PD control schemes provided relatively better responses in terms of the undershoot in LFC and overshoot percentage in AVR as compared to the NLTA-PID controller. From Tables 5 and 6 and Figure 11, it is concluded that the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD were effective for maintaining the frequency and voltage within the prescribed values with a satisfactory perfor-

**Settling** 

**Error**

MPSO-PI-PD 4.5407 0 −0.13 0 4.92 0 −0.135 0 LPBO-PI-PD 6.9478 0.005 −0.135 0 4.043 0 −0.17 0 AOA-PI-PD 6.6752 0 −0.115 0 4.69 0 −0.12 0

NLTA-PID [5] 0.1287 1.24 18.80 0 0.154 0.887 17.75 0 MPSO-PI-PD 0.6532 3.30 0 0 1.077 3.17 3.2971 × 10−4 0 LPBO-PI-PD 0.4546 1.22 0.28 0 0.464 1.381 0 0 AOA-PI-PD 0.610 1.23 0.27 0 0.435 1.499 0 0

**Table 5.** LFC responses (area-2).

**Table 6.** AVR responses (area-2).

**Time % Overshoot s-s**

**Settling** 

mance in a two-area IPS.

**Control Scheme Settling** 

**Control Scheme**

**Rise Time**


**Table 5.** LFC responses (area-2).



Figure 11 shows the graphical comparison of the performance parameters of NLTA-PID, MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a two-area IPS, respectively. It is very clear that the proposed PI-PD control schemes provided relatively better responses in terms of the undershoot in LFC and overshoot percentage in AVR as compared to the NLTA-PID controller. From Tables 5 and 6 and Figure 11, it is concluded that the proposed MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD were effective for maintaining the frequency and voltage within the prescribed values with a satisfactory performance in a two-area IPS.

#### *5.2. Three-Area, Three-Source System*

In this section, the proposed methodology is applied to a three-area IPS model with combined LFC-AVR. The model under study is presented in Figure 12, while the model parameters are provided in Appendix B.

The optimal values of MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD for a three-area IPS with combined LFC and AVR are given in Table 7. Figure 13 shows the frequency deviation response using MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a three-area IPS, respectively.




**Figure 11.** Graphical comparison of performance parameters. (**a**) ∆*f*1; (**b**) ∆*f*2; (**c**) Vt1; (**d**) Vt2. **Figure 11.** Graphical comparison of performance parameters. (**a**) ∆*f* <sup>1</sup> ; (**b**) ∆*f* <sup>2</sup> ; (**c**) Vt1; (**d**) Vt2.

In this section, the proposed methodology is applied to a three-area IPS model with combined LFC-AVR. The model under study is presented in Figure 12, while the model

*5.2. Three-Area, Three-Source System*

parameters are provided in Appendix B.

Area-1

Area-2

(**a**) (**b**)

NLTA-PID

Modified PSO-PIPD AOA-PIPD

LPBO-PIPD

NLTA-PID

Modified PSO-PIPD LPBO-PIPD

AOA-PIPD

NLTA-PID

Modified PSO-PIPD AOA-PIPD

LPBO-PIPD

NLTA-PID

Modified PSO-PIPD LPBO-PIPD

(**c**) (**d**)

AOA-PIPD <sup>0</sup>

*5.2. Three-Area, Three-Source System*

parameters are provided in Appendix B.

**Figure 11.** Graphical comparison of performance parameters. (**a**) ∆*f*1; (**b**) ∆*f*2; (**c**) Vt1; (**d**) Vt2.

In this section, the proposed methodology is applied to a three-area IPS model with combined LFC-AVR. The model under study is presented in Figure 12, while the model

**Figure 12.** Three−area IPS with combined LFC−AVR. **Figure 12.** Three−area IPS with combined LFC−AVR.

in a three-area IPS, respectively.

**Table 7.** Optimal values of controller parameters (area-3).

**Area Controller Parameters Proposed Control Schemes**

The optimal values of MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD for a three-area IPS with combined LFC and AVR are given in Table 7. Figure 13 shows the frequency deviation response using MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques

Kp2 1.2737 0.96 −0.38 Kd1 0.831 0.53 0.55 Kp3 1.5371 1.56 0.88 Ki2 1.965 1.62 1.91 Kp4 1.2543 0.85 1.13 Kd2 0.5936 0.56 0.5

Kp5 1.1106 0.77 0.86 Ki3 0.9076 0.61 0.71 Kp6 0.8639 1.48 1.55 Kd3 1.3118 1.03 0.86

**MPSO-PI-PD LPBO-PI-PD AOA-PI-PD**

Area-3

Kp7 1.7917 1.68 1.91 Ki4 1.8286 1.57 1.97 Kp8 0.9068 0.83 1.074 Kd4 0.6882 0.73 1.071

Kp9 0.7914 0.78 1.9 Ki5 1.0795 1.12 1.26 Kp10 1.2741 0.66 1.64 Kd5 0.8581 1.56 0.42 Kp11 1.2282 1.29 1.63 Ki6 1.4326 1.3 1.69 Kp12 0.9527 0.77 1.43 Kd6 0.5874 0.45 1.33 ITSE 0.3507 0.34485 0.4853

**Figure 13.** LFC response with PI−PD control scheme. (**a**) ∆*f*1; (**b**) ∆*f*2; (**c**) ∆*f*3. **Figure 13.** LFC response with PI−PD control scheme. (**a**) ∆*f* <sup>1</sup> ; (**b**) ∆*f* <sup>2</sup> ; (**c**) ∆*f* <sup>3</sup> .

Table 8 shows the summary of LFC responses of area-1, area-2, and area-3 using MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques, respectively. For area-1, LPBO-PI-PD provided 14% and 31% quick settling times as compared to the MPSO-PI-PD and AOA-PI-PD control schemes, respectively. The overshoot percentage and steady state error were zero in each case. Table <sup>8</sup> shows the summary of LFC responses of area-1, area-2, and area-3 usingMPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques, respectively. For area-1, LPBO-PI-PD provided 14% and 31% quick settling times as compared to the MPSO-PI-PDand AOA-PI-PD control schemes, respectively. The overshoot percentage and steady state error were zero in each case.


**Table 8.** LFC responses (area-3).

Further, MPSO-PI-PD exhibited 30% and 20% better undershoot responses as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. For area-2, LPBO-PI-PD yielded 3.3% and 11% quick settling times, as compared to the MPSO-PI-PD and AOA-PI-PD control schemes, respectively. The overshoot percentage and steady state

error were zero in each case. Further, AOA-PI-PD exhibited 4.16% and 5.74% better undershoot responses as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. For area-3, MPSO-PI-PD provided 25% and 16% quick settling times as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. The overshoot percentage and steady state error were again zero in each case. Further, AOA-PI-PD exhibited 22% and 34% better undershoot responses as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Figure 14 shows the terminal voltage responses of area-1, area-2, and area-3 using MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a three-area IPS, respectively. error were zero in each case. Further, AOA-PI-PD exhibited 4.16% and 5.74% better undershoot responses as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. For area-3, MPSO-PI-PD provided 25% and 16% quick settling times as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. The overshoot percentage and steady state error were again zero in each case. Further, AOA-PI-PD exhibited 22% and 34% better undershoot responses as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Figure 14 shows the terminal voltage responses of area-1, area-2, and area-3 using MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a three-area IPS, respectively.

Further, MPSO-PI-PD exhibited 30% and 20% better undershoot responses as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. For area-2, LPBO-PI-PD yielded 3.3% and 11% quick settling times, as compared to the MPSO-PI-PD and AOA-PI-PD control schemes, respectively. The overshoot percentage and steady state

**Time % Overshoot Undershoot s-s**

MPSO-PI-PD 5.43 0 −0.14 0 LPBO-PI-PD 4.65 0 −0.20 0 AOA-PI-PD 6.73 0 −0.175 0

MPSO-PI-PD 5.04 0 −0.120 0 LPBO-PI-PD 4.87 0 −0.122 0 AOA-PI-PD 5.46 0 −0.115 0

PSO-PI-PD 5.40 0 −0.122 0 LPBO-PI-PD 7.16 0 −0.143 0 AOA-PI-PD 6.40 0 −0.095 0

**Error**

*Sustainability* **2022**, *14*, x FOR PEER REVIEW 20 of 32

**Area Control Scheme Settling** 

**Table 8.** LFC responses (area-3).

Area-1

Area-2

Area-3

**Figure 14.** AVR response with PI-PD control scheme. (**a**) Vt1; (**b**) Vt2; (**c**) Vt3. **Figure 14.** AVR response with PI-PD control scheme. (**a**) Vt1; (**b**) Vt2; (**c**) Vt3.

Table 9 shows the summary of AVR responses of area-1, area-2, and area-3 using the MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control schemes, respectively. For area-1, AOA-PI-PD provided 26% and 2% quick rise times as compared to the MPSO-PI-PD and LPBO-PI-PD control techniques, respectively. Moreover, AOA-PI-PD yielded 38% and 29% fast settling times as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, it was observed that the percentage of overshoot and steady state error were almost zero with each tuning technique. For area-2, MPSO-PI-PD offered 3% and 13% quick rise times as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. Moreover, AOA-PI-PD provided 21% and 19% fast settling times as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, it Table 9 shows the summary of AVR responses of area-1, area-2, and area-3 using the MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control schemes, respectively. For area-1, AOA-PI-PD provided 26% and 2% quick rise times as compared to the MPSO-PI-PD and LPBO-PI-PD control techniques, respectively. Moreover, AOA-PI-PD yielded 38% and 29% fast settling times as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, it was observed that the percentage of overshoot and steady state error were almost zero with each tuning technique. For area-2, MPSO-PI-PD offered 3% and 13% quick rise times as compared to the LPBO-PI-PD and AOA-PI-PD control techniques, respectively. Moreover, AOA-PI-PD provided 21% and 19% fast settling times as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, it can be seen

can be seen that the overshoot percentage and steady state error were almost zero with each tuning technique. For area-3, LPBO-PI-PD produced 64% and 73% quick rise times

PI-PD and LPBO-PI-PD control schemes,respectively. Further, it can be seen that the overshoot percentage and steady state error were negligible with each tuning technique. Figure 15 shows the graphical comparison of the performance parameters of the MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a three-area interconnected sys-

> MPSO-PI-PD 1.53 3.48 5.8225 × 10−6 0 LPBO-PI-PD 1.15 3.01 4.5973 × 10−4 0 AOA-PI-PD 1.13 2.15 0.083 0

> MPSO-PI-PD 0.95 2.44 0 0 LPBO-PI-PD 0.98 2.37 0 0 AOA-PI-PD 1.09 1.92 0.37 0

> MPSO-PI-PD 1.32 3.30 0 0 LPBO-PI-PD 0.48 3.48 0.001 0 AOA-PI-PD 1.75 3.29 0 0

**Time % Overshoot s-s**

**Error**

**Area Control Scheme Rise Time Settling** 

tem.

**Table 9.** AVR responses (area-3).

Area-1

Area-2

Area-3

that the overshoot percentage and steady state error were almost zero with each tuning technique. For area-3, LPBO-PI-PD produced 64% and 73% quick rise times as compared to the MPSO-PI-PD and AOA-PI-PD control techniques, respectively. Moreover, AOA-PI-PD provided 0.3% and 5.45% fast settling times as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, it can be seen that the overshoot percentage and steady state error were negligible with each tuning technique. Figure 15 shows the graphical comparison of the performance parameters of the MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD control techniques in a three-area interconnected system.


**Table 9.** AVR responses (area-3).

Figure 16 shows the tie–line power responses of area-1, area-2, and area-3 using the MPSO-PI-PD, LPBO-PI-PD and AOA-PI-PD control schemes in a three-area IPS, respectively. It can be inferred that PI-PD-based control schemes including MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD yielded satisfactory tie–line powers responses with negligible undershoots and overshoot percentages in the three-area IPS.

#### *5.3. Sensitivity Analysis*

In this section, the robustness of the proposed nature-inspired computation-based PI-PD control techniques were tested with large variations in the system parameters of the three-area IPS with combined LFC-AVR. The generator time constant (Tg) and turbine time constant (Tt) were varied to Â ± 50% of their nominal values. The newer values of T<sup>g</sup> and T<sup>t</sup> after Â ± 50% variations are given in Appendix B. The optimum parameters of the PI-PD control scheme were the same as those used in Case 2. The AVR and LFC responses of the PI-PD control scheme with variations in T<sup>t</sup> and T<sup>g</sup> are depicted in Figures 17 and 18, respectively. Tables 10 and 11 show the summary of the performance parameters of LFC and AVR responses under parametric variations. From the obtained results, it is evident the overshoot percentages and steady state error were almost zero in each case. The AVR responses are almost indistinguishable to each other, despite the variation in system parameters. Figure 19 shows the graphical comparison of the performance parameters under this scenario. It is clearly observed that the system response under Â ± 50% variations was very identical to response with nominal values. This indicates that the proposed LPBO-PI-PD control technique was very realistic and robust under variations in the system parameters. These results clearly reveal that the re-tuning of the proposed controller is not necessary with large variations of at least Â ± 50%.

**Figure 15.** Graphical comparison of performance parameters. (**a**) ∆*f*1; (**b**) ∆*f*2; (**c**) ∆*f*3;(**d**) Vt1; (**e**) Vt2; (**f**) Vt3. **Figure 15.** Graphical comparison of performance parameters. (**a**) ∆*f* <sup>1</sup> ; (**b**) ∆*f* <sup>2</sup> ; (**c**) ∆*f* <sup>3</sup> ;(**d**) Vt1; (**e**) Vt2; (**f**) Vt3.

Figure 16 shows the tie–line power responses of area-1, area-2, and area-3 using the MPSO-PI-PD, LPBO-PI-PD and AOA-PI-PD control schemes in a three-area IPS, respectively. It can be inferred that PI-PD-based control schemes including MPSO-PI-PD, LPBO-PI-PD, and AOA-PI-PD yielded satisfactory tie–line powers responses with negligible un-

dershoots and overshoot percentages in the three-area IPS.

**Figure 16.** Tie−line power responses with PI−PD control scheme. (**a**) ∆*P*tie1; (**b**) ∆*P*tie2; (**c**) ∆*P*tie3. **Figure 16.** Tie−line power responses with PI−PD control scheme. (**a**) ∆*P*tie1; (**b**) ∆*P*tie2; (**c**) ∆*P*tie3.

*5.3. Sensitivity Analysis* **Table 10.** Settling time responses of PI-PD control scheme with variations in system parameters.


and AVR responses under parametric variations. From the obtained results, it is evident the overshoot percentages and steady state error were almost zero in each case. The AVR responses are almost indistinguishable to each other, despite the variation in system pa-**Table 11.** Overshoot and undershoot responses of PI-PD control scheme with variations in system parameters.


**Figure 17.** LFC and AVR responses with variations in Tt. (**a**) ∆*f*1 in area−1; (**b**) ∆*f*2 in area−2; (**c**) ∆*f*3 in area−3; (**d**) Vt1 in area−1; (**e**) Vt2 in area−2 and (**f**) Vt3 in area−3. **Figure 17.** LFC and AVR responses with variations in T<sup>t</sup> . (**a**) ∆*f* <sup>1</sup> in area−1; (**b**) ∆*f* <sup>2</sup> in area−2; (**c**) ∆*f* <sup>3</sup> in area−3; (**d**) Vt1 in area−1; (**e**) Vt2 in area−2 and (**f**) Vt3 in area−3.

**Figure 18.** LFC and AVR responses with variations in Tg. (**a**) ∆*f*1 in area−1; (**b**) ∆*f*2 in area−2; (**c**) ∆*f*3 in area−3; (**d**) Vt1 in area−1; (**e**) Vt2 in area−2 and (**f**) Vt3 in area−3. **Figure 18.** LFC and AVR responses with variations in Tg. (**a**) <sup>∆</sup>*<sup>f</sup>* <sup>1</sup> in area−1; (**b**) ∆*f* <sup>2</sup> in area−2; (**c**) ∆*f* <sup>3</sup> in area−3; (**d**) Vt1 in area−1; (**e**) Vt2 in area−2 and (**f**) Vt3 in area−3.

**Figure 19.** Graphical comparison of performance parameters in three-area IPS with variations in Tg and Tt. (**a**) Settling time; (**b**) undershoot. **Figure 19.** Graphical comparison of performance parameters in three-area IPS with variations in Tg and Tt . (**a**) Settling time; (**b**) undershoot.

#### **Table 10.** Settling time responses of PI-PD control scheme with variations in system parameters. **6. Conclusions and Future Work**

**Parameters/Variation Settling Time (LFC and AVR) ∆***f1* **∆***f2* **∆***f3 Vt1 Vt2 Vt3* Nominal Tg, Tt 4.65 4.87 7.16 3.01 2.37 3.48 Tg1, Tg2, Tg3/+50% 4.60 4.76 7.02 2.74 2.11 3.56 Tg1, Tg2, Tg3/−50% 4.71 4.95 7.32 3.25 2.59 3.56 Tt1, Tt2, Tt3/+50% 4.63 5.01 7.18 3.03 2.38 3.48 Tt1, Tt2, Tt3/−50% 4.60 4.71 7.11 2.99 2.36 3.48 The multi-area IPS included numerous control areas, which are connected through the AC tie–line. The transient and steady state performance of a multi-area IPS with AVR-LFC was thoroughly investigated in this research.Nature-inspired computation including MPSO-, LPBO-, and AOA-based PI-PD control technique was proposed for the optimization of the multi-area system. From the findings of Case 1, it is concluded that all proposed schemes provided relatively better undershoot responses as compared to the NLTA-PID controller [5] for LFC. Particularly, the AOA-PI-PD control scheme exhibited 60% and 56% better undershoots in the area-1 and area-2 LFC, respectively, as compared to the NLTA-PID controller. Similarly, NLTA-PID provided 18% and 17% overshoot in the area-1 and area-2 AVR, respectively, but the proposed PI-PD control scheme completely eliminated the overshoot percentage with each tuning algorithm. The results of Case 2 reveal that LPBO-PI-PD provided 14% and 31% quick settling times in area-1, whereas 3.3% and

11% quick settling times were provided in the area-2 LFC as compared to the MPSO-PI-PD and AOA-PI-PD control techniques, respectively. In the area-3 LFC, MPSO-PI-PD provided relatively lower settling times (25% and 16%) as compared to the LPBO-PI-PD and AOA-PI-PD control schemes, respectively. MPSO-PI-PD provided relatively better undershoot (30% and 20%) in the area-1 LFC, whereas the AOA-PI-PD control technique provided better undershoot in the area-1 (4.16% and 5.74%) and area-2 (22% and 34%) LFC, respectively. Moreover, the AOA-PI-PD control scheme provided 26% and 2% quick rise times, and 38% and 29% fast settling times in the area-1 AVR as compared to the MPSO-PI-PD and LPBO-PI-PD control schemes, respectively. Further, MPSO-PI-PD provided 3% and 13% fast rise times in the area-2 AVR as compared to the LPBO-PI-PD and AOA-PI-PD control schemes, respectively. AOA-PI-PD provided 21% and 19% quick settling times in the area-2 AVR, and 0.3% and 5.45% fast settling times in the area-3 AVR as compared to the MPSO-PI-PD and AOA-PI-PD control schemes, respectively. For the area-3 AVR, LPBO-PI-PD provided 64% and 73% fast rise times as compared to the MPSO-PI-PD and AOA-PI-PD control schemes, respectively. Finally, the resilience of the PI-PD control technique was assessed by varying the system parameters (Â ± 50%), and a comprehensive sensitivity analysis was carried out to confirm its robustness. The results confirm the superiority of the proposed PI-PD control scheme when applied to multi-area IPS with combined LFC and AVR. Keeping in mind the value of the present work, IPS with a combined LFC-AVR can be analyzed by incorporating multi-source and various energy storage devices to enhance the dynamic response of the power systems. Further, neurofuzzy and hybrid ANN controllers can also be utilized for multi-area multi-sources IPS. It will be worth employing PI-PD, neuro-fuzzy, or hybrid ANN to multi-area IPS under nonlinearity constraints. Moreover, very recently introduced nature-inspired computing techniques such as dandelion optimization, artificial rabbits optimization, and sea-horse optimization can be explored to find the optimal parameters of controllers in such types of application.

**Author Contributions:** Conceptualization, T.A., S.A.M., A.T.A.; Data curation, T.A., S.A.M., A.D., H.M.; Formal analysis, T.A., S.A.M., A.D., H.M., A.T.A., I.A.H.; Investigation, A.T.A., I.A.H.; Methodology, T.A., S.A.M., A.D., H.M., A.T.A., I.A.H.; Project administration, A.T.A., I.A.H.; Resources, T.A., S.A.M., A.D., H.M.; Supervision, A.T.A., I.A.H.; Validation, T.A., S.A.M., A.D., H.M.; Visualization, A.D., H.M., A.T.A., I.A.H.; Writing—original draft, T.A., S.A.M., A.D., H.M.; Writing—review & editing, T.A., S.A.M., A.D., H.M., A.T.A., I.A.H.; Funding Acquisition, I.A.H. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Norwegian University of Science and Technology.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Acknowledgments:** The authors would like to acknowledge the support of Norwegian University of Science and Technology for paying the Article Processing Charges (APC) of this publication. Special acknowledgment to Automated Systems & Soft Computing Lab (ASSCL), Prince Sultan University, Riyadh, Saudi Arabia. In addition, the authors wish to acknowledge the editor and anonymous.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**

**Table A1.** Area-2 [5].



**Table A1.** *Cont.*

#### **Appendix B**

**Table A2.** System parameters of area-1, area-2, and area-3.



**Table A3.** System parameters after Â ± 50% variations in T<sup>g</sup> and T<sup>t</sup> .

### **References**

