A New Robust Control Strategy for Parallel Operated Inverters in Green Energy Applications

Abstract

This research work puts forward a hybrid AC/DC microgrid with renewable energy sources pertaining to consumer’s residential area for meeting the demand. Currently, the power generation and consumption have experienced key transformations. One such tendency would be integration of microgrids into the distribution network that is characterized by high penetration of renewable energy resources as well as operations in parallel. Traditional droop control can be employed in order to get an accurate steady state averaged active power sharing amongst parallel inverters pertaining to hybrid AC/DC microgrid. It is presumed that there would be similar transient average power responses, and there would be no circulating current flowing between the units for identical inverters possessing the same droop gain. However, the instantaneous power could be affected by different line impedances considerably and thus resulting in variation in circulating power that flows amongst inverters, especially during unexpected disturbances like load changes. This power, if absorbed by the inverter, could result in sudden DC-link voltage rise and trip the inverter, which in turn causes performance degradation of the entire hybrid microgrid. When the hybrid generators act as unidirectional power source, the issue worsens further. In this research work, we have put forward a new distributed coordinated control pertaining to hybrid microgrid, which can be applied for both grid connected and islanded modes that include variable loads and hybrid energy resources. Also, in order to choose the most effective controller scheme, a participation factor analysis has been designed for binding the DC-link voltage as well as reducing the circulating power. Moreover, to both photovoltaic stations and wind turbines, maximum power point tracking (MPPT) techniques have been used in order to extract the maximum power from hybrid power system when there is discrepancy in environmental circumstances. Lastly, the feasibility and effectiveness pertaining to the introduced strategy for hybrid microgrid in various modes are confirmed via simulation results.

1. Introduction

The smart grid is regarded to be the next generation of traditional power system that relies on two-way power flow as well as information exchange competencies [1]. Conventional grid includes power generation, transmission, distribution, control centers, and end users [2]. Key benefits pertaining to smart grid include two-way communication services as well as its ability to further use the renewable energy resources [3]. Smart microgrids, as part of the bigger power system, could provide electricity to smaller areas operating in island or grid-connected mode. Smart microgrid also offer consumers with enhanced control, monitoring, and new service technologies [4]. With regards to grid connection mode, a common connection point (PCC) connects the microgrid to the utility network. In islanded mode, microgrid functions separately and meets the control objectives autonomously [5]. It is difficult to control microgrid in islanded mode as multilevel control techniques are needed to ensure that voltage and frequency security constraints as well as other system’s technical constraints are efficiently met.

A stable and appropriate power management strategy is needed to operate hybrid microgrid networks in order to control load sharing amongst DC and AC power sources. Controlling power-sharing in DC and AC microgrids has been a focus for various studies. However, until now, there have been less efforts in power management amongst AC and DC power sources concerning hybrid systems. The reliable power-sharing method, which has been motivated by droop control pertaining to conventional synchronous alternator governor systems, has been used for power-sharing control in microgrids [6]. The frequency droop features involve working of different converters together for determining the corresponding reference power. With regards to the drooping control theory, load balancing includes all drooping control converters based on drooping slope [7].

Few researchers have studied the power-sharing pertaining to hybrid microgrid independent mode [8,9,10]. For sharing realization pertaining to AC microgrid, commonly employed parameters are active power and frequency (P–F), reactive power and voltage (Q–V) features [11]. With regards to DC microgrid, it is recommended to employ both current and voltage (I–V), active power and voltage (P–V) droop features [12]. However, coordination between two microgrids, i.e., hybrid microgrid and power-voltage (P–V) droop features, is usually employed when it comes to DC microgrids being part of hybrid microgrid [13]. In order to extend active power sharing pertaining to sources of both microgrids, also referred as proportional active power sharing, proper amount of active power between two microgrids needs to be exchanged by interlinking converter. With regards to this, determination of reference active power pertaining to interlinking converter is done by considering AC microgrid frequency as well as DC microgrid voltage for ensuring proper exchange of active power and achieving proportional power-sharing [14,15].

However, there are some limitations pertaining to the droop-based sharing strategy, which are also referred as primary control in hierarchical control structures. For instance, there could be inaccurate sharing due to variation in voltage magnitude for microgrid bus as a result of line impedances when voltage forms a part of realization pertaining to droop features [16]. Various efforts have been put for eliminating or mitigating this error [17,18]. In [19,20], various corrective actions have been applied by injection of further harmonic signals into the line. However, such invasive methods include the drawback of line current distortion. In [21], enabling signals are activated at specified intervals in order to achieve reactive power-sharing. Nevertheless, during this time interval, load fluctuations could result in stability problems as well as reactive power sharing errors. Authors in [22] have given more accurate sharing by making use of low bandwidth communication links. However, implementation of such methods is key to secondary control. In [23], obtaining of precise reactive power-sharing is done amongst AC sources pertaining to microgrid, in which radial connection of sources to common bus has been established. In the same study, regulation of output impedance pertaining to source has been done in a manner that total source impedance and source capacity are inversely proportional. However, this control structure could also result in low voltage quality due to output impedance voltage droop. In [24], adjustment is done to the reactive power-voltage slope pertaining to AC featured of AC source in order to compensate voltage drop that is associated with source’s output impedance. In this work, it has been regarded that across the source output impedance, voltage droop is autonomous with regards to source voltage differences, which could result in errors pertaining to reactive power sharing. In [25], similar study has been proposed by employing current-voltage droop concept in order to offer accurate current sharing amongst DC sources, which is not a common thing for hybrid microgrid applications. In [26], authors have put forward power derivative integral terms pertaining to conventional droop controllers in order to enhance dynamic response as well as decrease circulation current amongst parallel inverters. In [27], addition of an adaptive derivative term is done for droop control that also allows minimizing overshoot with regards to the current as well as enhance stability. In [28], in order to accurately adjust drooping parameters when faced with various load conditions, central controller with low bandwidth communication has been employed.

The well-established droop control and virtual synchronous generator control schemes in the literature are proposed in [29] to allow power sharing through voltage source inverters. This means that the voltage source inverter and the voltage magnitudes can be associated with both actual and reactive output power through decentralized power sharing. In spite of its precise actual power sharing, it is noticed that due to the mismatched impedances of the grid, the reactive, unbalanced, and harmonic power share precision cannot be assured [30].

The voltage drop over distribution feeders has been calculated in [31] for the purpose of balancing load voltage and it is added to the voltage reference provided by the Q–V droop control method which leads to a balanced load voltage. An improved droop control strategy is proposed in [32] based on self-tuning droop curves to adjust each phase droop curve with the load variance in real time and to mitigate the voltage unbalance of PCC.

In its entirety, the following drawbacks were associated in the reviewed literature:

• The degree of issue mainly relies on the size of DC-link capacitance, which is chosen traditionally in order to deal with certain filter requirements like to decrease voltage fluctuation cause by switching current.
• Even though the above-mentioned studies stressed on enhancing transient dynamics pertaining to the average power control, not a single study was able to address the instantaneous transient power impact on parallel inverters’ stability. Moreover, not a single study regarded the impact cast by mismatched line impedances on the damping and circulated power of the microgrid system.

To deal with the aforementioned problems, this study puts forward a new robust control strategy pertaining to green energy application’s parallel operated inverters when faced with variable load as well as generation conditions. The key contributions of this work can be summarized as:

• Hierarchical control structure has been recommended for both frequency and voltage restoration pertaining to grid-connected microgrid as well as islanded microgrid by accounting for the complete nonlinear system model, regardless of the parametric disturbances as well as uncertainties.
• The impact pertaining to mismatched line impedances has been evaluated based on the performance of parallel inverters that were supplied by various energy sources, and it also investigates the instantaneous circulating power responses against microgrid’s stability. In addition, the system has been analyzed by employing microgrid’s small signal state space model that includes three inverters. In order to uphold microgrid’s stability, we recommend employing two controller schemes that account for supplementary phase as well as frequency loops. The most effective controller scheme possessing the least action linked with DC- link voltage have been employed for participation factor assessment. Simulation was applied to validate the proposed controller.
• Distributed consensus-based control scheme has been suggested to precisely carry out power sharing when there is frequency restoration. The proposed strategy can be regarded as fully distributed method that allows distributing communication and computational tasks amongst local controllers by working in parallel, which is also more scalable, flexible and insusceptible to single point failure.
• Hybrid microgrid has been studied by accounting for typical energy storage systems, unbalanced and balanced load switch and nonlinear and linear loads offered via description of model containing all necessary data to deal with the aforementioned studies.

The remainder of this paper is structured as follows. Section 2 presents the description of the proposed system. Section 3 describes the microgrid control methods. Section 4 presents the mathematical model of the distributed hybrid generation system, while Section 5 elaborates on the dynamic analysis of system. Section 6 presents the proposed control strategy for PV/wind/storage energy system. Section 7 presents the proposed controller to limit voltage of DC link. Section 8 presents the analysis of the participation factor, followed by the controller design procedure outlined in Section 9. Section 10 presents results of the proposed system and Section 11, concludes the paper.

2. Proposed System Description

The studied hybrid microgrid is composed of a photovoltaic array, together with batteries, wind turbine, and AC loads. These systems form a hybrid microgrid interconnected to a low voltage distribution. Figure 1 illustrates the configuration of a hybrid power system for the distribution microgrid. The photovoltaic system is connected to the DC-bus through a boost converter and the battery is connected via a bidirectional DC to DC converter to the charging regulator.

3. Microgrid Control Methods

This section defines the microgrid control layer which constitutes the hierarchical (primary, secondary, and tertiary) control method. The device level is controlled by the base layer, with the most rapid reaction, while the control level is controlled by higher layers with a slower reaction. To ensure the quality of power and the stable operation of the microgrid network, the main features of the multi-layer control scheme include filtering, harmonic current sharing, and reactive power, as well as control loops that regulate the flow of reactive and active power must be well studied. The droop method is mainly implemented on hierarchical controls in order to support parallel operation pertaining to multiple voltage sources that share network loads as well as to maintain power quality [33,34,35].

3.1. The Primary Control

The layer of primary control, which has been derived based on droop control, is used for managing power provided by each converter via regulation of voltage and frequency. Different from typical traditional power systems, electronically-coupled power generation units lack inertial characteristics that confer stability to system when present in synchronization phase. In its place, it offers full control over system dynamics as well as transients by offering rapid responses. Thus, to enhance microgrid stability as well as coordination amongst voltage sources that operate in parallel, electronic simulation of inertia characteristics pertaining to synchronous generator is done in network-forming converters via the droop method, thereby regulating frequency as well as voltage amplitude proportionally that result in reactive and active power components [36]. Equations (1) and (2) represent primary layer by accounting for the droop control method:

_ωi=_ωnom−_{mi Pi}

_Vi=_Vnom−_{ni Qi}

Here,_ωirepresents inverter angular which corresponds to measured active power_Pi,_ωnomsignifies the network’s nominal frequency and_midenotes the coefficient with regards to droop-based active power control method._Virepresents the output voltage amplitude pertaining to converter that is proportional to the measured reactive power_Qi,_Vnomdenotes nominal voltage amplitude and_niis defined as the coefficient associated with reactive power.

The frequency component can be defined as a global variable equally generated among constantly-state grid-forming converters, resulting in equal active power sharing where all inverters have constant and equal drop coefficient of_mi. By comparison, the voltage portion is considered to be a local variable; hence, even though we believe that the_ni droop coefficient is identical, the existence of an uneven voltage/amplitude in different microgrid nodes can not lead to a perfect reactive power sharing [37]. In order to demonstrate this, active and reactive power sharing relations for two grid-forming converters have been represented in Equations (3) and (4):

_{m1 P1}=_{m2 P2}

_{n1 Q1}=_{n2 Q2}+_V2−_V1

In the absence of hierarchical control, varying active powers are injected by each converter to grid, depending on detected line impedance at their nodes. To produce sinusoidal voltage reference for every converter, the voltage amplitude and angular frequency from droop equations are used, as shown in Equation (5):

_Vref=_Visin(_ωit)

The feed forward component must be incorporated as the active power in an angular frequency droop equation to improve dynamic transient response to reference voltage, as shown in Equation (6).

_ωi=_ωnom−_{mi Pi}−_mipdPdt

This aforementioned droop compensation component was shown to have compatibility with the proportional derived controller (PD). Here, the derived term refers to the feed forward signal whereas the_{mi Pi}gain contributes to a faster transient response to active variations in power. This parameter improves damping characteristics, improves the position of system’s closed loop poles and provides higher variability.

To reduce the effects of circulating currents, lessen distortion, and ensure harmonic current sharing under unbalanced and non-linear loads, an impedance at converter output is emulated by the droop method using an additional closed-loop control. Equation (7) shows that the default impedance is included as a new variable in voltage reference signal based on the output current.

_Vref=_Visin(ph)−(_{Rv io}+_Lvd_iodt)

wherephcorresponds to the integral over time of Equation (6),_iois output current of converter,_Lvand_Rvare the reactive inductive and resistive components of impedance_Zv, as shown in Equation (8):

_Zv=_Rv+j_Lv

The impedance variable undergoes modification to primarily provide an inductive network and ensure controllability of both the reactive and active powers. These are seen in droop Equation (9) and Equation (10). The low pass filter is used to obtain average power from instantaneous power.

P=_ωcs+_ωc×_pins

where_ωcis cut off frequency. The power angle and instantaneous output active power are related and are given by the following equation:

_pins=_Vo×_VPCCsin(_δO−_δPCC)_Xo

where_Voand_VPCCare the output voltages,_δOand_δPCCare inverter phase-angles andPCCnodes,_Xois reactance of output. By perturbing Equation (1) we obtain Equation (11):

Δω=−_mp×ΔP

By perturbing Equation (10) and assuming constant_Vo,_VPCCwe obtain Equation (12):

Δ_pins=_Hp×Δδ

where_HP=_Vo×_VPCCcos(_δeq)_Xo, andΔδ=Δ_δo−Δ_δPCC

The_δeqrefers to the phase difference’s equilibrium point where perturbation is conducted. From Equation (9), to Equations (11) and (12), the Low Pass Filter is ignored. Furthermore, the transfer function that provides a relationship between output powerΔPand bus frequencyΔ_ωPCCis given by Equation (13):

ΔP=−_Hps+_{mp Hp}×Δ_ωPCC

The droop control utilized in primary control provides an acceptable power balance among converters that serve as grid-forming units. However, it introduces errors in the steady-state voltage and frequency. Figure 2 illustrates the virtual output impedance loop regarding the droop control, voltage, and inner current loops.

3.2. Secondary Control Layer

One of the features of the secondary control layer is to compensate for voltage and frequency deviations resulting from constant droop control. These values are then returned to the specific reference while sustaining the power sharing that was obtained in the primary layer. Figure 3 illustrates the primary and secondary control procedures.

In the case a secondary layer is not implemented, both frequency and voltage of microgrid network depend on the deviations in the primary layer due to virtual inertia and droop control impedances. Two secondary control techniques are available:

(1) Centralized control: It is based on a central controller that requires one to all communication scheme. The central controller computes all components errors, using measurements from the PCC (Point of Common Coupling), and sends the corrective term to the other converters. This is a robust technique against communication constraints, however with low reliability (network is completely central dependent), tolerance and flexibility in case of fault of the converter, demanding controller duplication due to its master-slave configuration.

(2) Decentralized control: It is a hybrid system based on local and central controllers, that requires one to all communication schemes. As in the former technique, the signals errors are common to the inverters obtained by the central controller, but further properties as the corrective term are autonomously computed by each local controller. This technique is sensitive to communication issues, since the inverters can produce different corrective terms if receiving information in distinct time due to non-synchronized internal clocks between controllers, an effect known as clock drift.

3.3. Tertiary Control

The power flow between the main grid and the microgrid is enhanced in this control layer by the independent sensing and demand of renewable power sources. The management of energy exported (imported) to (from) utility networks depends mainly on technological and economic problems such as tolerance of network errors, harmonic mitigation, interruption, and energy prices. Therefore, effective two-way communication is needed to control power flow [38].

4. Distributed Hybrid Energy Generation System Simulink SimPowerSystems was used to develop, simulate, and test the system-level hybrid microgrid model made up of distributed generation systems. There are three key components of the hybrid microgrid: wind turbine system, photovoltaic system, and the battery storage system. Because of the network complexity, simulation of the hybrid microgrid model is done using the SimPowerSystem phasor mode in order to simulate the scenarios faster. 4.1. Modelling of Photovoltaic Cell

Figure 4 shows the equivalent circumference based on a single diode of the photovoltaic cell that can be characterized as a diode, current source, parallel resistance and resistance of series. The current and voltage characteristics of the photovoltaic cell are defined in the mathematical standard equation given below [39]:

I=_Iph,cell−_Io,cell[exp((q(V+IR_(s,cell) )/akT)−1]︸_Id,cell−V+I_{Rs,cellRp,cell}

where,_{Iph, cell}is the photocurrent (A) of the photovoltaic,_Io,cellis the saturation current or reversed leakage of the photovoltaic, k is constant of Boltzmann’s 1.38 × 10⁻²³ J/K, q is electron charge (1.602 ×¹⁰⁻¹⁹C),_Rp,cellis parallel resistance of photovoltaic,_Rs,cell, is series resistance of photovoltaic.

A photovoltaic cell model is considered an ideal one if the parallel and series resistances of the photovoltaic cells are not accounted for. Figure 5 illustrates the current and voltage curve for the ideal photovoltaic cell as derived from Equation (14) [40,41].

4.2. Modelling of Photovoltaic Module

A photovoltaic module is previously stated to consist of photovoltaic cells joint in parallel and in series. The basic mathematical equation is therefore obtained based on Equation (14) and the description for the voltage and current characteristics of photovoltaic module will be as follows [42]:

I=_Iph,cell−_I0[exp(V+I_RSa Vt)−1]−V+I_RSRp

where:_{Iph, cell}is the photocurrent (A) of the photovoltaic,_Iois photovoltaic reverse leakage current,_Vtis photovoltaic thermal voltage,_Rpis parallel resistance,_Rs is photovoltaic series resistance, Equation (15) generates the current and voltage curve as illustrated in Figure 6.

Photocurrent for a photovoltaic (_Iph ) is influenced by the quantity of the solar irradiance falling on module as well as temperature of the photovoltaic cell. This corresponds to Equation [43]:

_Iph=G_Gn(_Iph;n+_KiΔT)

where:_Iph;nrepresents the photocurrent when it is under the nominal condition (typically a 1000 w/m² irradiance and 25 °C temperature),ΔTrefers to the difference between the PV cell’s actual temperature (T) and the nominal temperature (_Tn) measured in °C,_Gnis the nominal irradiance (1000 w/m²),Gis the measured solar irradiance in w/m²_Kiis temperature coefficient. The open voltage circuit (_Voc) is affected by the cell temperature based on the following equation:

_Voc=_Voc;n+_KvΔT

where:_Kvis voltage of open circuit temperature coefficient and_Voc;nvoltage of open circuit as seen under nominal conditions. The following equation can be used to obtain the diode saturation current (_I0):

_Io=_Isc;n + _KiΔTexp(_Voc;n + _KvΔTa _Vt)−1

where_Isc;n is short circuit current under the nominal conditions [44].

In this paper, the electrical parameters of the SPR-305E-WHT-D solar is extracted and used to simulate this modular model. These parameters are listed in Table 1.

4.3. Battery Storage System

Battery storage systems are used for storing extra energy produced by renewable energy generation systems. However, if the renewable energy generation system does not generate enough energy, the load demand will be met by discharging battery bank. Battery model is as follows [45]:

SO_Cbat=100[1−(1_Qbat×∫0t_ibat(t)dt)]

_BAH=13600∫0t_ibat(t)dt

whereSO_Cbatrefers to the battery state ofcharge (%),_Qbatrefers to the maximum battery capacity (Ah),_ibatrepresents current of battery

4.4. Wind Turbine Model

One can design wind turbines (WT) as aerodynamic input torques that drive DFIG (Doubly Fed Induction Generator). Figure 7 illustrates the wind turbine’s characteristic output curves for varying wind speeds. One can then express mechanical power (_Pm ) extracted from wind turbines as [46].

_Pm=12ρ_{At Cp}(λ. β)_Vw3

where_Atis the area swept out by turbine blades,ρis air density in kg/m³,_Vwis wind speed,_Cp(λ. β)is performance coefficient of the wind turbine. The wind turbine performance factor (_Cp) depends on the blade aerodynamics and depicts the wind turbine’s efficiency. One can describe the performance factor as follows:

_Cp(λ.β)=0.5176(116_λi−0.4β−5)^e−(21_λi)+0.0068λ

1_λi=1λ+0.08β−0.035^β3+1

λ=ωRV

The PMSG wind turbine mathematical model is expressed in terms of the rotating reference frame d–q, and according to the following equations:

_Vd=−_{Rs Id}−_Ldd_Iddt+_{ωe φq}

_Vq=−_{Rs Iq}−_Lqd_Iqdt+_{ωe φd}

_ωe=p×_ωm

where:_Iqis q-axis current,_Idis d-axis current,_Rsrepresent resistant stator,_Lqrepresent inductance of stator in q-axis,_Ldrepresent inductance of stator in d-axis,prepresent polar number,_φqis flux in q-axis,_φdis flux in d-axis,

_φd=_{Ld Id}+_φf

_φq=_{Lq Iq}

The electromagnetic torque is:

_Tem=32×p[(_Ld−_Lq)_{Id Iq}+_{Iq φf}]

In this paper, the electrical parameters of the wind turbine are extracted and used to simulate this modular model. These parameters are listed in Table 2.

4.5. Concept of Stability

In general, stability denotes the power system ability to preserve synchronicity during disturbances. Stability issues are classified in terms of rotor angle stability, frequency stability and voltage stability. While stability of the rotor angle gives information about balance between the electromagnetic and mechanical torques required by a rotor generator to remain in synchronousness, the voltage and frequency stability refers to the power supply and demand balance of a network. Although the synchronous generator ‘s dynamic behavior plays a large role in preserving power system stability, rotor angle stability in a microgrid with converter-interacted distributed generators (DGs) is no longer a problem as DG inertia is inadequate, which allows major differences in the activity of a microgrid from the conventional grid [47]. Although the frequency and voltage stabilities of the microgrid are evaluated in an island mode, only the first is regulated by the grid-tied operation because of small impact on grid frequency, since voltage acts as a local variable. The stability of a microgrid is essentially dependent on its modes of operation, control topology, types of distributed generators, and network parameters [48], as well as small or large disturbance testing.

The small disturbance analysis involves evaluating the stability of a microgrid by linearizing the system around its operating points using an analysis of its own value, Nyquist or Routh Hurwitz stability criteria. Several small signal analyses that deal with load fluctuations, and influence of droop gains and variations in line impedance on frequency and voltage characteristics are conducted in [49].

A large disturbance stability analysis evaluates the capacity of a system to sustain the integrity of the network after a significant change including malfunction, service switching, or loss of generators or loads. Although it is widely used for conventional power systems, there is little information in the literature about a microgrid.

Microgrid stability can be improved in several ways, including by using control strategies, reactive power compensation, warehousing systems and load shedding. Although a high droop gain can enable a proper load sharing, it may cause instability in microgrid while a microgrid stability can be improved by optimizing droop gains and using a supplementary loop [50].

4.5.1. Microgrid Stability Classification

The main power supply in a microgrid is inverter-interfaced DGs. The dynamic method of microgrid stability is more complex than the stability of the conventional grid. The time and physical characteristics of the instability process are also fully taken into account. The classification of the microgrid is illustrated in Figure 8.

4.5.2. Stability of Islanded Microgrid

While microgrid is in islanded mode, the balance of power source and load variability should be maintained. A clear understanding of various types of instability is important for the satisfactory operation and design of the microgrids.

The properties upon which the stability of the microgrid depends, mainly include the configuration of the microgrid, the capacity of the DGs, and their related control strategy. Primary energy sources response speeds such as photovoltaic, and wind turbines are slow, but the power output is controlled by the interfaced inverters, which act very fast [51]. If load fluctuation occurs, controllable units could make corrective actions quickly and have a significant impact on stability.

In the same way as the grid-connected study of microgrid stability, studies on small signal stability in an islanded microgrid focuses primarily on effect of control gains, load fluctuations, and line impedance changes on the performance and frequency of the distributed generator. Islanded microgrid transient stability primarily consists of the effect on the operating phase of the microgrid under significant disruptions such as short-circuit fault, open-circuit failure, loss of distributed generators, and load. The ability of a microgrid to maintain stability can be addressed when major disturbances occur with the transient stability analysis process. 5. Dynamic Analysis of the System

This paper utilized the hybrid microgrid model to evaluate the system stability during abrupt changes in the load and with significant line impedance mismatch. The state-space model isx˙=ax , where a represents the state-space matrix. Therefore, each inverter can be modelled by its Thevenin equivalent circuit as shown in Figure 9. The equivalent circuit model consists of an ideal voltage source, V_O and an output impedance (_RO+jω_LO ). Figure 10 illustrates an islanded configuration of a microgrid. Each inverter is fitted with two cascaded space models. The first model shows the power-sharing controller while the second model represents the output impedance.

5.1. System Power-Sharing

Contribution of each inverter can be calculated with instantaneous output power:

_pins=32(_{vod iod}+_{voq ioq})_qins=32(_{vod ioq}−_{voq iod})

These calculated powers should pass through a low pass filter to attenuate the high frequency components. The average power is estimated:

P=_ωcs + _ωcpinsQ=_ωcs + _ωcqins

Since the d-component is chosen to be in line with the phase A voltage of the three-phase network, the q-component is zero, so each device has an optional droop control to share the power in the dq as:

ω=_ωo*−_mpP_vod=_Vo*−_nqQ_voq=0

where_ωo*and_Vo*are the nominal frequency and nominal voltage references,_mpis coefficient of frequency, and_nqis coefficient of voltage. The variables in Equation (33) should be disturbed around the equilibrium points when constructing the small signal model. The following equation represents the linearization of Equation (31):

Δ_pins=_Iod×Δ_vod+_Ioq×Δ_voq+_Vod×Δ_iod+_Voq×Δ_ioqΔ_qins=_Ioq×Δ_vod−_Iod×Δ_voq−_Voq×Δ_iod+_Vod×Δ_ioq

where_Vod, _Voq,_Iod,_Ioqare the equilibrium points for the output voltage and current. Symbol“Δ”denotes a signal variation in a state. Applying signal variations in Equations (32) and (33) gives:

sΔP=_ωcΔ_pins−_ωcΔPsΔQ=_ωcΔ_qins−_ωcΔQΔω=−_mpΔPΔ_vod=−_nqΔQΔ_voq=0

For any inverter the phase state calculated as:

sΔδ=Δω

The phase state at the PCC assumes a fixed frequency. This may be true in stiff grid connected mode. The rising frequency in islanded mode, however, may have minor variations in the system. To take that into account, Equation (36) is updated to:

sΔδ=Δω−Δ_ωcom

whereΔ_ωcomis a small variation of the common frequency of the system. It is used to convert any inverter state from its local context into a common framework so that all subsystems can be merged into one device.

By arranging the results of Equations (34)–(37), then the sub-state space model is derived by:

[Δδ˙ΔP˙ΔQ˙]=_AP[ΔδΔPΔQ]+_BPv[Δ_vodΔ_voq]+_BPi[Δ_iodΔ_ioq]+_BPωcom[Δ_ωcom]

where:

[Δω]=_CPω[ΔδΔPΔQ], [Δ_vodΔ_voq]=_CPv[ΔδΔPΔQ], _CP=[_CPωCPv], _AP=[0−_mp00−_ωc000−_ωc],

_BPv=[001.5_{ωc Iod}1.5_{ωc Ioq}1.5_{ωc Ioq}−1.5_ωc._Iod], _BPi=[001.5_{ωc Vod}1.5_{ωc Voq}−1.5_{ωc Voq}1.5_{ωc Vod}], _{BP ωcom}=[−100]

_CPω=[0−_mp0], _CPv=[00−_nq000]

5.2. Output Impedance Model

The output impedance model is a simple R-L network between two voltage nodes that assigns the output current to the inverter. The source voltage node is_voand the sink voltage node is_vn. The sub-state space model is provided in terms of the current state of the output:

[Δ_iod˙Δ_ioq˙]=_AZ[Δ_iodΔ_ioq]+_Bz1[ΔωΔ_vodΔ_voq]+_Bz2[Δ_vndΔ_vnq]

where:

_Az=[−_R0Loωo−_ωo−_RoLo], _Bz1=[_Ioq1_Lo0−_Iod01_Lo], _Bz2=[−1_Lo00−1_Lo]

5.3. Distribution Lines Model

The distribution lines in Figure 10 are another series of R-L networks between the node voltages_vn,ito_vn,h. In this paper, we call any source node for any line current as_vsourceand the sink node as_vsink. For each line state-space is:

[Δ_iLineD˙Δ_iLineQ˙]=_AN[Δ_iLineDΔ_iLineQ]+_BNω[Δω]+_BNvsource[Δ_vsourceDΔ_vsourceQ]+_BNvsink[Δ_vsinkDΔ_vsinkQ]

where:

_AN=[−_{RLineLLineωoωo}−_RLineLLine], _BNω=[_ILineQ−_ILineD], _BNvsource=[1_LLine001_LLine], _BNvsink=[−1_LLine00−1_LLine]

5.4. Loads Model

In the same way, the load small signal model is given by:

[Δ_iLoadD˙Δ_iLoadQ˙]=_AL[Δ_iLoadDΔ_iLoadQ]+_BLω[Δω]+_BLv[Δ_vnD,iΔ_vnQ,i]

where:

_AL=[−_{RLoadLLoadωoωo}−_RLoadLLoad], _BLω=[_ILoadQ−_ILoadD], _BLv=[1_LLoad001_LLoad]

5.5. Model of Inverters

Every inverter contains the state-space models of its power-sharing controller, Equation (38), and impedance model as shown in Equation (39). The inputs are node voltages,Δ_vnDQ, in common dq-frame. The output state is essentially the inverter’s output currentΔ_ioDQin the common dq-frame. The first inverter should give an extra output that is a frequency of the common dq-frame,Δ_ωcom. A state-space model of one inverter,i,can be built as in Equation (42) and the combined inverters models as in Equation (43):

[Δ_xinv,i˙]=_Ainv,i[Δ_xinv,i]+_Binv,i[Δ_vnDQ,i]+_Binvacom,i[Δ_ωcom,i]

where:

[Δ_ωiΔ_ioD,iΔ_ioϱ,i]=[_{Cinvω,iCinvc,i}][Δ_xinv,i], [Δ_xinv,i]=^{[Δ_δi Δ_Pi Δ_Qi Δ_iod,i Δ_ioq,i]T}

_Ainv,i=[_AP+_BPV⋅_CPvBPiBZ1⋅_CP+_BZ2⋅_CTvAZ], _Binv,i=[0_BZ2·_DTv], _{Binv ωcom,i}=[_{BP ωcom}00]

_Cinvω,i={[_CPω 0 0]i=1_[0]1×5i≠1}, _Cinvc.i=[_CTiDTi]

In the same way, the equations for all inverters and the combined model can be derived:

x=[ Δω Δδ ΔP ΔQΔ_VDC Δ_Iod Δ_Ioq︸PhotovoltaicΔω Δδ ΔP ΔQ Δ_VDC Δ_Iod Δ_Ioq︸Battery

Δω Δδ ΔP ΔQ Δ_VDC Δ_Iod Δ_Ioq︸Wind TurbineΔ_ILineD Δ_ILineQ︸Distribution linesΔ_ILoadD Δ_ILoadQ︸Loads ^]T

Figure 11a shows the microgrid dominant eigenvalues when three inverters generate a 1 mH output inductance given negligible line impedances (as shown in Figure 10), and when the value of_mPchanges from5×¹⁰⁻⁵to5×¹⁰⁻³ . The traces of the eigenvalues are identical and inverters are expected to behave in a similar manner. Figure 11b illustrates the same trace but with a different line inductance,L_line1=1×¹⁰⁻³HandL_line1=2×¹⁰⁻³H. It is apparent that when the band m_p is similar, the corresponding eigenvalues represent varying damping ratios, with poles observed to shift from their original locations. This can lead to an increase in overshoot and oscillations. In both figures, the eigenvalues at_mp=3×¹⁰⁻³are emphasized to signify that the same droop gain generates varying damping ratios of output power responses. It is desirable to increase_mPin order to obtain high sharing accuracy. This, however, degrades stability.

Therefore, when the load varies, even though equal design parameters are chosen, the paralleled inverters do not have identical transient replies. Compared to low-damped inverters that approach steady-state quicker, high damping inverters have a slower response to change. When the new steady-state value of active power is near to zero, power is transferred during transit from highly-damped inverters to lightly-damped inverters. It leads to a charging of the DC link condenser and pushing of DC voltage to higher values that can lead to a system trip.

Simulink was used to construct the hybrid microgrid model. When the load undergoes a sudden change from100% to 0%, the corresponding voltage and active power response of the DC-link are obtained. Figure 12 shows the response generated when using identical distribution lines,L_Line1= L_Line2=0 . The responses were observed to be quite identical and with only a slight effect on the DC-link voltages. In Figure 13, there is a different damping of the power response and a certain amount of power is imported. If DC to DC converter is bidirectional, similar to that in battery systems, power will sink from the DC link condenser. If the converter from DC to DC has a unidirectional system, similar to those in wind turbines and photovoltaic systems, one cannot bind DC voltage. This is shown in Figure 13b.

6. Control Strategy for PV/Wind/Storage Energy System

To ensure a stable response in the hybrid microgrid network during normal operation, Figure 14 and Figure 15 illustrates the proposed control methods. This paper conducted an investigation to determine which state has the greatest effect on DC-link voltage state. This can lessen the controller effort required to limit voltage of DC-link. As a result, small gain to controller is used, thus preserving stability.

6.1. Controller of the Converter Interfacing Photovoltaic

A single-phase boost is used to up the voltage from the panel and control maximum power point tracking. Input current (_IPV) sensing happens before the measurement of input capacitance (_Ci) in addition to the panel voltage (_VPV ). MPPT algorithm uses these two values. The MPPT algorithm determines a reference point, which, when maintained at the panel, provides maximum power to be extracted from the PV system. As depicted in Figure 14, MPPT is achieved using the voltage of the outer loop and the current flowing in the inner loop. Hence, the signs for the reference outer voltage compensator and feedback are reversed. The converter output is not adjusted. A voltage feedback mechanism provides input to the internal voltage comparators, which use pulse-width modulation (PWM) in an overvoltage situation to control the voltage from rising beyond the components’ rated voltage.

6.2. Controller of the Converter Interfacing WT

The wind turbine generator control scheme contains an MPP extractor for use in the independent variable wind turbine. It also features a DC bus voltage control and a PMSG (permanent synchronous magnet generator). The power converter makes correct adjustments so that it can deliver the maximum available wind turbine power through the rectified DC voltage and current drawn from that rectifier output. In addition, the loading station operates in standalone mode. Generation of the power reference is done to compare DC-link actual and reference values. PI controllers provide the power reference used to adjust the generated control signals. To keep the voltage of the microgrid at its desired value, each converter’s PWM signals are controlled regardless of loads and wind speed variations. A boost converter is used to control the DC voltage rectified to a higher voltage level to supply the provided power to the microgrid of the station. Figure 15 illustrates the converter’s control scheme. The parameters of PMSG are listed in Table 3.

6.3. Controller of Bidirectional Converter Interfacing Batteries

In the proposed setup, as depicted in Figure 16, the bidirectional converter has an output filter capacitor (_Cdc ), two switches (S1 and S2) and a high frequency inductor (L) that permit bidirectional flow of current. The power management system consists of a two voltage controller with appropriate limitation blocks to provide the necessary power flow under varying conditions. The controllers output a reference current for the storage of energy. The first aspect of control is DC-bus voltage regulation, while the other controllers regulate battery voltage. To facilitate better power management in the microgrid, backup energy storage is a part of the system, which comprises a battery connected to the DC bus using a bidirectional converter. The converter serves multiple purposes among which are: acting as a battery charge regulator during the grid-connected mode while serving as a boost converter to use battery energy and feeding it to the microgrid during insufficient power output from the photovoltaics and wind turbine to run the loads (both AC and DC) in islanded mode. During islanded operation, the optimal operating scenario is the equivalence between load energy requirement and PV and WT power generation, where there is no need for the converter to process energy. Figure 16 depicts the bidirectional control structure and a simplified phase of the power from the converter.

7. Proposed Controller to Limit DC-Link Voltage

To limit voltage of DC link during the normal operation of a microgrid, the DC-link voltage after exceeding a triggering value (_Vtr) will affect the frequency state which influences the active power.

To assure stable responses among hybrid microgrid networks in normal operations, Figure 17 describes two new control method proposals. When the voltage of DC-link surpasses a triggering value, the_Vtr, a loop is started. This supplementary loop influences either the (a) phase state or (b) the frequency. In this research, an investigation is conducted to determine the state that has key effect on DC-link voltages. In this way, the controller effort needed to control the DC link voltage and use smaller control gains, which in turn preserves stability.

Through analysis of the first control track, as shown in Figure 17a, wherein controller output signaling is used to manipulate frequency state, the small-signal output frequency described in Equation (1) is thereby derived:

Δω=−_mp×ΔP+_Pdc×Δ_vdc

where inΔ_Vdc represents the smaller signal state of voltage DC-link, which is derived as:

Δ_Vdc=−_kdcs×Δ_pins

where_Kdcis the linearizing factor of the DC voltage and the negative sign indicates the negative flow of power. In the assumption that instantaneous and average powers are equal and through substitution of Equation (45) for Equation (44), we derive Equation (46):

Δω=−(_mp+_Pdc×_Kdcs)×ΔP

Therefore, the phase state is obtained as Equation (47):

Δδ=−(_mps+_Pdc×_Kdc^s2)×ΔP

For the control system shown in Figure 17b, the frequency and phase-states are as follows:

Δω=−(_mps+_Pdc×_Kdc)×ΔP

Δδ=−[_mps+_Pdc×_Kdcs]×ΔP

An evaluation of Equation (46) with Equation (48) shows that the technique demonstrated in Figure 17a operates as the Proportional Integral (PI ) controller of the output power. However, this approach may not offer faster action against introduced power. Conversely, the technique demonstrated in Figure 17b operates as the proportional derivative (PD) controller. In this approach, the power derivative is not applied directly, although the action is realized in the proposed loop. The derivative term is well-known for providing the fastest response needed for limiting any introduced DC-link voltage and power rises. Additionally, when compared to Equation (49), an added pole is obtained at origin in Equation (47), thereby introducing a^90o phase lag that results in lower stability margins. In accord with the above discussions, the technique demonstrated in Figure 17b is shown to perform better overall.

8. Analysis of the Participation Factor

Participation analysis was performed to attain further insight into the excellence of the technique, as indicated in Figure 17b, over that demonstrated in Figure 17a. The participation factor_Pkiin Equation (50) of eigenvalues on a state represents a measure of the influence of the state regarding the eigenvalue. Conversely, the participation factor_πkiin Equation (51) of eigenvalues on a state represents a measure of whichever mode forms the majority of the response.

_Pki=_rki×_lki

_πki=^[Re(_lki)]2Re(_lki)×^[Re(_lki)]T

wherekdenotes the_kthstate,idenotes the_ith mode. The eigenvalues for the system are calculated from Equation (43) as in Table 4.

As the aim is to control the DC voltage state variable, we are looking for modes with key effects on the active power P state variable, given the direct relationship ofP with DC-link voltage. Table 5 describes participation factors for eigenvalues on the system states. Eigen values_{λ20 λ21 λ23 λ25} retain dominant influence on power responses. Operation of these modes can rework the active power response through transients, and thus, DC connection voltages can be reduced. Further analysis of the factor of participation is carried out to determine which states have key effects on these modes. Although the frequency is not a direct state, it includes a power output via a scalar function in relation to the power state as stated in Equation (11). Thus, Table 6 describes the phase and participation of the power state in the modes investigated, as was established early on. The results indicate that phase state has greater effect than power state on mode, which works to shape output power and consequently DC-link voltage. This validates similar conclusions made in the last section, in that the method displayed in Figure 17b performs better than that displayed in Figure 17a, for it results in smaller controller gains required to effect similar actions.

9. The Controller Design

After consideration of the method described in Figure 17b, and by perturbing Equation (1) and then computing the phase state variable, we derive:

Δδ=Δωs+_Pdc×Δ_vdc

SΔδ=−_mpΔP−_Pdc×_kdc×Δ_Pins

Equation (53) has been used in the Equation scheme state-space matrix Equation (43), whereinΔPandΔ_Pins are redefined. Figure 18 describes eigenvalue when proposed control gain varies as follows:0.0001>_Pdc>0.002.

From Figure 18, we can see following:

• The increase in the proposed gain of control (_Pdc)affects targeted eigenvalues_λ20,_λ21,_λ23 , in alignment with that displayed in Table 6. Nevertheless, certain other eigenvalues are similarly affected, namely_λ14,_λ15,_λ16,_λ17, as couplings between eigenvalues and state variables cannot be easily analysed. In addition, the targeted_λ25demonstrates negligible change, for it is mostly subject to determining filter bandwidth, in which_λ18,_λ19,_λ24, are not influenced.

• Increased_Pdcmovement targets eigenvalues farther to the left and therefore increases damping of system. Nevertheless, high frequency eigenvalue similarly move towards imaginary axis, whereby at certain higher gains, system instability ensues.

10. Simulation Results

The suggested hybrid microgrid network comprising three inverters, as displayed in Figure 1 in SimPowerSystem, is simulated in order to confirm the performance of the suggested control strategy. Inverters and converters are characterized in models of renewable energy source, as photovoltaics, wind turbines, and batteries. The parameters of the simulation are given in Table 7. DC-link voltage is acquired whenever loads instantaneously change from 100% of 6 kW down to 0%, with the resulting power response displayed in Figure 19a. Clearly, the response remains well-damped due to the performance of the proposed controller that activates automatically whenever DC voltage rises above_Vtr . This result is worthy of comparison as shown in Figure 13, which represents the same system except that the controller is not activated. Furthermore, DC-link voltage in Figure 19b is constrained and remains below trip levels, which validates the efficiency of the suggested approach.

Figure 20 illustrates average active energy and response of DC voltage when proposed controller’s gain is P_dc = 0.001. The system is unstable and oscillation frequency is 314 rad/second, this is consistent with the prediction in Figure 18 as the system becomes unstable and oscillation frequency is 312 rad/second.

To validate the recommended approach in view of the excellence of technique indicated in Figure 17b over that demonstrated in Figure 17a, Figure 21 describes response when control method, as displayed in Figure 17a, is implemented with_Pdc=0.005 . The gain was selected for its similar capacity to tuning DC-link voltage, as illustrate in Figure 19b. Clearly, the gain exceeds that used in the suggested method, as shown in Figure 17b, which therefore confirms that the suggested method requires less effort in regulating DC voltage. Furthermore, DC-link and power voltage responses can be highly oscillatory in comparison to that described in Figure 19.

11. Conclusions This research investigated the effect of mismatched line impedances on the stability of DC-link voltages under sudden load changes, as well as the performance of parallel inverters operating in hybrid microgrids. Studies showed that comparable droop gains do not assure similar power responses during transients in every case. If an appreciable line impedance mismatch exists, circulating power transients may degrade stability. The suggested technique was found to perform more efficiently than the frequency loop approach in stabilizing DC-link voltage states in microgrids. Any stable-state errors appeared solely during instances of zero-load transition. Nevertheless, this does not present a major issue when the load reconnects, as the level is sufficient to discharge the excess energy and thereby maintains zero error. The results of the simulation demonstrated also the improved performance and effectiveness of the suggested control approach. Additionally, the suggested control strategy exhibited much robustness against larger load variations. Furthermore, in an effort to avert the drawbacks of standard droop control, an enhanced droop control strategy derived from integrator current-sharing was developed. The MPPT method was used in both wind turbine and photovoltaic stations towards extracting maximum amounts of power from such hybrid power systems. The new strategy aimed to unify the control loops for all inverters in each mode of operation and for different energy sources to be DC-voltage regulators. This structure immunized the stability of the DC-link voltage states and the overall reliability of the microgrid against sudden unplanned islanding. The control scheme was analyzed and tested by simulation. As future works, the following points are suggested: (I) implementing the presented strategy on a real network and comparing the results (II) the proposed control scheme can be extended to address the network-related issues, such as data loss, packet jamming, link failure, cyber-attack, etc., using a cyber-physical systems framework.

Parameters	Value
Series connected modules	7
Parallel string	1
Voltage of open circuit(Voc)	6.42 (V)
Maximum voltage(Vmp)	54.7 (V)
Temperature coefficient of(Voc)	−0.27269 (%/°C)
Short-circuit current(Isc)	5.96 (A)
Maximum currentImp	5.58 (A)
Temperature coefficient of(Isc)	0.061745 (%/°C)
Shunt resistance(Rsh)	269.5934 Ω
Series resistance(Rs)	0.37152 Ω
Diode ideality factor	0.945
Diode saturation curentIo	6.3×10−1(A)
PV type	SPR-305E-WHT-D
Number of cells	96

Parameters	Value	Unit
Total inertia (J)	0.0119	Kg·m2
The state friction	5	N·m
Viscous friction coefficient (F)	0.001189	N·m·s
Inductance of d-axis	0.0082	H
Resistance of Stator	0.42	Ω
Inductance of q-axis	0.0082	H
Number of phases	3	N
Wind speed	12	m/s
Performance factor(Cp)	0.48	C
Pitch angle (β)	0	degree
air density(ρ)	1.225	Kg/m3
Base rotational speed	1.2	pu
ωe	376	rad/sec

Parameters	Values	Units
inductance of d-axis	0.0082	H
Resistance of Stator	0.42	Ω
inductance of q-axis	0.0082	H
Vpk/Krpm	98.7	-
Number of poles	8	-
Mechanical time	0.04	-
Inertia Moment	8 × 10−3	N–msec2

Eigenvalues	Locations	Eigenvalues	Locations
λ1	0.0	λ14	−151 + 318 i
λ2	0.0	λ15	−151 − 318 i
λ3	0.0	λ16	−90 + 318 i
λ4	−224,683,637 + 316 i	λ17	−90 − 318 i
λ5	−224,683,637 − 316 i	λ18	−4
λ6	−224,683,637 + 316 i	λ19	−6
λ7	−224,683,637 − 316 i	λ20	−20 + 7 i
λ8	−2,546,836 + 316 i	λ21	−20 − 7 i
λ9	−2,546,836 − 316 i	λ22	−23
λ10	−693,417 + 316 i	λ23	−25
λ11	−693,417 − 316 i	λ24	−31
λ12	−31,719 + 316 i	λ25	−31
λ13	−31,719 − 316 i	λ26	0

Mode/State	Δδ1	ΔP1	Δδ2	ΔP2	Δδ3	ΔP3
λ20,21	0.34	0	0.64	0	0.03	0
λ23	0.34	0	0.03	0	0.61	0
λ25	1.01	0	0	0	0	0

State	λ14, 15	λ16, 17	λ18	λ19	λ20, 21	λ22	λ23	λ24	λ25
Δδ1	0	0	0	0	0	0	0	0	0
ΔP1	0.002	0	0.07	0.1	0.2	0.1	0.6	0.1	0.5
ΔQ1	0.008	0	0	0.07	0.1	0.42	0.09	0.5	0.07
Δiod1	0.1	0	0	0	0	0	0	0	0
Δioq1	0.1	0	0	0	0	0	0	0	0
Δδ2	0.001	0	0.08	1.2	0.5	0	0.03	0	0
ΔP2	0	0	0.01	0.2	0.5	0	0.07	0.01	0.35
ΔQ2	0.014	0	0	0.14	0.2	0	0.01	0.3	0.07
Δiod2	0.2	00	0	0	0	0	0	0	0
Δioq2	0.2	0	0	0.009	0	0	0	0	0
Δδ3	0	0	1.3	0.083	0.04	0.15	0.45	0	0
ΔP3	0	0	0.16	0.017	0.04	0.26	1	0.1	0.5
ΔQ3	0.001	0	0	0.01	0.017	0.77	0.18	0.4	0.07
Δiod3	0.016	0.11	0	0	0	0	0	0	0
Δioq3	0.016	0.11	0	0	0	0	0	0	0
ΔiolineD1	0.1	0.07	0	0	0	0	0	0	0
ΔiolineQ1	0.1	0.2	0	0	0	0	0	0	0
ΔiolineD2	0.03	0.2	0	0	0	0.01	0	0	0
ΔiolineQ2	0.03	0	0	0	0	0.01	0	0	0
ΔioloadD1	0	0	0	0	0	0	0	0	0
ΔioloadQ1	0	0	0	0	0	0	0	0	0
ΔioloadD2	0	0	0	0	0	0	0	0	0
ΔioloadQ2	0	0	0	0	0	0	0	0	0
Δvdc1	0	0	0	0	0	0	0	0	0
Δvdc2	0	0	0	0	0	0	0	0	0
Δvdc3	0	0	0	0	0	0	0	0	0

Symbols	Descriptions	Values
Lo,iRo,i	Output impedance of inverter	1mH,1×10−1Ω
Lline,1,Rline,1	Impedance of line1	1 mH,2×10−3Ω
Lline,2,Rline,2	Impedance of line2	2 mH,3×10−3Ω
mp	Droop gain frequency	10-3 rad/s
nq	Droop gain voltage	10−3
Vo	Set point voltage	110 Vrms
fo	Set point frequency	50 Hz
ωc	Cut off frequency of filter	30 rad/s
VDC−Link	Voltage of DC link	200 V
Vtrip	Trip voltage of DC link	280 V
Vtr	Voltage of Triggering	215 V
Kdc	Linearization factor	2.4
Cdc_link	Capacitor of Dc-Link	2 mF
Pdc	Controller gain	5 ×10−4

Author Contributions

B.N.A.: writing-original draft, software, methodology, and validation; B.H.J.: formal analysis, investigation, resources, writing-review and editing; A.A.-M.: supervision, funding, writing-review and editing, F.B.: supervision, writing-review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.