Large-Signal Stability Analysis of Inverter-Based AC Microgrids: A Critical and Analytical Review

With the ever-increasing development of microgrids to improve power supply reliability and resilience, large-signal stability analysis of microgrids is needed to provide accurate insights into system dynamics during significant disturbances and to offer a robust and comprehensive assessment of microgrids’ ability to return to normal operation after contingencies. This paper presents a comprehensive analytical and critical review of large-signal stability analysis methods for inverter-based AC microgrids supported by both theoretical analyses and numerical case studies. These analyses and case studies are developed and implemented to analyze, compare, and identify the gaps in existing large-signal stability assessment methods, which have not been investigated in the literature. More specifically, the recent developments in large-signal stability analysis techniques for microgrids, including Lyapunov-based methods and energy function analysis, are reviewed and discussed. Also, as the accuracy of the dynamic models of inverter-based resources (IBRs) is a crucial factor for authentic and accurate stability assessment, the impacts of full-, reduced-, and second-order dynamic models of IBRs on the accuracy of large-signal stability assessment in microgrids are numerically scrutinized. The study concludes by identifying the applicability of existing stability analysis methods for microgrids (e.g., Krasovskii’s, Popov-Lure, and sum of squares (SOS)-based methods) and presenting their challenges for further investigation in the field of large-signal stability of microgrids using numerical case studies. This comparative assessment of large-signal stability assessment provides an informative analysis of the system’s capability to operate under contingencies and helps to set the protection system efficiently and prevent unnecessary outages and trips. The numerical assessment shows that SOS-based stability assessment can provide a more realistic and less conservative stability region. It is worth mentioning that other approaches are computationally more efficient and can be applied for online applications and control objectives.


I. INTRODUCTION
Microgrids are more vulnerable to becoming unstable due to faults and large disturbances than larger and interconnected power grids because microgrids are small-scale and localized power systems with low or no inertia distributed energy resources (DERs).Evaluating the robustness of microgrids The associate editor coordinating the review of this manuscript and approving it for publication was Zhilei Yao .
in the face of significant disturbances, the process referred to as large-signal stability assessment, helps to understand the dynamic behavior of DERs, ensures stable operation during islanded mode and generation shortages, addresses nonlinearities and control interactions, and ensures transient and steady-state stability.
The importance and distinction of large-signal stability analysis for microgrids can be attributed to their unique characteristics and operational requirements-DERs in microgrids exhibit fast dynamic responses due to their electronic converters and power electronics interfaces [1], [2], [3].Their response to large disturbances, such as sudden load changes or faults, can significantly impact the overall system stability.Moreover, microgrids can operate in an islanded mode when they are disconnected from the main utility grid.In this mode, the microgrid must maintain stable frequency and voltage without external support [4].Additionally, microgrids involve complex interactions between various control loops and power electronics converters.Nonlinearities, such as voltage and current limits and control interactions, can affect the system's stability.Furthermore, the high penetration of power electronics-based DERs in microgrids can introduce harmonic distortion and resonance phenomena [5], [6], [7], [8].Harmonics can affect system stability and lead to instability if not properly addressed.
Large-signal stability analysis considers both transient and steady-state stability aspects [5], [7].It assesses the system's response to large disturbances, such as faults or sudden changes in generation or load, and examines whether the system can recover and reach a stable operating point.This analysis ensures the robustness and reliability of microgrid operation under different operating conditions.
Several studies evaluate the stability of inverter-based DERs using various stability analysis methods and dynamic models of inverter-based resources (IBR).However, there are only a few comprehensive review studies that focus on the large-signal stability assessment of microgrids and summarize existing methods and approaches in a single study.Current review papers do not provide a comparative and analytical assessment of existing methods and approaches that are needed to examine the accuracy and authenticity of the large-signal stability analysis of microgrids.
In [9], authors provide a classification of the stability analysis based on the disturbance and duration, list the important features of microgrid and inverter-interfaced networks that must be considered in stability analysis, and process the factors that impact the large-signal stability of a converter-based microgrid.Moreover, the paper has tabulated the stability analysis tool and techniques.The work in [10] investigates the existing valuable studies about large-signal stability analysis on the distribution systems and IBRs as well as applications of Lyapunov-based stability assessment in detail.
The work in [11] presents one of the well-known reviews in the area of large-signal stability assessment.The authors have investigated the large-signal stability of the DC power system that includes motor drives and constant power load.Different direct methods for large-signal stability analysis of DC power systems are processed, including Takagi-Sugeno (T-S) multi-modeling, Brayton-Moser's mixed potential, block diagonalized quadratic Lyapunov function, and reverse trajectory tracking.The application of the T-S fuzzy modeling as one of the direct methods for stability analysis of a system with a DC-link connection is also proposed.Furthermore, the work in [11] proves the performance of the proposed method by providing numerical analysis and experimental results.The case study is an electric motor connected to a DC power source via an inadequately damped LC filter that is examined under the listed tools.
In [12], a comprehensive study on all aspects of stability challenges in microgrids is provided.The paper focuses on techniques for small-signal stability assessment and discusses the important contents of the stability assessment of microgrids under both small and large disturbances.Review studies [13] and [14] provide a review of the voltage stability of microgrids and categorize the voltage stability based on the operation modes (grid-connected or islanded), disturbance intensity, and controlling system.In [15], the stability of the microgrid is assorted based on the fault type, disturbance duration, and different energy resources.Also, the stability challenges of microgrids for both grid-connected and islanded modes are discussed.The work in [16] classifies the models and approaches of the stability analysis based on the large-and small-signal analysis.Moreover, the authors review the advantages and disadvantages of the models and approaches.
Although the aforementioned studies provide valuable reviews of existing research on the large-signal stability assessment of microgrids and IBRs, none of them analytically provides comparative scrutiny with numerical assessment and proofs in the large-signal assessment of AC microgrids.In other words, a critical review needs to compare and analyze different stability analysis methods and approaches to identify the accuracy and integrity of microgrid stability analysis.A numerical, analytical study and review of large-signal stability assessment of microgrids can help to determine the advantages and disadvantages of existing methods and approaches as well as introduce the best possible applications for each large-signal stability technique.
The compelling benefits, stability concerns associated with microgrids, and the research gap for a comparative investigation of different stability assessment techniques, as outlined previously, have inspired our concentration on their stability assessments.This involves a thorough exploration of pertinent definitions and a comprehensive, analytical scrutiny of the existing scholarly literature on the topic.
In this study, a comparative perusal of the large-signal stability assessment of AC microgrids and IBRs is provided such that the system's large-signal stability analysis is numerically scrutinized based on the different dynamic models and different mathematical methods for stability analysis.Based on the obtained results, the advantages and disadvantages of each method and dynamic model are discussed.This paper not only reviews and investigates the exciting techniques and methods for the microgrid stability assessment but also presents a numerical, comparative assessment of the large-signal stability analysis of IBRs, which has not been addressed in the other large-signal stability assessment surveys [9], [10], [13], [14].In this context, our work is focused on assessing the large-signal stability of the system without the implementation of any stabilization techniques.The study objective is to identify crucial parameters from a system stability standpoint.The contributions of this study are summarized as follows: • Provide a comprehensive literature review on the large-signal stability assessment of AC microgrids and classify existing stability analyses based on applied methods, microgrid dynamic models, and system configurations and characteristics.
• Propose a detailed mathematical dynamic state-space model and modify stability assessment methods for Lyapunov-based stability analysis accordingly.
• Conduct numerical assessments and comparative analytical review using the Lyapunov-based stability assessment including Krasovskii's, Popov-Lure, and sum of squares (SOS) method.
• Evaluate microgrid dynamic models including timedomain differential equations and synchronverter models in terms of accuracy, scalability, and computational time.The rest of the paper is organized as follows (classified in Figure 1).Section II presents the challenges of stability analysis in a microgrid.In this Section, classifications of small-and large-signal stability assessment are expanded and discussed, and the concept of the domain of attraction (DOA) and critical clearing time (CCT) is elaborated in detail.Section III represents the existing dynamic model of inverter-based DERs.In Section IV, the techniques and methods for stability analysis are developed and discussed.Section V develops a case study and numerically expatiates the approaches of existing stability analysis.The numerical analysis includes modeling, simulation, analysis, and comparison of the results.In Section VI, a discussion on the advantages and disadvantages of the simulated approaches is provided and the conclusion is presented.

II. STABILITY CHALLENGES IN MICROGRIDS
The stability analysis of microgrids is generally classified into two main categories (small-and large-signal stability analyses), which are explained in the following sections in detail.The small-signal stability is defined as steady-state stability for small disturbances, and the large-signal stability is defined as transient stability for large disturbances.Also, the large-signal stability assessment can present the microgrids' dynamic behavior during small disturbances.Although the main goal of this paper is to discuss and analyze the stability of AC microgrids to large disturbances, we briefly present the small-signal stability analysis for the sake of comparison and emphasize the importance of large-signal stability analysis for microgrids.

A. STEADY-STATE (SMALL-SIGNAL) STABILITY STUDY
Small-signal stability analysis of a microgrid involves studying the stability of the system under small disturbances and deviations from its steady-state operating conditions.The small-signal stability studies of inverter-based microgrids have been mainly conducted based on eigenvalue analysis and impedance-based approach [17], [18].
Using the eigenvalue technique, a state-space model of the system is linearized around the steady-state operating point to obtain the eigenvalues.The eigenvalues are then used as a stability criterion.For instance, when all eigenvalues of a continuous state-space model have negative real parts, the system is called stable.However, having at least one eigenvalue with a positive real part results in system instability.For instance, the eigenvalue assessment is applied to a system dynamic model to analyze the impact of changing controller parameters in inverters [19], filters, line, and load specifications on the system stability.
In addition, the sensitivity analysis of eigenvalues with respect to parameter changes holds significant importance within this field.Thus, the parameters that may affect the stability are realized and can be used as a viable solution for designing optimal controllers in the system.In [20], a novel Laplacian matrix eigenvalue is introduced to build a certified stability region.The computation burden of the proposed method is almost independent of the inverters' penetration level in the distribution grid.In [21], a smallsignal stability criterion based on the extended Gershgorin Theorem is formulated by reshaping the range of eigenvalue estimates.The test results suggest that the proposed method possesses the capability to evaluate the stability of high-order systems with minimal computational complexity.However, the method is applied for DC microgrids, not AC grids.
In [22], the stability analysis model of an islanded microgrid based on a linearized state-space model is broadened to encompass the dynamics of a smart load alongside converter-interfaced distributed generators.
Another technique used in small-signal studies is impedance-based analysis [23], [24], where the load and sources are represented by their input and output (I/O) impedance.This method evaluates the whole system instead of details of inner subsystems.It applies the Nyquist criterion to the ratio of output and input impedance to assess stability.The average linearization model is usually adopted for this analysis.
In [25], an efficient technique based on the perturbation signal quadratic residue binary (QRB) sequence is proposed for impedance measurement and stability improvement of microgrids with an adaptive control algorithm.The study given in [26] proposes a numerical-based technique to obtain impedance models of power-electronic dominated grids.The proposed technique using the finite-difference method shows accurate results compared with perturbation and analytical methods.In [27], authors address the small-signal assessment of AC islanded microgrid and show that active filtering function remarkably affects current-controlled converters (CCC), phase-locked loop (PLL), and the dynamic behavior of microgrid.They also offer useful guidelines for designing the PLL controller for multi-functional inverters.Microgrid modeling has also been an important topic in small-signal analysis.In [28], the limitation of different reduced-order models with virtual impedance is studied by evaluating the computation time, the parametric stability boundary, eigenvalues, and the relative error of the natural frequency and damping factor.In [29], the reduced-order model of the microgrid's distributed generator is developed, including virtual impedance and internal control.The authors also demonstrate that using quasi-stationary virtual impedance led to a larger stability region compared to the transient virtual impedances.
Small-signal stability assessment scrutinizes the stability of microgrids under small changes, which makes the results valid only around system operating points.Since microgrids are prone to both large and small disturbances, the results of small-signal analysis cannot be generalized.In contrast, large-signal stability assessment uses nonlinear models, enabling it to comprehensively evaluate system stability under both small and large disturbances.However, this can add computational complexity to the stability studies.Critically discussing the advantages and disadvantages of existing large-signal stability analysis methods and their associated computational challenges using theoretical and numerical analyses is the main motivation of this paper.

B. TRANSIENT (LARGE-SIGNAL) STABILITY ANALYSIS
Large-signal stability analysis for transmission power systems has been extensively studied in the literature [30], [31], and several methods have been commonly adopted by power companies.On the other hand, large-signal stability analysis for microgrids is still evolving especially with the high penetration of IBRs, where each method is applicable under certain system conditions.Applying the methods proposed for the stability assessment of the conventional power grids cannot present comprehensive and authentic results for AC microgrids that have complicated structures such as IBRs.The main differences between microgrids and conventional power grids are itemized as follows: • The microgrid includes inverters and other power electronics devices, making their dynamic response faster than that in transmission power systems equipped with mechanical components used in the generators [32].
• The complexity and nonlinearity of the dynamic model of the microgrids are elevated by the utilization of power electronic devices, rendering them arduous to investigate the stability [33].
• The IBRs in microgrids have lower inertia (or no inertia) compared to the synchronous generators and electrical machines in transmission grids [34], decreasing the system's stability and increasing its susceptibility to disturbances.
• The R/X ratio of lines in distribution grids and microgrids is larger than that in transmission systems, which can introduce convergence problems for the numerical solutions.
• Microgrids are considered small energy systems, and any changes and disturbances in microgrids can have an extreme effect on the whole system's operation.
• Microgrids are susceptible to experiencing unbalanced loads and power flow, leading to adding more complexity and nonlinearity to dynamic model [35], [36].
Consequently, the large-signal stability analysis for the inverter-based AC microgrid is inevitably complicated and needs close and extensive investigation.As the large-signal stability analysis of microgrids based on the time-domain simulation is an extremely time-consuming process, recently, research studies have focused on the mathematical and theoretical methods to appraise the large-signal stability in a computationally tractable manner.
Generally, the large-signal stability analysis can be classified based on the system dynamic model (discussed in Section III) and stability assessment method (discussed in Section IV), as shown in Figure 2.
In the realm of dynamic model-based classification for inverter-based microgrids, two main categories exist: the full-or reduced-order dynamic model and the synchronverter (second-order) model.The full-or reduced-order dynamic model is defined based on the time-domain differential equations of IBRs and circuit laws, including both controlling systems and circuits (discussed in Section V-B).For the synchronverter (second-order) model proposed by [37] and Responsibility of each control level in hierarchical control.[38], it has been assumed that an IBR mimics the dynamic behavior of a synchronous generator (SG).
Choosing an appropriate method for stability analysis as well as an accurate dynamic model for large-signal stability analysis of inverter-based microgrids are important factors, which vary depending on the application, microgrid configuration, and case study.For instance, in online applications where both fast and accurate responses are required, it is essential to utilize a dynamic model that is simple with reasonable approximations.Additionally, the selected stability analysis method must be fast and computationally tractable to meet real-time requirements.
On the other hand, for conducting a post-fault analysis to adjust the protection system and relay settings, the time-consuming computational process imposed by using an accurate and detailed dynamic model and the stability analysis method is not considered as an obstacle.
There are a few research studies addressing the complexity of the large-signal stability assessment of microgrids tabulated in Table 1.The applied method and investigated dynamic model are highlighted for each study.Moreover, there are a few survey studies on large-signal stability assessment of microgrids listed in Table 2. Generally, the large-signal stability assessment for microgrids has been conducted for a single IBR connected to an infinite bus.However, there are a few perusals that address the challenges of the transient stability of microgrids that are constructed with more than one IBR, such as [43], [50], and [54] (detailed in Table 1).It is worth noting that the main focus of this paper is on the theoretical and analytical large-signal stability analysis of AC microgrids rather than DC microgrids.However, there are a variety of studies on DC microgrids and the challenges of large-signal stability of the DC network [52], [55], [56], [57], [58], [59], [60], [61], [62], [63], [64], [65], [66], [67], [68].

C. CRITICAL CLEARING TIME AND DOMAIN OF ATTRACTION
The stability analysis has been carried out for two main purposes: determining critical clearing time (CCT) for temporary disturbances and the feasibility assessment of stable transitions for permanent disturbances [79], [83], [84], [85], [86], [87], [88], [89], [90], [91], [92], [93], [94].The CCT refers to the time duration following the occurrence of a contingency or disturbance in the power system.It is necessary to clear the contingency and restore the system to its pre-fault state within this duration in order to maintain stable operation.In other words, CCT represents the maximum duration during which a power system can tolerate a disturbance and remain stable without experiencing a system-wide collapse [95], [96], [97].
The CCT is an important factor in power system planning and operation since it serves as a crucial parameter by providing system operators with a time limit to take necessary corrective actions during events.It is also used as a criterion for designing protective relays and determining the necessary backup protection actions [54], [98], [99], [100], [101].This time is calculated based on the system parameters and network topology considering the fault location and type.Typically, precise CCTs for each contingency can be determined by conducting a suitable stability analysis of the grid, utilizing either time-domain simulation or theoretical evaluation.
In the case of temporary disturbances, the main objective of the stability analysis is finding the possible maximum time to clear the fault or contingencies without any blackouts or losing system operation-i.e., the fault or disturbance must be cleared before the determined clearing time duration to guarantee a stable operation.To determine the CCT using large-signal stability assessment [102], [103], [104], [105], [106], the evaluation and determination of the system parameters' critical boundaries and limits are undertaken to ensure system stability.The driven boundaries help in setting the protection system in a more realistic and less conservative manner to guarantee a stable operation and prevent unnecessary trips.As a consequential result, this approach allows for optimal utilization of the majority of the system's capabilities during transient faults and contingencies.
The other application for applying the large-signal stability assessment is for permanent and large disturbances or changes in the system.The assessment of large-signal stability during permanent disturbance is to determine whether the system can reach its new steady-state operational/equilibrium point (EP) after the system changes.Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.To ensure a stable transition regardless of permanent or temporary events, the disturbance trajectory must be inside the stability region of the system around the EP, representing the domain of attraction (DOA)-i.e., the domain of attraction, alternatively identified as the stability region, encompasses the set of all initial states that eventually converge to a specific equilibrium point of a system.Essentially, it demarcates the range of initial conditions that give rise to a stable state within the system.Each operational point or equilibrium point is confined to a specific domain of attraction [107], [108].For each DOA, a critical energy boundary is defined as determining the unstable boundary for any initial point and disturbance trajectories such that if the disturbance trajectory passes the critical energy boundary, the system cannot return to its original steady-state operation point or equilibrium point and experiences an unstable scenario or blackout [95], [96].
For a temporary disturbance, the critical energy boundary is exploited to identify CCT [109].For a permanent disturbance, if the pre-disturbance equilibrium point is not inside of the DOA of the new equilibrium point, the system experiences instability during the transient response.In other words, the previous operating point works as an initial point for the new equilibrium point to the steady-state operating point.According to the DOA definition [107], [108], if this initial point (previous steady-state operating point) is not inside of the DOA of the new steady-state operating point, trajectories will not converge to the new equilibrium point of the system, and cannot achieve a stable state within the system.
In this regard, for a specific transient stability assessment for permanent disturbances and changes in the system, the previous equilibrium point can be considered as a set of initial points and new trajectories for the new equilibrium point that is cleared at the time of starting permanent change or contingency.If the previous steady-state operation point is inside of the DOA of the new equilibrium point, it means that the energy of the system at clearing time is not more than the critical energy boundary of the new equilibrium point, and the system will successfully transfer to the new steady-state operational point.
For more clarification, consider a new equilibrium point in the middle of the light-green area shown in Figure 3; since the pre-disturbance equilibrium points B and C are outside of the new equilibrium point's DOA, their trajectories diverge from the new equilibrium point.On the other hand, the equilibrium point A is located inside of the new equilibrium point's DOA and has a successful transition to the new steady-state operational point (or new EP).
Generally, there are two ways to obtain the domain of attraction of a specific equilibrium point in a system, including time-domain simulation and direct method (e.g., Lyapunov-based methods) [46].As mentioned before, although the utilization of time-domain simulation is meticulous and authentic, it poses challenges in terms of computational burden and process complexity.On the other hand, direct methods offer not only time and cost efficiency but also achieve acceptable accuracy levels when compared to the accuracy of time-domain simulations [10], [110].
This study focuses on the large-signal stability analysis of inverter-based microgrids based on the direct methods and Lyapunov-based stability approaches.In the following sections, the challenges of large-signal stability assessment of IBRs are analytically discussed, and the studies in this area are numerically reviewed in detail.The numerical analysis and comparative study are presented for the existing techniques and methods.

III. MICROGRID COMPONENT DYNAMICS AND MODELING
The primary challenges associated with utilizing dynamic models for IBRs lie in achieving a balance between acceptable model accuracy and maintaining low computational costs during the stability analysis process.This is particularly important for online applications where the system needs to respond quickly to changes in system conditions.To address this problem, the reduced-order dynamic models of IBRs are introduced.They are designed to be computationally efficient and have reasonable approximations.This can be achieved by simplifying the model and reducing the number of states in the state-space model of the system [111], [112].
The reduced-order dynamic models of inverter-based DERs can provide a reasonable approximation of the system's dynamic behavior, satisfying the accuracy of the model.It is worth noting that the accuracy of the reduce-order dynamic models depends on the level of detail included in these models and the complexity of the system being analyzed.Therefore, to ensure the high accuracy of the dynamic model, it may be necessary to include more states and details in the state-space model of the system dynamics [46], [113].
The full-order dynamic model can be modified to a reduced-order model based on the approximation and desired outputs and achieved different numbers of order, including 11th, 9th, 7th, 5th, 3rd, and 1st order.The details of the reduced-order models can be found in [46].The full-order model is described and discussed in Section V-B.
The virtual synchronous generator (VSG) model can be deemed as a reduced-order dynamic model of IBRs.According to the VSG model, IBR's dynamics emulate the dynamic behavior of a synchronous generator in a power system.In fact, the VSG model is a second-order dynamic model of IBRs and is simple to apply for any online control application.The details of the VSG model are discussed in Section V-A.
Depending on the application and intended results, an appropriate dynamic model for inverter-based microgrids must be chosen.Nevertheless, the dynamic model is not the sole crucial element in the large-signal stability analysis of microgrids.The methodology for acquiring the system's stability region and defining stability criteria is another vital aspect of large-signal stability analysis for microgrids.In the subsequent section, various stability analysis methods are elucidated and formulated.

IV. TECHNIQUES FOR CONSTRUCTING DOA WITH THE LYAPUNOV-BASED METHOD
The Lyapunov-based method is a powerful tool used in control theory for analyzing the stability and domain of attraction of dynamic systems.The Lyapunov function search technique involves searching for a Lyapunov function that satisfies certain conditions, such as being positive definite and having a negative definite derivative [114], [115], [116].Once a Lyapunov function of a system is found, the DOA of the system around a specific equilibrium point can be estimated as the largest set of initial conditions for which the function is positive, and the derivative of the function is negative.
There are several techniques for constructing the DOA of a system around a specific equilibrium point using this method (e.g., LMI-based method, Krasovskii's, Popov-Lure, SOS-based).These techniques provide powerful tools for constructing the domain of attraction using the Lyapunovbased method.The choice of a technique depends on the specific system being analyzed and the complexity of the Lyapunov function needed to characterize the system's dynamic behavior.The existing well-known methods for stability analysis can be classified as follows.

A. KRASOVSKII'S METHOD
One of the simplest methods to construct the Lyapunov function for a nonlinear system is Krasovskii's method.In particular, Krasovskii's method involves constructing a Lyapunov function that satisfies certain conditions, such as being positive definite and having a negative definite derivative along the system trajectories.This method is often used in the analysis and design of nonlinear control systems, and it has been shown to be effective for a wide range of applications [42], [117].
Mathematically, for constructing a Lyapunov function of a nonlinear system based on Krasovskii's method, the state-space model of the system must be defined as follows: ẋ = f (x), and f (0) = 0. ( The Lyapunov function for this system can be defined as follows: where P is a symmetric positive definite matrix.The derivative of the Lyapunov function is as follows: where Jacobean matrix: Therefore, the derivative of the Lyapunov function can be modified as follows: where Q must be positive definite or V (x) must be negative definite to be asymptotically stable at the origin.

B. POPOV-LURE
If the state-space model of a system is defined as follows: the system will have two parts, linear and nonlinear parts determined by Ax, and bf (τ ), respectively.The transfer function of the linear part, Ax is: Lure theory and Popov criteria ( [118] and [119]) can provide conditions for stability analysis.According to Kalman [120], the state-space model of the system can be transferred to a reduced-order model as follows: where f (τ ) is a nonlinear function such that The aforementioned system has two parts, linear and nonlinear parts determined by Fx, and gf (τ ), respectively.The transfer function of the linear part, Fx is: According to the Lure-Lyapunov lemma [120], to construct the Lyapunov function, the two parameters in the following equation, real number γ and positive definite matrix K, need to be determined: To find the required parameters in (11), the next steps listed below need to be followed: where γ ≥ 0 and is a real number such that: where where K is positive definite.
By following the aforementioned steps in (12), the Lure-Lyapunov function in (11) can be obtained.To determine the region of stability, two following boundaries obtained by finding the local maximum of (11), are defined as follows:

C. LINEAR MATRIX INEQUALITIES (LMI)
The use of LMIs makes the Lyapunov function construction process tractable and computationally efficient, even for complex nonlinear systems [121], [122].Constructing a Lyapunov function for a nonlinear system utilizing an LMI-based method involves a rigorous set of steps to ensure stability analysis of the system with profound implications.The procedure can be summarized as follows: • Define the system dynamics: First, the system dynamics must be expressed in the form of a differential equation or a set of differential equations defined.Let us consider a nonlinear system described by the differential equation: where x represents the state vector and f (x) embodies the nonlinear dynamics of the system.
• Choose a candidate Lyapunov function: The next step is to choose a candidate Lyapunov function that satisfies certain properties.The Lyapunov function is a scalar function of the state variables that can help to determine the stability of the system.Choose a quadratic function of the form as follows: where P is a positive definite matrix.
• Compute the derivative of the Lyapunov function: Compute the derivative of the Lyapunov function with respect to time.Using the chain rule, V (x) is given by: Substituting ẋ with the nonlinear dynamics of the system f (x) and simplifying the expression, the derivative of the function can be modified as follows: • Choose a Lyapunov candidate: To guarantee the stability of the system, a positive definite matrix P must be found that satisfies the following linear matrix inequality: where A is the Jacobian matrix of f (x), Q is a positive definite matrix, and ''< 0'' denotes negative definiteness.
• Check the stability of the system: Once a positive definite matrix P is found that satisfies the LMI conditions, the Lyapunov function V (x) = x T Px can be utilized to analyze the stability of the system.If V (x) < 0 for all x in the state space, the system is stable.If V (x) > 0 for some x in the state space, the system is unstable.

D. TAKAGI-SUGENO FUZZY
The T-S (Takagi-Sugeno) fuzzy modeling technique is a popular approach for modeling nonlinear systems using a set of local linear models.This technique can also be used to construct Lyapunov functions for analyzing the stability of the system [121], [123].
The T-S fuzzy modeling technique involves partitioning the state space into a set of fuzzy regions, each of which is associated with a local linear model [124].The fuzzy membership function determines the degree of membership of a given state vector to each of the fuzzy regions.The T-S fuzzy model can be written as: 111476 VOLUME 11, 2023 Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.
where A i is the fuzzy set associated with the i th local model, and g i (x) is a linear function of x that describes the dynamics of the i th local model.
To construct a Lyapunov function for the T-S fuzzy model, we can use the following approach: • Define the system dynamics: Define a common quadratic Lyapunov function for each of the local models: where P i is a positive definite matrix associated with the i th local model.
• Choose a candidate Lyapunov function: Combine the local Lyapunov functions into a global Lyapunov function that satisfies the stability conditions for the entire system.This can be done using a weighted sum of the local Lyapunov functions as follows: where w i is a non-negative weight associated with the i th local model and K is the number of fuzzy regions.
• Compute the derivative of the Lyapunov function and check the stability of the system: Choose the weights w i such that the global Lyapunov function satisfies the conditions for the stability of the entire system.This can be done using LMI optimization techniques, similar to the approach used for constructing Lyapunov functions for nonlinear systems.The resulting Lyapunov function is a weighted sum of quadratic functions, which is a convex function that satisfies the necessary properties for analyzing the stability of the T-S fuzzy model.

E. SUM OF SQUARES
In Parrilo's thesis, the sum of squares method was first presented [125].In this method, numerous problems in system analysis that were previously challenging to answer have been addressed.One of the addressed challenges is the algorithmic stability analysis of nonlinear systems using Lyapunov techniques.In this paper, the application of the SOS-based Lyapunov function is suggested for stability analysis of the inverter-based microgrid to find a region of attraction with lower conservativeness.The SOSTOOL is applied to find the polynomial Lyapunov function.The sum of squares polynomial optimization programs can be created and solved using the free MATLAB toolbox SOSTOOLS [126].
Constructing a Lyapunov function for a nonlinear system using the sum of squares techniques involves the following steps [127], [128], [129]: • Define the system: Consider the following nonlinear system: where x is the state vector and f (x) represents the nonlinear dynamics of the system.
• Choose a candidate Lyapunov function: The next step is to choose a candidate Lyapunov function that satisfies certain properties.The Lyapunov function is a scalar function of the state variables that can help to determine the stability of the system.For the SOS-based Lyapunov function, choose a polynomial function of the form as follows: where p k is the coefficient of the polynomial; d is the degree of the polynomial; x n is the state; and n is the number of states in the state-space model of the system.For example, for a two-state state-space model with degree 2 of the polynomial, the following equation can be a candidate for the Lyapunov function: • Compute the derivative of the Lyapunov function: The derivative of the Lyapunov function with respect to time is calculated.Using the chain rule: where ∇V (x) is the gradient of the Lyapunov function and f (x) represents the nonlinear dynamics of the system.The derivative of the Lyapunov function in terms of the Lyapunov candidate can be computed using partial differentiation: • Choose a Lyapunov candidate: To guarantee the stability of the system, a Lyapunov candidate must be obtained that satisfies the following SOS constraint: where α is a positive scalar.
• Check the stability of the system: Once we find a Lyapunov candidate that satisfies the SOS constraint, we can use the Lyapunov function V (x) to analyze the stability of the system.If V (x) + αV (x) < 0 for all x in the state space, the system is stable.If V (x) + αV (x) > 0 for some x in the state space, the system is unstable.Constructing a Lyapunov function using SOS techniques involves choosing a candidate function, computing its derivative, and finding a Lyapunov candidate that satisfies an SOS constraint.Some numerical examples for simple systems can be found in [128] and [130].By analyzing the derivative of the Lyapunov function, the stability of the system and its domain of attraction can be determined.
The studies in other applications [90], [131] show that the SOS-based Lyapunov function can provide a larger and less conservative domain of attraction for a certain system compared to the other existing methods (Krososkii's, Popov-Lure, and T-S fuzzy methods), guaranteeing the accurate and less conservative stability assessment of the system.However, by increasing the order number of the state-space model of the system, the computational process gets extremely timeconsuming [131].
In the following sections, a numerical analysis of the large-signal stability of a microgrid with a single inverter-based DER is developed.A comparative study is presented by using different dynamic models and methods of Lyapunov function construction.

V. CASE STUDIES AND SIMULATION SETUPS
In the following sections, a numerical analysis of the existing methods and dynamic models is developed and assessed.The analytical, numerical assessment is provided for the second-order (VSG) and full-order dynamic models.The comparative assessment of the microgrid's stability is obtained based on the outcomes of the three Lyapunov-based stability methods named Krasoskii's, Popov-Lure, and SOS methods.

A. METHODS FOR CONSTRUCTING DOA FOR THE VSG DYNAMIC MODEL
The dynamic model employed for a system must be sufficiently accurate to predict the system's behavior reliably.However, using an excessively complex model can lead to computationally demanding processes.The dynamic model described in this section for inverter-based DERs is the second-order or virtual synchronous generator (VSG) model.This model treats the dynamic behavior of IBRs such as synchronous machines, resulting in defining an equivalent damping factor and inertia values.Consequently, the VSG model simplifies the modeling process while still capturing the essential dynamics of inverter-based DERs.

1) THE MATHEMATICAL MODEL FOR INVERTER-BASED MICROGRIDS BASED ON SECOND-ORDER MODEL -VSG
The active power can be effectively modeled based on the droop control mechanism employed by the inverters [132], [133], [134].The model can be expressed as follows: ) where m p is the droop gain of P/f loop.(.) d,ctrl and (.) out indicate the control and reference signals.(.) odq is the d-q frame voltage and current at the inverter connection node to the LC filter (shown in Figure 8).ω f , ω out , and ω d,ctrl are the low-pass filter and microgrid, control signal frequency, respectively.According to reference [37], microgrid droop control emulates the characteristics of synchronous machines related to their inertia and damping properties.This similarity is expressed through the following formulation [135]: According to (28a) and (28b) and transposing related terms, the aforementioned equation can be written as follows: To determine the virtual inertia and damping factor of IBRs, the equation obtained can be compared with the standard definition of the swing equation used for synchronous generators.This allows for the formulation of the following expressions: Numerous studies have confirmed that IBRs exhibit synchronous machine-like characteristics, with various definitions of inertia and damping factor proposed in [37], [38], and [135].Despite the differences in these proposed definitions, the resulting values are notably close to one another.For instance, the model developed in [38] treats the inverter-based source as an electrostatic machine and calculates the virtual inertia (M) and damping factor (D) using the inverter-based equations given below: where H is the inertia constant and is equal to 1 2 C dc V 2 dc ; C dc and V dc are DC-side capacitor and its voltage magnitude.Since the value of m p is significantly smaller than 1, the inertia M proposed in [38] can be considered practically equivalent to 1 m p ω f .
In the next step, the state-space dynamic model of the network is defined.The network includes an IBR connected to an infinite bus and the voltage at the point of common coupling (PCC) is assumed to remain constant.Consequently, the dynamic variables are limited to the angle and frequency.However, when applying this model to the islanded mode of operation, additional complexities arise due to voltage fluctuations.In such scenarios, the assumption of a constant voltage at PCC no longer holds true.
According to the obtained VSG model of IBR and the assumed network, the state-space dynamic model of the network can be formulated as follows: (33) 111478 VOLUME 11, 2023 Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.where θ Y , A, and B are defined based on the network and load parameters that can be found in Appendix A in [136].D and M are the same as obtained in (31).Then, the state-space model of the network is represented as: where x is defined as follows: Note that δ s is the equilibrium point that can be found by replacing P ref as desired produced power by IBR with P e in ( 33)-(*).As the system defined by ( 34) has a zero eigenvalue, a nonsingular transformation can be applied to convert it into the form utilized by Kalman [120] and can be transferred to the following reduced-order form, as demonstrated in the following expressions: where x is donated to reduced-order model from n orders to n − 1 orders.Once the appropriate dynamic model for the IBR is incorporated into the state-space dynamic model of the network, an appropriate method must be employed to analyze the network's stability.As detailed in Section IV-B, the parameters required to construct the Lyapunov function based on the Popov-Lure criteria can be determined using the following approach: Applying Popov Criterion: Here, γ is considered such that:

•
The real vector q is considered such that: Then, the Lur'e-type Lyapunov function is expressed in the following form: − cos(τ By replacing τ = δ−δ s and △ω with x 1 and x 2 , respectively, the Lur'e-type Lyapunov function is modified as follows: Using Lagrange multipliers, the minimum non-zero boundary indicated by R, can be obtained by substituting l = π −2(δ s +θ Y ) in (39).To determine the region of stability, the critical energy boundary is defined as follows:

2) NUMERICAL ANALYSIS FOR VSG MODEL
The state-space representation based on the VSG model is constructed for a grid-connected IBR detailed in Section V-A1.The Popov-Lur'e and sum of squares (SOS) methods, detailed in Sections IV-B and IV-E, respectively, are utilized to construct the Lyapunov function and the domain of attraction.Additionally, Table 3 provides the parameters for the IBR.

3) RESULTS AND DISCUSSION FOR VSG MODEL
The stability analysis results, which are based on the VSG model and SOS-based Lyapunov function, are presented in Figure 4 in a counter plot.Figure 4a and 4b illustrate the SOS-based Lyapunov function and its derivative, respectively.According to the Lyapunov stability theory, the intersection between the positive Lyapunov function and the negative derivative of the Lyapunov function represents the region of stability.By satisfying these two conditions, the domain of attraction is determined, as shown in Figure 4b.The comparison presented in Figure 7 between DOAs constructed using the SOS-based and Popov-Lyapunov methods reveals that the Popov-based DOA is overly conservative and smaller than the SOS-based DOA.This conservatism leads to unnecessary constraints on system parameters, which can have negative effects in practical applications such as protection systems and relay settings.Strict limitations may result in unnecessary triggering of the protection system, causing it to miss out on utilizing the full potential of the system.Therefore, an inaccurate and conservative DOA can impede the system's efficiency and effectiveness.

B. METHODS OF CONSTRUCTING DOA FOR THE FULL-ORDER DYNAMIC MODEL 1) THE MATHEMATICAL MODEL OF INVERTER-BASED MICROGRIDS BASED ON FULL-ORDER DIFFERENTIAL EQUATIONS
To create a further reliable dynamic model of the inverter-based microgrid, it is crucial to model all existing components, including DGs, loads, converters, and batteries, using appropriate approximations.
This section focuses on the IBR, its network, and control systems, as depicted in Figure 8.The network includes an IBR connected to the main grid and two loads.The dynamic model used for the IBR in this study is based on [50] and [137].The subsequent sections explain the theoretical model 111480 VOLUME 11, 2023 Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.   of the IBR, which covers power control, voltage control, current control, LCL filter, lines, connection, and loads.The state-space model of the microgrid is developed by combining and modeling all components according to the network configuration using the network mapping model.In this study, the DC-side components are excluded from the model to simplify the computation process.Additional information about the DC-side dynamic model can be found in [50].

a: POWER CONTROLLER
The power controller comprises several components, including the low-pass filter (ω f ), Q/V control loop, and P/f control loop.The nonlinear state-space model of the power controller is detailed below: In the rest of the paper, for simplification, P d,ctrl and Q d,ctrl are replaced with P and Q.
where n q and m p are the droop gain for Q/V and P/f loops, respectively.According to (41a) and (42b), the nonlinear dynamic model of the power controller is as follows: where where V n stands for d-axis nominal output voltage, ω n is the nominal frequency set in the inverter droop controller, ω indicates the local rotating reference frame of each inverter, and ω ref is the angular frequency of the global reference frame.The angle, δ, between the inverter's local d-q frame and the global D-Q reference frame is shown in Figure 9.

b: VOLTAGE CONTROLLER
The differential equation that represents the voltage controller dynamics is described as follows: where ϱ shows the variation of voltage in the d-q frame compared to the reference values.
The nonlinear dynamic model of the voltage controller is defined as follows: where In the given equations, K p v and K i v stand for the proportional-integral (PI) coefficients of the inverter's voltage controller.Moreover, F c denotes the voltage controller of the inverter's feed-forward coefficient.The grid frequency used in this study is set as the angular frequency of the global D-Q reference frame in the grid-connected microgrid.

c: CURRENT CONTROLLER
To formulate the state-space dynamic model of the current controller, the same approach as the voltage controller is applied.Equations in (48) describe the present controller's differential equation based on the circuit laws: where γ represents the variation of current in the d-q frame compared to the reference values.The state-space dynamic model of the current controller is obtained as follows: where In (50e), the PI coefficients for the current controller are K ic and K pc .

d: LCL FILTER
Using the circuit laws, the nonlinear differential equations of the LC filter are given by: Then, the state-space model of the LC filter can be described as follows: where The other parameter of (54b)(mentioned in bottom of next page) are defined as follows: where the output voltages for the Park frame are v bd and v bq .

e: STATE-SPACE MODEL OF AN IBR
The three-phase system is balanced free of harmonics, with a modulation ratio of 1 being assumed.To develop the state-space model for an IBR, the dynamic models of all components are consolidated, as described in [50] and [137].The combined state-space dynamic model is given by: where where C ′ P in (57), as shown at the bottom of the next page, is the rows 2 and 3 of C P in (44f).

f: NETWORK
Based on the circuit laws, the network's nonlinear dynamic differential equations at PCC are given by: di The state-space dynamic model can be obtained as below: where where (.) Ntwrk and (.) LINE are utilized to denote the network topology and line configuration, respectively.Furthermore, in accordance with Figure 8, v k and v b are defined as v b2 (or v 2 in Figure 8) and v b1 , respectively.The line current in the B LINE matrix is directed outward from node b (+1/L line ) and inward to node k (−1/L line ).All other elements within the B LINE matrix are 0.

g: LOAD
This analysis takes into account a common sort of constant resistive-inductive (RL) load.Below is a model of the load dynamics behavior at each node. di Then, the state-space dynamic model of all network loads is formulated as follows: where The transfer function for the global D-Q reference frame is represented by T DQ .

h: MICROGRID MODEL AND MAPPING MATRICES
To measure the voltage of individual nodes, a virtual resistor (r v ) is positioned between the ground and each node.In Figure 8, the virtual resistor is illustrated, with a value of 8000 being selected to minimize any adverse effects on the dynamic model's accuracy.The voltage of each node can be determined using the following expression: The load (L), line (Ntwrk), and inverter (IBR) models are developed and subsequently integrated to generate the voltage values, as expressed in (64).Rv matrix is determined based on the node configuration and the number of nodes, comprising diagonal components with a value of r v .The M IBR matrix is utilized to map the connection points of the inverters to the grid nodes, with the matching element being 1 or 0 depending on whether the inverter i is linked to node j.The relationship between loads and grid nodes is mapped using M L , with the element being 1 if no relationship exists between the load and the node, and 0 otherwise.The M Ntwrk matrix is responsible for mapping the lines to the grid nodes.
The matrix element corresponding to a node has a value of 1 when the current enters the node and a value of -1 when it exits [50], [137].The state-space model of the inverter-based microgrid, which incorporates all components, is elaborated in the following equation: where A Microgrid is presented in (67), as shown at the bottom of the next page, and additional information about the matrices of ( 67) is available in [137].

2) NUMERICAL ANALYSIS
The case study involves an IBR connected to the primary grid and two loads.Table 4 presents the network parameter.The stability analysis of Krasovoskii's approach and the SOS-based Lyapunov function are compared and elucidated.The sensitivity analysis of the SOS-based method is investigated for two situations: changes in load and control scenarios.

3) RESULTS AND DISCUSSION
The SOS-based technique has the potential to enhance the domain of attraction to a verifiable state space, supporting the accuracy claims mentioned above for the SOS-based Lyapunov function.Figure 10 displays a comparison of the DOA constructed using two techniques, namely Krasovskii's and the SOS-based methods.The DOA obtained via Krasovskii's method is exceedingly restricted compared to the SOS-based approach.The power system's response to contingencies in Krasovskii's method is rather conservative due to the limited domain of attraction constructed by 111484 VOLUME 11, 2023 Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.Krasovskii's method.The transient stability analysis, such as the domain of attraction, is typically employed to adjust the relay and protection systems.The domain of attraction can be used to determine the critical period required to clear faults and other contingencies, ensuring stable and reliable network operation.A limited domain of attraction results in a conservative adjustment of the protection mechanism, leading to an unnecessary outage or load shedding in the system.The Lyapunov function developed must meet the following criteria for stability region based on the Lyapunov-based stability: must be equal to zero at x n = 0. V (x n ) > 0 ⇒ the V must be positive in stability region.
V (x n ) ≤ 0 ⇒ the derivative of V must indicate a downward slope in stability region.(68) The network's state-space model is shifted to its equilibrium point or steady-state operation for the first condition, meaning that the network's equilibrium point is represented by the equation x n = x − x s 0 .Thus, V (x n ) is 0 for x n = 0. Figure 11 portrays the SOS-based Lyapunov function and its derivative.Negative contours on the Lyapunov function plot do not satisfy the second Lyapunov-based stability criterion.Hence, the areas in Figure 11a with negative contours are excluded from the stability region.The negative and zero contours in Figure 11b satisfy the third requirement; i.e., the negative derivative of the Lyapunov function, which excludes the light-colored regions (positive contours) from the stability region.
In the subsequent section, the impact of load variation events and control modifications on the stability region in inverter-based microgrids is examined.These scenarios are frequently used as an inevitable action of a smart system, they get modified during both temporary and permanent changes in the system to guarantee stable and reliable operation.The goal of analyzing these scenarios is to ensure that the microgrid can adapt to altering conditions while remaining stable, thereby avoiding any interruptions or system failures.By adapting to changing conditions and implementing appropriate control modifications, the stability region can be sustained, ensuring that the microgrid operates reliably and effectively.

a: LOAD CHANGING
Load reconfiguration is a solution that can enhance the reliability of a system during contingencies, thereby inducing a change in the stability region.Proper load reconfiguration or reduction can amplify the stability region, leading to an enhancement of the network's ability to withstand transient faults that may occur.By adjusting the load appropriately, the system's resilience to unexpected changes can be strengthened, thereby mitigating any disruptions or failures, and ultimately boosting its overall performance.
In addition, the impact of load increasing on the stability region is demonstrated in Figure 12.By increasing Load 1, the contour lines for the same V -slices become more condensed towards a smaller region of stability, indicating a reduction in the system's stability margin.This means that as the load increases, the system's ability to remain stable becomes more constrained, thereby increasing the risk of instability and potential system failure.Therefore, it is crucial to manage the load effectively to ensure that the system's stability is not  Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.
compromised, and the stability region is maintained or even expanded, thereby ensuring the network's reliability during both normal and contingency operations.

b: CONTROL MODIFICATION
Inverter switching systems implement various control techniques to ensure reliable and stable operation.A smart control system is designed to modify the control parameters during both normal and contingency operations.The control system employs several control objects to make appropriate decisions for the switching system, and the region of attraction can be an additional criterion for modifying the control parameters.Therefore, the control system can adjust the control parameters in response to changes in the system's behavior, ensuring that the system remains stable and reliable, ultimately improving its operational efficiency.
By increasing the K pv in the voltage controller loop, the stability region of the system is expanded.As depicted in Figure 13, the contour lines of the same value of the Lyapunov function are expanded.This results in an expansion of the stability region and an enhancement of the system's stability capability.Thus, by properly modifying the control parameters, the system's stability region can be improved, leading to a more robust and reliable operation of inverterbased microgrids.

VI. DISCUSSION AND CONCLUSION
This paper has provided a comprehensive critical review of large-signal stability analysis methods and models for microgrids.Numerical examples with critical analyses are provided to compare the performances of existing methods and models.Sections V-B and V-A have presented a comparative analysis of two different state-space models for an inverter-based microgrid.The study developed a nonlinear dynamic model to investigate a microgrid connected to the main grid, and the results demonstrated that (a summary is provided in Table 5): • Utilizing the Lyapunov-based stability method based on the sum of squares can produce a stability region that is less conservative than those generated by Krasovskii's and Popov-Lure's This precise stability region helps in accurately determining the critical clearing time for given contingencies.
• The type of chosen nonlinear dynamic model for the accuracy and authenticity of the stability assessment is crucial as an appropriate dynamic model allows for a more precise presentation of the microgrid's dynamic performance concerning the grid it connects to.
• The results emphasize the importance of choosing the appropriate stability analysis method, as it can significantly affect the accuracy of the predictions.The SOS-based Lyapunov method is superior in this regard and should be considered for future stability analyses of inverter-based microgrids.This can be particularly helpful in designing and optimizing the performance of such systems to improve the stability and reliability of power grids.
• It is worth mentioning that even though the Lyapunov function based on the SOS-based method has shown to be effective in analyzing the stability of microgrids, its computational process can be time-consuming, especially for higher-order state-space models.Additionally, in some cases, it may not converge to an acceptable result.Furthermore, for larger systems, using higher-degree Lyapunov functions with SOSTOOL can lead to a substantial rise in the number of optimization decision variables, presenting a challenge in terms of computational resources.
• To overcome the challenges associated with using the SOS-based method for stability analysis of large-scale systems, it is recommended to employ a simplified and smaller-order dynamic model for IBRs.This approach can reduce the computational demands of the SOS-based method, resulting in improved efficiency and effectiveness in analyzing the system's stability.However, it is crucial to ensure that the simplified models can still capture the critical dynamics and characteristics of the inverter-based microgrid being analyzed while having reasonable approximations.One model that has been shown to be both accurate and acceptable for inverters is the 5-order state-space model presented in [46].Another low-order state-space model, known as the second-order or virtual synchronous generator model [39], can only provide frequency and angle as the states of the dynamic model for grid-connected microgrid and IBRs.Therefore, selecting an appropriate simplified model for the inverter-based microgrid being analyzed is critical to balance the need for accuracy and computational efficiency.
• The selection of a model for stability analysis requires a trade-off between the number of states that can be monitored and the computational complexity of the model.Selecting the appropriate model should balance computational efficiency and the accuracy of the stability analysis results.The decision will depend on the specific needs and requirements of the stability analysis application.If accuracy is the top priority, then the 5-order state-space model may be the preferred option since it provides more states that can be monitored.However, if computational efficiency is the priority, then the second-order or VSG model may be the better choice for the stability analysis of inverter-based microgrids.Nevertheless, it is essential to recognize that simplified models may not fully capture all of the complex dynamics of the microgrid, and as such, the stability analysis results should be interpreted with caution.In conclusion, the selection of the model should consider the trade-offs between accuracy and computational efficiency, as well as the specific needs of the stability analysis.VOLUME 11, 2023 Authorized licensed use limited to the terms of the applicable license agreement with IEEE.Restrictions apply.

FIGURE 3 .
FIGURE 3. Domain of attraction with stable and unstable transients.
Figures 6a and 6b illustrate the 3-and 2-dimension presentations of the Lyapunov function constructed based on the Popov-Lure theory considering the VSG model of IBRs, respectively.As described in Section IV-B, the DOA is defined by the positive Lyapunov function and the negative derivative of the Lyapunov function around the original equilibrium point.The critical energy boundary is determined by Lagrange multipliers (l and R), described in (40), and is shown in Figure 6b by the white line surrounding the black surface.

FIGURE 6 .
FIGURE 6. DOA constructed by Popov-Lure method for VSG model of IBRs.

FIGURE 7 .
FIGURE 7. Comparison of DOAs constructed by SOS method and Popov-Lure method for VSG model of IBRs.

FIGURE 8 .
FIGURE 8. Network and the DC/AC droop controlled inverter.

FIGURE 9 .
FIGURE 9. Local d -q and global D-Q reference frames.

FIGURE 11 .
FIGURE 11.DOA constructed by SOS Method for the full-order dynamic model.

TABLE 2 .
Literature review studies on the large-signal stability of microgrids and IBRs.

TABLE 3 .
The parameters of the inverter control and network.

TABLE 4 .
The parameters of the inverter control and network.

TABLE 5 .
Advantages and disadvantages of the stability assessment methods according to the numerical assessments.