Grid-Connected Power Converters: An Overview of Control Strategies for Renewable Energy

: The move towards a greener energy mix to ﬁght climate change propels investments in converter-interfaced resources such as wind and photovoltaics, energy storage systems and electric vehicles. The ongoing evolution of the power system is occurring at a very fast pace, challenging transmission and distribution system operators to seek solutions that are not only adequate for this moment but also for future scenarios. Ongoing research in the ﬁelds of power electronics, power systems and control aims at developing control strategies that will help the energy transition to occur, while keeping a stable, secure and reliable power system. The objective of this paper is to present a critical review of the control strategies developed for grid-connected power converters found in renewable energy systems, energy storage systems and electric vehicles. The impact of grid-connected converters on the stability of power grids is also reviewed, highlighting the promising control strategies for enhancing system stability.


Introduction
Economic development and population growth are providing great opportunities for the early 21st century's technological development.The challenge of keeping economies growing and sustaining the needs of a growing number of inhabitants, whilst envisioning the near future, where the climate crisis will play an even larger role in everyday life, calls for research in topics that enable a sustainable future for generations to come.Today, much of the electricity is generated from fossil-fuel-based generation sources [1].Moreover, the transportation sector relies mainly on Internal Combustion Engines (ICE), which also contributes to the historically high levels of CO 2 emissions recently reported.
From the power generation sector's perspective, an environmentally friendly alternative power source is renewable energy, which is a clean energy source that does not contribute to any adverse effects on the environment [2].Another factor that led to the increase in the search for renewable energy was the decrease in the cost of implementing the energy generation system, such as solar photovoltaics (PV) or wind energy conversion systems [3].However, such forms of energy generation have the disadvantage of depending on environmental conditions, thus creating the need for a specific generation project for each region and its particularities [4].Such sources provide a fluctuating energy supply, which becomes a problem to supply critical demands; these problems have led to numerous scientific research studies in the area of renewable energy [5], especially regarding the joint operation with energy storage systems (ESSs) to alleviate the intermittency challenges.
In the transportation sector, the push for environmentally friendly solutions is to substitute the ICEs for electrical motors.The power source for the electric vehicles (EVs) is a battery ESS in place of fossil fuels.This is not, however, a simple transition.Without developments in the power sector in the direction of renewable energy, EVs are recharged with power from fossil-fuel-based generators, which would only displace the source of CO 2 emissions.Moreover, the power sector would have to provide the charging infrastructure for an increased number of EVs.
Renewable energy sources, ESSs and EVs can be classified as converter-interfaced resources (CIRs).They are connected to the main grid through power converters, responsible for the conversion of direct current (DC) into alternating current (AC), or vice versa.Another important role of power converters is the management of DC voltage and current levels that can be done by DC/DC converters.In PV systems, a voltage source converter (VSC) is employed as a DC/AC converter to inject the DC current generated at the PV panels into the AC electric grid.Moreover, a DC/DC buck or boost converter might be encountered in PV farms for coupling the voltage levels of the PV panel with the ones in the DC/AC converter.In wind power applications, a Back-to-Back (B2B) converter is found in Doubly Fed Induction Generators (DFIGs) and Permanent Magnet Synchronous Generators (PMSGs).The B2B converter is composed of two VSCs with a common DC link whose function is to control the machine and to manage the power flow into the grid by controlling the DC voltage.In ESS applications, the most common converters are DC/DC converters that allow the DC voltage level to be changed for charging or discharging of the ESS and a DC/AC converter when the charging station is connected into the main AC grid.
The increasing penetration levels of CIRs are translating into significant operational challenges for system operators across the globe as they exhibit a very different interaction with the electric grid than traditional synchronous generators [6,7].Hence, with growing shares of renewable energy and ESSs in the power system, more sophisticated control strategies are needed to maintain the reliability and stability of the power supply [8].The traditional control strategy of CIRs is based on well-established Proportional Integral (PI) controller schemes following the grid's phase angle for current injection.Recent advances in power electronics control, however, have explored strategies based on Model Predictive Control (MPC).
MPC is a family of control strategies that use a mathematical model of the system to predict future states and obtain the optimal control input to achieve the desired performance within a control horizon [9].The cost function to calculate the optimal input based on the converter's performance is usually tied with reference tracking and control efforts.Within the field of power electronics, Finite Control Set (FCS) MPC was proposed since the set of possible switching signals sent to the converter is finite [10,11], which simplifies the optimization problem.Moreover, MPC can also be combined with Repetitive Control to increase the robustness in non-ideal voltage conditions [12][13][14].
In recent years, the interest in control strategies capable of increasing the contribution of CIRs in power system stability has increased.For instance, modifications to traditional control structures to enable the participation of wind turbines have been studied in [15,16].Another rising interest of the industry is the Grid-Forming (GFM) converters.A GFM converter is an operation mode of power converters that enable renewable energy and energy storage resources to be controlled as a voltage source rather than as a Grid-Following (GFL) current source [17].Then, a GFM converter is able to generate its own phase angle, frequency and amplitude [18].As a consequence, GFM converters are able to improve the frequency stability of power systems [19][20][21].GFM converter technology is still at its early development stages and an increase in the number of works investigating it is expected in the coming years.
Considering that the development of control strategies for CIRs is an important task to be carried out to comply with the operational challenges of the near-future power system, the objective of this work is to present a review of the modeling and control of power converters in renewable energy applications.Within this context, we have included ESS and EVs, since, in order to fulfill their purpose of reducing carbon emissions, they need to be developed alongside renewable energy sources.In each section, we review the traditional PI-based control strategy and one MPC alternative.This paper is structured as follows: Section 2 presents the modeling and control of single-and two-stage PV farms.The single-stage PV farm is presented to illustrate the classic PI controller of a grid-connected VSC; the two-stage PV farm is studied with a focus on the DC/DC converter that regulates the voltage level for performing the MPPT algorithm.In Section 3, the wind turbine's DFIG is studied to illustrate the classic and predictive control of B2B converters.While the VSC tied to the grid controls the DC voltage, the one connected to the rotor side is responsible for the MPPT algorithm.In Section 4, we focus on the control of grid-tied ESSs.The storage considered can represent a grid-tied ESS or an EV.The models and MPC controllers presented complement the ones shown for the VSC of single-stage PV farms, since they present a similar model.Finally, Section 5 discusses the influence on the power system stability of the control mode chosen for the power converter.In this section, the control scheme for the GFM is presented and the motivations to research GFM control strategies are stated.Finally, Section 6 presents the conclusions.

Control of Solar Photovoltaic Farms
In this section, the control of solar PV farms is discussed.At first, a discussion on the general characteristics of PV farms, such as the development status of PV installed capacity and the main topologies for interconnecting PV farms to the grid, is presented; then, a review of the control strategy for DC/AC converters for PV farms is presented; finally, we review control strategies for the DC/DC converters applied in two-stage PV farms.

Characteristics and Connection Topologies of PV Farms
The global installed capacity of PV farms has soared in the last two decades.From nearly 1 GW in 2001, solar PV capacity has increased exponentially to 707.5 GW in 2020 [22].This sudden growth was made possible through incentive programs, declining costs of the mass production of PV panels and the investments in carbon-free power sources [23][24][25][26][27]. Integrating growing levels of PV farms to the electric grid, however, is not a straightforward task given the intermittent nature of the solar irradiance [28].
The solar irradiance is the primary resource of PV farms and its availability is subjected to cloud covering, the time of year and the time of the day [29].Moreover, due to the PV panel's construction, its current output, which translates into the generated power, is affected by the resource availability and the voltage applied to its terminals, as seen in Figure 1.In the early years of the deployment of PV energy, the main objective was to extract the maximum power available to maximize the financial returns.Hence, the main strategies for controlling a PV farm are derived from MPPT schemes [30][31][32].The MPPT for PV farms aims at producing a DC voltage reference value to be applied to the PV module that would result in maximum power extraction for the current solar irradiance condition, as given by the pink curve in Figure 1.
With the negligible share of PV power connected to the grid in the beginning of the development, the operation of PV farms concerning exclusively the injection of active power was possible.However, nowadays, the share of PV farms in some power systems cannot be neglected and changes are required for the control of PV farms [33].Therefore, grid codes have been updated and a greater contribution to power system stability is required from PV farms.As an example, the PV-STATCOM [34] proposed the use of the inverter to supply reactive power support to the power grid at nighttime when the converter was idle.Following work extended its application to daytime [35], where the remaining capacity of the PV farm could take part in the reactive power support of the grid if necessary.Despite the increase in power losses of this type of PV farm operation [36][37][38], this modern control philosophy allows PV farms to contribute actively to system stability, as seen in [39][40][41].
There are two possible topology connections for a PV farm with implications for its operational flexibility: single-stage and two-stage connections.In a single-stage PV farm, as shown in Figure 2, a DC/AC converter is employed to convert the DC power generated by the PV module into AC power to be injected into the grid [42][43][44].In a two-stage PV farm, as shown in Figure 3, the PV module interface with the electric grid is composed of a DC/DC converter and a DC/AC converter [43,45].Due to the extra power converter in a two-stage PV farm, single-stage PV farms are usually more efficient.However, in two-stage PV farms, the PV panels' output voltage associated with the MPPT algorithm is fully decoupled from the AC grid's peak voltage, which can increase operational flexibility [37].In Figures 2 and 3, v pv is the voltage applied to the PV module given by the MPPT algorithm; v dc is the DC voltage at the DC/AC converter's DC side in the two-stage connection; C and C pv are the capacitances connected to the DC/AC converter and to the PV module; i abc is the AC current injected into the grid; R f and L f are, respectively, the resistance and the inductance of the smoothing reactor.
The capacitor banks allocated in the DC circuit are responsible for filtering the DC current ripples and act as energy storage for the DC circuit.The value of the DC circuit capacitor banks is an important design parameter since it affects the dynamics of the DC voltage [46].The capacitor values can be chosen following the guideline (1) proposed by [47].
where v dc,min is the minimum converter DC voltage, S converter is the converter apparent power, f is the AC side frequency and ∆v dc is the maximum desired ripple in the DC circuit (typically considered as 2%).For the smoothing reactor L f , a typical value of 0.1 p.u. can be considered, whilst for R f , a typical parasitic resistance value of 0.01 p.u. can be considered.
Next, the control strategies of the converters associated with PV farms are discussed, starting with the control of DC/AC converters.

Control of DC/AC Converters for PV Applications
To develop a controller for the DC/AC converters shown in Figures 2 and 3, the starting point is a mathematical model for the main variables of the system.In this section, we present a model of the converter's interconnection to the electric grid and of the dynamics of its DC bus voltage.Further, we discuss control structures for PV farms.
A mathematical model for the DC/AC converter's currents can be obtained by applying the Kirchhoff voltage law to the AC side circuit: where L f and R f are the filter inductance and parasitic resistance, v c,abc is the voltage synthesized at the converter's terminals, v g,abc is the grid voltage at the point of connection and i abc is the converter current.By using the model given by (2) as a starting point for a controller, the designer is faced with a challenge.In the abc frame, the current references are time-varying as the converter is connected to an AC grid with sinusoidal currents.Therefore, one solution is to apply the Park Transform given by: Hence, the following model is obtained in the dq frame (assuming a balanced AC system): where ω is the AC grid angular frequency.Since the dq frame is rotating at a speed ω, when using a model developed in the dq frame, the controller is faced with the task of following constant reference values.This is a result of a constant projection of the currents in the rotating dq axis.Therefore, by employing a PI controller, it is possible to ensure zero steady state tracking error.
However, a cross-coupling between i d and i q poses a disturbance to the system.This disturbance can be dealt with considering the auxiliary inputs ∆v d and ∆v q that include a feedforward compensation as [48]: Then, substituting (6) and ( 7) into (4) and ( 5), the following model is obtained: The models given by ( 8) and ( 9) can be written as first-order transfer functions G(s) given by: The current dynamics can be controlled using a PI controller G i (s) as: with k p = L f /T s and k i = R f /T s [49], where T s is the desired response time for the current loop.Now, the task at hand is to obtain the reference i * d and i * q that the current loops should track.
Given the cascaded control structure briefly discussed at the beginning of this section, the i * d reference value is obtained from the DC voltage control.Regarding the single-stage PV farm, the DC voltage reference is generated by the MPPT and it is applied to the PV panel to ensure maximum power operation.As for the two-stage PV farm, the DC voltage is a constant value that can be chosen by the PV farm operator according to its operation strategy.
The mathematical model of the DC voltage dynamics is obtained by observing the power balance in the DC side capacitance.The following model can be obtained for the DC voltage for both topologies [48,50]: where i v represents the current injected by the PV panel in the single-stage PV farm or the current injected by the DC/DC converter in the two-stage PV farm.When considering i v as a disturbance, a transfer function model can be written for (12), relating v dc to the reference i * d as: A PI controller G dc (s) with the same structure as (11) can be used as an outer loop to regulate the DC voltage by changing the value of i * d .The tuning of the outer loop PI controller can be done by following the methods proposed in [48,49].
Finally, considering that the system is balanced, the reactive power reference can be translated into a reference value i * q as: The control block diagram of the strategy using this cascaded PI controller structure is represented in Figure 4. From the contributions of ∆v d and ∆v q , a voltage reference v * abc is obtained and a Pulse Width Modulation (PWM) is responsible for the generation of the gating signals sent to the converter's switches.
The simulation of a 250 kW single-stage PV farm is shown in Figure 5.The connection of the PV farm to the electric grid follows the one presented in Figure 2, with parameters given by Table 1.The considered solar irradiance profile is shown in Figure 5a.In Figure 5b, the DC voltage tracking is shown, where the reference value v * pv is given by the MPPT, following the pink-crossed line in Figure 1.The cascaded PI controller scheme is able to accurately track the value for v pv with the irradiance variations shown.Finally, Figure 5c shows the power injected by the PV farm into the electric grid.The average power follows the solar irradiance profile of Figure 5a, reinforcing that the DC voltage controller can quickly follow the MPPT reference.

Control of DC/DC Converters for PV Applications
DC/DC converters are widely used in PV farms.Their main application is to step up the voltage of the PV array to the voltage level of a DC/AC converter in a two-stage PV farm [51][52][53] or to interface the PV panel with a DC microgrid [54,55].In both applications, the DC/DC converter is responsible for controlling the DC voltage applied to the PV panel to operate in the reference value provided by an MPPT algorithm.In this work, we focus on the control of the DC/DC boost converter for two-stage PV farms.
A DC/DC boost converter for PV systems is shown in Figure 6.The DC/AC converter representing the second stage is assumed to be an ideal voltage source.The modeling of boost DC/DC converters is not a straightforward task [56,57].The system changes states as the switch turns on and off [58,59].An average model of the system given by ( 15) and ( 16) can be obtained by using averaging theory [60].
where v irr is the voltage of the PV panel, v pv is the voltage at the PV panel's terminals, V dc is the voltage of the second stage of the PV farm, R is the internal resistance of the PV panel, L is the inductance of the DC/DC converter and d is the converter duty cycle.
The structure of the model in ( 15), (16) suggests that the control of the boost converter can be achieved in a cascaded structure, where the outer loop is responsible for generating the reference for the current i L to steer the PV panel voltage v pv to its reference value given by the MPPT.Then, the inner current control loop is responsible for tracking the reference of i L by manipulating the duty cycle.The duty cycle is the input to a PWM generator responsible for determining the switching signals S sw for the switch sw.This control strategy is represented in the block diagram of Figure 7.The outer and inner control loops can be designed with PI controllers using the structure given in (11).A procedure for tuning the PI controllers by using the Ziegler-Nichols method is presented in [61].More sophisticated techniques for tuning the controllers of boost converters applied to PV systems based on optimization algorithms are presented in [62][63][64].
To avoid using an average model of the converter, more advanced control strategies were proposed for the DC/DC boost converters used in PV systems.The two most popular control strategies are based on sliding mode control [58,[65][66][67] and FCS [68][69][70].In sliding mode control approaches, the boost converter modeling takes into consideration both states of the circuit in Figure 6, depending on the position of the switch, sliding between them as the switch changes its position.
On the other hand, FCS-based control strategies, as a member of the MPC family, take advantage of a model of the system and the limited number of possible states to solve an optimization problem to choose if the switch must be on or off.Moreover, the design of FCS controllers follows comprehensible steps and their design is fairly simple.As an example, the FCS control strategy proposed by [70] is shown in Figure 8.In Figure 8, ( 18) and ( 17) are given by: where the argument k explicitly indicates the discrete time nature of the control strategy and T sw is the switching frequency.As seen in Figure 8, the MPPT provides v pv (k + 1), the desired voltage to be applied to the terminals of the PV model, at the next time step to ensure that the maximum power is extracted.In this sense, v pv (k + 1) can be interpreted as a reference value for the voltage v pv .Then, ( 17) is used to predict the value of the inductor current in the next time step if the PV panel is operating at the v pv (k + 1) voltage given by the MPPT.Finally, the duty cycle d can be obtained from (18), driving the DC/DC converter towards the desired value for the inductor current i L (k + 1).The desired duty cycle is then used as an input to the PWM generator to obtain the switching signals S sw for the switch sw.At each sample, as new measurements become available, the process is repeated.
To illustrate the behavior of the FCS controller shown in Figure 8, the experimental setup used in [70] was subjected to step changes in the inductor current and in the voltage applied to the PV panel side.In Figure 9, the behavior of the inductor current of the DC/DC converter is shown when a step change is applied to the reference value i * L .The controller is able to track the reference change within 2 ms.When the step change is applied to the reference v * pv , the FCS controller is able to regulate the voltage applied to the PV panel v pv within 10 ms, as seen in Figure 10.The voltage controller is five times slower than the current controller, explicitly showing the time-scale separation between the two cascaded loops.Step response of the voltage v PV using the predictive controller.

Control of Doubly Fed Induction Generator-Based Wind Turbines
This is a dynamic model; therefore, it serves to analyze and study the transient behavior of the Doubly Fed Induction Generator (DFIG) and, more importantly, to lay the groundwork for the subsequent development of machine control.Through this dynamic model, it is also possible to know the behavior of the machine in steady state and, in general, it is possible to determine the continuous performance of the machine and the behavior of its variables, such as torque, currents and magnetic fluxes, under certain supply voltage conditions.
Likewise, the use of this information makes it possible to know how the transition from one state to another will take place, being able to detect unsafe behaviors, such as instabilities and high transient currents, in addition to providing information on dynamic oscillations, torque or current ripples, among others.

DFIG's Mathematical Model in the Synchronous Frame (dq)
According to the AC machine models developed and discussed by different authors, such as [71], a model for the DFIG can be obtained by assuming that the three-phase system is symmetrical, the permeability of iron is infinite, the machine is magnetically linear, the resistances and reactances are constant, the voltages and currents are sinusoidal in steady state and, finally, that the effect of the harmonics in the iron is not considered.The model is obtained by writing Kirchhoff laws for the equivalent circuits in the abc frame and can subsequently obtained in the dq (or synchronous) reference frame by applying a Park transformation.
where the subscripts s, r indicate variables corresponding to the stator and the rotor, respectively; the resistance and inductance of windings is given by R and L; the mutual inductance between the rotor and the stator is given by L m ; the magnetic flux is represented by ψ; the electromagnetic torque is T m ; p is the number of pairs of poles; and ω r , ω s and ω m are, respectively, the rotor angular frequency, the synchronous speed and the rotor mechanical speed.

Stator Flux-Oriented System
The Stator Flux-Oriented Vector System is a control strategy that uses the DFIG's mathematical model, shown in the previous section, that aligns the "dq" reference frame with the stator flux.Hence, it is possible to consider that | ψ dq,s | = ψ d,s and consequently ψ q,s = 0, allowing us to rewrite the expressions ( 19) and ( 21) as ( 24) and ( 25), respectively.
According to (25) and separating this equation for the two axes "d" and "q", the stator currents can be written as ( 26) and (27) for each axis.
When the machine is in steady state, the voltage drop across the stator resistance can be neglected and the stator voltage can be considered perpendicular and proportional to the stator flux, resulting in: Using this information, the expressions of the active power P and the reactive power Q of the stator in the stator flux-oriented synchronous reference frame are given by Equations ( 29) and (30) [71]: This strategy offers a decoupling between the "dq" axes for the stator currents and powers, allowing the possibility of individually controlling the currents or powers.This strategy is very useful, as the DFIG's stator is directly connected to the electrical system and the standard control strategy is to regulate the stator variables through the rotor currents i dq,r .The stator flux angle can be computed using a PLL [72,73].

State Space Model
The intention is to take advantage of the state space model to apply control over the rotor currents using the rotor voltages as the manipulated variables.Thus, by substituting ( 26) and ( 27) into (22) and further inserting the result into (20), the expression (31) is obtained.
where ω r is the slip speed and the angular frequency of the rotor, given by ω r = ω s − ω m and σ = 1 − (L 2 m /L 1 L 2 ).Now, from (31), it is possible to deduce the equations of the state space model, in which the current components i d,r and i q,r are the state variables (and also the output variables by considering the output matrix C c as an identity matrix), the voltages v d,r and v q,r are the input variables or control signals, and ψ s is considered as a disturbance.
In this way, the model represented in (32) is obtained. where Discretizing the model (32) and considering T as the sampling period and k as a sampling instant, using a Zero-Order-Hold (ZOH) [74,75] without delay, the discrete model from (32) can be written as (34). where

Proportional Integral Control for DFIG
For the direct power control of the DFIG, a cascaded control scheme is used to control the power produced by the DFIG, being injected or absorbed by the grid.Following the cascade scheme through Equations ( 29) and ( 30) from the active and reactive power references provides a reference value for the currents i d , r and i q , r .
Finally, calculating the error between the reference currents and the measurement, a PI processes the error for both the direct axis and quadrature axis current to provide the voltage values that must be applied to the rotor through an inverter controlled by PWM.
V r,q = ( i r,qRe f The diagram presented in Figure 11 represents the full scheme for the DFIG's direct power control using PIs.

Model Predictive Control for the DFIG
An MPC controller uses a model of the system to predict its future actions and calculate the optimal inputs based on the minimization of a cost function.Additionally, it allows for the easy handling of the system's constraints.
From the discrete state space model obtained in the previous section, the prediction of the components of the rotor current vector can be performed as [76]: where In the previous equations, n y is the prediction horizon of the grid current, and I g is the prediction of the rotor current vector.
The quadratic cost function J in matrix representation is given by: where I rRe f ∈ R (n y •q)×1 represents the rotor current vector future references, M y ∈ R (n y •q)×(n y •q) weighs the controlled grid currents (outputs) and their predictions, M u ∈ R n u ×n u weighs the control efforts of inverter voltage minus rotor voltage vectors, U ∈ R n u ×n u is the input, and n u is the control horizon.
The expression of the minimized cost function (Equation ( 45)) can be represented algebraically from ∂J ∂U = 0 and isolating U. Hence, the representation is given by: The rotor current vector control using MPC calculates the voltage vector of the inverter in the synchronous reference frame using Equation (46), considering the value of the rotor voltage vector components in the fundamental component frequency.This voltage vector in the rotor frame is the input of the PWM algorithm as space vector modulation.In this way, the rotor current will reach their references.
The MPC controller needs the information about the angle and the fundamental frequency component of the rotor flux to compute the appropriate voltage vector of the inverter to operate during rotor voltages.
A block diagram of the MPC control of the DFIG is shown in Figure 12.
V r,q Figure 13 presents an experimental result of the MPC control applied to DFIG.Table 2 presents all the DFIG parameters present in the experimental bench.We show the voltage and current signal in the stationary phase, together with the current values in the synchronous reference, where it can be observed that when the current reference i q,r is changed from 1 to 3 A and the i d,r reference is kept; there is an increase in the amplitude value of the current in the stator of the machine and thus more power is injected into the grid consequently.In addition to the MPC presented in this section in a system with normal operation, in the literature, the same type of controller has been applied in a system with disturbances.In [77][78][79][80], the authors present predictive control in a system where there is an imbalance in the voltages of the electrical grid.Another situation found in practical applications is the distortion present in the main voltage due to harmonics generated by non-linear loads, as presented in articles [14,78,81].

Control of Grid-Connected Energy Storage Devices and Electric Vehicles
Nowadays, electric vehicles (EV) are a sustainable transport option that generates less noise and pollution when compared to vehicles with internal combustion engines.EVs use electric motors and batteries to propel the vehicle, operating more efficiently and at lower operating costs.The advancements in lithium-ion batteries along with fast charging technology are factors that have helped in the launch of EVs [82,83].
When an EV is connected to the electric grid, it can behave as a load (when it is charging) or it can provide a greater contribution to the electric grid by behaving as a grid-connected battery at the expense of additional control loops [84].In this latter case, the EV is considered a grid-connected ESS.When controlling the DC/AC converter of an EV, the objective of the control system is to directly regulate the apparent power elements (S = P + jQ) exchanged with the electrical grid.The active and reactive power references are generated by a higher control level called an energy management system (EMS).In this paper, we have considered a simple EMS in the form of a PI controller that generates the active power reference P g,Re f based on the the desired DC voltage V Re f Bat to be applied in the battery and a fixed reference Q g,Re f that can be chosen by the distribution system operator.However, more complex schemes based on MPC to obtain such references can be adopted (see [85,86]) that can help to manage the lifetime and the state of charge of the battery.The signal of the active power reference will determine if the EV behaves as load or a battery.Moreover, the possibility to inject/absorb reactive power allows for the participation in the support of the grid.
For the required additional control strategy, it is necessary to measure the voltages and currents on the grid side.To facilitate the use of classical vector control approaches, the grid currents must be in a specified referential state, (α, β) or (d, q) Hence, the phase currents i g,a , i g,b and i g,c are transformed into vectors in the desired frame of reference.On the other hand, the grid voltage is measured and its angle detected through the Phase-Locked Loop (PLL) [87], in order to obtain the voltages in the referential state v g,d and v g,q when using the synchronous reference.This angle is used to transform the variables in the frame (a, b, c) to (d, q) and vice versa.
When aligning the phase a with the synchronous reference frame, the direct axis current i d directly controls the active power exchanged with the grid P g ; in the same way, the quadrature axis current ig, q directly controls the reactive power Q g exchanged with the grid.The advantage is that the control of these two components i g,d and i g,q is independent and provides an effective means for controlling active and reactive power on the grid side [88].
Therefore, by decomposition into components, the equations of the model of the converter and filter connected to the electrical grid are obtained by: where v and i are the voltage and current spatial vectors, respectively.R g and L g are the resistance and inductance of the filter, respectively; the subscripts g and inv denote the electrical grid and the converter, and ω g is the angular speed of the electrical grid.
From the voltage and current vectors, one can calculate the apparent power expression: From ( 49), the active and reactive power of the grid can be calculated.
Obtaining the grid model and calculating the values of active and reactive power makes possible the design of control techniques such as the Proportional Integral controller and the Predictive Control.The purpose of the control is to maintain the DC bus voltage from a reference provided by a battery management system (BMS) at a higher level of supervision to determine the charging and discharging of battery EVs.

Proportional Integral Control for Storage
Figure 14 demonstrates in a block diagram the entire closed loop using proportional integral (PI) controllers [89].From the mains voltage, the system changes to the frame (d, q), thus making the active and reactive power directly proportional to the currents i g,d and i g,q , since the system voltage oriented v d = 0 [49].To carry out the Voltage-Oriented Control (VOC) strategy, the direct axis of the system in the synchronous reference frame (d, q) is aligned with the grid voltage vector | v g,dq |, i.e., the direct axis voltage v g,d has the same magnitude as the grid voltage vector | v g,dq | = v g,d , so the component on the quadrature axis v g,q is equal to zero [88].
According to Equation ( 52), and with expressions ( 50) and ( 51), it is possible to calculate the active and reactive power that flow through the converter connected to the electrical grid and depend directly on the current vector i g,dq , as: From ( 53) and ( 54), it can be concluded that the active and reactive power can be controlled by controlling the direct axis and quadrature currents, respectively.
The direct axis current i g,d represents the active power P g on the grid side; this control loop regulates the power injected into the grid, when the converter works in steady state.The active power in the circuit is kept constant at the predetermined value P g,Re f as a reference; this reference is provided from a PI that calculates the error between V Re f Bat (sent from BMS) and V dc presented in (55).
Thus, the PI controller processes the active power error P g and obtains the reference current i dRe f through the following expression: The quadrature axis current i g,q represents the reference reactive power Q gRe f .This reference value can be adjusted in the following ways: positive to obtain an inductive power factor, negative to obtain a capacitive power factor and zero to obtain a unity power factor.It can be calculated through the following expression: With algebraic manipulation from ( 47) and ( 48), it is possible to obtain the following equation [49]: where ω g is the speed of the synchronous reference (grid angular frequency).Equations ( 58) and (59) show that the derivatives of the direct axis and quadrature axis component of the current depend on variables both on the quadrature axis and on the direct axis, respectively.This indicates that the control system is cross-coupled, which can lead to difficulties in the dynamic performance of the controller.
The coupling terms e g,d and e g,q are the voltages induced due to the three-phase transformation of the reference inductance L g to the synchronous reference (d, q), and are used to improve the behavior of the system; they can be calculated as follows: e g,d = −ω g L g i g,q (60) To improve the dynamic performance of the PI controllers, the voltages U g,d and U g,q of the controller decoupled from the mains can be defined as: This solution shows the independent behavior of the two control loops.Decoupled control makes the designing of PI controllers more convenient and the system is stabilized easily.
Finally, in order to be able to calculate the coordinate transformations for the voltages and currents of the electrical grid, the grid voltage angle θ g is needed.With this angle, the converter can be connected to the electrical grid, and the frequency of the converter output signal must be synchronized with the frequency of the electrical grid.A technique to synchronize with the electrical grid is the phase-locked loop (PLL), which needs the three phase voltages (a, b, c) of the grid, and the output of the PLL is the phase angle θ g [49].
The main purpose of the PLL operation is to control the phase angle through the processing of a PI regulator; thus, the three-phase sinusoidal voltages have a transformation in the referential state (d, q), through the same angle initially estimated.This quadrature axis voltage is processed by a PI controller in which the angular velocity of the electrical grid is obtained, until this signal approaches zero, thus ensuring the alignment of the direct axis d with the electrical grid vector v g,dq , which allows us to infer that there is a synchronism between the grid and the converter.
The angular speed of the electrical grid ω g can be calculated as follows: where k vqP and k vqi are the gains of the PLL controller PI.The mains angle is calculated using: There are PLLs that provide greater robustness to the system when it comes to harmonics distortions, as presented in [90][91][92].

Finite Control Set Model Predictive Control (FCS-MPC) for Energy Storage Systems
FCS is a controller that falls under the umbrella of MPC.Hence, it is a control strategy that relies on a model of the system to predict its future behavior, and through minimization of a cost function, it obtains the optimal input to be applied to the system.In the case of converters based on power electronics, as is the case of the DC/AC converters found in grid-connected ESSs, FCS relies on the finite set of possible control inputs (the state of the converter switches) to evaluate every possibility to find the input that minimizes the cost function [93].Figure 15 shows a block diagram of the inverter connected to the grid with predictive control of the type finite control set, which summarizes this control strategy for grid-connected ESSs.

PI Predictive
Model eq 68 Calculation of P (k+1) and Q (k+1) eq 69 and 70 Minimization of the cost function g eq 71 An FCS-MPC controller takes advantage of the fact that the switches of the two-level DC/AC converter shown in Figure 15 can only assume a finite number of states that generate the output voltages shown in Table 3.The notation with the complex number j implies in this scenario that the system is modeled in the stationary reference frame.

S a S b S c
Vector Voltage v inv (j) The converter currents and voltages are obtained in the stationary frame following [11]: where x denotes the grid current i g , the grid voltage v g or the inverter output voltage v inv .
In discrete time, the difference equation model of the converter can be used to predict the current of the next step as [93]: From ( 68), the active and reactive power exchanged by the converter with the electric grid can be predicted by: From the predicted active and reactive power values ( 69) and ( 70), the following cost function g is formulated: where P g,Re f and Q g,Re f are the active and reactive power references, respectively.Finally, since, for a short sampling time, it can be considered that v g (k + 1) ≈ v g (k), the cost function ( 71) can be minimized by investigating the effect that the finite set of control inputs shown in Table 3 has on the cost function g by evaluating the predicted currents using (68).The vector v inv,αβ that results in the lowest value of g is chosen and the corresponding switching is applied to the converter.Figure 16 summarizes the how the FCS-MPC algorithm works in a flowchart.
Figure 17 shows the transient response of an active power step of the implemented controller in an experimental bench using the DSP TMS320F28335, with electronic boards to condition the voltage (v) and current (i) through Hall effect magnetic sensors and a back-to-back inverter of the model "Semikron SKS20F (B6CI)2P + E1CIF + B6U 14V12" connected to the grid.To connect the inverter to the grid, there is an electrical contactor to integrate the experimental bench into the grid.
The active power step varies from 500 W to 750 W and has a settling time of around 1.8ms.The next test presents a step in the V Re f Bat voltage and the variation of the active power provided by the PI controller of the voltage loop.
The results in Figures 17 and 18 support the idea of the controller being able to follow a requested power reference when there is a variation in the bus voltage.The parameters of the experimental bench and electrical grid components are given in Table 4.With the voltage signal coming from an external battery voltage management loop, the control proposed in the literature is able to maintain the voltage level with a fast dynamic response, since the voltage reference is reached with a time of 0.1 s.Moreover, once the bus voltage is reached, the power goes to zero until a new voltage value is sent by the battery management system.

Stability of Power Systems with High Penetration of Grid-Connected Converters
The dynamic characteristics of power grids have significantly altered with the increasing penetration of power-electronic converter (PEC)-interfaced generation sources, transmission and distribution systems and loads [94].More specifically, these technologies include wind and photovoltaic generation sources, energy storage systems (e.g., battery energy storage systems, flexible ac transmission systems (FACTS), compensation devices (static-synchronous compensator (STATCOM) and static-var compensators (SVC)), high-voltage direct current (HVDC) converter stations and power electronic-interfaced loads [18,95].Therefore, the PEC-interfaced stability issues are not only limited to generation sources, but could also originate from other PEC-interfaced devices in the power grid, such as HVDC converters and STATCOMs [94,96].
The PEC topology used in these PEC-interfaced sources plays a significant role in power grid stability.The PECs could be broadly categorized into two major types based on their topology [97,98]: (1) Grid-Forming (GFM) converters, and (2) Grid-Following (GFL) converters.The GFM and GFL converter topologies are shown in Figure 19.
Generation sources interfaced with GFM converters act as a voltage source (more precisely, as a controlled voltage source) and they can maintain the internal voltage vector constant by maintaining a constant voltage magnitude and frequency under disturbances [97].Thus, they have a number of superior features, such as the black start capability, enhanced synchronization performance in weak grids, inertia and frequency support capability (to improve the rate-of-change-of-frequency (RoCoF) and frequency nadir [18,99,100].However, generation sources interfaced with the GFL converters [97] are designed to follow the phase angle/frequency of the voltage waveform imposed by the synchronous generators in the power system, and inject the power into the grid according to a fixed power reference or a power reference set by the maximum power-point tracking (MPPT) system.Therefore, typically, the GFL converter-interfaced sources are not typically designed to support system inertia and reactive power requirements, and also have a small fault current contribution [101].These features have created multifold operational and stability challenges in power grids, and the impact depends on the degree of the penetration level of GFL converters.A small penetration of GFL converters will not influence the operation of the power system, and hence their impacts can be ignored.However, when power systems begin to integrate a large number of GFL converters with high capacities, their impacts will not be restricted to the local region/area, but will begin to influence the entire power network.Subsequently, the overall dynamic performance of the power system is significantly affected.On the contrary, the GFL topology can be improved with auxiliary control schemes, such as frequency and voltage controllers, to provide enhanced support to the power grid [97].However, their fundamental topology (i.e., GFL) remains unchanged despite these additional support services provided to the power grid.
The stability issues associated with the PEC-interfaced sources are well reported in the published literature [96].For example, PEC-interfaced sources have influenced the transient stability [102], small-signal stability [103], voltage stability [104,105] and frequency stability [8,106].However, the present-day stability issues originating from the PEC-interfaced sources can be broadly discussed under the following three themes: 1.
Low short-circuit strength and weak grid issues; 2.
Inertia reduction and frequency stability; 3.
Converter-associated interactions and resonance.
The above stability issues are discussed in detail in the following sections.

Low Short-Circuit Strength and Weak Grid Issues
Short-circuit strength is an imperative factor for power grid stability which assists in maintaining the power grid voltage following disturbances in the power grid [101].The short-circuit strength or short-circuit capacity of a given busbar depends on the Thévenin equivalent impedance, Z th , at that location.Hence, the short-circuit capacity (SCC) at a given location is given by: where V t is the rated voltage at the busbar.Short-circuit strength is high closer to synchronous generators and lower at remote locations (due to high Thévenin equivalent impedance).The Short-Circuit Ratio (SCR) is defined as the short-circuit capacity in MVA at the point of connection (POC) divided by the rated apparent power output in MVA of the generation source/or PEC; The effect of the short-circuit ratio was investigated considering the test system shown in Figure 20.The test system consisted of a 1000 MVA synchronous generator connected through a 110 kV, 50 km transmission line.The rated power factor of the synchronous generator was 0.8; hence, the synchronous generator was rated at 800 MW.The synchronous generator was represented by the simple excitation system model and simple governor model.Typical parameters were considered for the synchronous generator, the excitation and the governor.The short-circuit ratio was calculated between the external grid busbar short-circuit MVA and the generator rated MVA; hence, the external grid busbar short-circuit MVA was varied to analyze the performance under different SCRs.
When the SCR falls between 3 and 2, the system is considered to be weak, and when the SCR falls below 2, the system is considered very weak.Under weak grid conditions, any small disturbance at the generation source terminal could lead to instability.Figure 21 illustrates the generation source terminal voltage under three different short-circuit ratios (SCRs).According to Figure 21, when the SCR was at 1.5, even a 50 ms 0.12 p.u. voltage drop at the generator terminal caused the terminal voltage to oscillate and become unstable.However, in the other two instances, the terminal voltage managed to recover following the voltage drop.Other than the voltage stability issues, additional issues, such as power transfer issues and converter synchronisation issues exist, under weak grid conditions.
GFL converters contribute to weak grids from two aspects: (1) connecting at remote locations in the power grid with high Thévenin equivalent impedance; (2) low short-circuit current contribution during faults.Typically, renewable energy source and required geographical characteristics are available at remote locations; hence, GFL-based generation sources inevitably need to be situated at these remote locations.Therefore, when the penetration level of the GFL converters is significantly increased, a power system may experience weak grid issues [107].Although the first factor is unavoidable, the second factor could be mitigated by deploying a GFM converter topology and synchronous condensers [18].Therefore, power system operators are trying to keep the minimum required shortcircuit strength to maintain power system security by keeping a minimum level of conventional synchronous generators online in the system.This is one of the limiting factors of the grid GFL converter-interfaced generator integration [108].The rating of the GFL converter is usually determined to keep the SCR at 3 at the point of connection (POC) when the minimum number of synchronous machines are online in the system under N1 critical contingency events [109].

Inertia Reduction and Frequency Stability Issues
Conventional synchronous machines have an inherent inertial response capability (response proportional to RoCoF and stored energy capacity), which is a very important attribute that helps power systems to arrest the frequency drop before reaching unstable frequencies that could ultimately lead to a blackout [106].Compared with synchronous generators, the inertia of the GFL converter is very low, due to the decoupling of the electrical dynamics of the grid and the electro-mechanical dynamics of the source by the PEC [8].The solar PV systems do not even have stored kinetic energy as they do not have any rotating parts [110].Thus, the future power system with a high share of PEC-interfaced generation sources no longer has the traditional inertial response characteristics based on rotational kinetic energy.
To substitute the primary frequency response capability of the synchronous generator, battery energy storage systems (BESSs) are employed in power grids [111].However, the BESSs still have the limitation that their capacity and lifetime incrementally deteriorate when they are subjected to frequent transient power fluctuations and deep charge-discharge operation cycles [111].Moreover, different inertial response control strategies have been developed [112,113]; however, their capability is limited due to the lack of response duration and dependency on the renewable energy source.
The system frequency response with wind generation was investigated using the test system shown in Figure 22.The two synchronous generators (i.e., SG-1 and SG-2) were rated at 500 MVA and their active power rating was 400 MW.The two synchronous generators were represented by the steam-turbine governor model and simple excitation system model.Typical parameters were considered for the synchronous generator, the simple excitation and the steam-turbine governor.Prior to the frequency event, a 50 MW load was not connected to the network and both SG-1 and SG-2 were operating at 250 MW active power output.Then, in a 50% wind generation scenario, SG-2 was disconnected and connected a full-converter wind generator (FCWG)-based wind farm.The FCWG-based wind farm had a rated capacity of 500 MW (100 × 5 MW) and operated at 250 MW power output.The frequency response from the FCWG was obtained based on the ∆f-based droop approach [18].
Figure 23 illustrates the system frequency response for a low-frequency event under three (3) scenarios: a) with synchronous generation, (b) 50% wind generation with no frequency response, (c) 50% wind generation with frequency response.The wind generation was represented by a full-converter wind generator (FCWG) and the frequency response was obtained based on the ∆f approach.According to Figure 23, when 50% of the system synchronous generation was replaced by the FCWG-based wind generation, the frequency nadir dropped to 49.69 Hz, which is a 0.11 Hz drop from the synchronous generator only case.However, when the ∆f-based frequency response strategy was adopted, the frequency nadir was improved by 0.04 Hz, and further improvement could be achieved with a high gain for the ∆f-based frequency response scheme.

Other Converter-Associated Stability Issues
Among the other stability challenges, low-and high-frequency interactions and resonance issues can be considered as emerging stability challenges in the PEC-interfaced sources.Among these converter-associated stability issues, certain stability issues occur during weak grid conditions [101].The following are the major converter-associated stability issues:

•
Weak grid-driven converter-associated stability issues-The phase-locked loop and inner current control loop can destabilize the PEC under weak grid conditions, while leading to sustained oscillations in the power grid [114,115].The PEC-interfaced generation sources and voltage source converter (VSC)-based HVDC converter systems have reported PLL-based stability issues under weak grid scenarios.High PLL gain is identified as the root cause of this instability issue [114].Various alternative PLL schemes have been proposed to overcome these weak grid-related stability issues [116].Moreover, inner current loop gains are also reported to be the cause of instability in PEC-interfaced sources under weak grid conditions [117].

•
Fast interactions-When the PEC-interfaced sources are involved with oscillation frequencies higher than the nominal power frequency, these interactions are categorized as fast interactions.These oscillation frequencies could be high as the kHz range.There are several mechanisms by which these fast interactions could occur: (1) between PEC controllers and the passive power network elements, (2) between one PEC controller and another PEC controller or a PEC cluster, (3) between PEC controllers and other network devices, such as synchronous generators.These fast interactions are not limited to the PEC-interfaced generation sources, but also present in other PEC-interfaced devices, such as STATCOMs and HVDC systems [118,119].• Slow interactions-The slow interactions can also originate from the same mechanisms outlined in the fast interactions category.However, the frequency range is less than the power frequency (typically less than 10 Hz); therefore, the frequency range of these interactions depends on the controller parameters and network elements.Other than these mechanisms, weak grid-related oscillations also occur at this low-frequency range, but should be analyzed under the weak grid category.

•
Resonance issues-The resonance issues could arise when the network elements (series capacitors) together with the source-end impedances form a resonance circuit at sub-synchronous frequencies, leading to effective negative resistance [120].Ultimately, this would lead to high current and voltage oscillations while damaging power grid equipment.This type of resonance condition has been reported in a number of power networks (e.g.,China and ERCOT-Texas) [120], and researchers have mainly referred to these resonance issues as sub-synchronous control interactions (SSCI).Compared with the conventional resonance issues, in this particular resonance issue, PEC converters play a major role by affecting the negative resistance by PEC controllers.Therefore, the resonance condition is governed by the PEC controller.

Conclusions
This article provided a review of challenges and control methods in renewable energy systems.Much of the challenge lies in the intermittent nature of the energy generation, which can be partially solved with control techniques or by integrating more than one generation method into the power system.In addition to the above, some devices, such as batteries or diesel engines, can be added to the system to improve the robustness under fluctuating energy production.The control techniques presented in this work provide the necessary methodologies for the implementation of the proposed techniques and the results support the validity of the technique.The stability issues are posing threats to power grids with high renewable power generation, and they are not limited to power generation sources, but can also be applicable to other devices with a power electronic converter interface at the grid connection end.In connected systems, autonomous power generation presents a major problem in terms of power quality due to voltage and frequency fluctuations and harmonics generated by the weak grid.The article presented classic controls in generation systems focused on renewable energy based on different studies in the literature for wind, photovoltaic and grid-connected systems, in addition to the stability issues.Furthermore, instability mitigation strategies were also presented in this paper.

Figure 1 .
Figure 1.PV panel power characteristics as a function of the DC voltage and solar irradiance.

Figure 2 .
Figure 2. Single-stage connection of a PV module to the electric grid.

Figure 3 .
Figure 3. Two-stage connection of a PV module to the electric grid.

Figure 4 .
Figure 4. Block diagram of the DC/AC converter control in the dq frame.

Figure 5 .
Figure 5. Simulation of a single-stage PV farm following a change in solar irradiance: (a) irradiance profile; (b) DC voltage response; (c) Power response.

Figure 6 .
Figure 6. Circuit model of a DC/DC boost converter for applications in PV systems.

Figure 7 .
Figure 7. Block diagram of the cascaded PI controller scheme for the DC/DC boost converter.

Figure 8 .
Figure 8. Block diagram of the FCS predictive controller scheme for the DC/DC boost converter proposed in [70].

Figure 9 .
Figure 9.Step response of the current i L using the predictive controller.
Figure 10.Step response of the voltage v PV using the predictive controller.

Figure 11 .
Figure 11.Block diagram of the current control based on power references for DIFG.

Figure 13 .
Figure 13.MPC dynamic response for DIFG: (a) current and voltage responses α components; (b) current response in the stationary frame dq.

Figure 14 .
Figure 14.Block diagram of the grid current control in dq frame using VOC.

Figure 15 .
Figure 15.Finite control set diagram connected to the grid.

Figure 18 .
Figure 18.Dynamic response for voltage variation: (a) DC voltage response; (b) active power response.

Figure 20 .
Figure 20.Single-machine test system for dynamic performance analysis with different short-circuit ratios.

Figure 21 .
Figure 21.Generation source terminal voltage under different SCRs.

Figure 22 .
Figure 22.Single-machine test system for dynamic performance analysis with different short-circuit ratios.

Figure 23 .
Figure 23.System frequency response with wind generation.

Table 1 .
Single-stage PV farm system parameters.

Table 2 .
Parameters of the DFIG and the experimental bench.

Table 3 .
Switching states and voltage vectors of the inverter.