Transient Stability Analysis of a Multi-Machine Power System Integrated with Renewables

: The impact on the stability of power systems is rising as the penetration level of renewable energy with sporadic natures rises rapidly on the grid. However, the impact of different types of renewable energy sources (wind, solar) and their combination on system stability varies even with the same penetration level. This paper concentrates mainly on the stability analysis of multi-machine systems connected to various types of renewable energy sources. The study presents a simple and novel control technique named automatic reactive power support (ARS) for both single and combinations of renewable sources by injecting the available reactive power into the system during fault through converters to enhance system stability. The permanent magnet synchronous generator (PMSG) and doubly fed induction generator (DFIG) are both considered as wind generators in this paper for comparison. In addition, transient stability enhancement is carried out by improving critical clearing time of a three-phase fault in the power system. With the creation of a 3-phase fault at various buses, stability analysis is carried out on the 9-bus WSCC test bus system and also on the 68-bus IEEE test system. Comparative analysis of six test case conditions is provided and the considered cases are without renewable source, with DFIG as a wind generator, PMSG as a wind generator, solar PV farm, wind farm with DFIG and solar PV in combination and the combination of wind farm with PMSG and solar PV. Moreover, the improvement in critical clearing time of the system is compared using conventional and proposed controls with all the aforementioned renewable sources. Comparative results show that the proposed control technique improves system stability and also that the combination of renewable energy sources ought to enhance the critical clearing time of system.


Introduction
The integration of non-conventional energy sources, specifically solar PV power and wind power generation sources, with the grid has risen drastically in the past decade [1].As advanced power electronics develop, solar and wind energy become the most attractive and promising types of renewable energy sources [2,3].However, the challenges, such as the impact of a high penetration level of renewable energy sources (RES), and their possible solutions are focused and addressed in the literature [4,5].
In [6], studies were conducted to provide electric energy to a local community with 100 households and a health center by using a wind-solar hybrid power generation system in Ethiopia.The effect of partial shading in the PV module is investigated in [7][8][9].The benefits and drawbacks of various maximum power point tracking (MPPT) strategies in solar PVs and their importance is explained in [10][11][12] to assure maximum power is delivered to the system.However, the difficulties in the control techniques of MPPT with proportional-integral-derivative controllers are minimized by using the model predictive Energies 2022, 15, 4824 2 of 18 control [13] and the sliding mode control [14].The workings of wind energy systems with an implementation of MPPT are investigated in [15][16][17][18][19]. Studies of different types of wind generators are carried out in [20][21][22][23].The impact of grid-integrated wind farms, with a DFIG as wind generator, on the transient stability of power systems is investigated in [24,25].Stability analysis of grid-integrated solar PV systems is seen in [26][27][28].The decoupled control approach is a well-adapted control method which allows for controlling real and reactive power independently and is considered in [27].To enhance stability, the grid-side converter (GSC) behaves as a STATCOM while the DFIG operates at a constant speed during the fault.Solar system inverters can be configured to regulate renewable power generation in response to changes in system frequency and voltage [28].
Recently, studies have been more intense regarding the integration of combinations of renewable sources which ensure the optimum utilization of resources and hence improve overall efficiency as compared to single modes of operation.The hybridization of solar and wind systems exploits the advantages of both while at the same time easing limitations.The power generated by the hybrid renewables can be extracted more efficiently by implementing new technologies in power electronics.The implementation and working of hybrid systems of wind and solar energy, their advantages and disadvantages.and different topologies of hybridization are explained in [29], and the power system equivalents for dynamic and transient studies are created and analyzed in [30].The effect on transient stability with hybrid models of hydro power, solar, and wind power injected to a single machine connected to an infinite bus system is investigated in [31].The improvement of system stability connected to hybrid renewable sources with arbitrary amounts of reactive power support is mentioned in [32].
The main aim of the proposed research work is introducing a simple, new technique to improve multi-machine system (MMS) transient stability coupled with both individual and combinations of renewable sources.The concern about system stability is increasing as the share of renewable energy injection with variable power natures is expanding.Objectives of the present work are as follows: • Modifying the controllers of the converters such that all RES connected to the system will be injecting the maximum available reactive power into the grid during faults.Normally, the converters are configured to inject only real power into the grid (with unity power factor).The injection of reactive power during faults improves the bus voltage profile and hence enhances the overall stability of system.

•
Performing a comparative analysis of critical clearing time for a MMS coupled with different renewable sources.
The proposed controller's performance is validated with an individual RES and as well as hybrid renewables along with different fault locations at different fault instants.The study is performed on a WSCC 9-bus system and also IEEE 68-bus system.
The rest of the article is organised as follows.Section 2 presents the mathematical modelling of a complete power system including the RES.Both the existing control techniques and the proposed control techniques used in the RES are explained in Section 3. Section 4 is carried out with analysis of system results followed by the conclusion in Section 5.

Mathematical Model of the System
Modelling of power system components, such as synchronous machines, a power system network, solar PVs, DFIGs, and PMSGs, along with the turbine, rotating mass, and also the integration of these RES to the power system network are presented in this section.

Modelling of the Multi-Machine System
The synchronous machine model considered for study is the flux decay model with a static exciter.As the time constant of the speed governor system (electro-mechanical phenomena) is in the order of seconds and the time constant of machine dynamics (electromagnetic phenomena) followed by sudden disturbance is in the range of milliseconds to a second [33], the dynamics of the speed governor are neglected.The machine equations (for Energies 2022, 15, 4824 3 of 18 n machine systems with i = 1 through n) including a static exciter with one time constant and one gain can be summarized as [24,34] where δ represents the rotor angle, ω s and ∆ω i are the synchronous speed and deviation in rotor speed (per unit), respectively, K D is damping coefficient, H represents the inertia constant, P e and P m are the output electrical power and input mechanical power, respectively, the q-axis and d-axis components of synchronous reactance are represented by x q and x d , respectively, x q and x d are the q-axis and d-axis components of transient reactance, E qi is the q-axis component of voltages behind the transient reactance of the ith generator, T doi is the d-axis open circuit time constant, T A and K A are the time constant and gain of the exciter, E fd represents the exciter voltage, V is the per unit terminal voltage of the machine, and angle is represented by θ.

Modelling of Solar PV
Photovoltaic (PV) cells are used to convert solar energy into electricity through a photoelectric effect.An ideal representation of a PV cell is a current source in parallel with a diode.However, two resistors, one in series and the other in parallel, are also included in practice.Here, series resistors serve as internal resistance to current source and shunt resistors include the leakage current.
The equation to describe characteristics I-V of a PV cell is given as [7]: where I L and I D are photocurrent and reverse saturation current of diode, respectively.K denotes Boltzmann's constant, T is the temperature, 'α' is the diode ideality factor, and 'q' represents charge of electron.R SH and R S represent the equivalent shunt and series resistance of the PV cell.A connection diagram of a solar PV farm to the grid is shown in Figure 1.The power generated by the solar PV is injected to the grid through a boost converter, voltage source converter (VSC) filter, and transformer.Equations of solar PVs are modelled assuming that power generated is always maximum.

Modelling of Turbine and Rotating Mass
The wind farm contains turbines and generators which have rotating masses, unlike the solar PV farm.These turbines and rotating masses have their own dynamics and the same are modelled in this section.A two-mass model representation is used to model both the wind generator and the wind turbine.The equations are as follows [24,35].

Modelling of Turbine and Rotating Mass
The wind farm contains turbines and generators which have rotating masses, unlike the solar PV farm.These turbines and rotating masses have their own dynamics and the same are modelled in this section.A two-mass model representation is used to model both the wind generator and the wind turbine.The equations are as follows [24,35].
Here, ω t and ω r are the mechanical speed of the turbine and rotor, respectively, H g and H t are the generator and the turbine inertia, respectively, θ tw represents the torsional angle of shaft, and C sh and k sh are the damping coefficient and shaft stiffness, respectively.T m and T e are the mechanical and electrical torque, respectively.

Modelling of the DFIG
The connection diagram of a wind farm with a DFIG is shown in Figure 2. The DFIG is connected to the system through a transformer, and the rotor of the DFIG is also connected to the system through a back-to-back converter.The total power injected into the grid is the algebraic sum of powers delivered from the DFIG and the grid-side converter (GSC).The mathematical model of the DFIG [24] is represented in state space form with d-q representation and the following are the equations. where where 'v' and 'i' are voltage and current of the machine, respectively.The suffixes 'qs', 'qr', 'ds', and 'dr' are q-axis and d-axis components of stator and rotor quantities, respectively.'ω eB ' is base electrical speed and 'ω s ' is the synchronous speed.'L' and 'R' are the inductance and resistance of the machine, respectively.The suffixes 'ss', 'rr', and 'm' are self and mutual components, respectively.In addition, the equivalent q-axis and d-axis source voltages behind transient reactance "e ' q " and "e ' d " are as follows 'ωeB' is base electrical speed and 'ωs' is the synchronous speed.'L' and 'R' are the inductance and resistance of the machine, respectively.The suffixes 'ss', 'rr', and 'm' are self and mutual components, respectively.In addition, the equivalent q-axis and d-axis source voltages behind transient reactance "e ' q" and "e ' d" are as follows (10)

Modelling of the PMSG
The PMSG (Permanent Magnet Synchronous Generator) is mathematically modeled from the equivalent d-axis and q-axis circuit of the machine.Here, dq0-axis modelling of the PMSG is used, in which the alignment of the d-axis is along the magnet axis.The daxis and q-axis voltage equations of the PMSG in terms of current are given by [36,37]  The electromagnetic torque is as follows ω s e q i qs + e d i ds . (10)

Modelling of the PMSG
The PMSG (Permanent Magnet Synchronous Generator) is mathematically modeled from the equivalent d-axis and q-axis circuit of the machine.Here, dq0-axis modelling of the PMSG is used, in which the alignment of the d-axis is along the magnet axis.The d-axis and q-axis voltage equations of the PMSG in terms of current are given by [36,37] where V qs , V ds , i qs , and i ds are q-axis and d-axis components of stator voltages and currents, respectively.L q and L d are the q-axis and d-axis inductances, respectively.The stator resistance is represented by R s , ω rot is the electrical rotor speed in rad/s, and λ M represents the rotor magnetic flux produced by the permanent magnet.The connection diagram of a wind farm with a PMSG is shown in Figure 3.The PMSG is connected to the system through the back-to-back converter and the transformer.
Energies 2022, 15, x FOR PEER REVIEW 6 of 19 where Vqs, Vds, iqs, and ids are q-axis and d-axis components of stator voltages and currents, respectively.Lq and Ld are the q-axis and d-axis inductances, respectively.The stator resistance is represented by Rs, ωrot is the electrical rotor speed in rad/s, and λM represents the rotor magnetic flux produced by the permanent magnet.The connection diagram of a wind farm with a PMSG is shown in Figure 3.The PMSG is connected to the system through the back-to-back converter and the transformer.

Integration of RES with the Grid
Integration of renewable energy sources to the grid is a great task with an increased penetration level.The integration of hybrid renewables to the existing power system creates technical challenges including harmonic distortion, flicker, voltage regulation, etc.The challenges involved in the integration of RES with the grid and the advanced techniques to overcome those challenges are presented in [38,39].The power generated by the RES is intermittent in nature and can cause stability issues grid-side.The synchronization of voltage magnitude and the frequency of RES with the grid is the most important aspect in integrating RES to the grid.The injection of power must also occur with less harmonics.

Integration of RES with the Grid
Integration of renewable energy sources to the grid is a great task with an increased penetration level.The integration of hybrid renewables to the existing power system creates technical challenges including harmonic distortion, flicker, voltage regulation, etc.The challenges involved in the integration of RES with the grid and the advanced techniques to overcome those challenges are presented in [38,39].The power generated by the RES is intermittent in nature and can cause stability issues grid-side.The synchronization of Energies 2022, 15, 4824 6 of 18 voltage magnitude and the frequency of RES with the grid is the most important aspect in integrating RES to the grid.The injection of power must also occur with less harmonics.

Controllers of Renewable Energy Sources
Various types of controllers used in wind (DFIG, PMSG) and solar PV systems are explained in this section.Controllers are a vital part of any system.Three system variables must be controlled strictly in a wind generation system.They are: (1) the optimal power generated; (2) the injected power (active and reactive) to the utility grid; and (3) the DC-link voltage.In this study, the conventional controllers are modified such that the transient stability of the system can be enhanced.Details of conventional and modified controllers are presented in this section.Controllers of all converters of RES are modelled using the dq0 reference frame.

Conventional Controllers in the DFIG
Wind farms with a DFIG basically have two converters.One is a rotor side converter (RSC), and the other is a grid-side converter (GSC).The conventional controllers used in these two converters are presented here.

Controllers in the Rotor Side Converter (RSC)
Orientation of the d-axis is considered along with the voltage for the RSC of a DFIG control.The expressions for electromagnetic torque and reactive power in terms of d-axis and q-axis rotor current are as follows [24] The reference values of both q-axis and d-axis components of rotor currents are calculated by re-arranging ( 13) and ( 14) as follows 3.1.2.Controllers in the Grid-Side Converter (GSC) The main aim of the GSC of both PMSG and DFIG controllers is to uphold the DC-link voltage at a constant level and to regulate the flow of both reactive and active power independently between the grid and the inverter.In the three-phase balanced system, the instantaneous active (P) and reactive power (Q) outputs are described by the following equations where, V q , V d , I q , and I d are q-axis and d-axis components of grid voltages and currents, respectively.For the easy control of P and Q, a control approach is applied based on voltage orientation.It is assumed that the alignment of the d-axis of the reference frame is aligned with the space vector of grid voltage.This makes the q-axis component of space vector for the grid voltage become zero (i.e., V q = 0).Thus, the expression for active and reactive power of Equations ( 17) and ( 18) can be modified as follows: From ( 19) and ( 20) we can control P and Q independently by d-axis and q-axis components of current (I d and I q ), respectively.The reference for the q-axis component of current is chosen to be zero (i.e., I qpref = 0) in order to maintain unity power factor for the power injected into the grid.Voltage equations between the grid and inverter are as follows [24] Here, R g and L g are the resistance and inductance of line, respectively.The conventional control block of the GSC is shown in Figure 4.

Conventional Controllers in the PMSG
Wind farms with a PMSG also have two converters.One is a machine side converter (MSC), and the other is a grid-side converter (GSC).The conventional controllers used in these two converters are presented here.

Controllers in the Machine Side Converter (MSC)
Extraction of maximum power is the aim of the MSC controller for the PMSG by controlling the speed of the PMSG rotor.The MPPT controller starts to operate when the speed of the wind is greater than the cut in speed, and it will stop when the wind speed surpasses the rated value.The torque (Te) equation of the permanent magnet synchronous generator is given by [37] ( ) where p is the number of pole pairs, λM is the magnetic flux produced by the permanent magnets in the PMSG rotor, and remaining quantities hold the same definition except that these quantities belong to the PMSG.The electromagnetic torque Te can be controlled independently by q-axis current iq alone by assuming that the d-axis current id is equal to zero.Thus, the torque expression will be simplified to

Conventional Controllers in the PMSG
Wind farms with a PMSG also have two converters.One is a machine side converter (MSC), and the other is a grid-side converter (GSC).The conventional controllers used in these two converters are presented here.

Controllers in the Machine Side Converter (MSC)
Extraction of maximum power is the aim of the MSC controller for the PMSG by controlling the speed of the PMSG rotor.The MPPT controller starts to operate when the speed of the wind is greater than the cut in speed, and it will stop when the wind speed surpasses the rated value.The torque (T e ) equation of the permanent magnet synchronous generator is given by [37] where p is the number of pole pairs, λ M is the magnetic flux produced by the permanent magnets in the PMSG rotor, and remaining quantities hold the same definition except that these quantities belong to the PMSG.The electromagnetic torque T e can be controlled independently by q-axis current i q alone by assuming that the d-axis current i d is equal to zero.Thus, the torque expression will be simplified to It is clear from (24) that torque can be controlled solely by i q .The control structure of the MSC is derived from the voltage-current equation and torque equation as shown in Figure 5. Reference torque is calculated from optimal speed of rotation.I qref can be calculated from reference torque which is kept to an optimum value.The error signals are generated by comparing actual currents with the reference values.Error signal through PI of the controller gives d-axis and q-axis reference rotor voltages (V dref and V qref ).Reference phase voltages in the abc frame will be obtained by converting dq0 voltages.With the use of these signals in the abc frame and PWM, switching pulses of the MSC are then generated.Extraction of optimal wind power delivers the needed torque or power reference, which is given by where Here, ωtB and SB are the base speed and the base power of the wind turbine, respectively.CPmax is the maximum value of CP, the wind turbine's coefficient of performance, which is obtained when pitch angle (β) = 0°; λopt is the tip speed ratio when CP = CPmax.

Controllers in the GSC
The GSC of the wind farm with the PMSG is the same as the wind farm with the DFIG.The control strategy, as illustrated in Figure 4, and the equations to model the controller are also the same.

Pitch Angle Control
If the speed of the wind exceeds its rated value, a mechanical method is used to protect the blade from being damaged by controlling blade angle.For wind speed below the rated value, the maximum power point tracking technique is implemented such that the maximum amount of power can be extracted from wind velocity.If the wind speed is very high (more than the rated velocity), the power extraction is limited by pitch angle control.This is achieved by turning the blades away from the wind.

Conventional Controllers in Solar PVs
Solar PVs and the grid are connected through back-to-back converters consisting of a DC-DC chopper circuit followed by a DC-AC voltage source converter (VSC).Here, ω tB and S B are the base speed and the base power of the wind turbine, respectively.C Pmax is the maximum value of C P , the wind turbine's coefficient of performance, which is obtained when pitch angle (β) = 0 • ; λ opt is the tip speed ratio when C P = C Pmax .

Controllers in the GSC
The GSC of the wind farm with the PMSG is the same as the wind farm with the DFIG.The control strategy, as illustrated in Figure 4, and the equations to model the controller are also the same.

Pitch Angle Control
If the speed of the wind exceeds its rated value, a mechanical method is used to protect the blade from being damaged by controlling blade angle.For wind speed below the rated value, the maximum power point tracking technique is implemented such that the maximum amount of power can be extracted from wind velocity.If the wind speed is very high (more than the rated velocity), the power extraction is limited by pitch angle control.This is achieved by turning the blades away from the wind.

Conventional Controllers in Solar PVs
Solar PVs and the grid are connected through back-to-back converters consisting of a DC-DC chopper circuit followed by a DC-AC voltage source converter (VSC).

Controllers in the DC-DC Converter
DC-DC converters mainly use the MPPT technique to transfer maximum power into the system.The MPPT controller is a completely electronics-based control system used to extract the maximum power available at the PV module.In this work it is assumed that the power injected into the grid is always the maximum.

Controllers in the VSC
A VSC is present next to the DC-DC chopper circuit which converts DC voltage to 3-phase AC.VSC controllers aim to maintain DC-link voltage as a constant value and to regulate the flow of reactive and active power between the grid and inverters independently.The same control strategy as in Figure 4 is implemented for the GSC.The equations to model the controller are also identical to ( 19)- (22).Keeping the reactive power reference to zero is also followed here such that the power injected into the is maintained at unity power factor.

Proposed Controllers
The integration of hybrid renewable energy sources to the grid is increasing in capacity and also affects the stability of the system.However, if the injected power (both active and reactive) is controlled tactically during a small duration of time of disturbance (such as a fault), system stability can be enhanced.Bus voltage magnitude will decrease during a fault.Injection of reactive power helps to improve the voltage profile and hence improves system stability.
Here, the controllers of the VSC in solar PVs and the GSC in both the DFIG and the PMSG are adjusted in a way so that both the RES inject maximum available reactive power throughout fault duration to achieve the enhancement in stability.This automated control of converters is termed as Automatic Reactive-power Support (ARS) in this paper.The procedure of ARS is the same for the VSC of solar PVs and the GSC of both DFIGs and PMSGs.Reactive power injected into the grid will be zero when there is no disturbance in the system.In contrast, during fault periods maximum available reactive power is injected to grid.This is achieved by changing the reference value of q-axis current to the maximum available current (i.e., I qref = I drated ) and I dref = 0.This modification in the controller is shown in Figure 6.The first block in Figure 6 is the switch, which has three inputs and one output.The detection of faults is carried out by comparing the RMS voltage of the bus (in pu) with 0.6pu.If V rms ≤ 0.6, then the switch gives I qref = I drated as output, otherwise the switch gives I qref = 0 as output.The combined effect improves terminal voltage and hence ensures stability enhancement.The real power injected will be zero in the case of solar systems and PMSG-based wind farms since these systems are connected to the grid through back-to-back converters, whereas in the DFIG-based wind farm the real power transferred through the back-to-back converter only kept to zero to achieve the proposed modification in the controller.A VSC is present next to the DC-DC chopper circuit which converts DC voltage to 3phase AC.VSC controllers aim to maintain DC-link voltage as a constant value and to regulate the flow of reactive and active power between the grid and inverters independently.The same control strategy as in Figure 4 is implemented for the GSC.The equations to model the controller are also identical to ( 19)- (22).Keeping the reactive power reference to zero is also followed here such that the power injected into the grid is maintained at unity power factor.

Proposed Controllers
The integration of hybrid renewable energy sources to the grid is increasing in capacity and also affects the stability of the system.However, if the injected power (both active and reactive) is controlled tactically during a small duration of time of disturbance (such as a fault), system stability can be enhanced.Bus voltage magnitude will decrease during a fault.Injection of reactive power helps to improve the voltage profile and hence improves system stability.
Here, the controllers of the VSC in solar PVs and the GSC in both the DFIG and the PMSG are adjusted in a way so that both the RES inject maximum available reactive power throughout fault duration to achieve the enhancement in stability.This automated control of converters is termed as Automatic Reactive-power Support (ARS) in this paper.The procedure of ARS is the same for the VSC of solar PVs and the GSC of both DFIGs and PMSGs.Reactive power injected into the grid will be zero when there is no disturbance in the system.In contrast, during fault periods maximum available reactive power is injected to grid.This is achieved by changing the reference value of q-axis current to the maximum available current (i.e., Iqref = Idrated) and Idref = 0.This modification in the controller is shown in Figure 6.The first block in Figure 6 is the switch, which has three inputs and one output.The detection of faults is carried out by comparing the RMS voltage of the bus (in pu) with 0.6pu.If Vrms ≤ 0.6, then the switch gives Iqref = Idrated as output, otherwise the switch gives Iqref = 0 as output.The combined effect improves terminal voltage and hence ensures stability enhancement.The real power injected will be zero in the case of solar systems and PMSG-based wind farms since these systems are connected to the grid through back-toback converters, whereas in the DFIG-based wind farm the real power transferred through the back-to-back converter only kept to zero to achieve the proposed modification in the controller.

Analysis of System Results
Two test bus systems, WSCC 9-bus and IEEE 68-bus, are considered for simulation study.In order to make the results comparable, a local additional load is connected at the same bus where the renewable source is connected such that the steady state solution of

Analysis of System Results
Two test bus systems, WSCC 9-bus and IEEE 68-bus, are considered for simulation study.In order to make the results comparable, a local additional load is connected at the same bus where the renewable source is connected such that the steady state solution of the test bus system is the same with and without a renewable source.Transient stability analysis is carried out by finding the critical clearing time (CCT) of the generators in the system following a 3-phase fault near the bus.CCT of generators is calculated and compared for both the systems with no renewable source, with a DFIG in a wind farm, a PMSG in a wind farm, a solar PV farm, and hybrid renewables (combination of wind farm with a DFIG and solar farm, as well as the combination of wind farm with a PMSG and solar farm) with both conventional and the proposed automatic reactive-power support.A total of 100 MW is considered as base power for the whole system.

WSCC 9-Bus Test System
The system under study, WSCC 9-bus, is shown in Figure 7 with renewable energy sources (RES).An additional load is also connected at bus 8.Additional load is used to nullify effects on the initial state of the system.The total power from RES injected into the system is 80 MW.However, for both the combinations of RES, it is considered that the power injected through the wind farm is 60 MW, and 20 MW power is injected through the solar PV farm so that the total injected power into the system would remain 80 MW to make the results comparable.A 3-phase fault is created at bus 5 and the rotor angle of both generators is observed.Study of the system was carried out in two steps and involved constant power injection with the assumption that wind speed and solar irradiation were constant with time, and the injection of varying power while considering that wind speed and solar irradiation changes with time.
Energies 2022, 15, x FOR PEER REVIEW 11 of 19 farm) with both conventional and the proposed automatic reactive-power support.A total of 100 MW is considered as base power for the whole system.

WSCC 9-Bus Test System
The system under study, WSCC 9-bus, is shown in Figure 7 with renewable energy sources (RES).An additional load is also connected at bus 8.Additional load is used to nullify effects on the initial state of the system.The total power from RES injected into the system is 80 MW.However, for both the combinations of RES, it is considered that the power injected through the wind farm is 60 MW, and 20 MW power is injected through the solar PV farm so that the total injected power into the system would remain 80 MW to make the results comparable.A 3-phase fault is created at bus 5 and the rotor angle of both generators is observed.Study of the system was carried out in two steps and involved constant power injection with the assumption that wind speed and solar irradiation were constant with time, and the injection of varying power while considering that wind speed and solar irradiation changes with time.

9-Bus Test System with Constant Power Injection
Initially, the wind speed is considered as 13.95 m/s and the solar irradiation as 1000 W/m 2 .A 3-phase fault is created at 0.2 s and cleared after 431 ms.Rotor angle of generator-2 is shown here in Figure 8a with and without different RES. Figure 8a shows that stability is improved with the power injection of RES with conventional control when compared to the system without RES.Rotor angle variations of generator-2 with ARS are shown in Figure 8b when the fault is created at 0.2 s and cleared after 484 ms, since the lowest CCT of generator-2 is 483 ms (when connected with a DFIG wind farm because of low power rating of back-to-back converters).The system voltage profile at bus 8 when connected to different types of RES is shown in Figure 9 for both conventional and proposed control techniques with the fault duration of 440 ms.It is clear from Figure 9a-e that the voltage profile is improved with the proposed control technique.The improvement in the voltage profile caused by the fault is because of the reactive power supply during the time period of the fault and hence improves the CCT of the system.

9-Bus Test System with Constant Power Injection
Initially, the wind speed is considered as 13.95 m/s and the solar irradiation as 1000 W/m 2 .A 3-phase fault is created at 0.2 s and cleared after 431 ms.Rotor angle of generator-2 is shown here in Figure 8a with and without different RES. Figure 8a shows that stability is improved with the power injection of RES with conventional control when compared to the system without RES.Rotor angle variations of generator-2 with ARS are shown in Figure 8b when the fault is created at 0.2 s and cleared after 484 ms, since the lowest CCT of generator-2 is 483 ms (when connected with a DFIG wind farm because of low power rating of back-to-back converters).The system voltage profile at bus 8 when connected to different types of RES is shown in Figure 9 for both conventional and proposed control techniques with the fault duration of 440 ms.It is clear from Figure 9a-e that the voltage profile is improved with the proposed control technique.The improvement in the voltage profile caused by the fault is because of the reactive power supply during the time period of the fault and hence improves the CCT of the system.The reactive power support given by the DFIG wind farm during the fault period is less compared to the PMSG wind farm with proposed control because of the low rating of the DFIG grid-side converter which is shown in Figure 10. Figure 10 clarifies that the reactive power supplied to the system during fault is larger when the combination of the The reactive power support given by the DFIG wind farm during the fault period is less compared to the PMSG wind farm with proposed control because of the low rating of the DFIG grid-side converter which is shown in Figure 10. Figure 10 clarifies that the reactive power supplied to the system during fault is larger when the combination of the The reactive power support given by the DFIG wind farm during the fault period is less compared to the PMSG wind farm with proposed control because of the low rating of the DFIG grid-side converter which is shown in Figure 10. Figure 10 clarifies that the reactive power supplied to the system during fault is larger when the combination of the PMSG wind farm and solar PVs is connected to the system since the rating of converters connected to the RES are the same as the rating of the RES.PMSG wind farm and solar PVs is connected to the system since the rating of converters connected to the RES are the same the rating of the RES.CCT of generator-2 (δ21) with no renewable source is 430 ms.CCT values of generator-2 for different renewable sources with conventional and proposed control techniques are tabulated in Table 1.Hybrid renewable sources (the combination of the DFIG wind farm and solar PVs) give a CCT of 465 ms with the conventional control which is an improvement of 35 ms, whereas the hybrid renewable source (PMSG wind farm and solar PV) gives a CCT of 551 ms with the proposed control, an improvement of 121 ms compared to when no renewable source (No RES) is connected.It is clear from Table 1 that the improvement in CCT is better when the combination of PMSG wind farm and solar PVs is connected to the system with proposed control ARS when compared to remaining cases.The profiles of both wind speed and solar irradiation are shown in Figure 11.Four different instances of time, namely A-D, are considered here to support different combinations of variations in both wind speed and solar irradiation such that one profile is in a positive slope and the other is in a negative slope, or vice versa, or one profile is in local minimum and the other is with some slope.The fault instants (time instants) are simulated by creating the fault at that particular time instant, and these fault instants represent time along the profiles of wind speed and solar irradiation at which the fault was created in the system for study.For the case of fault instant 'C', the fault is created at 40 s with wind speed 15.8 m/s and solar irradiation 800 W/m 2 , which can also be seen in Figure 11.CCT values of generator-2 for different RES and with variable power injection into the system for both conventional and proposed control techniques are tabulated in Table 2 with faults CCT of generator-2 (δ 21 ) with no renewable source is 430 ms.CCT values of generator-2 for different renewable sources with conventional and proposed control techniques are tabulated in Table 1.Hybrid renewable sources (the combination of the DFIG wind farm and solar PVs) give a CCT of 465 ms with the conventional control which is an improvement of 35 ms, whereas the hybrid renewable source (PMSG wind farm and solar PV) gives a CCT of 551 ms with the proposed control, an improvement of 121 ms compared to when no renewable source (No RES) is connected.It is clear from Table 1 that the improvement in CCT is better when the combination of PMSG wind farm and solar PVs is connected to the system with proposed control ARS when compared to remaining cases.

9-Bus Test System with Variable Power Injection
The profiles of both wind speed and solar irradiation are shown in Figure 11.Four different instances of time, namely A-D, are considered here to support different combinations of variations in both wind speed and solar irradiation such that one profile is in a positive slope and the other is in a negative slope, or vice versa, or one profile is in local minimum and the other is with some slope.The fault instants (time instants) are simulated by creating the fault at that particular time instant, and these fault instants represent time along the profiles of wind speed and solar irradiation at which the fault was created in the system for study.For the case of fault instant 'C', the fault is created at 40 s with wind speed 15.8 m/s and solar irradiation 800 W/m 2 , which can also be seen in Figure 11.CCT values of generator-2 for different RES and with variable power injection into the system for both conventional and proposed control techniques are tabulated in Table 2 with faults at different instances which elucidates that the combination of PMSG and solar with ARS gives better improvements in CCT even with variable power injection into the system.at different instances which elucidates that the combination of PMSG and solar with ARS gives better improvements in CCT even with variable power injection into the system.Moreover, Table 1 specifically shows the CCT of generator-2 when the system is connected to RES and the power injected by the RES is maintained at a constant throughout the study.Table 2 shows the same CCT with power injected into the system varying with respect to time.Therefore, the CCT values of particular renewable energy sources are different in Tables 1 and 2.Moreover, Table 1 specifically shows the CCT of generator-2 when the system is connected to RES and the power injected by the RES is maintained at a constant throughout the study.Table 2 shows the same CCT with power injected into the system varying with respect to time.Therefore, the CCT values of particular renewable energy sources are different in Tables 1 and 2.

IEEE 68-Bus Test System
The system under study, IEE 16-machine 68-bus, is shown in Figure 12.A renewable source is connected at different buses considering five different zones [24].The total power from RES injected into the system is 400 MW.However, for both combinations of RES, it is considered as 300 MW through the wind farm and 100MW through the solar farm in order to make the results comparable.The buses in the different zones are as follows:

IEEE 68-Bus Test System
The system under study, IEE 16-machine 68-bus, is shown in Figure 12.A source is connected at different buses considering five different zones [24 power from RES injected into the system is 400 MW.However, for both com RES, it is considered as 300 MW through the wind farm and 100MW throu farm in order to make the results comparable.The buses in the different zone lows:

68-Bus Test System with Constant Power Injection
A 3-phase fault is created at bus 53 when the renewable energy source i at bus 28 in zone-II and the rotor angle of all generators is observed.Rotor a least stable generators are shown in the Figure 13a for both with and witho RES.CCT of the system with no renewable source is 664 ms.Rotor angle var ARS are shown in Figure 13b.It is clear from Figure 13b that the variations in are better when the combination of wind farm with a PMSG and solar PV f nected to the system with proposed control ARS.CCT values for differen sources with conventional and proposed control techniques are tabulated in T brid renewable sources (combination of PMSG wind farm and solar PV) giv 684 ms with conventional control which is an improvement of 20 ms, where combination gives a CCT of 751 ms with the proposed control, an improvem compared to when no renewable source is connected.In addition, the impr CCT with ARS is better in the combination of wind farm with PMSG and sol

68-Bus Test System with Constant Power Injection
A 3-phase fault is created at bus 53 when the renewable energy source is connected at bus 28 in zone-II and the rotor angle of all generators is observed.Rotor angles of the least stable generators are shown in the Figure 13a for both with and without different RES.CCT of the system with no renewable source is 664 ms.Rotor angle variations with ARS are shown in Figure 13b.It is clear from Figure 13b that the variations in rotor angle are better when the combination of wind farm with a PMSG and solar PV farm is connected to the system with proposed control ARS.CCT values for different renewable sources with conventional and proposed control techniques are tabulated in Table 3. Hybrid renewable sources (combination of PMSG wind farm and solar PV) gives a CCT of 684 ms with conventional control which is an improvement of 20 ms, whereas the same combination gives a CCT of 751 ms with the proposed control, an improvement of 87 ms compared to when no renewable source is connected.In addition, the improvement in CCT with ARS is

Energies 2022 , 19 Figure 1 .
Figure 1.Connection diagram of a solar PV farm to the system.

Figure 1 .
Figure 1.Connection diagram of a solar PV farm to the system.

Figure 2 .
Figure 2. Connection diagram of a wind farm with a DFIG to the system.The electromagnetic torque is as follows

Figure 2 .
Figure 2. Connection diagram of a wind farm with a DFIG to the system.

Figure 3 .
Figure 3. Connection diagram of a wind farm with a PMSG to the system.

Figure 3 .
Figure 3. Connection diagram of a wind farm with a PMSG to the system.
Energies 2022, 15, x FOR PEER REVIEW 8 of 19Here, Rg and Lg are the resistance and inductance of line, respectively.The conventional control block of the GSC is shown in Figure4.

Figure 6 .
Figure 6.Proposed controller (ARS) for the q-axis control block.

Figure 6 .
Figure 6.Proposed controller (ARS) for the q-axis control block.

Figure 8 .Figure 9 .
Figure 8. Variations in rotor angle (δ21) with: (a) conventional control and (b) proposed control (ARS) when the fault is at bus 5.

Figure 9 .
Figure 9. Variations in system voltage at bus 8 when connected to: (a) a wind farm with a DFIG; (b) a wind farm with a PMSG; (c) a solar PV farm; (d) a combination of a wind farm with a DFIG and a solar PV farm; and (e) a combination of a wind farm with a DFIG and a solar PV farm.

Figure 10 .
Figure 10.Reactive power support given different RES.

Figure 10 .
Figure 10.Reactive power support given different RES.

Figure 11 .
Figure 11.Profiles of wind speed and solar irradiation.

Figure 11 .
Figure 11.Profiles of wind speed and solar irradiation.

Table 1 .
CCT of generator-2 with conventional and proposed control when connected to different RES.

Table 1 .
CCT of generator-2 with conventional and proposed control when connected to different RES.

Type of RES Connected CCT of δ 21 in ms Conventional Control Proposed Control Improvement
Energies 2022, 15, x FOR PEER REVIEW 14 of 19

Table 2 .
CCT of generator-2 when connected to different RES at different fault instants.

Table 2 .
CCT of generator-2 when connected to different RES at different fault instants.