Control of a Three-Phase Current Source Rectiﬁer for H 2 Storage Applications in AC Microgrids

: The share of electrical energy from renewable sources has increased considerably in recent years in an attempt to reduce greenhouse gas emissions. To mitigate the uncertainties of these sources and to balance energy production with consumption, an energy storage system (ESS) based on water electrolysis to produce hydrogen is studied. It can be applied to AC microgrids, where several renewable energy sources and several loads may be connected, which is the focus of the study. When excess electricity production is converted into hydrogen via water electrolysis, low DC voltages and high currents are applied, which needs speciﬁc power converters. The use of a three-phase, buck-type current source converter, in a single conversion stage, allows for an adjustable DC voltage to be obtained at the terminals of the electrolyzer from a three-phase AC microgrid. The voltage control is preferred to the current control in order to improve the durability of the system. The classical control of the buck-type rectiﬁer is generally done using two loops that correspond only to the control of its output variables. The lack of control of the input variables may generate oscillations of the grid current. Our contribution in this article is to propose a new control for the buck-type rectiﬁer that controls both the input and output variables of the converter to avoid these grid current oscillations, without the use of active damping methods. The suggested control method is based on an approach using the ﬂatness properties of differential systems: it ensures the large-signal stability of the converter. The proposed control shows better results than the classical control, especially in oscillation mitigation and dynamic performances with respect to the rejection of disturbances caused by a load step.


Introduction
In recent years, the development of renewable energies has accelerated with the aim of reducing the share of electricity production that is carried out by means that emit high levels of greenhouse gases, such as gas-fired power stations or, even worse, those using coal [1,2]. In this study, we will consider the context of AC microgrids, in which there are several sources of renewable energy and several loads ( Figure 1). The AC network has been used as the standard choice for energy transmission systems since the 19th century for its many advantages such as [3]: • Power transmission over a long distance; • AC voltage magnitude conversion.
The renewable energy sources are, by nature, intermittent and it is, therefore, difficult to balance energy production with consumption when the share of renewables in the energy mix is significant. To compensate for the intermittency of these sources and the fact that their production cannot be controlled, ESSs are a true solution [4]. Among the different types of ESSs, it is possible to highlight two categories: those using batteries [5] and those using hydrogen production by water electrolysis [6]. In the following, we will consider the This solution allows for the maintenance of the grid frequency, even in the case of low power demand. Hydrogen production based on renewable energy or low-carbon power is called green hydrogen. It has significantly lower carbon emissions than grey hydrogen, which is produced by the steam-reforming of natural gas, which constitutes the bulk of the hydrogen market [8]. Green hydrogen appears to be an attractive energy vector to reduce the carbon footprint of the transportation sector and industrial processes. Unfortunately, due to the high cost of production, green hydrogen represents less than 0.1% of total hydrogen production [9].
To convert the excess electricity production into hydrogen, it is necessary to use a converter connected to this AC microgrid to feed the electrolyzer. A state-of-the-art representation of the conversion structures is presented in [10]. Most topologies can be classified into four main categories, which are presented in Figure 2, among which three of them consider controlled converters.
Currently, thyristor-based solutions are the most widely used in industrial applications [11,12], especially in electrolysis, due to the necessity of high current generation to supply high-power electrolyzers. However, the use of thyristors implies the generation of reactive power, which dramatically reduces the power factor. In addition, the output current ripple is quite high, which leads to an increase in energy consumption [13,14].
The second category is a combination of the rectifier with a DC chopper. It avoids the use of passive filters as it can improve the power quality. It is important to combine new emerging DC-DC converters with power management strategies to improve the energy efficiency and the hydrogen flow rate [15].
The third category is the pulse-width modulation (PWM) current-source rectifiers (CSR). They appear as attractive topologies regarding the power quality and their dynamic behavior, even if they seem to be limited to medium-voltage applications to maintain good efficiency, which is in the scope of the application [16].
In this study, the CSR is used to interface the electrolyzer with the grid. To control an electrolyzer, two ways are possible: either impose the voltage or the current. Since the current control, in the case of a cell malfunction, leads to an increase of the electric potential at the cell terminals which may irreversibly damage the cell, the choice was made to use a rectifier that allows one to control the voltage at the electrolyzer terminals. It improves the durability of the systems, especially in the case of cell aging.  As the hydrogen production rate is directly linked to the current flowing through the electrolyzer, it is necessary to use a converter that provides high currents for (relatively) low voltages. Since the typical voltage of a cell is about 1.8 V, and as a classical stack contains 100 of these cells in series, the voltage is about 180 V, whereas the current is about 1000 A. To link the electrolyzer to a classical US grid AC network, the use of a three-phase buck is suggested here, since it is a suitable option for applications that require a lower DC output voltage than the AC mains. The topology has already attracted the attention for different applications such as wind energy conversion systems [17] and on-board electric vehicle chargers [18,19]. This topology is becoming more attractive for motor-drive applications. The classical topology is used in [20], while a novel one is presented in [21] to improve the reliability. Commutation methods in current source inverters (CSI) that use four-quadrant switches are presented in [22] to improve the efficiency compared to conventional CSI with series power diodes.
The classical topology of a buck-type rectifier is studied here because it is a fair compromise between good performance and a low number of semi-conductors to achieve the step-down function. In a single conversion stage, this converter allows one to adjust the output DC voltage, while ensuring a sinusoidal grid current in phase with the main phase voltage [23].
The classical control deals only with output variables of the converter at the price of grid current oscillations, which are undesirable for AC microgrids. The presence of oscillations inside the grid current is explained by the resonance of the L s C filter excited by the harmonics generated by the switching of the transistors. These oscillations are traditionally reduced with damping methods [24]. The proposed synthesis naturally integrates the presence of the input filter and guarantees the stability of the system for all operating points.
The improvement of the control that is proposed here is a stategy based on the control of input and output variables at the same time, in order to directly control grid current oscillations without the need of damping methods. As the proposed control method is based on an approach using the flatness properties of differential systems, it also ensures the large-signal stability of the converter.
The paper is organized as follows: Section 2 describes the operation principle and the modeling of this converter; Section 3 details the proposed control to improve the converter performance; Section 4 deals with the sensibility analysis; and, finally, Sections 5 and 6 display a simulation and experimental results to validate the proposed solution.

Instantaneous Model
The topology is based on an unidirectional CSR with six transistors in series with six diodes to ensure a bidirectional voltage blocking capability. The harmonics injected by the buck rectifier into the microgrid are reduced by an input differential mode L s C filter.
One of the advantages of this topology is the absence of a high start-up current and the naturally induced protection under legs short circuit [25]. Due to the presence of an inductance L d at the output, the current is approximately constant.
At each time, only one high-side SiH and one low-side SiL switch have to be turned on. If, simultaneously, more than one switch is conducting at the upper side or at the lower side, the state of their series diode will be defined by the sign of the phase-to-phase voltage associated with the two bridge legs. Due to voltage ripple and the fact that at each sector crossover, the line voltage is equal to zero, this can lead to the series diodes turning on simultaneously. This generates a distortion of the grid current [26].
A freewheeling diode, D f , is added: without any command on the switches, the DC output voltage is zero. This functionality could be performed by short-circuiting a leg, but at the price of higher switching losses (the pattern would contain a higher number of vectors), and higher conduction losses (the current would flow in a higher number of semiconductors) [27]. The freewheeling diode, D f , also ensures the continuity of the current inside the output inductance, L d . It is conducting when all switches are open. If the switches S1H and S2L (shown in Figure 3) are conducting, the output current, I dc , flows from phase a and returns back into phase b.
In this case, the line-to-line voltage U ab is applied to the output. The maximum output voltage of this converter is [27]: where cos(Φ) is the power factor and U N the grid nominal line voltage. To get the maximum achievable output voltage, the rectifier has to operate at a unity power factor; that is to say that grid current and phase voltage have to be in phase. In this case, the reactive power generated by the L s C input filter is totally compensated. To establish the time varying differential equations of the converter, the boolean switching variables SiH and SiL are considered with i = {1, 2, 3} (1 value for a closed switch, 0 when open). They can be put in a vector form: The equations of this converter can be written as follows (where the exponent T stands for the transpose): As this model uses boolean switching variables, it requires the Space Vector Modulation (SVM) part to determine the dynamics of the converter. Unfortunately, adding this part makes the model more difficult to implement and leads to an increase of the computer simulation time; the characteristic time of the switches is far lower than the time constants of the other variables. This is why it is preferable to work with an average time model in which all variables are replaced by their average values computed over a switching period T.

Space Vector Modulation
The SVM technique consists of controlling the reference vector I re f in the αβ frame in which the converter generates only six active current vectors, with each vector being a combination of different active switches. With six switches, and two states for each switch (open or closed), we have 2 6 = 64 theoretical possibilities; however, due to the operating conditions stated above, there are six active vectors and four zero vectors. The zero-state vector will be performed here by the opening of all the switches, thanks to the freewheeling diode D f .
The importance of the modulation scheme has been studied in [28] and an overview has been presented in [29]. Since the study of the modulation scheme is not the scope of this article, the choice was made to use the six-sector Modified Full wave Symmetrical Modulation (MFSM) scheme ( Figure 4) to facilitate its practical implementation into the FPGA. The change of sector is eased since the zero-state vector is used on both neighboring sectors. Next, the FPGA will generate the command for the drivers of the switches. An SVM is decomposed in three steps:

1.
Identify in which sector the reference vector I re f is; 2.
Compute the application times of adjacent vectors and the zero-state vector during a switching period T; 3.
Choose the modulation sequence during a switching period.
The six vectors shown in Figure 5 correspond to an unique combination of the six switches. For example, the vector − → I 1 corresponds to the following combination: If we directly use these vectors, we are in the full-wave command. To control the reference vector − → I αβ in phase and magnitude, a temporal aspect has to be added: this is the Pulse Width Modulation (PWM). To apply the reference vector, we build an average vector that results from the association of two adjacent vectors, ( − → I j , −→ I j+1 ), and one zero-state vector, − → I 0 , whose order of application depends on the selected pattern. In Figure 5, the area of application of the reference vector that corresponds to the circle inscribed inside the hexagon can be seen. Each active vector has a modulus equal to: The maximum modulus for the reference vector is: The modulation index, M, is defined by the ratio between the reference vector modulus and its maximum: Its value is between 0 and 1.

Average Model
The source phase voltages are assumed to form a direct balanced system: To get almost constant variables, Concordia and Park transformations are performed to move from an abc to a dq frame. The transformations of the variables X are done as follows: The phase voltages are in the Park's reference frame V d = √ 3 V rms and V q = 0. Using (4a) and (9), the grid current equation becomes: The instantaneous expression averaged during one switching period of the input current of the rectifier in α β frame can be obtained using: This allows us to define: Using (12) and (13), the AC capacitor voltage equation becomes: The governing equation of the output current, I dc , is: Moreover, we can write: Therefore, it is possible to write the equation as a power balance between the input and the output of the converter, where I dqre f × V cdq corresponds to the reference of active power consumed by the rectifier.
This leads to the following system of differential equations: This formulation allows one to work with the reference vector I dqre f = [I dre f , I qre f ], where its components I dre f and I qre f are the control variables of the active power and the reactive power consumed by the converter, respectively.
As the model only needs values for I dqre f in a state-space model, its simulation is faster than the whole model with the SVM part.

The Principle of Flatness of Differential Systems
A non-linear systemẋ = f (x, u), composed of a state vector x ∈ R n and an input vector u ∈ R m , is considered as flat if, and only if, there exists a flat output vector y composed of m elements and the functions ϕ : (R m ) r+1 → R n , Ψ : (R m ) r+2 → R m and φ : R n × (R m ) s+1 → R m , such that it is possible to write [30]: where r and s are integers. When a system is flat, the dynamics of the state variables are described by the flat outputs and a finite number of their derivatives. Therefore, if the trajectory of the flat output is well controlled, the trajectory of the state and control variables are also well known, even during transient. The trajectory planning of the flat outputs allows one to generate the command of the system in the case where the system is well known. If the model is exact, the applied commands allow one to reach the desired equilibrium point. Due to parasitic elements and modeling errors, an open-loop operation without status feedback will result in a static error between the desired, and the actual, operating point. To ensure a zero static error, a linearising regulator is used with feedback control as described in [31,32].

The Case of the Three-Phase Buck-Type Rectifier
In this work, the control of the converter will be realized using two loops. Indeed, the model (18) will be divided into two subsets of equations: the first four equations on one side (i.e. the dq components of grid currents and filter capacitor voltages) and the last two equations on the other side (governing the DC output current and voltage).
The link between the AC system and the DC system is the power consumed by the rectifier I dqre f × V cdq : it is taken from the grid and injected in the DC system. The DC variables do not explicitly appear in the AC system; therefore, the AC variables can be treated separately. Two different control loops can be considered.
An internal loop will control the grid current with fast dynamics to reject oscillations based on the principle of flatness control applied to the reduced system defined by (18a-d).
An external loop with a slow dynamic will control the output voltage by regulating the total energy contained in the output inductor and capacitor as shown in Figure 6.

Internal Loop
In the following, we will verify the flatness properties of the AC subsystem described by Equations (1)-(4) of (18). This system contains the following state vector x and control vector u: As this system contains two control variables, it can be considered as flat if, and only if, it exists as a flat output vector y = [y 1 , y 2 ] following (19). We choose the grid currents as the flat outputs associated to the control inputs I dre f and I qre f .
V cd and V cq can be expressed in terms of the flat outputs y 1 , y 2 , and their derivatives: The second derivative of the flat outputs leads to: Finally, we can express the control vector u as a function of y 1 , y 2 , and their derivatives: As the state variables can be expressed as a function of the flat outputs and their derivatives, and the control variables as a function of the flat outputs and their first and second derivatives, then this system can be defined as flat.
A feedback linearization will be used to ensure the control of the flat outputs to their reference trajectories. This technique first introduces a fictitious control, ν = [ν d ν q ], defined as following:ÿ Next, the calculation of the variables ν d and ν q are generally synthesized by imposing a control law that allows setting the dynamics of errors thanks to the following second order differential equation, where the variables y 1re f and y 2re f correspond to I Ldre f and I Lqre f : Finally, it is possible to express the control variables I dqre f based on (22) by replacing y 1 andÿ 2 by the defined fictitious control variables ν d and ν q : Coefficients are chosen in such a way that the polynomials of its characteristic equation have the poles with negative real parts. Therefore, the error converges exponentially to 0 [33].
An identification with a third order polynomial is done to find the coefficients ∀i ∈ {1, 2}: K 1,i , K 2,i and K 3,i .
The pole of the first-order polynomial is placed on the same real part line as poles of the second-order polynomial in order to benefit from good dynamics.
The values of the parameters ξ and ω i (equal to ω 1 and ω 2 , respectively, for the d and q axis of the grid current) define the desired dynamic behavior. Indeed, the values of ω 1 and ω 2 correspond to the desired bandwidth of the current loops. The bandwidth of the controllers must be high enough to attenuate the oscillation,s but a trade-off has to be done with the noise sensitivity.

External Loop
The external loop is based on a power regulator. It consists of controlling the total energy (y) stored inside the output inductance, L d , and the output capacitor, C dc . It is possible, based on the variation of the total electromagnetic energy storedẏ, to bring out a power balance between the input and the output of the rectifier.
If the same power balance is made on the AC side: The left-hand side of the equation is the variation of energy of L s and C. The term V cd I dre f + V cq I qre f is the coupling term between AC and DC (power sent from the AC system to the DC system), and it vanishes if the whole system AC + DC ((28) + (29)) is taken into account. If we assume that both dynamics of variations of the magnetic energy stored in the input filter, L s C, (slow variations) and the serial resistance, r s , are neglected, thenẏ can also be written as:ẏ It is then possible to express the reference of active power: The bandwidth selected for the external loop (which corresponds directly to the value of ω BF ) must be relatively low when compared to the one of the internal loop (ω 1 and ω 2 ). The values of the control parameters are defined in Table 1. Similarly, the following control law will be used to computeẏ, where the integral component allows one to reject modeling errors and y re f corresponds to the total reference of energy obtained by replacing in (28) V dc by V dcre f .
y =ẏ re f + K 1 (y re f − y) + K 2 (y re f − y) (33) This time, the coefficients will be classically identified with a second-order polynomial:

BF
In this section, we proposed a new control for the three-phase buck-type rectifier. While the classical control of this converter is based only on the control of the output variables, the proposed control deals with input and output variables. It is based on the flatness properties of differential systems. It uses two loops with large differences of bandwidth. The internal loop directly controls the grid current to damp the oscillations and controls the input variables of the rectifier. The external loop with a low dynamic controls the total energy stored in the output inductor, L d , and capacitor, C dc . The proposed solution does not depend on the operating point compared to the classical one using cascaded PI controllers. In addition, the tuning of the proposed controllers are simpler than the classical solution because the desired bandwidth is selected directly by the variables ω 1 , ω 2 , and ω BF .

Sensitivity Analysis
The aim of this analysis is to study the stability of the system in a closed loop under parameter variations by observing the location of eigenvalues. This analysis is based on the average model of the converter added with the control regulators presented before. To achieve this, the enhanced system is written in the following form and the eigenvalues of the state matrix are determined after using the jacobian at the operating point (Test I in Table 2), where y, y 1 , and y 2 are related to the control part.
We will study the stability of the system with variations of the AC filter parameters (C and L s ). Indeed, these values can change during operation and time. To realize this study, control laws will use constant values of C = 39 µF and L s = 225 µH, while system values will be increased, respectively, from 20 µF to 60 µF and from 110 µH to 330 µH. The other parameters used to realize this test are listed in Tables 1 and 3. Figure 7 presents the stability analysis of the system in a closed loop with parameter variations. Errors on AC filter parameters are introduced and the sign of eigenvalues are studied to make conclusions on the stability of each point.
We found that even with ±50% errors on both values of the AC filter, the system remains stable because all the eigenvalues have a negative real part.

Stability analysis under parameters variations
Case when control and converter parameters are identical Stable cases  Tables 1 and 3, and the operating point corresponds to Test I in Table 2.

Simulation Results
Tables 1 and 4 display the control parameters which will be used for the simulation and experimental tests. Table 3 shows the parameters of the buck-type current-source rectifier used in the simulation and the experimental tests. In Figure 8, the simulation of the instantaneous model with SVM and the one of the average models are carried out, when the system is in an open loop (meaning that the values for both control variables I dre f and I qre f are set directly). Both simulations show the same dynamic behavior, which validates the approximation of the average model.
To ensure good continuity of service of the electrolyzer, the DC-output voltage of the rectifier has to be maintained as much as possible during a sudden voltage drop of the network. The simulation (using the average model) of the system during a drop of the rms grid voltage (reduction from 55 V to 40 V) is displayed in Figures 9 and 10 for an electrolyzer with 25 cells in series (2 V at the terminals of each cell). Figure 11 shows the system behavior for grid voltage reduction from 110 V to 90 V for an electrolyzer with 50 cells in series (2 V at the terminals of each cell). The system remains stable during the disturbance: this one is removed quickly after 50 ms. There is a transient voltage drop on the DC system of about 5 V in both cases, followed by a minor overvoltage. The DC current oscillation is about 1A in both cases.  Table 3 when the voltage reference is increased from 100 V to 150 V on a 10 Ω resistive load with a voltage grid at 110 V.   Figure 10 displays the influence of the value of the control parameter ω BF on the output voltage V dc . It can be seen that the higher the value of ω BF is, the faster the disturbance is removed. In addition, it is important to notice that a high value of ω BF generates a lower voltage drop.  Tables 1 and 3. A simulation is performed to compare the classical control using two cascaded PI controllers with and without the presence of an active damping method with the proposed solution. PI controllers can be written as follows for the current and voltage loops: The active damping method used consists of emulating a damping resistor, without any additional power loss by making the CSR produce an additional current added to the reference current of the d axis (I dre f ). It is based on the voltage measurement across the AC capacitor. After using a high pass filter to keep only the alternative part, the obtained voltage is divided by the virtual resistor to generate the damping current [34].
The points of comparison will be: the magnitude of the oscillations on the grid current in steady state and the performances under disturbances generated by a resistive load whose value is changed from 10 Ω to 20 Ω. The two controllers are tuned to obtain the same bandwidth in a closed loop in each case. Two cases are presented with each control corresponding to different settings: one is relatively slow (cutoff frequency of the voltage loop is 85 rad·s −1 ) and one is relatively fast (cutoff frequency of the voltage loop is 175 rad·s −1 ).
In Figure 12, it is possible to see the grid current in steady state with: the classical control without an active damping method (top), the classical control with an active damping method (middle), and the proposed control (bottom). Converter parameters are listed in Table 3 with a voltage grid at 110 V and a resistive load of 10 Ω. Control parameters are listed in Tables 1 and 4. Controllers are tuned to be slow on the left side and fast on the right side. It is possible to highlight that the presence of an active damping method reduces the magnitude of oscillations (from 3A to 1A) compared to the solution without. However, the proposed control shows a greater reduction of the magnitude of oscillations (from 3 A to 0.15 A). In addition, to highlight this result, the Total Harmonics Distorsion (THD) has been computed in the different cases and presented in Figure 13. It appears that the proposed solution has the lowest THD value of the grid current.
The behavior of the system subjected to a sudden change of the load has been tested. A comparison between the classical control with and without the active damping method is first presented in Figure 14. This figure confirms that the introduction of the damping method does not affect the dynamic behavior of the control. Next, a comparison between the classical and the proposed control is given in Figure 15. The classical and the proposed solutions are compared in two cases depending on the tuning of the control parameters. One case is considered as slow and the other as fast. In the first case where the controllers are slow, the overvoltage on V dc is reduced from 50 V to 30 V with the proposed solution. Additionally, the voltage drop with the classical control is about 10 V when there is no voltage drop with the proposed control. Regarding the output current, I dc , there is a current drop of about 1.5 A with the classical control, while it is equal to zero with the proposed control.
In the second case where the controllers are faster, with the proposed control, the overvoltage on V dc is reduced from 32 V to 22 V. Additionally, while there is no voltage drop with the proposed control, the voltage drop with the classical control is about 5 V. In this case, with the proposed control, the current drop on I dc is not eliminated, but reduced from 1.5 A to 0.8 A.
From these results (resumed in Table 5), it is then possible to see that in both cases (slow and fast), the proposed control shows better dynamic performances compared to the classical solution. Furthermore, it allows one to reduce the magnitude of the grid current oscillations, which leads to an improvement in the value of the grid current THD.  Tables 3 and 4.

Experimental Results
In order to verify the proposed solution, a test bench has been set up as shown in Figure 16. It is first composed of a three-phase controllable AC source Chroma 61702. Second, a MicroLabBox DS1202 is used with Matlab/Simulink and ControlDesk to monitor the state variables of the system and compute the reference variables. Next, these variables are sent to a FPGA DE0-Nano, which is responsible for sending command signals to the semiconductors using optical fibers. The FPGA has been programmed using Quartus software in VHDL language. Finally, the system also includes an active load (EA Elektro-Automatik 9750-22) that can be controlled with analog inputs directly coming from the DAC port (Digital Analog Converter) of the MicroLabBox.
The quasi-static characteristic of an electrolyzer has been implemented in the active load. This characteristic allows us to emulate the behavior of a selected electrolyzer and establish the operating point according to the voltage imposed on its terminals. A typical characteristic of one cell is presented in Figure 17: it displays the current density (A/cm 2 ) as a function of the voltage (V) across the cell. A stack corresponds to an association between serial and parallel cells. In this study, the surface is selected at around 7 cm 2 with 50 cells in series for Test I and Test II, and 25 cells in series for Test III. The power extracted by the electrolyzer as a function of the voltage at its terminals is shown in Figure 17. In order to extract the maximum power while maintaining the electrolyzer in good conditions, it was decided to apply 2 V to the terminals of each cell as an operating point for Test I and Test III. Another operating point has been tested with Test II, corresponding to 1.8 V to the terminal of each cell.

Operating Point
Test I

Operating Point
Test II Figure 17. Typical characteristic of a PEM cell (left) and power as a function of voltage across a stack of 50 cells in series (right). Figures 18 and 19 present the behavior of the system in steady-state with the proposed control during Test I and Test II. The parameters during these tests are indicated in Tables 1-3. In these figures, it is possible to see the grid voltage V a , the grid current I La , the output current I dc , and the output voltage V dc . It can be seen that grid voltage, V a , and current, I La , are in phase. It means that the reactive power, coming in particular from the differential mode input, L s C, filter is totally compensated. The proposed control makes it possible to operate with a unit power factor. In addition, the oscillations on the grid current are damped using the proposed control without the need of any damping element or virtual active damping control.
Idc(5A/div) ILa(5A/div) Va(50V/div) Vdc(100V/div) Figure 18. Experimental results in a steady state when the voltage at the terminals of each cell (50 cells in series) is set to 2 V. Test I ( Table 2) with parameters in Tables 1 and 3.
Idc(5A/div) ILa(5A/div) Va(100V/div) Vdc(100V/div) Figure 19. Experimental results in a steady state when the voltage at the terminals of each cell (50 cells in series) is set to 1.8 V. Test II (Table 2) with parameters in Tables 1 and 3. Figure 20 displays the behavior of the system in a steady state with the proposed control during Test III. The parameters during this test are shown in Tables 1 and 2. We obtain the same results as in Test I and Test II: no oscillations on the grid current, grid voltage and current are in phase, and there are low ripples on the output current. In Figure 21, it is possible to see the value of the THD of the grid current during the test. The definition chosen for the THD is the ratio between the amplitude of the harmonics and the fundamental value of the signal (Total Harmonic Distortion compared to the fundamental). Here, it is less than 5%, which indicates a grid current with very few disturbances in the network induced by the rectifier.  It is important, in order to guarantee a good continuity of service of the electrolyzer, to maintain the desired voltage at its terminals, for example during a micro-cut of the electrical network, i.e., during a sudden variation of its rms voltage. In order to study the response of the system during a sudden drop in the rms voltage of the network, Test III was carried out. In Figure 22, it is possible to see the grid voltages, V abc , the output current, I dc , and voltage, V dc , during a sudden drop of the rms grid voltage from 55 V to 40 V with two different values of ω BF . The system remains stable and the desired output voltage is maintained in both cases. Further, as in simulation (Figure 9), it appears that a high value of ω BF , allows one to reduce the voltage drop (7 V instead of 9 V) and removes the perturbation rapidly (50 ms instead of 100 ms).
The simulation corresponding to Test III was displayed in Figure 9. The behavior of the system is similar between the simulation and experimental tests; for example, with ω BF = 120 rad·s −1 , the DC transient voltage drop is about 7 V (5 V by simulation) and the current drop 1.25 A (1.2A by simulation). The DC current exceeding the set value is less than that in the simulation: it appears to be damped. It may be explained by modeling errors or a misidentification of parameters, but also the way in which the rms voltage drop is generated.

Conclusions
AC microgrids, composed of several renewable energy sources and several loads, need solutions to regulate the power produced as a function of the power consumed. To cope with the intermittency of these sources, ESS based on the electrolysis of water to produce hydrogen is a popular option. To improve the connection of the electrolyzer to the AC grid, the use of a single-stage conversion, three-phase buck-type rectifier has been evaluated.
First, an average model to model the rectifier has been described. Simulation times are faster than those of the instantaneous model (since it does not require the SVM part) without altering the system dynamics.
Next, we proposed a new control of this converter based on the flatness properties of the differential systems. It allows one to control the oscillations of the L s C input filter without using active damping method and guarantees the large signal stability of the converter. The proposed solution is easy to tune, does not depend on the operating point, and presents better results regarding oscillation mitigation and dynamic performance than the classical control with active damping.
A simulation and experimental tests have been carried out to prove the efficiency of the solution. An active load has been used to emulate the electrolyzer, described by its quasi-static characteristic. In a steady state, the unitary power factor was achieved. Further, the THD of the grid current is kept low (4%) during operation.
Finally, the behavior of the converter has been studied under a sudden drop of the rms grid voltage. We noticed that the converter remains stable and removes the perturbation quickly. This allows one to maintain the operating condition of the electrolyzer and guarantees a good continuity of service.
The proposed solution is a suitable one-stage conversion option to use the excess of electricity production to store hydrogen based on water electrolysis, to maintain a constant voltage at the terminals of the electrolyzer in the case of a grid disturbance, all while the grid current is stabilized with a unity power factor.