On Small Signal Frequency Stability under Virtual Inertia and the Role of PLLs

This paper presents a methodology that aims at identifying virtual inertia (VI) gain limitations from virtual synchronous generators (VSGs) while maintaining the frequency stability considering the delay associated with the frequency measurement process. The phase-locked loop (PLL) is typically used for frequency estimation that is used to calculate the rate of change of frequency (RoCoF) and it drives the VI loop. The PLL is generally accompanied by a low-pass filter that aims to suppress the impact of harmonics. This filter introduces a delay that when used with the VI control loop causes stability issues for high values of VI gain. A comparison of various PLL approaches suggests that certain variants tend to permit higher value of cut-off frequencies which can be utilized to increase the VI gain limit from VSG. This study presents a method by which the upper limit on VI gain can be quantified and related to the cut-off frequency of the PLL low pass filter that is indirectly representing the delay. It is performed using small signal frequency stability analysis on the frequency domain model of the grid with virtual inertia emulating VSG. The effective maximum VI gain from VSG is explored while satisfying the frequency measurement accuracy specification considering harmonics. The results show that the requirements of reaching a stable operation with sufficient stability margins can still be met with a faster PLL-based system and the potential increases in VI support from VSG can be quantified using the proposed method. The study has been first performed on a single machine single inverter bus (SMSIB) system and is generalized to the multi-machine and multi-inverter system.


Introduction
The climate change impacts have necessitated the need of transition from traditional combustion based generation to environmentally friendly sources of energy.In relation, recent decades have experienced an increasing trend towards distributed generation from wind and solar energy [1].This trend has important consequences on system stability.The power system is stable if mechanical power provided to synchronous generators is in balance with the electrical demand.An imbalance causes a shift in the frequency as rotor of the generator tend to slow down (speed up) depending on the torque experienced by the rotor resulting in frequency decrease (increase).The inertia response of the machines acts instantaneously and the kinetic energy in the rotor is converted to electrical power, limiting the RoCoF.The phenomenon is termed as inertia response from the generator.In addition, the frequency variation is sensed by the governor and the mechanical input of the generator is controlled accordingly.This action is termed as the primary frequency response.The mechanical inertia in the rotor of synchronous generators has been instrumental in limiting the RoCoF leading to reduced oscillations and peaks in the frequency termed as "frequency nadir".The frequency settles to a new steady-state value with the generator output power determined by the droop gain of the power-frequency governor loop.As synchronous generators are replaced with inverter interfaced generation, the mechanical inertia is replaced by a limited amount of energy available at DC link capacitors in case of PV, and in the rotors of wind turbines [2].Additionally, in the inverter dominated generation, there is a delay associated with the frequency measurement process.The decline in inertia and delayed response has direct consequence on the frequency stability resulting in increased RoCoFs and more extreme frequency nadirs that can cause undesired load-shedding, cascading failures or large-scale blackouts.
Recently, the impact on frequency stability due to a decline in system inertia has become a matter of significant importance as reflected by different reports and studies [3].The report of electricity reliability council of Texas (ERCOT) [4] discusses the reduction in inertia caused by the economic displacement of synchronous generators in meeting the inertia requirements when wind generation is serving the highest share of load.In the European network of transmission system operators (ENTSO-E) framework, the primary reserve needs to be activated when the frequency deviation exceeds ±10 mHz from the nominal frequency of 50 Hz.In recent years, ENTSO-E has experienced increasing levels of deterministic frequency deviations resulting in the activation of a significant share of primary frequency control reserves that have been initially intended and dimensioned for large generation and load outages [5] The synchronous area continental Europe has experienced RoCoF of up to 1 Hz/s in some control areas having the highest risk of imbalance accounting to a 20% difference between supply and demand.This level is expected to rise to 2 Hz/s with a 40% imbalance ratio [6].Along with the high dynamics of the frequency transients, a significant drop in frequency can lead to triggering of under frequency relays that can cause load shedding.Some case studies are an exception, such as a study of the RoCoF in Austrian grid reported in [7] that suggests that no critical increase in RoCoF has been observed in the first quarter of 2018 despite 74% share of instantaneous non-synchronous generation.The result cannot be generalized and is possibly due to the conventional generation running at partial loading or in idle mode, thus the inertia is still available.However, when the conventional generation is dominantly replaced by non-synchronous generation, the RoCoF is expected to increase.
Various steps have been taken to counteract the declining inertia situation.In several countries, the inertia response from wind farms has been made mandatory [8,9].In addition a new product called enhanced frequency response (EFR) has been introduced in the UK [10] that requires a response time of under 1 s.Among market participants, energy storage systems (ESS) have emerged as a key player [11].The control system of ESS provides the inertia support by emulating the response of a synchronous generator.However, the main limiting factors of the inertia support from ESS are their stored energy and the capability to respond to the fast frequency transients.In this context, the estimation and processing of frequency from the time domain voltage signal is the most critical part during the inertia response [12].
In this paper, the inertia emulating source under current control is termed as virtual synchronous generator (VSG).It operates in grid following mode in terms of frequency control and can represent an energy storage system, PV, wind turbine or other inverter interfaced generation.The virtual inertia emulation from a VSG requires calculation of RoCoF, for which the frequency measurement is central.The best frequency measurement method should ideally be accurate, easy to implement, converge quickly to the expected value, reject noise and harmonics and require low computational power [13].Several frequency measurement methods have been proposed in the literature.Some of them without being exhaustive include zero-crossing [14], least squares [15], demodulation methods [13], Kalman filtering [16,17] and the phase-locked loop (PLL) [18].Each method either applied to single phase or preferably to three phases must deal with the harmonics, voltage imbalance, phase angle jump and non-stationary environments [16].This paper focuses on the frequency estimation using PLL due to its wide spread application in practical inverters.From the application perspective, authors in [19] have explored the fast frequency response potential from a offshore multi-terminal high voltage direct current scheme where a real time digital simulator based PLL module has been used that introduced a delay of 120 ms in detection of voltage variations.Similarly in [20], the synthetic inertia contribution from electric vehicles is explored.The frequency estimation introduces a delay of 200 ms during signal propagation and filtering in the PLL.This delay plays an important role in limiting the virtual inertia support while studying the stability of the power system.In literature, the quantification of virtual inertia support and its relevance with the delay in frequency estimation has not received significant attention.Authors in [3] reviews new stability analysis methods that consider delays inherently using delay differential-algebraic equations (DDAEs) representing converters in low inertia studies.This approach aims to reduce the large scale DDAEs into low-frequency and critical eigen-values whose influence on the stability margin can be studied [21,22].This approach needs further extension to consider multiple and time varying delays while considering the virtual inertia gains of the VSGs.
This motivates the presented study that models the influence of delays in frequency measurement as the cut-off frequency of the PLL low-pass filter and relates it to the eigen-value based stability analysis of the system.The damping criterion determines the virtual inertia support limitations from a VSG.The contributions from this work can be classified into two categories: firstly, a parametric study is performed to investigate the influence of low-pass filter cut-off frequency ( f c ) on the frequency estimation accuracy using various PLL approaches.A comparative study of PLL variants suggests that the second order generalized integrator frequency locked loop (SOGI-FLL) approach with the low-pass filter provides a good estimate of frequency.A parametric study identifies the maximum value of f c that satisfies the frequency estimation accuracy requirements while considering harmonics in reference to EN50160 guidelines [23].SOGI-FLL allows to increase the f c and thus increasing the bandwidth of the VI loop.The second part of paper focus on the quantification of maximum provision of the VI gain given the higher value of the f c .The limits are identified using small signal stability analysis of the equivalent frequency domain model of the grid having VSGs.The most relevant paper to this study has been [24], where authors have performed the small signal analysis of a microgrid using aggregate models.The impact of f c and the effective virtual inertia gain are separately studied.This paper extends the study in quantifying the increase in maximum virtual inertia gain from the VSG for a given range of f c .The result is crucial to quantify the increase in inertia support that can realized if a fast and robust frequency estimation and processing technique is used.The frequency stability is assessed with a conservative damping ratio taken as a reference.Analysis of the proposed approach is performed for a SMSIB system and is generalized to a system having multiple synchronous machines and VSGs.
Following this introduction, the paper is organized as follows: Section 2.1 presents the frequency estimation methods.A parametric study of finding the f c limits and a harmonic sweep analysis is performed to explore the impact of filtering on harmonic reduction in Section 2.2.It is followed by the discussion on inertia from VSG and small signal frequency stability model of the VSG and grid in Section 2.3.The models are used in the case studies for single and multi-inverters providing virtual inertia in Sections 3.1 and 3.2, respectively.The paper ends with an conclusion and outlook in Section 4.

Modeling and Analysis of PLL and Virtual Inertia in Grid
This section provides a comparison of the various PLL-based frequency estimation approaches under a step event.The delay introduced by filtering is studied and subsequently a parametric study is performed to study the influence of PLL low-pass filter cutoff frequency on the frequency measurement accuracy.The objective is to maximize the overall bandwidth by increasing the cut-off frequency in order to support higher levels of the VI gain while limiting the impact of harmonics in the frequency measurement process.This is followed by the discussion about inertia concept and virtual inertia provision from VSG.A small signal stability model of the aggregate system is derived combining frequency domain model of VSG with its PLL, dynamic loads and the grid frequency response.

Frequency Estimation Methods
A typical process of frequency estimation using PLL starts with the monitoring of time domain voltage signal.Since the voltage measurement is accompanied by harmonics and noise, a filtering process is necessary.Typically a low-pass filter is used at the output of PLL in order to suppress the impact of harmonics on the frequency estimation [24,25].This filtering introduces some delay, which plays, as shown in the next sections, a decisive role in the system stability.In order to get an insight in the delay introduced by the filtering process, a frequency step function is considered.It is to note that a characteristic frequency response function is not used here as the objective is not to test the behavior of PLL across the spectrum of frequency measurement but rather is to study the limiting factors of obtaining accurate frequency measurements while considering harmonics.A step function also provides basis for comparison of various PLL techniques (with and without filter) in order to observe the delay introduced due to the filtering process.Three variants of PLL have been used in this work to assess the performance of frequency measurement from PLL:

•
Single phase PLL with automatic gain control as discussed in [26] with its implementation referred here as Matlab function [27].• SOGI-FLL introducing a control loop to auto-adapt the center frequency of SOGI to the input frequency [28].It operates on the single phase voltage signal.• SOGI-FLL with a sixth order low-pass Butterworth filter at its output.
The comparison of three models is shown in Figure 1, where SOGI-FLL based model with the low-pass filter gives comparatively improved results in terms of settling time while satisfying the 1 mHz accuracy requirement [5,29].SOGI-FLL is selected and referred to as PLL in rest of the paper.

Impact of the PLL Cut-Off Frequency on the Robustness against Harmonics
In this section, harmonics are added to the time domain voltage signal and PLL is used to track the frequency, moreover the influence of f c is explored.In order to investigate the impact of disturbances, a range of harmonics are added to the frequency signal.The standard EN 50160 provides guidance on the harmonic levels that are not expected to be exceeded (basically 95% percentiles) in Europe [23].The harmonic order and levels are given in Table 1.The instantaneous value of the voltage signal for the subsequent examination is generated for the given target frequency signal.The phase voltage signal is calculated for a given frequency and harmonics level as, where, T s is the sampling time while K H i and A H i are the per-unit harmonic order and amplitude of respectively.Time domain signal is shifted by 120 and 240 degrees to make three phase voltage signals.
The f c of PLL impacts the delay in frequency measurements and hence the achievable accuracy.It should be sufficiently low to suppress the harmonics in the acceptable range (minimal measurement accuracy) and should be as high as possible in order to decrease the delay in the frequency measurement process.This leads to a parametric study in which f c is sequentially increased, while considering the harmonics.The harmonics according to the levels from Table 1 are added to the time domain voltage signal collectively.Subsequently, the PLL is used to measure the frequency from this time domain signal.Figure 2 shows the measurement error (base to peak amplitude of frequency signal) after the frequency step event.The measurement error increases as f c is increased.It can be observed that a value of 40.5 Hz leads to the violation of the maximum error limit (±1 mHz).Based on the result, the maximum permissible value f c is taken as 40 Hz for the stability studies.In order to study the impact of PLL input signal harmonics on its output and how filter can suppress them, a frequency sweep of harmonic content is performed.The frequency of harmonic signal is varied between 50 and 200 Hz.During this process, the amplitude of harmonic is set to 8% corresponding to maximum value mentioned in Table 1 and in par with other standards.The procedure for harmonic sweep and subsequent analysis of frequency measurement using PLL is outlined in Algorithm 1.This analysis is performed with the step change in the input frequency and the output frequency is analyzed as the response reaches a steady-state.Short time Fourier transform (STFT) of the output frequency signal is used to study its harmonics contents.The analysis yields to the results shown in Figure 3.It can be observed that harmonic content in range of 60 to 120 Hz impacts the PLL frequency measurement most significantly, while the low-pass filter effectively suppresses the harmonics content in the PLL output if the input signal harmonic order is more than 2.
Algorithm 1: Algorithm for PLL performance analysis considering harmonics.
1 Specify the harmonics order with amplitude and f c range. 2 for f c ∈ F c do (For each cut-off frequency value) 3 for h i ∈ H do (For each harmonics order h i )

4
Synthesize time domain current signals for each phase and add harmonics.

5
Define the f c of Butterworth filter and zero-crossing detection signal as reference.

6
STFT of the PLL output for instance using 100 fast Fourier transform points leads to time resolution of 1.6 s and a frequency resolution of 0.15 Hz.

7
Find the amplitude (of the frequency oscillations) and frequency (Hz) of harmonics in the output frequency of PLL.In order to understand the concept of inertia, it is fundamental to discuss the determining factors for frequency stability.The frequency is established by the generator with primary energy source rotating the turbine.The dynamics of frequency are governed by the classical swing equation [30], where, P m and P e are the mechanical input and electrical output power of the generator, D g is the generator damping torque coefficient due to windage and friction, ω is the frequency and M is the angular momentum.The angular momentum of the rotor (M = Jω m ) is related to the inertia constant (H) as, The inverters that aim to generate similar behavior as synchronous machines emulates the inertia analogous to angular momentum of the generator rotor.This is achieved by various inertial control topologies, a summary of which can be referred to Table 1 of [31].Among them, the most straightforward approach is the frequency-power response based inertia emulation method.The inverters employing such a method are generally referred to as VSGs.In this method, the VSG emulates the behavior of a synchronous generator by releasing/absorbing energy while responding to the magnitude and dynamics of the frequency given as, One can observe a striking similarity between Equation ( 4) and ( 2).Here, P VSG 0 denotes the primary output power of VSG, K D acts like frequency droop and emulates the damping of synchronous machine, while K I is the inertia emulating gain.The VSG measures frequency deviation (∆ω) and acts as a dispatch-able current source regulating its output energy.The aggregate small signal stability model of virtual inertia emulating VSG integrated in the power grid requires modeling of grid dynamics, frequency dependent loads and VSG.The following sub-sections formulate the aggregate analytical model of the frequency response of the system.

Frequency Dependent Load Model
A portion of the system load consists of rotating machines whose energy consumption is influenced by the input power frequency.Such that, if the frequency decreases due to an increase of load in the system, it reduces the load torque and thus their power consumption.As a result the load is reduced by the factor ∆P f .This counteracts naturally the load increase and is therefore referred in literature as self-regulating effect of loads.It is a common way to model frequency dependency of loads with a reference parameter D. It is expressed as ratio of normalized change in load to the frequency, The swing equation of synchronous generator given by Equation ( 2) for a small perturbation while neglecting the machine damping torque is given as, On substituting the inertia constant (H) from Equation (3) in ( 5) and converting to per-unit values of ω, P m and P e , the swing equation transforms to,

Frequency Domain Model of Power Grid
If the system is considered in the steady-state condition and change in the reference mechanical power (∆P m ) is zero, then a change in the frequency (∆ω) can be directly related to a change in load (∆P L ) as shown in Figure 4.The frequency domain model of the power grid from Equation ( 7) is given as, The system is influenced by the generator whose droop gain (R) is given as, The primary load frequency control (LFC) loop of the generator is modeled as, Here, τ g and τ t are the governor and turbine time constants.LFC changes the active power output of the generator with respect to the change in measured frequency, where power output of a generator is modeled by the time constants influencing the output of governor droop.The overall grid transfer function without G VSG is given as, Frequency domain block diagram representing LFC of a single machine connected to an infinite bus.

Frequency Domain Model of Virtual Inertia and Primary Response from VSG
The frequency domain model for the VSG is shown in Figure 5 where frequency measurement is carried out using PLL.Typically a rate limiter and saturation block are added to the output of frequency domain block.However, they are neglected while deriving the aggregate frequency domain model in order to assess the effective limit of virtual inertia gain.The evaluated value of K I is termed as maximum effective virtual inertia gain K I .The small signal model for PLL is defined by phase tracking dynamics and is typically expressed by a second order transfer function [24], Here, the K PLL p and K PLL i represents the proportional and integral gains of PLL, while ρ indicates an exponential delay term associated with the signal propagation delay.
The low pass filter associated with the PLL is taken as a sixth order low-pass Butterworth filter as discussed in Section 2.1 to remove the disturbances caused by harmonics.The control variable for the filter is its cut-off frequency (ω c = 2π f c ) and transfer function is given as, The dynamic response of PLL involves the measurement delay due to essential preprocessing of frequency identification and filtering, therefore these delays are to be considered properly for the small signal stability analysis.For this purpose, a Padé equivalent [32] has been used that combines both the signal propagation delay of PLL (ρ) and the low pass filter associated delay.In the Padé model, the time delay is represented by the transfer function of the nth order, i.e., The overall transfer function of PLL is therefore given as, Once the frequency is measured, it is fed to the LFC loop given as, The VSG time delay can be modeled as a cumulative time constant as, This leads to the overall transfer function of VSG, The VSG transfer function (G VSG ) and of the grid (G grid ) are combined to formulate the aggregate transfer function,

Single Machine Single Inverter Bus System
The single machine single inverter bus system (SMSIB) consists of a VSG with virtual inertia provision, generator and a load.The line diagram is shown in Figure 6 and the system parameters are given in Table 2.The frequency stability is evaluated for the step increase in the load from the initial value of 1 MW to 1.2 MW.The small signal stability analysis are performed for the ∆P ref mentioned in Table 2, however the stability results obtained are based on the system pole locations and are independent of ∆P ref .The frequency domain model of the network in reference to Equation ( 11) is given as, where, The actual and per unit quantities are represented by the capital and small letters, respectively.The per-unit quantities are calculated as, The VSG transfer function is evaluated according to Equation (18) where the parameters K PLL p and K PLL i are set to 10.The transfer function of the LFC loop is given as, A frequency response analysis is performed on the aggregate frequency domain model of the test case.During parametric study, the f c is varied between 20 to 40 Hz and K I is increased from 5 to 40.Meanwhile, the most critical pole is identified as the one having minimum damping ratio.The minimum damping ratio of 0.05 is required by many utilities as criteria for stability [33,34]; this value has been therefore taken as target.During the process, it is assumed that the VSG internal time constant is taken as 0.02 s.The results of this study are shown in Figure 7, where, for the cut-off frequency of 40 Hz the virtual inertia gain can be increased to 36 s (without violating the stability limit) from the value of 27 s calculated for f c of 20 Hz.The higher value of the upper limit of virtual inertia corresponds to more inertia capability that can be derived from the VSG.Alternatively, if the virtual inertia gain is set to 27 and 36, one can find the upper limit of f c , the result of which is shown in Figure 8.Here the objective is to find the minimum value of the cut-off frequency satisfying the 0.05 damping criterion, the lower f c is, the better harmonic rejection becomes.The minimum level of f c is identified as 20 and 40 Hz respectively, verifying the results from Figure 7.In the next step, we observe the influence of generator nominal power (P nom gen ), VSG nominal power (P nom VSG ) and f c on the maximum effective virtual inertia gain of VSG for the reference event.If the generator nominal power is increased, it is anticipated that more virtual inertia gain could be used for a given cut-off frequency of PLL as the generator is able to stabilize the system due to increased inertia impact (H is expressed normalized to S n ).Similarly, greater the nominal power of virtual inertia source implies lesser inertia gain limit as shown in the Figure 9.
Figure 10a shows the influence of f c on the virtual inertia gain limit of the VSG.From this comparison, we observe that K I can be increased by 33% if we increase f c from 20 to 40 Hz, due to reduction of delay.Further increase in the upper limit of K I is possible by increasing f c but as discussed in Section 2.1, it cannot be increased above 40 Hz in order to meet the harmonics rejection requirements.The maximum inertia contribution from VSG increases as the generator nominal power is increased as shown in Figure 10b.However, this trend saturates when the VSG capacity is less than 10% of generation capacity.

Multi-Machine System and Stability Analysis
The frequency response analysis of multi-machine systems aim at generalizing the results obtained for the SMSIB system.The multi-machine system considered consists of two VSG, two synchronous generators and a load.The single diagram is shown in Figure 11 and system parameters are given in Table 3.The equivalent droop effect of the generators is derived as, While, the system level transfer function is given as, where, r i is the ith generator droop gain, τ i is the first order governor delay, d is the damping introduced by the frequency dependent load and h i is the inertia of ith generator.The P final L is the final load.The per-unit quantities are calculated as, For the multi-machine system, if f c is set to 20 Hz, the corresponding maximum cumulative inertia gain is K I being the sum of K I of all VSGs has an upper limit of 13.6 s as shown in Figure 12.During the analysis, a load step of 1 MW has been applied to the system.If f c is increased to 40 Hz, then it can be observed that the gain can be increased to 18 s, i.e., 32% increase.The analysis of maximum cumulative virtual inertia gains with increasing share of VSGs and nominal power of generators is shown in Figure 13.The impact of f c on K I can be observed in Figure 14a.Here, the results are similar to Figure 10 and we observe a linear relationship between K I and f c .The combined inertia contribution from all VSGs (MWs) as function of the system inertia (MWs) is shown in Figure 14b where, all the bold symbols are the sum of respective parameters for generators and VSGs.In this analysis, the total generation capacity increases from 1 to 20 MW while total nominal power of VSGs is fixed at 1 MW.We can observe that for smaller microgrids, where the cumulative capacity of VSGs is at least 10% of the generation capacity, the upper limit of the inertia support has an increasing trend.Once the VSGs capacity is less than 5% of the generation capacity, the inertia support saturates.In light of the results from this study, we conclude that the maximum cumulative virtual inertia gain K I from VSGs is influenced by its nominal power and delay in the frequency measurement process.Moreover, the response to the reference incident during stability analysis is governed by the nature of frequency perturbation, overall inertia, cumulative droop gain of the power system and droop gains of participating VSGs.On the one hand, we observe that the maximum acceptable cut-off frequency of the PLL is determined by the desired robustness against harmonics.In general, the higher cut-off frequency leads to the smaller delay and hence greater is K I .

Conclusions and Outlook
This work addresses the key question in the industry about virtual inertia support from VSG in improving the small signal frequency stability.In this paper, a methodology is presented to identify the maximum permissible value of virtual inertia gain from VSG and identify the influencing factors.It is reported that by deploying an effective and faster PLL, the f c can be increased and more inertia can be emulated with a given VSG while offering significant level of robustness against harmonics.Due to a range of participating factors in the system, that depends on the load, inertia availability from generators and their droop gains, the upper limit of K I can be experimentally calculated using the presented methodology.In conclusion, with regards to supporting the inertia provision from VSG, the frequency measurement process plays a central role and the potential advantage of its improvement can be quantified by the frequency response analysis as presented in this paper.
Besides the academic interest, the work is also related to research and industry working in standardization committees in particular for the group "CENELEC WG7 Power frequency measurement for DER management" in TC8X [35] and the "RG-CE System Protection and Dynamics Sub Group" of ENTSO-E [36].New standards are necessary to normalize the behavior of generators aiming at providing virtual inertia to ensure a stable operation of the system.Testing methods are necessary to determine the accuracy, robustness, delays from different implementations of virtual inertia in non-rotating generators.This paper provides a method for evaluating the virtual inertia limitations from stability perspective due to delays in frequency measurement and processing.However, in order to reach to more generic conclusions, the multi-dimensional nature of factors influencing the maximum virtual inertia gain of a VSG requires further study of the influences of load type, inertia availability and its distribution in the network, droop gains of the generators and the contingency specific evaluation.This study can be extended to hardware in the loop platforms that consider the detailed dynamics of inverters and increase the insight on further limiting factors on the virtual inertia.Moreover, the approach can be useful to compare virtual inertia provided by grid following inverters to the grid forming control in the future [37].Recently, a comparison of various PLL approaches on the basis of accuracy and sensitivity to noise suggests Low-Pass Filter (LPF) PLL

Figure 1 .
Figure 1.Frequency estimation performance analysis for various PLL configurations.The zero crossing (ZC) detection is for reference while SOGI (SOGI-FLL) with-out and with low-pass filter and PLL implementation (1Ph-PLL) from Matlab are compared.The f c has been set to 20 Hz.

Figure 2 .
Figure 2. Base to peak amplitude of PLL output measured frequency around steady-state value as function of f c of PLL, when harmonics are considered collectively.

3 .
Virtual Inertia from VSG and Small Signal Stability Model

Figure 3 .
Figure 3.The figure shows impact of increasing harmonic order in the input of PLL to the harmonic contents of frequency estimated at the output.Each sub-figure has on x-axis the harmonic frequency order (K H ) in per-unit to the nominal frequency of 50 Hz.As the harmonic order is increased, the output harmonic frequency ( f H ) and the amplitude (A H ) are shown on the y-axes, where the dominant suggests amplitude and (L, R) are referring to left and right axes.The figures are arranged as, without filtering on left-side (a,c,e) and considering PLL low-pass filter on the right-side (b,d,f), respectively.The f c of PLL is set to 20 Hz.

Figure 5 .
Figure 5. Frequency domain block diagram of VSG shows the PLL stage followed by the control loops and G τ VSG representing VSG time-constant.The aggregate transfer function is given by G VSG .

Figure 6 .
Figure 6.System overview diagram of SMSIB system in test case.

Figure 7 .
Figure 7. Stability locus of most critical pole showing the maximum virtual inertia gain for given values of f c in SMSIB system.Damping ratio of 0.05 and 0.1 are shown in the background.

Figure 8 .
Figure 8. Stability locus of most critical pole showing the minimum value of f c for different values of virtual inertia gains.

Figure 9 .
Figure 9. Maximum virtual inertia gain limit for satisfying damping criteria as function of nominal power of the generator in (a) and VSG in (b) ( f c = 40 Hz).

Figure 10 .
Figure10.(a) Maximum virtual inertia gain limit as function of f c ; (b) Inertia contribution from VSG to the system in MW seconds as function of system inertia for increasing level of generator nominal power and while keeping H fixed at 5 s.

Figure 11 .
Figure 11.Single line diagram of the multi-machine system.

Figure 12 .
Figure 12.Stability locus of most critical pole indicating maximum value of cumulative virtual inertia gain of all VSGs for given value of f c in the multi-machine system.

Figure 13 .
Figure 13.Cumulative virtual inertia gain as function of percentage share of VSG nominal power in (a) and of generator nominal power is shown in (b).The percentage share is to the total generation and VSG nominal power capacity.The f c of PLL is set to 40 Hz.

Figure 14 .
Figure 14.(a) Cumulative virtual inertia gain as function of cut-off frequency of PLL and (b) Maximum inertia support from VSGs compared to inertia in grid for increasing level of generation capacity.

Table 1 .
Harmonic frequency contents from EN 50160.

Table 3 .
Parameters for multi-machine system.