Power System Stability Improvement of FACTS Controller and PSS Design: A Time-Delay Approach

: The existence of low-frequency oscillations in power systems is the cause of power angle instability, limiting the transmission of maximum tie-line power. One of the effective ways to improve the stability limits is by installing a power system stabilizer and supplementary excitation control to augment with an automatic voltage regulator (AVR) supplemental feedback stabilizing signal. This paper proposes a new strategy for simultaneously tuning the power system stabilizer (PSS) and FACTS controller, considering time delays. The design of the proposed controller is modeled as an optimization problem, and the parameters of the controller are optimized through the grasshopper optimization algorithm (GOA). The suggested controller’s efﬁcacy is evaluated for both single-machine inﬁnite bus systems and multi-machine power systems under various disturbances. It also investigated the performance of the proposed controller with variations in signal transmission delays. The results obtained from GOA optimized proposed controller are compared with those obtained from the differential evolution algorithm, genetic algorithm, and whale optimization algorithm. In this context, the proposed GOA optimized controller reduced the objective function value by 16.32%, 14.56%, and 13.72%, respectively, in the SMIB system and 1.41%, 9.98%, and 13.31%, respectively, for the multi-machine system compared with the recently published WOA, and the well-established GA and DE. Further, the proposed controller is found to be stable and effectively increases stability even under small disturbances.


Introduction
A sustainable power system accommodates higher penetration of renewable generations via effective planning and operational approaches to ensure system reliability and stability generally undergoes structural and technological transformations [1][2][3][4][5][6]. System stability, being one of the essential aspects [7][8][9][10][11][12], numerous stability issues may be caused by heavily loaded long tie-lines in an interconnected power system. There have been different approaches to mitigating negative impacts of inreased renewable penetration into the grid [13][14][15][16][17]. On this note, the design of a suitable power system stabilizer (PSS) is currently a research interest. With the increasing usage of FACTS devices for the enhancement of power system oscillations damping, proper synergism of FACTS controller and PSS is essential. As the power system is non-linear in nature, the control parameters of FACTS devices and PSS are estimated, taking the system's non-linearity into consideration [18,19]. When the FACTS controller and PSS are utilized simultaneously in the power system, investigating the coordination strategy becomes as essential as presenting an appropriate function to prepare this procedure. Generally, it can prove to be an optimization process to estimate the values of their controller parameters at the same time.
When a damping controller based on a PSS and a FACTS device is designed and operated in a coordinated manner, it leads to the enhancement of power system performance. This coordinated process works well and is an important topic of research and study these days [20][21][22][23][24]. Installation of the PSS at a location close to the generator and that of a FACTS device away from the generator have to be taken into consideration. Therefore, the transmission delay that is caused by each sensor and signal transmission in the power system should be considered during the design of the damping controller based on FACTS-PSS [25]. Moreover, different algorithms like the linear matrix inequality technique, residue method, and multi-model adaptive approach for control have been applied to enhance the damping of the coordinated controller. The verification of the robustness of the design approach in the coordinated technique is the main topic of discussion, as the high efficiency and effectiveness of the techniques ensure a robust controller design that varies with the system configuration.
In recent times, a lot of research is going on, which are heuristic in nature that uses likenesses of a social system [26]. Because of their ability to identify optimal solutions to multi-modal, non-differentiable, and complex objective functions, these techniques can be promising problem solvers when applied in the research community.
In recent times, various algorithms have been used for modeling damping controllers that are based on PSS and FACTS. From the present literature analysis, it is perceived that the combined coordinated structure of PSS and SSSC have significant relevance for the leadlag controller. Nowadays, various intelligent methods have been extensively used in this coordinated design. PSS has been designed in combination with other FACTS controllers in the literature, such as PSS and SSSC by GA [27], and PSS and SVC by improved teachinglearning techniques [28]. The author in [29] discussed the coordinated control of two FACTS controllers, such as TCSC and SSSC, using various disturbances. The author investigates the time delay method for the coordinated SSSC and PSS controller using a hybrid PSO-GSA technique [30]. The author analyses the coordinated structure of SSSC and PSS using the ACO technique [31]. The author investigates the coordinated SSSC and PSS design using a hybrid technique [32]. In [33,34], the author discussed simultaneous tuning of PSS and SSSC controllers using modified SCA. Similarly, in [35], PSO is used to tune SSSC and PSS controllers. In [36], authors suggested a coordinated structure of PSS and SSSC controllers using modified WOA. The author investigated the PSS and STATCOM using the learning BAT algorithm in [37]. In [38] the PSS and SVC parameters were tuned by using the BFPSO algorithm. The coordinated STATCOM and PSS design using SOA discussed in [39]. The MWOA-NM algorithm was proposed in [40] to co-ordinately tune the PSS and SSSC parameters. Furthermore, in [41], a PSO with an exponential delay algorithm is proposed to co-ordinately tunning AVR and PSS parameters. In [42], the author discussed the simultaneous tuning of SSSC-POD with PSS using the mayfly optimization technique. The various intelligent techniques are extensively used by several researchers in modeling the coordinated structure of FACTS controller together with PSS.
In [43], Saremi et al. proposed a grasshopper optimization algorithm (GOA) based on the social activity of grasshoppers. The GOA has gained popularity over the years because of its simple structure, gradient-free mechanism, and convergence to an optimum global value instead of an optimum local value. The GOA processed the optimization problem in two stages: exploration and exploitation. This technique randomly varies the solutions in the exploration stage, leading to a more extensive search space and discovering its promising areas. After exploration, the GOA exploits the search area to get the optimum solution. However, an improper balance between exploration and exploitation leads to converging the solution to a local one instead of tracking for the global solution. The obtained results for the optimization problem with GOA are compared with PSO, FPA, BA, FA, SMS, and

System Modeling
It is necessary to do the evaluation of the performance of any damping controller based on PSS and SSSC along with their proper design. Therefore, a SMIB system is considered, as shown in Figure 1, comprising a synchronous generator that is linked to an infinite bus with the help of a transformer and an SSSC. Moreover, the hydraulic turbine and governor (HTG) equipped in the generator signifies a nonlinear hydraulic turbine model, a servomotor, and a PID governor system. Apart from this, some other equipment existing within the generator includes the excitation system consisting of a DC exciter along with a voltage regulator, in the absence of saturation function of the exciter and a PSS.
researchers are more interested in simultaneously tuning the FACTS controller and the PSS when it comes to improving stability. In this present work, the following objectives are accomplished.

•
This paper investigates the simultaneous tuning of the PSS and FACTS controller considering the potential time delays. Its goal is to improve power system stability in the presence of potential time delay. • GOA is utilized to fine-tune control parameters for the proposed fuzzy lead-lag controller design . To illustrate the robustness of the proposed design methodology,  simulation results for both single-machine infinite-bus and multi-machine power  systems are presented under different disturbances and

System Modeling
It is necessary to do the evaluation of the performance of any damping controller based on PSS and SSSC along with their proper design. Therefore, a SMIB system is considered, as shown in Figure 1, comprising a synchronous generator that is linked to an infinite bus with the help of a transformer and an SSSC. Moreover, the hydraulic turbine and governor (HTG) equipped in the generator signifies a nonlinear hydraulic turbine model, a servomotor, and a PID governor system. Apart from this, some other equipment existing within the generator includes the excitation system consisting of a DC exciter along with a voltage regulator, in the absence of saturation function of the exciter and a PSS.
In Figure 1, the transformer is represented by T; the infinite bus and generator terminal voltages are represented by VB and VT, respectively; V1 and V2 are the bus voltages, the output voltage of the SSSC converter and the DC voltage source is represented by Vcnv and VDC, respectively; the line current is I and the total flow of real power in the transmission lines and that of a single line are each represented by PL and PL1, respectively.  In Figure 1, the transformer is represented by T; the infinite bus and generator terminal voltages are represented by V B and V T , respectively; V 1 and V 2 are the bus voltages, the output voltage of the SSSC converter and the DC voltage source is represented by V cnv and V DC , respectively; the line current is I and the total flow of real power in the transmission lines and that of a single line are each represented by P L and P L1 , respectively.

Model of SSSC
An AC voltage can be generated and controlled with the help of a solid-state voltagesourced converter (VSC) called SSSC. The line impedance is virtually compensated by the controllable voltage (Vq) injection by SSSC, which is connected in series with the transmission lines of the power system. The quadrature component of the line current Vq influences the flow of power in the transmission lines independent of the line current magnitude by reproducing the reactance of either an inductor or a capacitor. Moreover, a VSC that is present on the transformer's secondary side is used to alter Vq. Additionally, Sustainability 2022, 14, 14649 4 of 22 a capacitor linked on the direct current side of the VSC usually acts as a DC voltage source. This DC voltage source then provides AC voltage using forced commutated power electronics devices. The dynamic control of polarity and magnitude of Vq determines the extent to which the transmission lines are compensated as the device operates in both capacitive and inductive modes [44]. However, charging of capacitors and losses occurring in the transformer and VSC draw some amount of active power from the line. The PWM inverters, which are based on IGBT, have been utilized in the current study [45][46][47].

Proposed Fuzzy Structured FACTS Controller Structure with PSS
The ability of a fuzzy logic controller (FLC) to boost the performance of the lead-lag controller has been recognized by many authors. It deals with structure parameters and updates the online parameters of the controller [48]. Membership function (MF) tuning, as well as rule tuning, being the two types of fuzzy control tuning methods, play a significant role in fuzzy control performance enhancement. The strong rule base and standard MFs techniques are used by focusing on the applications. Scaling factors are to be tuned in ideal fuzzy lead-lag control [49]. In real-time applications, designing the controller is done with the help of triangular MFs, owing to the cost-effectiveness of its practical parametric illustration in comparison with others. In this paper, we selected triangular MFs to design the controller. On the other hand, a steady characterization of MFs is normally considered to make efficient use of memory, computational efficiency, and necessity of performance analysis [50]. As a result, for obtaining the output of the FLC, MFs' error derivative is chosen. The MFs are based on five fuzzy linguistic factors: positive big, negative big, positive small, negative small, and zero. The MFs considering the output of FLC, error derivative, and error is shown in Figure 2.
An AC voltage can be generated and controlled with the help of a solid-state voltagesourced converter (VSC) called SSSC. The line impedance is virtually compensated by the controllable voltage (Vq) injection by SSSC, which is connected in series with the transmission lines of the power system. The quadrature component of the line current Vq influences the flow of power in the transmission lines independent of the line current magnitude by reproducing the reactance of either an inductor or a capacitor. Moreover, a VSC that is present on the transformer's secondary side is used to alter Vq. Additionally, a capacitor linked on the direct current side of the VSC usually acts as a DC voltage source. This DC voltage source then provides AC voltage using forced commutated power electronics devices. The dynamic control of polarity and magnitude of Vq determines the extent to which the transmission lines are compensated as the device operates in both capacitive and inductive modes [44]. However, charging of capacitors and losses occurring in the transformer and VSC draw some amount of active power from the line. The PWM inverters, which are based on IGBT, have been utilized in the current study [45][46][47].

Proposed Fuzzy Structured FACTS Controller Structure with PSS
The ability of a fuzzy logic controller (FLC) to boost the performance of the lead-lag controller has been recognized by many authors. It deals with structure parameters and updates the online parameters of the controller [48]. Membership function (MF) tuning, as well as rule tuning, being the two types of fuzzy control tuning methods, play a significant role in fuzzy control performance enhancement. The strong rule base and standard MFs techniques are used by focusing on the applications. Scaling factors are to be tuned in ideal fuzzy lead-lag control [49]. In real-time applications, designing the controller is done with the help of triangular MFs, owing to the cost-effectiveness of its practical parametric illustration in comparison with others. In this paper, we selected triangular MFs to design the controller. On the other hand, a steady characterization of MFs is normally considered to make efficient use of memory, computational efficiency, and necessity of performance analysis [50]. As a result, for obtaining the output of the FLC, MFs' error derivative is chosen. The MFs are based on five fuzzy linguistic factors: positive big, negative big, positive small, negative small, and zero. The MFs considering the output of FLC, error derivative, and error is shown in Figure 2. The SSSC injected voltage (Vq) is governed by the damping controller, which is based on the PSS and SSSC. This Vq is used for damping and hence can be considered the controller output while inputting the speed deviation (∆ω). Here the proposed controller is comprised of two lead-lag components, as detailed in Figure 3. The structure comprises a gain block with gain KPS, a signal washout block serving as a high-pass filter, and a phase compensation block to provide appropriate phase-lead characteristics between input and output signals and which are two-staged, as shown in the figure. The output of the PSS structure shown in Figure 4, Vs cumulate with the reference voltage of the excitation The SSSC injected voltage (Vq) is governed by the damping controller, which is based on the PSS and SSSC. This Vq is used for damping and hence can be considered the controller output while inputting the speed deviation (∆ω). Here the proposed controller is comprised of two lead-lag components, as detailed in Figure 3. The structure comprises a gain block with gain K PS , a signal washout block serving as a high-pass filter, and a phase compensation block to provide appropriate phase-lead characteristics between input and output signals and which are two-staged, as shown in the figure. The output of the PSS structure shown in Figure 4, Vs cumulate with the reference voltage of the excitation system, V ref . From the washout function's perspective, the time constant value T WS = T WP may lie somewhere between 1-20 s and is not critical. In the present analysis, K PS and K PP are the controller gains, k 1s , k 2s , k 3s and k 1p , k 2p , k 3p are the scaling factors and time constants (T 1s , T 2s , T 3S, and T 4S ) to be computed.  In the present study, time delay means communication delay that is inherently present in any control system. A communication delay in a power grid is defined as the time between sending a signal from the source device to the receiving of the signal at the destination device. Typically, SSSC devices are installed at locations that are far away from the generating units. The speed deviation of the generator is taken as an input signal to the proposed fuzzy lead-lag structured SSSC controller, as shown in Figure 3. Transmitting the speed deviation signal from the generator to the SSSC controller will take some time; therefore, the time delay is included for the SSSC controller, as shown in Figure 3. Owing to the advancement in optical fiber communication, synchronous phasor measurement can be done by a wide-ranging measuring system that can then be transmitted simultaneously to control centers. Thus, considering generator angle and speed deviation as remote signals, they can be put to use for designing purposes. However, the delay in the transmission network poses a major problem while using these types of signals as the delay results in significant degradation of that particular transmission line performance. Even in the worst scenarios, a dedicated channel of communication should avoid exceeding a delay of 50 ms during signal communication [45]. So, the remedial solution is to take care of the time delays of this order during the designing of the controller, as it is done in most cases these days. If improvement of the power system stability is taken into consideration, the remote speed deviation signal is preferred over the local tie-line power and line current signals for damping controllers  In the present study, time delay means communication delay that is inherently present in any control system. A communication delay in a power grid is defined as the time between sending a signal from the source device to the receiving of the signal at the destination device. Typically, SSSC devices are installed at locations that are far away from the generating units. The speed deviation of the generator is taken as an input signal to the proposed fuzzy lead-lag structured SSSC controller, as shown in Figure 3. Transmitting the speed deviation signal from the generator to the SSSC controller will take some time; therefore, the time delay is included for the SSSC controller, as shown in Figure 3. Owing to the advancement in optical fiber communication, synchronous phasor measurement can be done by a wide-ranging measuring system that can then be transmitted simultaneously to control centers. Thus, considering generator angle and speed deviation as remote signals, they can be put to use for designing purposes. However, the delay in the transmission network poses a major problem while using these types of signals as the delay results in significant degradation of that particular transmission line performance. Even in the worst scenarios, a dedicated channel of communication should avoid exceeding a delay of 50 ms during signal communication [45]. So, the remedial solution is to take care of the time delays of this order during the designing of the controller, as it is done in most cases these days. If improvement of the power system stability is taken into consideration, the remote speed deviation signal is preferred over the local tie-line power and line current signals for damping controllers In the present study, time delay means communication delay that is inherently present in any control system. A communication delay in a power grid is defined as the time between sending a signal from the source device to the receiving of the signal at the destination device. Typically, SSSC devices are installed at locations that are far away from the generating units. The speed deviation of the generator is taken as an input signal to the proposed fuzzy lead-lag structured SSSC controller, as shown in Figure 3. Transmitting the speed deviation signal from the generator to the SSSC controller will take some time; therefore, the time delay is included for the SSSC controller, as shown in Figure 3. Owing to the advancement in optical fiber communication, synchronous phasor measurement can be done by a wide-ranging measuring system that can then be transmitted simultaneously to control centers. Thus, considering generator angle and speed deviation as remote signals, they can be put to use for designing purposes. However, the delay in the transmission network poses a major problem while using these types of signals as the delay results in significant degradation of that particular transmission line performance. Even in the worst scenarios, a dedicated channel of communication should avoid exceeding a delay of 50 ms during signal communication [45]. So, the remedial solution is to take care of the time delays of this order during the designing of the controller, as it is done in most cases these days. If improvement of the power system stability is taken into consideration, the remote speed deviation signal is preferred over the local tie-line power and line current signals for damping controllers based on FACTS. Thus, a 15 ms sensor time constant and 50 ms time delayed signal transmission are each considered for PSSs and SSSC-based controllers, respectively.

Optimization Problem
T WS =T WP = 10 s is used in the present work and is pre-specified [20]. The objective is to determine the time constants along with gains associated with the controller. As ∆V q is zero during steady-state conditions, V qref becomes constant. However, during dynamic conditions, the modulation of the injected series voltage V q is done by applying a certain algorithm in order to damp out the system oscillations. V qref can be assumed to be constant due to the slow operation of the power flow loop during the steady state. Thus, under dynamic conditions, the effective value of V q can be formulated as With the help of any one of the remote signals like speed-deviation, tie-line power, or the deviation in the power angle, the system oscillation can be observed, and minimizing any one of such deviations could be the objectives of researchers. In the proposed work, the two objective functions that are considered are an integral time absolute error of speed deviation for the SMIB power system and that of speed signals corresponding to the inter-area and local oscillation modes. Both the objective functions are expressed below: For the SMIB power system: For power systems consisting of multiple machines: where, t = time range of simulation.
The time domain simulation of the above power system model is performed for a specified time period for the purpose of the objective function. The range of the damping controller and PSS lie within the prescribed bounds. The aim of improving the system response in terms of settling time and overshoots is achieved by minimization of the objective function. Therefore, the optimization problem is expressed from the discussed design approach as follows: It is to be noted that two gains, eight-time constant parameters, and six scaling factors are needed to be optimized for a SMIB system consisting of one damping controller and one PSS. Similarly, for a power system composed of multiple machines such as multiple PSSs and damping controllers, which is equivalent to the no. of generators, optimization is carried out on all the parameters.

GOA Method
It has been demonstrated that the Grasshopper Optimization Algorithm (GOA) [43] benefits from high exploration while exhibiting very fast convergence speed. Exploration and exploitation balance smoothly with the unique adaptive mechanism in this algorithm. These features possibly enable the GOA algorithm to deal with greater exploration of search space and outperform other methods. Furthermore, computational complexity is superior to that of many techniques, which is displayed in the literature. These powerful characteristics of the GOA technique motivated us to propose this present study.

Overview of GOA
Grasshoppers are an insect. Because of their harm to agriculture and crop production, they are treated as a pest. A real grasshopper is shown in Figure 5. The life cycle of grasshoppers appears in Figure 6. Grasshoppers show some substantial behaviors during the food search phase. Typically, grasshoppers engage themselves in groups and can form some of the largest swarms in the animal kingdom. The discovery of swarming activities in the same nymph and adult period characterizes the grasshopper group. As mentioned above, the nature-inspired methods reasonably split the search space into two tendencies: exploration and exploitation. These two factors, as well as the seeking of food, are carried out by grasshoppers. Therefore, a mathematical model of grasshopper behavior can be designed by considering the above-discussed aspects.
they are treated as a pest. A real grasshopper is shown in Figure 5. The life cycle of grasshoppers appears in Figure 6. Grasshoppers show some substantial behaviors during the food search phase. Typically, grasshoppers engage themselves in groups and can form some of the largest swarms in the animal kingdom. The discovery of swarming activities in the same nymph and adult period characterizes the grasshopper group. As mentioned above, the nature-inspired methods reasonably split the search space into two tendencies: exploration and exploitation. These two factors, as well as the seeking of food, are carried out by grasshoppers. Therefore, a mathematical model of grasshopper behavior can be designed by considering the above-discussed aspects.  they are treated as a pest. A real grasshopper is shown in Figure 5. The life cycle of grasshoppers appears in Figure 6. Grasshoppers show some substantial behaviors during the food search phase. Typically, grasshoppers engage themselves in groups and can form some of the largest swarms in the animal kingdom. The discovery of swarming activities in the same nymph and adult period characterizes the grasshopper group. As mentioned above, the nature-inspired methods reasonably split the search space into two tendencies: exploration and exploitation. These two factors, as well as the seeking of food, are carried out by grasshoppers. Therefore, a mathematical model of grasshopper behavior can be designed by considering the above-discussed aspects.  The mathematical model designed in this paper for the swarming behavior of grasshoppers is presented below: where X i symbolizes the i-th grasshopper's position, social interaction denotes S i , F i is the nth grasshopper's gravitational force, and W i denotes wind speed factor. The above circumstances can be constructed as X i = r 1 I i + r 2 F i + r 3 W i . The above-mentioned r 1 , r 2 , and r 3 are randomized values selected between [0, 1].
where d kl signifies the distance of k-th from l-th grasshopper and may be calculated as d kl = |X l − X k | and s is the force that describes the potency of social interaction forces and is a unit vector from k-th grasshopper to l-th grasshopper. This is known as social interaction forces (s function), and is calculated as: where a signifies intensity of attraction and l indicates attractive length scale. This function s describes the effect on grasshopper communal interaction gathered from [39]. We consider the distance of function s between 0 and 15 as in [43]. It is worth noting that the repulsion occurs when the distance lies between 0 and 2.079. No attraction or repulsion occurs when the distance between two grasshoppers is 2.079. It is called the comfort zone. According to [43], the attraction steadily increases from 2.079 to 4 distance units before slowly decreasing beyond 4 distance units. The change in parameter values l and a in Equation (8) leads to various social behaviors of artificial grasshoppers. Figure 7 depicts a pictorial representation of the comfort zone and social interaction of grasshoppers. noting that the repulsion occurs when the distance lies between 0 and 2.079. No attraction or repulsion occurs when the distance between two grasshoppers is 2.079. It is called the comfort zone. According to [43], the attraction steadily increases from 2.079 to 4 distance units before slowly decreasing beyond 4 distance units. The change in parameter values l and a in Equation (8) leads to various social behaviors of artificial grasshoppers. Figure 7 depicts a pictorial representation of the comfort zone and social interaction of grasshoppers.  The function 's' segregates the space between two grasshoppers into attraction, repulsion, and comfort zones. However, this function returns zero value if the separation between two grasshoppers is more than 10. Hence, with larger separation, this function fails to apply forces between two grasshoppers. Therefore, to address this issue, this paper mapped the separation between the grasshoppers to a value between 1 to 4. In Equation (6), the F component is calculated as: where the gravitational constant is denoted by g and the unity vector directed towards the earth's center is represented by > e g . The component W in Equation (6) is estimated as where the drift constant is represented by u and the unity vector in the wind direction is denoted by > e w . In the nymph stage, a grasshopper does not have wings. Hence, its movement is predominantly influenced by wind speed. By replacing Equation (7), (9), and (10) in Equation (6), solving the grasshopper swarming behavior can be represented as follows: where N is the number of grasshoppers. As nymphs primarily stay on the ground, their location should lie below the threshold value. However, we cannot use this equation in this optimization algorithm because it avoids the search space condition such as exploration and exploitation. Nevertheless, this model cannot be utilized specifically to explain optimization issues, mostly because the grasshoppers quickly attain the comfort zone and the swarm does not converge to a predetermined point. To resolve the optimization problems, a modified form of the above equation is considered as follows.
where ub d and lb d represent upper and lower bounds, respectively, in the dimension D-th dimension, The target value denotes > T d and c is a decreasing constant that reduces the attraction, repulsion, and comfort zones. In Equation (12), the term S represents social interaction, as in Equation (6). Gravity is ignored in Equation (11) and the wind direction is assumed to always be towards the target.
The first term of Equation (12) gives information about the position of the other grasshopper and the interaction between them in nature. The second term counts the tendency of them to move toward the food. In Equation (12), the parameter c represents the grasshopper's deceleration as it approaches and takes food. To increase the applicability of the optimization technique random behavior of the grasshopper is considered. To incorporate randomness into grasshopper interaction, each term in Equation (12) is multiplied by random variables. Exploiting and exploring the search space requires mathematical formulas and tuning the level of exploitation and exploration requires a search agent. In the nymph stage, the grasshopper search for food in the local search space due to the absence of wings. However, in the adult stage, they seek a larger search space region. In exploration, the grasshopper discovers a favorable place with a higher chance of finding food. Exploration arises first in this strategy, as it does in stochastic optimization because the first attention is on locating promising search space. Following exploration, exploitation is used to locate food in potential local locations in order to obtain the best global value. As the number of iterations increases, c has a significant impact on the grasshopper's behavior. Outside c reduces the search space around the target maximum, whereas inner c reduces grasshopper attraction or repulsion forces.
The value of c should decrease relative to the number of iterations to maintain the proper balance between exploration and exploitation. This phenomenon encourages the algorithm towards exploitation as the number of iterations increases. The decreasing coefficient c diminishes the size of the comfort zone in proportion to iteration counts, as in Equation (13).
where c m and c n are the maximum and minimum values, the current iteration is denoted by I, while the maximum number of iterations is denoted by L. The value of c m and c n in this technique are 1 and 0.00001, respectively. The proposed flowchart of the GOA technique is represented in Figure 8.

Results and Discussion
For the calculations and design of the damping controller, a toolbox named Sim Power Systems (SPS) was used extensively. Engineers may easily create and build models in the SIMULINK environment with the MATLAB-based design tool SPS. Models of typical power equipment such as transformers, transmission lines, power electronics, and machines, are all present in its libraries. Load flow and initialization of the three-phase machines can be performed with the help of the 'Powergui' block. Thus, the analysis of the developed models can be done by the graphical user interface (GUI) tools that are provided by it. Figure 9 depicts the SMIB power system model designed in MATLAB consisting of a generating unit in connection with a double-circuited parallel line transmission line, generally. The above connection is established by connecting a step-up transformer and an SSSC in between. For tuning the proposed controller, the GOA algorithm is employed using an m file. Simulation of the entire model is done along with simultaneous calculation of the objective function owing to the occurrence of any disturbance in the system using the MATLAB R2016a environment. Further, the controller parameters can be found by minimization of the fitness values of Equation (2), whose optimization is carried out using the GOA algorithm.  Table 2. In this proposed simulation work, different parameters have been initialized for the application of GOA, GA, DE, and WOA. Maximum no. of generation and population size is the common control parameters of the algorithms. In this paper, the comparison between the optimization techniques is carried out by considering the same search agents (30 search agents) and total iterations (500 iterations). The selection of various parameter values for the GOA algorithm must be made carefully for implementation and efficient performance. The following are the many cases that were chosen here: calculation of the objective function owing to the occurrence of any disturbance in the system using the MATLAB R2016a environment. Further, the controller parameters can be found by minimization of the fitness values of Equation (2), whose optimization is carried out using the GOA algorithm.  Table 2. In this proposed simulation work, different parameters have been initialized for the application of GOA, GA, DE, and WOA. Maximum no. of generation and population size is the common control parameters of the algorithms. In this paper, the comparison between the optimization techniques is carried out by considering the same search agents (30 search agents) and total iterations (500 iterations). The selection of various parameter values for the GOA algorithm must be made carefully for implementation and efficient performance. The following are the many cases that were chosen here:   The performance of the suggested controller is demonstrated at P e = 0.85 pu and δ 0 = 52.3 deg, which are the nominal loading conditions in terms of the occurrence of a severe disturbance in the system. A 3-cycles, 3-phase fault is imposed at the mid-section of the transmission line linking bus-2 and bus-3 at time t = 1 s, and the system returns to its previous state after it is cleared. The system's various responses are depicted in Figures 10-13, including speed deviation in pu, tie-line power P L in MW, power angle in δ (degree), and SSSC injected voltage V q in pu. From the different responses, we reached the conclusion that the GOA-optimized suggested controller provides improved dynamic response when compared with the well-established GA, DE algorithm, and recently published WOAoptimized controller. It can also be shown that the recommended GOA-optimized controller has good low-frequency oscillation damping capabilities and can easily stabilize the device by adjusting the SSSC-injected voltage. In comparison with recently published WOA, DE, and GA-optimized controllers, the suggested technique offers improved dynamic response in terms of minimal overshoot, minimum undershoot, and settling time, as shown in Figures 10-13. As a result, the suggested controller increases the restriction on power system stability and capacity. Figure 14 displays various transmission delays considering speed deviation. Here it is observed that transport delays have a significant inverse effect on the system responses.

Case A: Nominal Loading Condition
The performance of the suggested controller is demonstrated at Pe = 0.85 pu and 0  = 52.3 deg, which are the nominal loading conditions in terms of the occurrence of a severe disturbance in the system. A 3-cycles, 3-phase fault is imposed at the mid-section of the transmission line linking bus-2 and bus-3 at time t = 1 s, and the system returns to its previous state after it is cleared. The system's various responses are depicted in Figures 10-13, including speed deviation in pu, tie-line power PL in MW, power angle in  (degree), and SSSC injected voltage Vq in pu. From the different responses, we reached the conclusion that the GOA-optimized suggested controller provides improved dynamic response when compared with the well-established GA, DE algorithm, and recently published WOA-optimized controller. It can also be shown that the recommended GOA-optimized controller has good low-frequency oscillation damping capabilities and can easily stabilize the device by adjusting the SSSC-injected voltage. In comparison with recently published WOA, DE, and GA-optimized controllers, the suggested technique offers improved dynamic response in terms of minimal overshoot, minimum undershoot, and settling time, as shown in Figures 10-13. As a result, the suggested controller increases the restriction on power system stability and capacity. Figure 14 displays various transmission delays considering speed deviation. Here it is observed that transport delays have a significant inverse effect on the system responses.

Case B: Light Loading Condition
To test the superiority of the suggested controller, the generator's loading state is changed to light loading (Pe = 0.5 pu and δ 0 = 33.23 deg). At the midpoint of the transmission line, a 5-cycle 3-phase fault arises, followed by load removal at t = 1.0 s at bus-1. Figures 15-17 show the system's response to this possibility, which clearly illustrates the quality of the suggested controller for changes in working conditions and types of disturbance. In addition, when compared with the well-known GA, DE algorithm, and the recently published WOA-optimized controller, the proposed GOA approach delivers improved transient response. It can be seen that the GOA-optimized coordinated design of the PSS and SSSC controller substantially suppresses the rotor angle swing, which provides good damping characteristics for electromechanical oscillation modes.  Figures 15-17 show the system's response to this possibility, which clearly illustrates the quality of the suggested controller for changes in working conditions and types of disturbance. In addition, when compared with the well-known GA, DE algorithm, and the recently published WOA-optimized controller, the proposed GOA approach delivers improved transient response. It can be seen that the GOA-optimized coordinated design of the PSS and SSSC controller substantially suppresses the rotor angle swing, which provides good damping characteristics for electromechanical oscillation modes.

Case C: Heavy Loading Condition
The robustness of this controller is carried out in this condition (Pe = 1.0 pu and δ 0 = 60.73 deg). In this situation, the load is disconnected near bus 1 at t = 1.0 s for 300 ms. Figure 18 depicts the speed deviation response under heavy loading conditions. In comparison to the well-established GA, DE algorithm, and recently published WOA-optimized controller, the effectiveness of the GOA-optimized proposed controller provides more reliable performance. It is evident from the figure that system stability is retained, and power system oscillations with the proposed controllers are effectively damped out. The ITAE values for the three situations mentioned above are plotted in Figure 19 to demonstrate the enhancement by the proposed GOA technique. The robustness of this controller is carried out in this condition (Pe = 1.0 pu and 0  = 60.73 deg). In this situation, the load is disconnected near bus 1 at t = 1.0 s for 300 ms. Figure 18 depicts the speed deviation response under heavy loading conditions. In comparison to the well-established GA, DE algorithm, and recently published WOAoptimized controller, the effectiveness of the GOA-optimized proposed controller provides more reliable performance. It is evident from the figure that system stability is retained, and power system oscillations with the proposed controllers are effectively damped out. The ITAE values for the three situations mentioned above are plotted in Figure 19 to demonstrate the enhancement by the proposed GOA technique.    The SMIB system's transient response parameters using proposed controller optimized with DE/GA/WOA/GOA techniques are displayed in Table 3.   Table 3.

Extension to Multi-Machine Power System
The power system involving multiple machines and tie lines is represented in Figure 20. In between bus 5 and bus 6, the SSSC controller is placed as shown in Figure 20.

Extension to Multi-Machine Power System
The power system involving multiple machines and tie lines is represented in Figure 20. In between bus 5 and bus 6, the SSSC controller is placed as shown in Figure 20. The various scenarios selected are as follows.

Scenario 1: Three-Phase Fault Disturbance
At t = 1 s, in one of the line sections joining bus 1 and bus 6, a 3-phase self-clearing fault is implemented. The duration of fault is taken as a 3-cycle. At bus 6, after clearing it, the system gains back to its original state. The response of the system is exposed in Figures 21 and 22. From these figures, the local mode and inter-area modes of oscillations are observed to be greatly oscillatory in the absence of controllers. So, the suggested controllers aim at improving the stability of the power system significantly by modulation of the SSSC injected voltage and stabilization of the signals of PSSs and thereby suppress these The various scenarios selected are as follows.

Scenario 1: Three-Phase Fault Disturbance
At t = 1 s, in one of the line sections joining bus 1 and bus 6, a 3-phase self-clearing fault is implemented. The duration of fault is taken as a 3-cycle. At bus 6, after clearing it, the system gains back to its original state. The response of the system is exposed in Figures 21 and 22. From these figures, the local mode and inter-area modes of oscillations are observed to be greatly oscillatory in the absence of controllers. So, the suggested controllers aim at improving the stability of the power system significantly by modulation of the SSSC injected voltage and stabilization of the signals of PSSs and thereby suppress these oscillations. The effectiveness of the proposed coordinated design approach is examined by varying the transport delay signals. At t = 1 s, in one of the line sections joining bus 1 and bus 6, a 3-phase self-clearing fault is implemented. The duration of fault is taken as a 3-cycle. At bus 6, after clearing it, the system gains back to its original state. The response of the system is exposed in Figures 21 and 22. From these figures, the local mode and inter-area modes of oscillations are observed to be greatly oscillatory in the absence of controllers. So, the suggested controllers aim at improving the stability of the power system significantly by modulation of the SSSC injected voltage and stabilization of the signals of PSSs and thereby suppress these oscillations. The effectiveness of the proposed coordinated design approach is examined by varying the transport delay signals.    Initially, a 50 ms time delay for signal transmission is considered. The effectiveness of the proposed control approach is examined by varying the delay time. Figures 23 and 24 illustrate the response of the system, which shows that the variation in time delay doesn't have a remarkable effect on the proposed controller performance.

Scenario 2: Small Disturbance
Further, a small or minor disturbance in the system is considered to ensure the complete examination of the suggested controller's performance. The small disturbance is simulated by detaching the load at bus 4 at t = 1.0 s for a duration of 100 ms. Figures 25 and 26 illustrate the system response for the above-mentioned situation and from which the robustness and efficient damping ability of the proposed controller even under small disturbance conditions is clearly justified. The figure shows the efficacy of PSS and SSSC in the damping of low-frequency oscillations. From the figures, it is obvious that the GOA-optimized coordinated design is stable and ensures stability for large variations in loading conditions. As can be seen from Figures 25 and 26, the proposed GOA-optimized controller significantly outperforms the WOA, DE, and GA-optimized proposed controller in terms of dampening power system oscillations and reducing settling time. The ITAE values for the above-mentioned two scenarios of the multi-machine system are displayed in a bar plot, as displayed in Figure 27, to better illustrate the improvement by the proposed approach. Further, a small or minor disturbance in the system is considered to ensure the complete examination of the suggested controller's performance. The small disturbance is simulated by detaching the load at bus 4 at t = 1.0 s for a duration of 100 ms. Figures 25 and 26 illustrate the system response for the above-mentioned situation and from which the robustness and efficient damping ability of the proposed controller even under small disturbance conditions is clearly justified. The figure shows the efficacy of PSS and SSSC in the damping of low-frequency oscillations. From the figures, it is obvious that the GOAoptimized coordinated design is stable and ensures stability for large variations in loading conditions. As can be seen from Figures 25 and 26, the proposed GOA-optimized controller significantly outperforms the WOA, DE, and GA-optimized proposed controller in terms of dampening power system oscillations and reducing settling time. The ITAE values for the above-mentioned two scenarios of the multi-machine system are displayed in a bar plot, as displayed in Figure 27, to better illustrate the improvement by the proposed approach.      Table 4 shows the optimum parameters of the proposed controller for multi-machine systems. Table 5 illustrates that GOA methods that take into consideration multiple instances had the lowest ITAE values when compared with DE, GA, and WOA methods. From the simulation the following observations are made:  Table 4 shows the optimum parameters of the proposed controller for multi-machine systems. Table 5 illustrates that GOA methods that take into consideration multiple instances had the lowest ITAE values when compared with DE, GA, and WOA methods.  From the simulation the following observations are made: 1.
Initially, a SMIB system with a fuzzy lead-lag structured controller is considered, and the dominance of GOA as related to WOA, DE, and GA is demonstrated.

2.
The percentage decrease in J value with the GOA technique compared with recently published WOA, well-established GA, and DE are 16.32%, 14.56%, and 13.72%, respectively, in the SMIB system. The percentage decrease in J value with the proposed GOA technique compared with recently published WOA, well-established GA, and DE are 1.41%, 9.98%, and 13.31%, respectively, for the multi-machine system. 3.
The effectiveness of the GOA is validated with various load conditions and results confirmed that GOA performed better compared with the recently published WOA, and well-established GA and DE optimization algorithms.
The wide adoption of the GOA algorithm in recent optimization problems is due to its exploration ability in the initial stage of optimization. Furthermore, this assists in exploring the promising regions in the search space. Due to large attraction forces between the grasshoppers, exploitation is high in the last step. The GOA maintains a smooth balance between exploration and exploitation to avoid convergence to local optima points. This behavior is due to the proposal of the adaptive comfort coefficient. The application of GOA in the electrical engineering domain enhances its reliability and accuracy by its capability of solving many unimodal as well as multi-modal complex problems.
The trends toward the growth of meta-heuristic algorithms enhanced in recent years due to their potential benefits, including (1) dynamic condition flexibility, (2) fast convergence speed, and (3) solving convex problems.
Computer engineering and electrical engineering fields dominate the application of GOA from the analysis of the literature.

Conclusions
In the present work, the investigation of the enhancement of the power system steadiness by the coordinated application of numerous damping controllers is thoroughly achieved. A time-domain objective function is utilized to diminish the oscillations of the power system for the controller design problem that has been proposed. Then a GOA technique is employed to achieve proper tuning of the controller parameters optimally and coordinately. Further, the effectiveness of the suggested coordinated strategy is tested through the simulation results obtained at various load conditions and disturbances. It is observed that the suggested controllers are robust enough to manage fault sites and variations in operating conditions. Thereby, the stability of the system is enhanced by generating the stabilizing output signals. Finally, when the proposed approach is applied to a multi-machine system and simulation is carried out, the consequences lead to the decision that the GOA-optimized intended controller becomes superior compared with the well-established GA, DE algorithm, and recently published WOA-optimized controller. The percentage decrease in J value with GOA technique compared with recent published WOA, well-established GA and DE are 16.32%, 14.56%, and 13.72%, respectively, in the SMIB system. The percentage decrease in J value with proposed GOA technique compared with recent published WOA, well established GA and DE are 1.41%, 9.98%, and 13.31%, respectively, for a multi-machine system. The simulation outcomes substantiate that the intended controller is particularly versatile for the real-world application to power systems.