Comprehensive Dynamic Interaction Studies in Inverter-Penetrated Power Systems

: In a renewable-energy-penetrated power system (RPPS), inverter-based resources (IBRs) pose serious challenges to power system stability due to their completely different dynamic characteristics compared with conventional generators; thus, it is necessary to study the dynamic interactions between IBRs and power systems. Although many research efforts have been dedicated to this topic from both power electronics and power system researchers, some research from the power electronics field treats the external power system as a voltage source with an impedance, therefore ignoring the dynamic characteristics of a power system, while most of the research from the power system field applies simulation-based methods, for which it is difficult to directly interpret the interaction mechanism of IBRs and external system dynamics. Thus, none of these studies can explore the accurate dynamic interaction mechanism between IBRs and power systems, leading to performance degradation of IBR-integrated power systems. Our study takes into account the dynamic characteristics of both IBRs and the external power system, resulting in the development of a new open-loop transfer function for RPPSs. Based on this formulation, it is observed that under certain operating conditions, the dynamic interactions between the inverter and the power system help enhance IBR-penetrated power system stability compared with the case for which the external power system is controlled as a voltage source. The study also reveals how the inverter (phase-locked loop, control parameters, etc.), external power system (network strength) and penetration ratio in an IBR-penetrated power system affect the dynamic interactions between IBRs and the external power system using the proposed quantified interaction indices.


Introduction
Out of concern about the quick depletion of fossil fuels and their negative environmental impacts, the past decade has witnessed a fast increase in renewable energy generations [1][2][3].However, the integration of renewable energy into the existing power grid poses new stability challenges due to the different dynamic characteristics of these IBRs from those of synchronous generators [4,5], which, in turn, often leads to curtailing renewable generation in order to maintain system stability [6], causing the waste of clean, renewable generating resources.
To address these challenges, one of the most urgent research topics for power system planners and operators is to explore and understand the dynamic interaction mechanisms between IBRs and power systems in order to help develop efficient methods to improve power system stability and better accommodate IBRs into existing systems.Meanwhile, a competent understanding of the interaction mechanisms also benefits power electronics engineers by providing them with more precise inverter design guidance, thus allowing them to take full advantage of inverter capability while maintaining the stability of the power system.Although grid-forming inverter (GFMI) technology is currently a prominent research area due to its capability to support grid voltage and frequency, the majority of inverters deployed today function in grid-following inverter (GFLI) mode.Consequently, this paper explores the dynamic interactions between GFLI and the grid to better understand their operational mechanism.
Recent research on the interaction mechanism between IBRs and power systems can be divided into three categories: (1) In the field of power electronics, the external power grid is typically considered as an ideal voltage source with an impedance in renewable energy systems.This is referred to as a voltage source inverter system (VSIS), which takes into account various grid strengths and fully models renewable energy generators, particularly inverters.Based on this model, numerous studies have investigated the influences of inverter control on system stability as well as controller designs, with a focus on the impact analysis arising from phase angle tracking.The authors of [7] proposed a hybrid phase-locked loop (PLL) to improve the inverter's dynamic performance.The authors of [8] investigated the small-signal stability in low-inertia power systems and compared various control strategies through diverse bifurcation studies.The authors of [9][10][11] explored why a PLL deteriorates the stability of the inverter system, the extent of a PLL's influence, and how to improve the stability of a weak power system.While these researchers thoroughly account for the dynamic characteristics of PLLs, overlooking the dynamic characteristics of power systems prevents revealing the dynamic interactions between the IBRs and power systems.(2) In contrast to the aforementioned research, some stability studies on RPPSs adopt more detailed power system models, but they have simplified inverter models.The authors of [12][13][14] investigate how a high penetration of PVs affects a power system's static performance and transient stability.However, simplified inverter models are adopted in these studies.The authors of [15] developed a generic wind generator model for power system stability studies, but the inverter's PLL dynamics are neglected, which is, nevertheless, crucial to the stability of a weak power system.In this case, the incompleteness of the inverter model may lead to the final conclusions deviating from reality.
(3) Among all the renewable-energy-penetrated power system dynamics studies, simulations are the most widely used approach.The major challenge of simulation-based approaches is that analysis and control are heavily related to pre-designed scenarios, making the research conclusions hard to generalize due to the inability to emulate all the possible operational scenarios [5,16].Although simulations can analyze a system's dynamic characteristics at some operating points, a process of extensive and repeated simulations is required to study these characteristics across all crucial operating points.Furthermore, simulation-based approaches also face difficulties with revealing dynamic interaction mechanisms.
This paper employs a dq-impedance-based stability approach in the frequency domain to analyze the interaction dynamics and stability of RPPSs by considering the dynamics of both the IBR and the external power system.The model effectively captures their interactions by integrating the dynamics of both elements into a new transfer function matrix.Subsequently, multivariable control theory [17] is used to elucidate the dynamic interaction mechanisms between the IBR and the power system.In addition, the dynamic interactions are demonstrated through an RPPS under different factors, such as different PLLs, control parameters, SCRs and operating points.Furthermore, we also provide quantified indices to assess the impacts of interactions on stability, which, to our best knowledge, has not been done before.The major contributions of this paper are threefold: (1) introduce a transfer function matrix for the first time that simultaneously contains dynamic information from both the inverter and the power system, (2) reveal the interaction mechanism of an IBR and an external power system through dq-impedance models with the influence of different factors, and (3) develop quantitative indices for the complex interactions between an IBR and a power system.The research outputs of this work are of value not only to power system planners and operators, as the work sheds insights on the dynamic interaction mechanisms between IBRs and power systems, but also to power electronics engineers, as it provides useful guidance for inverter controller design with power system stability being considered.
The rest of the paper is structured as follows.Section 2 provides detailed modeling of the dynamic interactions of the RPPS and deduces the open-loop transfer function matrix.The dynamic interaction mechanism of the system is then explored using multivariable frequency domain control theory in Section 3. The consistency of dynamic interaction characteristics between theoretical analysis and simulation is then verified in Section 4. Section 5 concludes the whole paper.

Modeling Dynamic Interactions between an IBR and a Power System
For the sake of convenience, this paper analyzes a three-phase power system under the synchronous rotating reference frame (SRF), where the variables of the power system are described in the dq-vector form: x abc → x dq := (x d , x q ), while the dq components of the inverter are under the inverter's SRF.Therefore, the state and control variables of the IBR, SG, transmission line and load can be summarized in (1), ( 2) and (3), respectively.A coordinate transformation between the system reference and the inverter reference is required [18].
Figure 1 shows an inverter-connected power system with a detailed control diagram of the GFLI.The transfer function model of the inverter has been widely studied in power electronics research papers and will be cited below while referring interested readers to [19][20][21]: where G PLL is the transfer function of the PLL control loop; V f is the amplitude of the inverter output voltage; I max is the maximum output current of the inverter; and ω p , ω v and ω i are the bandwidths of the active power loop, voltage loop and inner current loops, respectively.Figure 2 shows the interaction model used in our work.The major modeling difference between other works and ours is demonstrated by comparing Figures 1 and 2. Other research works model the external power system as an ideal voltage source with series impedance, while herein we are going to establish a transfer function model of the external power system, including transmission lines, load and generator dynamics.Thus, in their models, v godq and Z g (i.e., Z gequ in Figure 1) are real vectors/matrices, but in our models, they are all transfer functions of the complex frequency (s), with the closed-loop transfer function model of the whole system as: (s) . (5) From ( 5), it can be seen clearly that Y c (s)Z g (s) is the open-loop transfer function of the system and is a perfect candidate to explore the interaction mechanisms between IBR and the external power system, as Y c (s) represents the dynamics of the inverter, while Z g (s) is the dynamics of the external power system.

Transfer Function on External Power System Dynamics
Without losing generality, the third-order SG model is adopted here as follows [22].Undoubtedly, higher-order generator models can likewise be incorporated.
where P m and P e are the mechanical and electromagnetic power, respectively, of the generator.The interface of SG and the network is modeled as follows: Here, we apply a Laplace transform on Equations ( 6) and ( 7); then, combining the results through algebraic manipulation in the frequency domain yields: where the details of G x ′ d , G r and G ∆P m are shown in Appendix A. Furthermore, the transmission lines and loads can be modeled as (9).The nodal current balance constraint in the network is modeled as (10).
Apply again the Laplace transform on Equations ( 9) and ( 10) and then eliminate the algebraic variables i lcdq , i l1dq , i ldq , i lgdq and i l2dq .The transfer function matrix of the external power system dynamics, which takes into account the dynamics of generator and transmission lines as well as loads, can be succinctly formulated as follows: where ∆v gdq , ∆i gdq , ∆v f dq and ∆i cdq are complex variables; the expressions of Y cj , Y gk and Y cjgk can be found in Appendix A.

Modeling Dynamic Interactions
Following the arguments below Equation ( 5), the equivalent impedance of the power system Z g (s) needs to be solved in order to use Y c (s)Z g (s) for modeling the interactions between the IBR and the external dynamic system.For ease of understanding, pure resistive loads are utilized in this study.Nevertheless, it is noteworthy that the methods are also applicable to any other load models.Combining ( 8) with (11) gives: where the details about Z j X j , Y j , B j1 and B j2 can be examined in Appendix A.
Then, the impedance of the external power system can be expressed as: It should be noted that Z g (s) in ( 13) retreats into Z gs of ( 14) when the dynamics of the external power system are neglected.Furthermore, the system's short-circuit ratio (SCR), which is widely used to assess the strength of a power system, is actually SCR = 1 ∥Zgs∥ [19].Therefore, both Z gs and SCR are special cases of the general expression of Z g (s) proposed here.
From Equation ( 14), if the external power system is represented by a voltage source with internal impedance, the impedance can be simplified into a matrix form: where R g = Re[Z gs ] and L g = Im[Z gs ].

Dynamic Interactions Analysis
By applying the dynamic interaction model proposed in Section 2, this section analyzes the impact of various factors (SCR, PLL, control parameters and penetration ratio) on the stability of the entire system resulting from the dynamic interaction between the IBR and the external system.Case study parameters are shown in Table 1 [19].Since the transfer function models of the interaction have been established, Bode plots become a handy instrument to study how various control, network and operational factors influence the dynamic interactions between the IBR and the external power system.While the complicated interactions in the time domain can become much clearer when they are observed in the frequency domain, it is still challenging to quantify the interaction level.This issue becomes severe when considering a multi-input-multi-output (MIMO) dynamic system, as the proposed transfer function model matrix reflects the input/output characteristics of the entire system, while Bode plots can only reflect the channel-to-channel frequency characteristics.
To further quantify the interaction, the generalized Nyquist criterion (GNC) [21,23] is utilized by extending the classical Nyquist diagram into MIMO dynamic systems.It should be noted here that although GNC has been widely used in power electronics studies, it is rare to see its applications in power systems.Based on (4), ( 13) and ( 15), the generalized Nyquist diagram of the open-loop transfer function Y c Z g and Y c Z gequ is shown in Figure 3 for a weak power system with SCR = 1, where λ 1 and λ 2 are the eigenvalues of the respective transfer function matrix.The term (margin) can then be defined as 1 + L, with L being the value of the leftmost point in the generalized Nyquist diagram.The motivation to use this margin to reflect the dynamic interaction between the IBR and the system is two-fold: (1) it is a quantified value with clear and important physical meaning, as it measures the distance of the whole dynamic system towards being unstable; (2) it is an immediate consequence of the interactions between various IBR (Y c (s)) and external system dynamics (Z g (s)).This section will utilize the mentioned methods to analyze the dynamic interactions between the IBR and the external system.The analysis will be conducted under various factors such as grid strength, PLL, controller parameters and penetration ratio.To clearly and concisely describe system characteristics, henceforth we use pllmode = 1 and 2 to indicate the classical PLL controller design and the ideal PLL (which can always achieve phase tracking), respectively, and use statemode = 1 and 0 to indicate RPPS and VSIS, respectively.

RPPS with Different SCRs
SCR is a widely used grid strength index.Calculating SCR does not consider the system dynamics and is based on the grid impedance Z gequ in (14).Typically, the power system is considered strong when SCR ≥ 3 and weak when SCR < 3 [24].The Bode plots of Z gequ and Z g at SCR = 1 (weak system) and 3 (strong system) are shown in Figure 4. Here, to reflect that both Y c and Z g are matrices, the following symbols are used: }, where DD and QQ channels are decoupling channels, while DQ and QD channels are coupling channels.
Figure 4 shows that when the angular frequency ω > 2 rad/s, the amplitudes of Z g3(DD) and Z g3(QQ) are smaller than those of Z g1(DD) and Z g1(QQ) , respectively.This results in a reduction in the amplitudes of all channels of Y c Z g3 , which can help to decrease the error in the control variable (∆i codq ) and output variable (∆i cdq ), as per Equation (5).Then, when ω < 15, 000 rad/s, the amplitudes of Z g3(DQ) and Z g3(QD) are smaller than those of Z g1(DQ) and Z g3(QD) , respectively.This means that in a stronger power system, there is smaller dq-axis coupling, which is important for power electronics engineers for controller design as it reduces the requests for compensation in the dq-coupling in controller design and facilitates achieving the desired control objectives in each axis.Additionally, when systems with different SCRs have the same Y c , the difference between Y c Z g1(DD) (Y c Z g(DD) =Y c(DD) Z g(DD) +Y c(DQ) Z g(QD) ) and Y c Z g3(DD) mainly comes from the difference between Z g1(DD) and Z g3(DD) because the coupling terms have less effect.Specifically, at a frequency of approximately 1770 rad/s, the magnitude of Y c Z g3(QQ) is much smaller than that of Y c Z g1(QQ) , and both of them are below 0 dB, indicating that Y c Z g3 has greater stability margins.This is in line with the results shown in Figure 5, which shows how the proposed interaction index changes with respect to changes in the grid strength.Frequency (rad/s) label Therefore, a system with higher SCR has lower coupling between the d-axis and q-axis, lower control error and higher interaction index.On the other hand, for power electronics engineers, a power system with a higher SCR also means the controller is more easily designed and more robust.

RPPS vs. VSIS
As shown in Figure 4, in the middle-frequency range, under the same SCR, there is a noticeable similarity in the frequency responses of both Z g and Z gequ in all channels.However, over almost the entire frequency range, Z g(DD) and Z g(QQ) have a lower magnitude than Z gequ(DD) and Z gequ(QQ) , respectively.This indicates that the system has a lower control error regarding any changes in the input if the system's dynamic is considered compared with the VSIS model.In the high-frequency range, the magnitudes of all channels of Z gequ are higher than those of Z g , suggesting that the system is more prone to high-frequency oscillation in VSIS.
In the middle-frequency domain, Bode plots of Y c Z g and Y c Z gequ have similar shapes in all channels, which means the two systems have similar dynamics within this frequency range.Moreover, Figure 4 shows that, compared to Y c Z gequ , the magnitudes of all channels of Y c Z g attenuate extremely fast at high frequencies, indicating that there is much lower possibility of high-frequency oscillations in RPPS.Y c Z g(QQ) and Y c Z gequ(QQ) have almost the same phase Bode plots when the external power systems have the same SCR.However, when ω = 1770 rad/s (268 Hz), Y c Z gequ1(QQ) has almost no stability margin, while Y c Z g1(QQ) still has a stability margin.This indicates that RPPS is more stable than VSIS under these conditions as dynamic interactions between IBR and the power system are taken into account, which is consistent with Figure 3.
Figure 5 shows that the dynamic interaction index increases with increasing SCR.It can be seen clearly for a very weak power system that the interactions between IBR and the external system have discernible differences for different external system dynamics (RPPS and VSIS).The increasing SCR can help reduce the impact of external system dynamics on its interactions with IBR.It can also be observed that RPPS always lies above VSIS, indicating that the interactions between the external system dynamics and the IBR actually improve the system stability level.This can be attributed to the frequency spectrum dynamic performances introduced by the coherent inertia of the generators, as can be observed through the Bode plot studies in Figure 4.

RPPS with different PLLs
Figure 6 demonstrates how the PLL influences the characteristics of the IBR and then affects the interaction between IBRs and the external power system when SCR = 1.Y ca f (a new Y c when the PLL of the IBR is ideal) still has the same characteristics as Y c in DD and QD channels, but when 3 rad/s < ω < 2000 rad/s, the magnitude of Y ca f (DQ) is smaller than that of Y c(DQ) , which can reduce the coupling between the d-axis and q-axis of Y ca f .When ω < 300 rad/s, the magnitude of Y ca f (QQ) is smaller than that of Y c(QQ) , which is beneficial for Y ca f to improve system stability.When ω = 1770 rad/s and the phase equals −180 degrees, the magnitude margin of Y c Z g(QQ) is very small, which can potentially cause the system to become unstable.However, at a frequency of 300 rad/s and a phase angle of −180 degrees, Y ca f Z g(QQ) still has a significant magnitude margin.
Several researchers have examined the effects of PLLs on the system to improve PLL performance to increase system stability in weak power systems [7,9,10,19].Different from such studies, this paper provides an intuitive demonstration of the PLL's impacts on the system dynamics using the proposed interaction index, which can provide guidance for power electronic engineers to enhance PLL design.Figure 7 shows how the proposed interaction index changes under the influence of the PLL and the change in I cq , which represents different operating points.When pllmode = 2, the system has an improved interaction index, demonstrating that PLL has a negative impact on the interactions between IBR dynamics and the external power system.

RPPS vs. VSIS
Based on Figure 7, it is evident that the interaction indices of RPPS and VSIS remain similar across different operating points when utilizing an ideal PLL.For the system with a traditional PLL, the index of RPPS is above that of VSIS for most operation points, which demonstrates that the interactions between external system dynamics and IBR dynamics increase the whole system's stability level, which is consistent with the observations from Figure 6.This also indicates that the dynamic interactions between the IBR and the external system hinge heavily on the interactions between the PLL and the external dynamics.

Control Parameters and Dynamic Interaction
Selecting the bandwidths for different control blocks in an IBR is critical in inverter controller design.Five bandwidths need to be decided in Figure 1: namely, ω i , ω v , ω p , ω pll and ω LPF , corresponding to the bandwidths of inner current loop, voltage loop, active power control loop, PLL controller and the low-pass filter, respectively.Among them, ω i and ω pll have the most significant influence on the interaction.Based on Figure 8, varying combinations of ω i and ω pll result in different frequency responses for Y c and thus different effects on Y c Z g .Specifically, (1) Y c1 and Y c4 , as well as Y c2 and Y c3 , have the same ω i and exhibit the same frequency response in the DD and QD channels, indicating that ω pll has no effect on these two channels; (2) Y c1 and Y c2 , as well as Y c3 and Y c4 , have the same ω pll and the same frequency response in the DQ and QQ channels, suggesting that ω i has little effect on these two channels; (3) optimal stability occurs at ω i = 3000 and ω pll = 50, where the phase angle equals −180 degrees and the amplitude margin of Y c3 Z g(QQ) is the largest, which is consistent with the peak point in Figure 9.
The proposed index can be used to evaluate the dynamic sensitivity of the system stability with respect to these key control parameters.Figure 9 illustrates how the index changes regarding the different combinations of ω i and ω pll .The figure of contours shows the region where the choice of ω i and ω pll can make the whole system unstable as well as how to tune the key parameters of the IBR to increase the whole system's stability by leveraging the interaction between the IBR and the external system dynamics.Both, we believe, are very beneficial to power electronics engineers for IBR controller design.Therefore, this paper can form a handy tool for inverter controller designers for finding reasonable parameters with system stability being considered.It should also be noted that from Figure 9b, in the region where ω i > 1200 rad/s and ω pll > 50 rad/s, the dynamic interaction between the IBR and the generator dynamics contributes to improvement in the stability level compared to that between the IBR and the static external system.

Penetration Ratio and Dynamic Interaction
The ratio between the IBR's real power output and the total generation in the whole system is widely used as a measure of the IBR's penetration level.As shown in Figure 10, when the penetration ratio equals 0.25 or 0.75 and the angular frequency ω > 100, Y c1 and Y c2 , as well as Z g1 and Z g2 , have very similar frequency responses in DD, DQ and QD channels.However, the magnitude of Y c1(QQ) is smaller than that of Y c2(QQ) , which leads to the magnitude of Y c1 Z g1(QQ) being much lower than 0 dB compared to Y c2 Z g1(QQ) when the phase angle equals −180 degrees.Therefore, a lower penetration ratio can result in better stability for the RPPS.
From Figure 11, it can be seen that the interaction index decreases with the increase in the penetration ratio regardless of the power system model adopted.As this interaction index aligns with the stability level of the whole system, the decreasing index with the increase in IBR penetration indicates that the GFLI reduces the stability level of the system, and its negative impacts increase along with the penetration level.
Furthermore, when using the VSIS model, the interaction index declines more quickly with respect to the increase in IBR penetration.Specifically, the critical value of the penetration ratio using the VSIS model is approximately 0.6, whereas, with the RPPS model, it is approximately 0.8, which means that the interactions between the IBR and VSIS have more severe adverse impacts on the stability of the system.This confirms the earlier observation that interaction between the IBR and generator dynamics might improve the system stability under certain operational scenarios as outlined in Section 3.1.2.label

Simulation Results
Numeric simulation studies are carried out to validate the developed IBR-integrated power system model and the correctness of the proposed index.Detailed Simulink models corresponding to the analytical models in Section 2 are developed.The steady-state values are shown in Table 2, and the controller parameters of the IBR are identical to those in Table 1.Based on the system's designed parameters and the operation point, the interaction index of RPPS is 0.22, indicating that the system is stable.However, the interaction index of VSIS is around 0.077, which is at the edge of stability.The trajectories of IBR voltage and active and reactive power are shown in Figure 12 for two different dynamically interacting systems in the paper.The disturbance is created by adding a 3 MW load at the inverter terminal bus between 5 s and 5.1 s.The simulation results fully validate the theoretical analysis results: (i) the RPPS shown in Figure 12a remains stable after experiencing the numerical disturbances; (ii) the VSIS shown in Figure 12b, although still stable initially due to the positive interaction index, finally loses its stability because of the disturbance, which confirms it is at the edge of stability.Compared with the frequency domain analysis in the previous sections, it should be noted that the time domain simulation reveals whether a system is stable or unstable, but it cannot clearly provide the margin to instability and unveil how the IBR interacts with the external system dynamics to affect system stability.As described in Section 3.1, the interaction index of RPPS is 0.22 when SCR = 1, whereas the interaction index is 0.76 when SCR = 3.This means the RPPS with SCR = 3 has a larger stability margin and very strong anti-interference capability owing to its larger interaction index compared to the RPPS with SCR = 1.To verify that systems with various SCRs have different interaction indices, we build two other simulations for which these two systems with different SCRs remain stable before 5 s and then suffer from a disturbance, as shown in Figure 13.The disturbance is created by adding a 5 MW load at the terminal bus of the inverter between 5 s and 5.1 s.The simulation results show: (i) the RPPS with SCR = 1 shown in Figure 13a loses stability after the disturbance; (ii) the RPPS with SCR = 3 shown in Figure 12b remains stable after the disturbance, which is attributable to its higher positive interaction index.
Combining Figure 12 and Figure 13, it can be observed that the RPPS with SCR = 1 can withstand a 3 MW disturbance but loses stability under a 5 MW disturbance.However, the RPPS with SCR = 3 is capable of maintaining stability when suffering from the 5 MW disturbance, which is a result of its larger interaction index compared to that of SCR = 1.Therefore, the interaction index provides an intuitive way to reflect all factors that impact system interactions.

Conclusions
This paper derives a new transfer function matrix of a system by considering both the dynamic characteristics of the IBR and the external power system simultaneously for the first time.Not only can the respective dynamic characteristics of the IBR and the external system be revealed by analyzing the newly developed transfer function matrix from a frequency domain perspective, but we can also determine how they interact with each other.This paper proposes a numerical index to quantify the interaction impacts of inverter-based resources (IBRs) on overall system stability.Through a detailed examination of the dynamics between the IBR and two distinct external power system models, we find significant differences in stability outcomes, highlighting the critical need for in-depth studies of these interactions.Furthermore, the study conducts an extensive sensitivity analysis of the system's stability in relation to network strength, control components, controller parameters and penetration levels.These insights are invaluable for power system planners, operators and inverter design engineers and offer targeted guidance for enhancing system robustness.Extending research to a grid-forming inverter and a multi-machine power system is undergoing.
Funding: This research received no external funding.

Figure 3 .
Figure 3. Nyquist plot of Y c Z gequ and Y c Z g with SCR = 1.(a) Nyquist plot of Y c Z gequ and (b) Nyquist plot of Y c Z g .

Figure 4 .
Figure 4. Bode plots of Z g , Z gequ , Y c Z g and Y c Z gequ with different SCRs.The numbers in the subscript represent the values of SCRs.

Figure 5 .
Figure 5. Margins of RPPS and VSIS with different SCRs.The margin value comes from the GNC of the system's open-loop transfer function.

Figure 6 .Figure 7 .
Figure 6.Bode plots of Y c , Y ca f , Y c Z g , Y ca f Z g under the effect of PLL.

Figure 9 .
Figure 9. Impact of controller parameters (ω i and ω pll ) on the stability of the system in RPPS and VSIS.In (b), the solid line represents RPPS and the dashed line represents VSIS.(a) Change in index with respect to current loop bandwidth and PLL bandwidth.(b) Contours of margin.

Table 1 .
Parameters of the inverter.

Table 2 .
Design parameters of the simulation.