Comparative Performance Assessment of Different Energy Storage Devices in Combined LFC and AVR Analysis of Multi-Area Power System

: This paper made an attempt to put forward the comparative performance analysis of different energy storage devices (ESDs), such as redox ﬂow batteries (RFBs), superconducting magnetic energy storage (SMES) device and ultra-capacitors (UCs), in the combined frequency and voltage stabilization of a multi-area interconnected power system (MAIPS). The investigative power system model comprises two areas, and each area consists of the power-generating sources of thermal, hydro and gas units. The intelligent control mechanism of fuzzy PID was used as a secondary controller optimized with a hybridized approach of the artiﬁcial electric ﬁeld algorithm (HAEFA) subjected to the minimization of integral time absolute error (ITAE) objective function. However, the superiority of fuzzy PID in dampening the deviations of combined load frequency control (LFC) and automatic voltage regulator (AVR) responses was revealed upon comparison with conventional PI and PID. Further, the LFC-AVR combined analysis was extended to incorporate different ESDs one after the other. The simulation results reveal the efﬁcacy of incorporating ESDs with the LFC-AVR system and the supremacy of RFBs in damping out the ﬂuctuations in frequency and voltage.


Introduction
One of the central issues that practicing engineers face in electrical power systems is the simultaneous control of terminal voltage and area frequency. The impairment of any of these parameters strongly impairs the life expectancy and performance of other operating equipment associated with the power system. Controlling devices that have been installed in large complex power systems are intended to deal with small load disturbances PID fine-tuned with the HAEFA mechanism was implemented as a secondary regulatory technique in both the loops of LFC and AVR.
However, the design and implementation of the secondary regulator mechanism alone is not sufficient for the power system network to maintain stability during load perturbations [25]. These secondary regulators are only capable of restoring stability during small perturbations. So, the complex power system definitely needs an ancillary strategy as a territorial control approach to overcome the hunting mechanism of the power system apparatus during large load disturbances. Ancillary strategies include overlying the DC line with the existing AC transmission line (AC/DC links) [26] in parallel and the incorporation of ESDs [27,28] with the system. In this paper, different ESDs were incorporated with the system one after the other to assess their performance in obtaining the simultaneous stabilization of voltage and frequency. A comparison of different ESDs' performances in a multi-area combined LFC-AVR system has so far not been reported in the literature, and the intelligent fuzzy PID regulator for a combined system is also a maiden attempt.
In light of the abovementioned aspects, this paper contributes the following: a.
An LFC-AVR combined model of a two-area system was designed for simultaneous voltage and frequency stabilization. b.
Intelligent-based fuzzy PID was implemented as a secondary regulator in the loops of LFC and AVR. c.
The implementation of fuzzy PID for a combined system has not been reported in the literature so far. d. The HAEFA algorithm was utilized to tune the parameters of the secondary controller. e.
The supremacy of fuzzy PID was demonstrated with controllers of traditional PI/PID. f.
The necessity of coupling AVR with LFC loop was clearly demonstrated and justified. g.
The performance assessment of different ESDs in a combined effect was analyzed.

Methodology
A HAEFA-based fuzzy PID control approach is presented in this paper for a combined LFC-AVR analysis. To carry out the analysis, the LFC-AVR system was designed and developed in MATLAB/Simulink version of (R2016a). Disturbance of 10%SLP was injected in area-1 of test system to analyze fuzzy PID control performance. Efficacy of fuzzy PID was revealed with traditional PI/PID controllers. Moreover, the efficacy of the controller depends on the selection of optimization approach. The newest technique of HAEFA was implemented in this paper after testing it on different benchmark standard functions. Further, various ESDs were considered individually with the LFC-AVR system as territorial controllers to boost the performance. Simulation results reveal the supremacy of RFBs over other ESDs.

Modeling of LFC Loop
The power system displayed in Figure 1 considered for investigation consists of two areas of equal generation capacity. Each area comprises identical generation sources of thermal, hydro and gas units. The dynamical analysis of the power system was carried out upon laying area-1 with 10%SLP. The parameters of generation units such as gain and time constants in the investigated power system are provided in [29]. The modeling of the power-generating utilities employed in this work is given as follows: time constants in the investigated power system are provided in [29]. The modeling of the power-generating utilities employed in this work is given as follows:   (1) Hydro system: Thermal system: Hydro system: Energies 2022, 15, 629 5 of 22 Gas system:

Modeling of AVR Loop
AVR safeguards the voltage of the synchronous generator to prevent the violation of predetermined limits. The model of AVR with coupling coefficients is displayed in Figure 2. The AVR system mainly consists of a sensing unit, generator field, amplifier unit and excitation units. The sensor continuously senses the terminal unit and generates the error voltage (∆V) signal upon comparing it with the reference voltage signal. The error voltage signal can then be amplified and fed to the excitation unit to change the synchronous generator field excitation. With this the terminal voltage deviation is quickly compensated in order to make the system stable. The detailed mathematical modeling of AVR unit is as follows [30]:

Amplifier
Exciter Gen.Field

Modeling of LFC and AVR Coupling
The main aim of analyzing the effect of the combined LFC and AVR system control frequency and terminal voltage simultaneously. Frequency can be controll regulating real power generation, and terminal voltage is maintained by controllin generator field excitation system. However, as the actions taken by the AVR system in 'V' variation, there is a considerable effect of terminal voltage on real power gene [22], given by So, the automatic voltage regulator has an instantaneous considerable impact on frequency control loop. From [1] it was found that under dynamic load ) P (Δ d condi there will be deviation in the system frequency ( Δf ), resulting in a shift in the gen load angle δ) (Δ . On incorporating these effects, the change in real power generat mathematically expressed as where ' E is EMF induced in the generator stator winding, 1 K is the deviation in ' a small deviation in ' E , and the field output voltage of the generator, S P , is synchr ing the power coefficient chartered as power angle characteristics slope at initial oper Amplifier: Exciter: Generator: Sensor:

Modeling of LFC and AVR Coupling
The main aim of analyzing the effect of the combined LFC and AVR system is to control frequency and terminal voltage simultaneously. Frequency can be controlled by regulating real power generation, and terminal voltage is maintained by controlling the generator field excitation system. However, as the actions taken by the AVR system result in 'V' variation, there is a considerable effect of terminal voltage on real power generation [22], given by So, the automatic voltage regulator has an instantaneous considerable impact on load frequency control loop. From [1] it was found that under dynamic load (∆ P d ) conditions, there will be deviation in the system frequency (∆f), resulting in a shift in the generator load angle (∆ δ). On incorporating these effects, the change in real power generation is mathematically expressed as ∆P e = P S ∆δ + K 1 E where E is EMF induced in the generator stator winding, K 1 is the deviation in 'V' for a small deviation in E , and the field output voltage of the generator, P S , is synchronizing the power coefficient chartered as power angle characteristics slope at initial operating angles. The variation in voltage with respect to the alteration in the rotor angle is termed as The transfer function of the generator field considering the change in load angle is approximated by K 1 , K 3 , K 4 Coefficients are positive, but K 2 may be negative depending on total system reactance (X S ).
Where X S = X 1 + X 2 + X 12 , The K coefficients are mentioned in [30].

Controller and Objective Function
In general, classical PID regulators are widely accepted and utilized for the monitoring and regulation of power system frequency, owing to the benefits of design simplicity. However, these traditional controllers are not effective or competent for a system with uncertainties. Intelligent-based fuzzy regulators are more effective in dealing with complex realistic power system networks. Despite that, the performance of fuzzy controllers strongly relies on the design of appropriate fuzzy rule base interface system and membership functions (MFs), which is a complex job. The structure of the fuzzy-assisted PID regulator employed for the combined LFC and AVR study in this paper is displayed in Figure 3 [31]. The fuzzy logic controller (FLC) is provided with the input as area control error (ACE) and derivative of ACE. Usually, triangular, bell-shaped and trapezoidal MFs are utilized for FLC systems because of their simplicity in real-time execution and requirement of low-level memory. Thus, triangular MFs were utilized in this paper for the FLC as both the input and output. Five fuzzy linguistic variables that were considered for the FLC system were zero (Z), small negative (SN), big negative (BN), small positive (SP) and big positive (BP), as displayed in Figure 4 [32], and the rules are noted in Table 1. The center of gravity method of the defuzzification procedure was used in this study to calculate the fuzzy output. The output of the FLC system was then fed to the traditional PID to present the final output. However, the parameters of the traditional PID needed to be optimally rendered and were thus subjected to the time-domain-based performance index. Out of all the performance indices, ITAE is meant to be the most effective in reducing response overshoots and settling time to be implemented in this work, as given in Equation (12) [33]. center of gravity method of the defuzzification procedure was used in this study to calculate the fuzzy output. The output of the FLC system was then fed to the traditional PID to present the final output. However, the parameters of the traditional PID needed to be optimally rendered and were thus subjected to the time-domain-based performance index. Out of all the performance indices, ITAE is meant to be the most effective in reducing response overshoots and settling time to be implemented in this work, as given in Equation (12) [33].   center of gravity method of the defuzzification procedure was used in this study to calcu-late the fuzzy output. The output of the FLC system was then fed to the traditional PID to present the final output. However, the parameters of the traditional PID needed to be optimally rendered and were thus subjected to the time-domain-based performance index. Out of all the performance indices, ITAE is meant to be the most effective in reducing response overshoots and settling time to be implemented in this work, as given in Equation (12) [33].

Ultra Capacitors
UCs, usually called super capacitors, possess the tendency to store bulk power by charge separation between the metallic plates of a large surface area. Though the UCs are more costly compared to other ESDs, they have been extensively utilized with the benefits of a large power density and extravagant specific power. With a longer cycle time and lower and easy maintenance, they are suitable for peak demand situations. However, the UC is very clean and eco-friendly with zero emissions, which facilitates its extensive application. The transfer function of UCs utilized in this paper is provided in Equation (13).

Superconducting Magnetic Energy Storage
SMES devices are one of the most prominent ESDs, as their static nature leads to higher efficiency. SMES involves no chemical reactions, and thus no harmful gases are exposed to atmosphere. This device comprises a magnetic circuit enclosed in a closed chamber, which is maintained at a cryogenic temperature to attain super conductivity. SMES devices also comprise power converter unit and a step-down transformer. The magnetic coil charging and discharging can be achieved through the transformer, and the interface between the SMES and grid is facilitated by the power converter unit. The charging of coil in SMES will be done during off-peak load durations and discharging at peak demands. The compensator model of SMES utilized in this paper for stabilizing frequency is displayed in Figure 5 [34]. The SMES device parameters are stabilization gain K SMES , and T 1 , T 2 , T 3 and T 4 are time constants.

Superconducting Magnetic Energy Storage
SMES devices are one of the most prominent ESDs, as their static nature leads to higher efficiency. SMES involves no chemical reactions, and thus no harmful gases are exposed to atmosphere. This device comprises a magnetic circuit enclosed in a closed chamber, which is maintained at a cryogenic temperature to attain super conductivity. SMES devices also comprise power converter unit and a step-down transformer. The magnetic coil charging and discharging can be achieved through the transformer, and the interface between the SMES and grid is facilitated by the power converter unit. The charging of coil in SMES will be done during off-peak load durations and discharging at peak demands. The compensator model of SMES utilized in this paper for stabilizing frequency is displayed in Figure 5 [34]. The SMES device parameters are stabilization gain KSMES, and T1, T2, T3 and T4 are time constants.

Redoxflow Batteries
RFBs are one of the most effective ESDs based on the electrochemical reaction theory. Due to the benefits of a quick response nature and load-leveling characteristics, they have an extensive range of applications. Sulfuric acid solution with ions of vanadium acts as an

Redoxflow Batteries
RFBs are one of the most effective ESDs based on the electrochemical reaction theory. Due to the benefits of a quick response nature and load-leveling characteristics, they have an extensive range of applications. Sulfuric acid solution with ions of vanadium acts as an electrolytic solution in RFBs for chemical reactions. RFBs comprising of two isolation chambers eliminate the problem of self-discharge. The oxidation and reduction reactions of the electrolytic solution in RFBs are responsible for the charging and discharging process. Its operation range depends upon the quantity of the electrolyte solution employed. In this paper, the workings of RFBs for the LFC study are given in Equation (14) [35].

HAEFA Algorithm
HAEFA is one of the newest meta-heuristic searching strategies developed by the researchers in [25]. Based on Coulomb's law of electric field theory, the AEFA approach was implemented by the authors in [36]. This AEFA algorithm was implemented by the researchers especially for the control and operation of modern power system analysis and LFC studies and detected a drawback of updating the step size, which affects the effective utilization of search space. To overcome the above and to make the AEFA more effective in obtaining solutions for engineering optimization problems, the researchers in [25] put forward the hybridization mechanism and termed it the HAEFA algorithm. The implementation of this HAEFA approach is explained as follows: First, initially randomize the gains of the controller in both LFC-AVR loops of the considered power system model as The particles of each of the populations are updated in simulation. After performing the simulation for a given time period, the objective function values are calculated for each of the populations and treat them as local best solutions. The objective function conceived in this paper was ITAE, and the corresponding fitness function was evaluated as Among the local best solutions, select the population that results in the lowest objective function value with the highest fitness value as the global best population. The respective solution is considered to be the global best solution, as the formulated objective function handles frequency and power deviations pertaining to multiple areas, which increases the problem complexity. Due to this the diversity of the solutions is increased. In order to bring the local best solutions towards the global best solution a novel velocity equation is formulated. In this equation the chaotic constants C 1 and C 2 are generated dynamically to change the particles' velocity to favor the local and global best solutions. Using these local and global best solutions, the velocity of each of the populations can be calculated as Later, the particles position will be updated using where V (t+1) i and X (t+1) i represent velocity and the position of ith particle in the (t + 1) iteration, respectively. The chaotic parameters are varied dynamically to limit the chances of solution divergence using The acceleration (a) and electric field (E) of ith particle at the (t + 1) iteration will be calculated as shown below.
The total force exerting on ith particle in (t + 1) iteration as The force magnitude exerted on ith due to the jth particle in the (t + 1) iteration is where 'K' represents Coulomb's constant, which can be calculated as where K 0 and α are constant parameters. In order to maximize the force exerted by the electric field on the particles, half of the equivalent Euclidian distance is considered, as follows.
The charge of the ith and jth particles as indicated with Q (t+1) i and Q (t+1) j at the (t + 1) iteration is as follows.
The problem formulated in this paper is a minimization one, and hence the best and worst fitness values will be evaluated as This procedure is repeated for three iterations to train the populations to increase the searching capability. This is considered to be the first stage of the proposed algorithm. Here, after starting the iterative process, it takes a minimum of three iterations to obtain the decreased variance of deviation of the solutions. Using this statistical analysis, the developed algorithm generated first-stage solutions after three iterations. After this, the pair-wise comparison process was initiated to divide the total populations in half. The populations related to these solutions were forwarded to the second stage of the proposed methodology. Due to this, the remaining iterative process was performed with maximum efficacy to obtain the best solution within fewer iterations. The process of pair-wise comparison is as follows below: Let f itness 1 , f itness 2 , f itness 3 . . . . . . . . . .. f itness n be the fitness values of the solution for 'n' particles; If f itness 1 > f itness 2 , then send f itness 1 to the next stage; If f itness 3 > f itness 4 , then send f itness 3 to the next stage; If f itness 5 > f itness 6 , then send f itness 5 to the next stage; If f itness (n−1) > f itness (n) , then send f itness (n−1) to the next stage. This process of pair-wise comparison was performed on all populations to obtain the best half of the entire population. The remaining algorithmic operations were performed on these populations.
After reviewing the literature carefully, it was noticed that the hybridization of multiple algorithms increases the algorithm's performance to the maximum extent. From this, it was determined that in this paper a new hybridization approach was established by taking the advantage of crossover operation. As the considered objective function makes the problem more complex, it was necessary to obtain a converged solution instead of a diverged solution. From the fourth iteration onwards the existing crossover operation was performed on second-stage populations to minimize the solution divergence with updated velocity. Crossover can also be called recombination, which works to create offspring from two parents by combining the genetic information in a specified manner. It is the only way of generating new offspring stochastically from an existing population, which is analogous to the biological human reproduction system. There are various types of crossover operations available in the literature; however, after the conclusion of an extensive literature survey, it was determined that uniform crossover operation may yield better results, and it was thus implemented. (30) where λ is the random number between [0, 1] and Table 2, given below, indicates the uniform crossover operation.  After crossover operation the positions are updated; this whole procedure will be repeated till the completion of the maximum number of iterations, and then the best values are displayed. HAEFA flowchart is displayed in Figure 6. After crossover operation the positions are updated; this whole procedure will be repeated till the completion of the maximum number of iterations, and then the best values are displayed. HAEFA flowchart is displayed in Figure 6. Before implementing this strategy for power system optimization, the efficacy of hybridization was tested on the benchmark function of sphere function given in Equation (31). Both the algorithms were implemented with a population of 100 for 100 trials. The efficacy of the algorithm was assessed in terms of initial and final function values, as displayed in Figure 7. From Figure 7, it can be concluded that the function values at the initial and final instants were good, but a better solution was found using HAEFA as compared to the AEFA approach. Moreover, the final values of the function with HAEFA were almost below the mean value, which shows its effectiveness and tendency to maintain exploration and exploitation. Before implementing this strategy for power system optimization, the efficacy of hybridization was tested on the benchmark function of sphere function given in Equation (31). Both the algorithms were implemented with a population of 100 for 100 trials. The efficacy of the algorithm was assessed in terms of initial and final function values, as displayed in Figure 7. From Figure 7, it can be concluded that the function values at the initial and final instants were good, but a better solution was found using HAEFA as compared to the AEFA approach. Moreover, the final values of the function with HAEFA were almost below the mean value, which shows its effectiveness and tendency to maintain exploration and exploitation.

Simulation Results
The simulation of the combined LFC and AVR system of MAIPS was designed in the domain of MATLAB/Simulink version (2016a), and the HAEFA optimization algorithm was developed in a (.m file) format. In order to assess the dynamical analysis, the combined LFC and AVR system was targeted with 10%SLP on area-1. The HAEFA algorithm was implemented with a population size of 50 for 100 iterations.

Case-1: Analysis of Combined LFC and AVR System with HAEFA-Tuned Controllers
The combined LFC and AVR model of MAIPS were analyzed with different controllers such as traditional PI/PID and intelligent-based fuzzy PID in both the loops of LFC and AVR one at a time by injecting area-1 with a load of 10% SLP. The parameters of these controllers were rendered optimally with a hybridized algorithm of the HAEFA approach. Figure 8 depicts the responses of the combined LFC-AVR system under various controllers fine-tuned with the HAEFA approach, and these responses were interpreted in terms of settling time, as provided in Table 3. After observing the responses rendered in Figure 8, we came to the decision that the undershooting and overshooting of the responses were very comfortably mitigated by the fuzzy PID compared to traditional PI/PID controllers. Moreover, the terminal voltages quickly reached specified values with the HAEFA-tuned fuzzy PID regulator. From Table 3, it is obvious that the responses reached the steady condition in much less time with the fuzzy PID controller than traditional PI/PID. Thus, the intelligent fuzzy PID showed its predominance in regulating the behavior of complex MAIPS of the LFC-AVR combined model. The parameters of PI/PID/fuzzy PID controllers that were retrieved optimally with the HAEFA technique are noted in Table 4. Further, the ITAE index value of the HAEFA-tuned fuzzy PID controller was improved by 54.19% with the HAEFA-based PI and 41.66% with HAEFA-based PID controllers.

Simulation Results
The simulation of the combined LFC and AVR system of MAIPS was designed in the domain of MATLAB/Simulink version (2016a), and the HAEFA optimization algorithm was developed in a (.m file) format. In order to assess the dynamical analysis, the combined LFC and AVR system was targeted with 10%SLP on area-1. The HAEFA algorithm was implemented with a population size of 50 for 100 iterations.

Case-1: Analysis of Combined LFC and AVR System with HAEFA-Tuned Controllers
The combined LFC and AVR model of MAIPS were analyzed with different controllers such as traditional PI/PID and intelligent-based fuzzy PID in both the loops of LFC and AVR one at a time by injecting area-1 with a load of 10% SLP. The parameters of these controllers were rendered optimally with a hybridized algorithm of the HAEFA approach. Figure 8 depicts the responses of the combined LFC-AVR system under various controllers fine-tuned with the HAEFA approach, and these responses were interpreted in terms of settling time, as provided in Table 3. After observing the responses rendered in Figure 8, we came to the decision that the undershooting and overshooting of the responses were very comfortably mitigated by the fuzzy PID compared to traditional PI/PID controllers. Moreover, the terminal voltages quickly reached specified values with the HAEFA-tuned fuzzy PID regulator. From Table 3, it is obvious that the responses reached the steady condition in much less time with the fuzzy PID controller than traditional PI/PID. Thus, the intelligent fuzzy PID showed its predominance in regulating the behavior of complex MAIPS of the LFC-AVR combined model. The parameters of PI/PID/fuzzy PID controllers that were retrieved optimally with the HAEFA technique are noted in Table 4. Further, the ITAE index value of the HAEFA-tuned fuzzy PID controller was improved by 54.19% with the HAEFA-based PI and 41.66% with HAEFA-based PID controllers.    In this sub section, a comparison has been made between the responses of the considered MAIPS without and with considering AVR coupling under the supervision of the HAEFA-tuned fuzzy PID controller. In order to obtain the analysis in more comparative approach, the power system with and without AVR coupling was analyzed for the same load disturbance, 10%SLP on area-1, and the responses are displayed in Figure 9. From the responses rendered in Figure 9, it can be noticed that the deviations in system behavior were slightly heightened with the consideration of AVR coupling, which can be justified by Equation (8). The stability of the power system depends on both terminal voltage and area frequency. The monitoring and regulation of voltage and frequency of the power system must be carried out with intense care. The mathematical analysis given in Equation (8) demonstrates the cumulative and simultaneous impact of variations in the AVR loop leading to the fluctuations in the LFC loop and vice versa. Considering all these analyses the ITAE objective function provided in Equation (12) is formulated for simultaneous voltage and frequency stabilization.

Case-3: Analysis of Combined LFC and AVR System with Different ESDs
The fuzzy PID regulator based on the HAEFA mechanism showed the supreme performance; as revealed in aforementioned cases, further analysis was performed under its supervision. The design of the secondary regulator alone was not sufficient to maintain the stability of the large complex interconnected power system networks. Thus, territorial control strategy was needed. Different ESDs such as RFBs, UCs and SMES were placed in both the areas of the LFC loop one after the other for the same load disturbances. The responses of the combined LFC-AVR system with different ESDs are compared in Figure 10 and were assessed numerically from a settling time point of view, as noted in Table 5, and they are further compared in the bar chart displayed in Figure 11 for easy observation. Analyzing the results depicted in Figure 11, we came to know that the deviations in frequency and line power flow were mitigated further, and the respective area terminal voltages were enhanced to a great extent. We also observed that for the different ESDs, such as RFBs, SMES and UCs, the RFBs showed more predominance in bringing the deviations in system responses to a steady condition more effectively. This was possible only because the quick response nature and lower cycle time of RFBs facilitates the complex combined LFC-AVR system to prevent hunting during load fluctuations compared to other ESDs.  (8) demonstrates the cumulative and simultaneous impact of variations in the AVR loop leading to the fluctuations in the LFC loop and vice versa. Considering all these analyses the ITAE objective function provided in Equation (12) is formulated for simultaneous voltage and frequency stabilization.

Case-3: Analysis of Combined LFC and AVR System with Different ESDs
The fuzzy PID regulator based on the HAEFA mechanism showed the supreme performance; as revealed in aforementioned cases, further analysis was performed under its supervision. The design of the secondary regulator alone was not sufficient to maintain

Conclusions
The stability of frequency and voltage of the interconnected power system can be addressed simultaneously by connecting the AVR loop with the LFC using coupling coefficients. In this paper a HAEFA-tuned fuzzy PID was employed as a secondary regulator in both the LFC and AVR loops, and the dynamical analysis was carried out by laying area-1 with 10% SLP. The responses of LFC and AVR showcase the supremacy of fuzzy PID over classical PI and PID. Moreover, with fuzzy PID the ITAE function value was

Conclusions
The stability of frequency and voltage of the interconnected power system can be addressed simultaneously by connecting the AVR loop with the LFC using coupling coefficients. In this paper a HAEFA-tuned fuzzy PID was employed as a secondary regulator

Conclusions
The stability of frequency and voltage of the interconnected power system can be addressed simultaneously by connecting the AVR loop with the LFC using coupling coefficients. In this paper a HAEFA-tuned fuzzy PID was employed as a secondary regulator in both the LFC and AVR loops, and the dynamical analysis was carried out by laying area-1 with 10% SLP. The responses of LFC and AVR showcase the supremacy of fuzzy PID over classical PI and PID. Moreover, with fuzzy PID the ITAE function value was enhanced by 54.19% with PI and 41.67% with PID. However, the necessity of considering AVR coupling with the LFC loop was demonstrated and justified. Further, the combined system was incorporated with different ESDs such as RFBs, SMES and UCs one after the other in both the areas, and the system dynamical responses were compared to reveal the best one. The comparative results reveal that there as an enhancement in system dynamical behavior when incorporating ESDs in the view of the reduction in peak undershoot, overshoot and settling time. Moreover, out of all the ESDs the RFBs showed more predominance in the enhancement of system deviations and is recommended as the preferred territorial controller for the stability of interconnected power systems over the other ESDs.