Proton Exchange Membrane Electrolyzer Modeling for Power Electronics Control: A Short Review

: The main purpose of this article is to provide a short review of proton exchange membrane electrolyzer (PEMEL) modeling used for power electronics control. So far, three types of PEMEL modeling have been adopted in the literature: resistive load, static load (including an equivalent resistance series-connected with a DC voltage generator representing the reversible voltage), and dynamic load (taking into consideration the dynamics both at the anode and the cathode). The modeling of the load is crucial for control purposes since it may have an impact on the performance of the system. This article aims at providing essential information and comparing the di ﬀ erent load modeling.


Introduction
Hydrogen is the most innumerable and simplest element on Earth. It can store and deliver usable energy. However, it does not exist by itself in nature and must be made from different elements that contain it. For instance, it can be combined with carbon (e.g., oil, natural gas) and with oxygen in water (H 2 O) [1]. Hydrogen has the highest specific energy per kilogram of all fuels (i.e., 120-140 MJ/kg), but its energy density is less suitable for storage (i.e., 2.8-10 MJ/L) depending on the physical-based storage (e.g., compressed (350-700 bar), liquid) [2]. On one hand, the global hydrogen production from natural gas, coal, and oil by using the reformation process represents approximately 96%. On the other hand, the use of the water electrolysis process to split the deionized water into hydrogen and oxygen represents around 4% of the global hydrogen production [3]. Albeit hydrogen being an intrinsically clean energy vector, it requires energy to be produced; the kind of adopted source makes the difference. Hydrogen produced by fossil fuels is known as grey hydrogen due to indirect pollution.
To supply the water electrolysis process, renewable energy sources (RES) (e.g., wind turbines, photovoltaic) are the most suitable since they can limit the environmental impact. In this way, the so-called green hydrogen is obtained. Blending this kind of hydrogen into the existing natural gas pipeline network has been proposed as a means of increasing the output of renewable energy systems. Delivering blends of hydrogen and methane by pipeline also has a long history; recently, the rapid growth in installed wind power capacity and interest in the near-term market readiness of fuel cell electric vehicles has increased the stakeholder interests [4,5]. The water electrolysis process is carried out by an electrolyzer. Nowadays, there are three main types of electrolyzers: alkaline electrolyzer (AEL), proton exchange membrane electrolyzer (PEMEL), and solid oxide electrolyzer (SOEL) [3,6]. Currently, AEL and PEMEL are commercially available, and SOEL is still under research and development. AEL is the oldest technology widespread around the world, commonly used for large-scale systems. The advantage of this technology is its long lifespan and lower capital cost due to the cheaper catalysts based on Nickel material. However, AELs have several drawbacks such as low current density and operational pressures [6]. Besides, their response time in case of dynamic operations is slower compared to PEMELs. Indeed, since AEL is based on a liquid electrolyte leading up to more inertia, the ion transportation is slower compared to PEMEL [7]. This issue is particularly crucial when an electrolyzer has to be operated under fast dynamics, usually met when RES are adopted. Indeed, a RES depends on weather conditions (as solar irradiance or wind speed) that can vary abruptly during the operation. For this reason, this technology is not generally considered to be coupled with RES; conversely, it is employed for high-power stationery operated plants supplied by the grid where the harmonics impact on the grid has to be considered [8][9][10].
In comparison, PEMELs are generally used for small-scale hydrogen production. Anyway, PEMELs can be connected in series to form a multi-stack electrolyzer for large-scale systems. Furthermore, PEMELs offer some advantages compared to AELs, such as high current densities, high cell voltage efficiency, and fast system response when operating dynamically [3,11], since PEMEL has faster ion transportation due to the solid and thin membrane. As a result, this review article is focused on PEMEL technology since it can be connected with RES to produce hydrogen.
On the other side, power converters are mandatory for the hydrogen production system supplied by RES to interface the RES with PEMEL. In these systems, DC-DC converters play a crucial role since most of RESs generate high DC voltage output and PEMELs need to be supplied with low DC voltage input [12]. However, the system needs an AC-DC converter if a wind turbine conversion system or power grid is used [10,13]. Dealing with PEMEL supplied by power converters, it is clear that two systems with their dynamics have to interact with each other. For this reason, the PEMEL model needs a deep investigation; it becomes of the utmost importance when RES is employed since the delivered power can vary with a dynamic that solicits the whole conversion chain.
To simulate the hydrogen production system by using PEMEL supplied by RES, it is necessary to know the PEMEL electrical characteristics. The model must replicate the real behavior of PEMEL (i.e., static, dynamic operations) to investigate its interaction with power converters. In the literature, three types of models have been reported (e.g., resistive load, static load, and dynamic load) [14][15][16][17][18][19][20][21][22][23][24][25]. The last modeling comes from the need to investigate the dynamic of the conversion chain with RES.
The main objective of this article is to present each model used for power electronics control and how these models have been developed. Besides, the validity of these models is discussed. Finally, some guidelines are provided to determine the parameters of the proposed models (i.e., static and dynamic).
This article is divided into three sections. After the introduction presenting the current state-of-the-art, and the reasons to carry out such review work, Section 2 is focused on the proton exchange membrane electrolyzer system for hydrogen production including the DC-DC converter. Finally, a review of proton exchange membrane electrolyzer modeling for control propose and a comparison of the static and dynamic models from the control point of view are provided in Section 3.

Proton Exchange Membrane Stack
The first proton exchange membrane electrolyzer has been developed by the General Electric (GE) Company in the 1960s. This technology aims at overcoming the weaknesses of AEL technology [3]. The PEMEL converts the electrical energy directly into chemical energy. The principle of PEMEL is shown in Figure 1. The water is oxidized electrochemically within the anode catalyst layer combine with electricity at the anode side (W ele , including Gibbs energy and energy losses); then water is broken down into oxygen gas, protons, and electrons, and this reaction called oxidation reaction is given by (1). The protons transport across the proton conductive membrane to the cathode side while the electrons move through the outer circuit and also reach the cathode side. The protons and the electrons electrochemically react within the cathode catalyst layer and produce hydrogen gases. This reaction, called reduction reaction, and the overall reaction, called redox reaction, are given by (2) and (3), respectively [3]. It has to be remarked that the redox reaction (3) requires energy to be performed, whereas the reduction (2) is almost spontaneous; it justifies the need to supply the PEMEL and the different dynamics expected by the two reactions.
with electricity at the anode side (Wele, including Gibbs energy and energy losses); then water is broken down into oxygen gas, protons, and electrons, and this reaction called oxidation reaction is given by (1). The protons transport across the proton conductive membrane to the cathode side while the electrons move through the outer circuit and also reach the cathode side. The protons and the electrons electrochemically react within the cathode catalyst layer and produce hydrogen gases. This reaction, called reduction reaction, and the overall reaction, called redox reaction, are given by (2) and (3), respectively [3]. It has to be remarked that the redox reaction (3) requires energy to be performed, whereas the reduction (2) is almost spontaneous; it justifies the need to supply the PEMEL and the different dynamics expected by the two reactions.
Oxidation reaction Reduction reaction In PEMEL, the key components are bipolar plates with flow channels, current collectors, and membrane electrode assembly (MEA) [26]. These components are shown in Figure 2 giving a crosssection of mass transport in a PEMEL. The bipolar plates are mandatory and have crucial functions in PEMEL operation, such as conducting electrons, connecting single cells to realize a stack, arranging a flow path for pure water-sharing over the current collectors, isolating hydrogen and oxygen, supporting the membrane and electrodes, and bringing thermal conduction to handle the PEMEL temperature [26]. To meet these functions, bipolar plates must feature specific properties such as high thermal and electrical conduction, low gas permeability, and high mechanical and corrosion resistance. Generally, bipolar plates based on Titanium (Ti) material are employed since Ti has a high mechanical and corrosion resistance. Oxidation reaction In PEMEL, the key components are bipolar plates with flow channels, current collectors, and membrane electrode assembly (MEA) [26]. These components are shown in Figure 2 giving a cross-section of mass transport in a PEMEL. The bipolar plates are mandatory and have crucial functions in PEMEL operation, such as conducting electrons, connecting single cells to realize a stack, arranging a flow path for pure water-sharing over the current collectors, isolating hydrogen and oxygen, supporting the membrane and electrodes, and bringing thermal conduction to handle the PEMEL temperature [26]. To meet these functions, bipolar plates must feature specific properties such as high thermal and electrical conduction, low gas permeability, and high mechanical and corrosion resistance. Generally, bipolar plates based on Titanium (Ti) material are employed since Ti has a high mechanical and corrosion resistance. be converted into hydrogen. For this reason, the converter needs a maximum power point tracking. Furthermore, to increase the input voltage of PEMELs, multi-stack PEMELs can be used, consisting of connecting several PEMELs in series to meet the requirements (i.e., voltage ratio, hydrogen flow rate, energy efficiency). Finally, power converters slightly decrease the overall efficiency. Since the efficiency of power converters is generally higher than 92%-95%, the expected lowering of the whole conversion chain (including the PEMEL), is around 2%-6%.

AC-DC and DC-DC Converters
Power converters can be classified based on the type of RES or bus systems configuration (i.e., DC or AC). From the literature, four types of hydrogen production systems including PEMELs supplied by solar energy, wind energy, DC microgrid system, and the power grid have been reported. In any case, step-down DC-DC converters are needed to supply PEMELs with a low DC voltage (i.e., less than 10 V) [25,34]. In comparison, current collectors are used to move the electrons from the catalyst layer to bipolar plates and to eliminate the gases (i.e., hydrogen, oxygen) from the catalyst layer [26]. It is important to point out that current collectors must operate in the same environment conditions (i.e., acid, high overvoltage) as bipolar plates. Therefore, porosity and electrical conductivity are important requirements in choosing the materials of current collectors. On one hand, in PEM fuel cells, carbon (C) paper or cloth is employed as current collectors (commonly known as gas diffusion layer (GDL)) at both sides of the electrodes. On the other hand, in PEMELs, carbon paper or cloth is not suitable for either the catalyst layer or current collector of the anode side since the anode potential is higher than the cathode potential [26]. This high potential leads to the corrosion of the carbon material during the water electrolysis process. As a result, at the anode side, a material based on Ti is preferred since it features low corrosion even under high anode potential and acid environment. Regarding the cathode side, its potential is smaller than the anode side and close to that for PEM fuel cell operation. Therefore, a material based on carbon is perfectly fit not only for the catalyst layer but also for the current collector.
Finally, MEA (i.e., PEM in Figure 2) is essential to move the protons from the anode to the cathode side and drive the electrons to migrate around a flow path to the cathode. Moreover, the MEA allows the electrical insulation between the anode and cathode while acting as a reactant barrier against gas crossover [27]. Currently, fluoropolymer (PFSA) Nafion membranes from DuPont Company are the most widespread in PEMELs, since they feature high thermal stability and proton conductivity (i.e., 0.1 S cm −1 at 100 • C) and thin membrane (25-254 µm). The choice of the membrane thickness results in a compromise between the expected operating pressures across the membrane, mechanical resistance, low gas crossover, and ohmic resistance [27]. Different types of Nafion membranes are available in the market, such as Nafion TM 211, Nafion TM 212, Nafion TM 115, Nafion TM 117, and Nafion TM 1110, depending on their thickness [28]. Based on previous works reported in the literature [27,29], it has been emphasized that very thin membranes allow reducing ohmic losses and operating the membrane at high pressure due to its high mechanical resistance. However, the higher the operating pressure, the higher the equivalent current of hydrogen crossover [30]. As a result, faradaic losses due to the hydrogen and oxygen crossover increase. Faradaic losses related to Faraday's efficiency are particularly noticeable at low current densities. Since the energy efficiency of the PEMEL is linked to the Faraday's efficiency, it leads to a decrease of energy efficiency. This important issue has been discussed in [31]. Furthermore, equivalent currents of hydrogen and oxygen crossover are inversely proportional to the membrane thickness [30]. The higher the membrane thickness, the smaller the currents of hydrogen and oxygen crossover. High currents of hydrogen and oxygen crossover may damage the membrane and cause failures of the PEMEL stack [32]. For this reason, membrane thickness and operating pressures are key issues to enhance the performance of the PEMEL [31]. Currently, PFSA Nafion membranes C 2020, 6, 29 5 of 20 suffer from having a high cost, and their performances are strongly dependent on the relative humidity (RH) and temperature (T), which impede their development at large scale and market penetration [26]. Indeed, these two operating parameters play a key role in drying or flooding the membrane linked to its water content. The membrane must be perfectly hydrated to enhance the move of protons from the anode and cathode and to reduce the gas crossover. To cope with this important issue from the cost and performance point of view, membranes based on hydrocarbon material have much to offer. Indeed, this type of membrane is a promising solution to replace Nafion membranes, since it presents several advantages such as low-cost, fit proton conductivity (>80 mS cm −1 ), low gas permeation, and able to operate at low RH (<25%) and large temperature range (from −40 • C to 120 • C) [33].
In PEMEL system, the stack is not only the key component. Indeed, PEMEL system includes also several ancillaries such as power conditioning system (i.e., AC-DC and/or DC-DC converter), pure water circulation system (i.e., tanks, circulation pump, and valves), hydrogen processing system (i.e., storage tanks, hydrogen separator, and valves), and cooling system (i.e., cooling pumps, valves) [26]. In Figure 3, a PEMEL system is provided showing only the stack and power electronics part. Given that PEMELs operate at low DC voltages, it is not feasible to directly couple them with the power grid or RES in the hydrogen production system. For this reason, power converters are mandatory to meet the requirements of PEMELs, particularly in terms of low conversion gain. The control of the converters differs depending on the source. For grid-connected applications, the available power significantly overcomes the power required by the electrolyzer; hence, it can be regulated based on the quantity of hydrogen to be produced. Differently, a RES requires that all the available power must be converted into hydrogen. For this reason, the converter needs a maximum power point tracking. Furthermore, to increase the input voltage of PEMELs, multi-stack PEMELs can be used, consisting of connecting several PEMELs in series to meet the requirements (i.e., voltage ratio, hydrogen flow rate, energy efficiency). Finally, power converters slightly decrease the overall efficiency. Since the efficiency of power converters is generally higher than 92-95%, the expected lowering of the whole conversion chain (including the PEMEL), is around 2-6%. C 2020, 6, x FOR PEER REVIEW 5 of 21 meet the requirements of PEMELs, particularly in terms of low conversion gain. The control of the converters differs depending on the source. For grid-connected applications, the available power significantly overcomes the power required by the electrolyzer; hence, it can be regulated based on the quantity of hydrogen to be produced. Differently, a RES requires that all the available power must be converted into hydrogen. For this reason, the converter needs a maximum power point tracking. Furthermore, to increase the input voltage of PEMELs, multi-stack PEMELs can be used, consisting of connecting several PEMELs in series to meet the requirements (i.e., voltage ratio, hydrogen flow rate, energy efficiency). Finally, power converters slightly decrease the overall efficiency. Since the efficiency of power converters is generally higher than 92%-95%, the expected lowering of the whole conversion chain (including the PEMEL), is around 2%-6%.

AC-DC and DC-DC Converters
Power converters can be classified based on the type of RES or bus systems configuration (i.e., DC or AC). From the literature, four types of hydrogen production systems including PEMELs supplied by solar energy, wind energy, DC microgrid system, and the power grid have been reported. In any case, step-down DC-DC converters are needed to supply PEMELs with a low DC voltage (i.e., less than 10 V) [25,34].

AC-DC and DC-DC Converters
Power converters can be classified based on the type of RES or bus systems configuration (i.e., DC or AC). From the literature, four types of hydrogen production systems including PEMELs supplied by solar energy, wind energy, DC microgrid system, and the power grid have been reported. In any case, step-down DC-DC converters are needed to supply PEMELs with a low DC voltage (i.e., less than 10 V) [25,34]. First, the hydrogen production systems supplied by solar energy consist of using photovoltaic (PV) panels, step-down DC-DC converter, and PEMEL. The illustration diagram of this system is presented in Figure 4. The photovoltaic enables harvesting energy from the sun to generate DC voltage and then coupling with the buck (step-down) converter to supply DC current to the PEMEL to produce hydrogen gas. It should be remarked that operating with PV panels, the available power depends on solar irradiance. To maintain a high efficiency of the conversion chain, all this power has to be used to supply the PEMEL. It is guaranteed by the maximum power point (MPPT) algorithm that both controls the step-down converter and supplies the suitable current to the PEMEL. C 2020, 6, x FOR PEER REVIEW 6 of 21 First, the hydrogen production systems supplied by solar energy consist of using photovoltaic (PV) panels, step-down DC-DC converter, and PEMEL. The illustration diagram of this system is presented in Figure 4. The photovoltaic enables harvesting energy from the sun to generate DC voltage and then coupling with the buck (step-down) converter to supply DC current to the PEMEL to produce hydrogen gas. It should be remarked that operating with PV panels, the available power depends on solar irradiance. To maintain a high efficiency of the conversion chain, all this power has to be used to supply the PEMEL. It is guaranteed by the maximum power point (MPPT) algorithm that both controls the step-down converter and supplies the suitable current to the PEMEL. Second, wind energy systems are composed of the wind turbine, step-down DC-DC converter, and PEMEL. However, wind turbines allow exploiting the kinetic energy of the wind to be transformed into AC voltage at the output of a three-phase generator. The system needs a three-phase rectifier AC-DC converter interfacing the wind turbine, and the buck converter. The schematic diagram of this system is shown in Figure 5. The wind generation system shows a fast dynamic due to gusts. It is partially mitigated by the DC bus capacitor C1, but it must operate with the power transfer to the PEMEL. On one hand, if the available power rises abruptly and it is not transferred to the PEMEL, the voltage of C1 can overcome its rated voltage with risk of damaging. On the other hand, the power flowing through the PEMEL has to be varied without voltage overshoots that could damage the same PEMEL. The architecture shown in Figure 5 is widely used for permanent magnet synchronous generator (PMSG)-based wind turbines. The simple diode rectifier on the generator side gives a cost-efficient solution since active power flows unidirectionally, and no reactive power is required. Unfortunately, the uncontrolled rectifier might cause low-frequency torque pulsation exciting shaft resonance. The availability of cheap controlled power devices and microprocessor platforms allows conceiving more complicated schemes with better performance. The uncontrolled rectifier can be substituted by a PWM active rectifier, as in Figure 6; it enables full power controllability (four-quadrant operation). The simple structure (available as an integrated power module) is robust and reliable. Besides many Second, wind energy systems are composed of the wind turbine, step-down DC-DC converter, and PEMEL. However, wind turbines allow exploiting the kinetic energy of the wind to be transformed into AC voltage at the output of a three-phase generator. The system needs a three-phase rectifier AC-DC converter interfacing the wind turbine, and the buck converter. The schematic diagram of this system is shown in Figure 5. The wind generation system shows a fast dynamic due to gusts. It is partially mitigated by the DC bus capacitor C 1 , but it must operate with the power transfer to the PEMEL. On one hand, if the available power rises abruptly and it is not transferred to the PEMEL, the voltage of C 1 can overcome its rated voltage with risk of damaging. On the other hand, the power flowing through the PEMEL has to be varied without voltage overshoots that could damage the same PEMEL. First, the hydrogen production systems supplied by solar energy consist of using photovoltaic (PV) panels, step-down DC-DC converter, and PEMEL. The illustration diagram of this system is presented in Figure 4. The photovoltaic enables harvesting energy from the sun to generate DC voltage and then coupling with the buck (step-down) converter to supply DC current to the PEMEL to produce hydrogen gas. It should be remarked that operating with PV panels, the available power depends on solar irradiance. To maintain a high efficiency of the conversion chain, all this power has to be used to supply the PEMEL. It is guaranteed by the maximum power point (MPPT) algorithm that both controls the step-down converter and supplies the suitable current to the PEMEL. Second, wind energy systems are composed of the wind turbine, step-down DC-DC converter, and PEMEL. However, wind turbines allow exploiting the kinetic energy of the wind to be transformed into AC voltage at the output of a three-phase generator. The system needs a three-phase rectifier AC-DC converter interfacing the wind turbine, and the buck converter. The schematic diagram of this system is shown in Figure 5. The wind generation system shows a fast dynamic due to gusts. It is partially mitigated by the DC bus capacitor C1, but it must operate with the power transfer to the PEMEL. On one hand, if the available power rises abruptly and it is not transferred to the PEMEL, the voltage of C1 can overcome its rated voltage with risk of damaging. On the other hand, the power flowing through the PEMEL has to be varied without voltage overshoots that could damage the same PEMEL. The architecture shown in Figure 5 is widely used for permanent magnet synchronous generator (PMSG)-based wind turbines. The simple diode rectifier on the generator side gives a cost-efficient solution since active power flows unidirectionally, and no reactive power is required. Unfortunately, the uncontrolled rectifier might cause low-frequency torque pulsation exciting shaft resonance. The availability of cheap controlled power devices and microprocessor platforms allows conceiving more complicated schemes with better performance. The uncontrolled rectifier can be substituted by a PWM active rectifier, as in Figure 6; it enables full power controllability (four-quadrant operation). The simple structure (available as an integrated power module) is robust and reliable. Besides many The architecture shown in Figure 5 is widely used for permanent magnet synchronous generator (PMSG)-based wind turbines. The simple diode rectifier on the generator side gives a cost-efficient solution since active power flows unidirectionally, and no reactive power is required. Unfortunately, the uncontrolled rectifier might cause low-frequency torque pulsation exciting shaft resonance. The availability of cheap controlled power devices and microprocessor platforms allows conceiving more complicated schemes with better performance. The uncontrolled rectifier can be substituted by a PWM active rectifier, as in Figure 6; it enables full power controllability (four-quadrant operation). The simple C 2020, 6, 29 7 of 20 structure (available as an integrated power module) is robust and reliable. Besides many turbines where the electrical generator is a squirrel cage induction generator (SCIG), a double-fed induction generator (DFIG) is adopted in this scheme in the back-to-back configuration for grid connection [35][36][37].
C 2020, 6, x FOR PEER REVIEW 7 of 21 turbines where the electrical generator is a squirrel cage induction generator (SCIG), a double-fed induction generator (DFIG) is adopted in this scheme in the back-to-back configuration for grid connection [35][36][37]. In low-power wind turbines with a PMSG, the three-phase high-frequency semicontrolled rectifier, shown in Figure 7, is attractive, since it is simple and robust. Moreover, all active switches are connected to a common point, and a short-circuit through a leg is not possible. A higher distortion of the generator currents results [38]. A three-level neutral-point diode clamped front-end topology (3L-NPC) gives two DC outputs that can be connected to two electrolyzers; each DC level has one half of the voltage obtained by the two levels of the active rectifier. It is sketched in Figure 8. The 3L-NPC topology is one of the most commercialized multilevel converters on the market, and it is usually proposed as a back-to-back topology in wind turbines. A potential drawback is the midpoint voltage fluctuation of the DC bus; however, this problem is minimized by the controlling of the redundant switching status [35,36,39].
Third, like hydrogen production systems supplied by solar energy, the DC microgrid system is shown in Figure 9 (the input filter is not drawn for the sake of clarity). The input of the buck converter in this system interfaces with the DC bus in DC microgrid and the output coupling with PEMEL [13]. In this case, the dynamic behavior depends on the need to use the hydrogen as storage depending on the constraints of the microgrids.
Finally, the hydrogen production system connected with the utility grid is present in Figure 10. This system consists of a delta-star transformer and an AC-DC rectifier interfacing the grid system and PEMEL. Since the current of the electrolyzer must be controlled, thyristors-based rectifiers are generally preferred for this purpose [8][9][10]. In Figure 10, a 6-pulse thyristor bridge rectifier is shown. Currently, thyristor-based rectifiers are particularly employed for industrial and power-to-gas applications where a high voltage and current are required to supply electrolyzers. In [8][9][10], 6-pulse and 12-pulse thyristors-based rectifiers connected to alkaline electrolyzers have been investigated In low-power wind turbines with a PMSG, the three-phase high-frequency semicontrolled rectifier, shown in Figure 7, is attractive, since it is simple and robust. Moreover, all active switches are connected to a common point, and a short-circuit through a leg is not possible. A higher distortion of the generator currents results [38]. turbines where the electrical generator is a squirrel cage induction generator (SCIG), a double-fed induction generator (DFIG) is adopted in this scheme in the back-to-back configuration for grid connection [35][36][37]. In low-power wind turbines with a PMSG, the three-phase high-frequency semicontrolled rectifier, shown in Figure 7, is attractive, since it is simple and robust. Moreover, all active switches are connected to a common point, and a short-circuit through a leg is not possible. A higher distortion of the generator currents results [38]. A three-level neutral-point diode clamped front-end topology (3L-NPC) gives two DC outputs that can be connected to two electrolyzers; each DC level has one half of the voltage obtained by the two levels of the active rectifier. It is sketched in Figure 8. The 3L-NPC topology is one of the most commercialized multilevel converters on the market, and it is usually proposed as a back-to-back topology in wind turbines. A potential drawback is the midpoint voltage fluctuation of the DC bus; however, this problem is minimized by the controlling of the redundant switching status [35,36,39].
Third, like hydrogen production systems supplied by solar energy, the DC microgrid system is shown in Figure 9 (the input filter is not drawn for the sake of clarity). The input of the buck converter in this system interfaces with the DC bus in DC microgrid and the output coupling with PEMEL [13]. In this case, the dynamic behavior depends on the need to use the hydrogen as storage depending on the constraints of the microgrids.
Finally, the hydrogen production system connected with the utility grid is present in Figure 10. This system consists of a delta-star transformer and an AC-DC rectifier interfacing the grid system and PEMEL. Since the current of the electrolyzer must be controlled, thyristors-based rectifiers are generally preferred for this purpose [8][9][10]. In Figure 10, a 6-pulse thyristor bridge rectifier is shown. Currently, thyristor-based rectifiers are particularly employed for industrial and power-to-gas applications where a high voltage and current are required to supply electrolyzers. In [8][9][10], 6-pulse and 12-pulse thyristors-based rectifiers connected to alkaline electrolyzers have been investigated A three-level neutral-point diode clamped front-end topology (3L-NPC) gives two DC outputs that can be connected to two electrolyzers; each DC level has one half of the voltage obtained by the two levels of the active rectifier. It is sketched in Figure 8. The 3L-NPC topology is one of the most commercialized multilevel converters on the market, and it is usually proposed as a back-to-back topology in wind turbines. A potential drawback is the midpoint voltage fluctuation of the DC bus; however, this problem is minimized by the controlling of the redundant switching status [35,36,39].
Third, like hydrogen production systems supplied by solar energy, the DC microgrid system is shown in Figure 9 (the input filter is not drawn for the sake of clarity). The input of the buck converter in this system interfaces with the DC bus in DC microgrid and the output coupling with PEMEL [13]. In this case, the dynamic behavior depends on the need to use the hydrogen as storage depending on the constraints of the microgrids.
Finally, the hydrogen production system connected with the utility grid is present in Figure 10. This system consists of a delta-star transformer and an AC-DC rectifier interfacing the grid system and PEMEL. Since the current of the electrolyzer must be controlled, thyristors-based rectifiers are generally preferred for this purpose [8][9][10]. In Figure 10, a 6-pulse thyristor bridge rectifier is shown. Currently, thyristor-based rectifiers are particularly employed for industrial and power-to-gas applications where a high voltage and current are required to supply electrolyzers. In [8][9][10], 6-pulse and 12-pulse thyristors-based rectifiers connected to alkaline electrolyzers have been investigated from the point of view of power quality, gas quality, and output current ripple. It has been demonstrated that the use of a 12-pulse thyristor-based rectifier enables enhancing the power quality and the gas quality while reducing the output current ripple.
C 2020, 6, x FOR PEER REVIEW 8 of 21 from the point of view of power quality, gas quality, and output current ripple. It has been demonstrated that the use of a 12-pulse thyristor-based rectifier enables enhancing the power quality and the gas quality while reducing the output current ripple.   The four types of hydrogen production need a suitable controller to control the PEMEL current or voltage to manage its hydrogen flow rate and/or its energy efficiency. Because of the intermittent behavior of solar and wind energy systems, the input voltage of the DC-DC converter may change quickly [23]. As a result, without a suitable and dedicated controller, output hydrogen production and energy efficiency as well can be impacted.
Different control techniques can be applied to control the DC-DC converter for providing a constant current or voltage to PEMEL such as the maximum power point tracking (MPPT) to  from the point of view of power quality, gas quality, and output current ripple. It has been demonstrated that the use of a 12-pulse thyristor-based rectifier enables enhancing the power quality and the gas quality while reducing the output current ripple.   The four types of hydrogen production need a suitable controller to control the PEMEL current or voltage to manage its hydrogen flow rate and/or its energy efficiency. Because of the intermittent behavior of solar and wind energy systems, the input voltage of the DC-DC converter may change quickly [23]. As a result, without a suitable and dedicated controller, output hydrogen production and energy efficiency as well can be impacted.
Different control techniques can be applied to control the DC-DC converter for providing a constant current or voltage to PEMEL such as the maximum power point tracking (MPPT) to  from the point of view of power quality, gas quality, and output current ripple. It has been demonstrated that the use of a 12-pulse thyristor-based rectifier enables enhancing the power quality and the gas quality while reducing the output current ripple.   The four types of hydrogen production need a suitable controller to control the PEMEL current or voltage to manage its hydrogen flow rate and/or its energy efficiency. Because of the intermittent behavior of solar and wind energy systems, the input voltage of the DC-DC converter may change quickly [23]. As a result, without a suitable and dedicated controller, output hydrogen production and energy efficiency as well can be impacted.
Different control techniques can be applied to control the DC-DC converter for providing a constant current or voltage to PEMEL such as the maximum power point tracking (MPPT) to The four types of hydrogen production need a suitable controller to control the PEMEL current or voltage to manage its hydrogen flow rate and/or its energy efficiency. Because of the intermittent behavior of solar and wind energy systems, the input voltage of the DC-DC converter may change quickly [23]. As a result, without a suitable and dedicated controller, output hydrogen production and energy efficiency as well can be impacted. Different control techniques can be applied to control the DC-DC converter for providing a constant current or voltage to PEMEL such as the maximum power point tracking (MPPT) to overcome the intermittent issues in these systems. Besides, the modeling of the PEMEL is important for control purposes, since it might have an impact on the performance of the system [25].

Proton Exchange Membrane Electrolyzer Modeling for Control Purpose
The proton exchange membrane electrolyzer modeling is mandatory for control purposes. Indeed, to develop efficient controllers, the PEMEL model has to be coupled with the model of the DC-DC converter. The objective is to get the overall model (i.e., DC-DC converter and the load) under the form of a transfer function or a state-space model. The developed controllers can be used to supply the PEMEL with a constant input energy (to optimize hydrogen flow rate and energy efficiency) or variable input energy to manage the hydrogen flow rate based on requirements (e.g., state-of-charge of hydrogen tanks) [25]. Based on the current literature, three types of modeling have been reported and are investigated in the following section.

Resistive Model
First of all, the PEMEL can be modeled as a resistor as shown in Figure 11. The resistor is the simplest model to represent the electric power transferred to the PEMEL. C 2020, 6, x FOR PEER REVIEW 9 of 21 overcome the intermittent issues in these systems. Besides, the modeling of the PEMEL is important for control purposes, since it might have an impact on the performance of the system [25].

Proton Exchange Membrane Electrolyzer Modeling for Control Purpose
The proton exchange membrane electrolyzer modeling is mandatory for control purposes. Indeed, to develop efficient controllers, the PEMEL model has to be coupled with the model of the DC-DC converter. The objective is to get the overall model (i.e., DC-DC converter and the load) under the form of a transfer function or a state-space model. The developed controllers can be used to supply the PEMEL with a constant input energy (to optimize hydrogen flow rate and energy efficiency) or variable input energy to manage the hydrogen flow rate based on requirements (e.g., state-of-charge of hydrogen tanks) [25]. Based on the current literature, three types of modeling have been reported and are investigated in the following section.

Resistive Model
First of all, the PEMEL can be modeled as a resistor as shown in Figure 11. The resistor is the simplest model to represent the electric power transferred to the PEMEL. This simple model has been used to study the power quality of an alkaline electrolyzer [10] and to develop controllers applied to a synchronous [14] and classic buck converter [17], full-bridge multi resonant converter [18], soft-switching full-bridge converter [15,16], and three-level interleaved buck converter [23]. However, only in [10,23], details are provided to determine the value of the resistive load. In [10], the authors have expressed the equivalent resistance of the electrolyzer based on the static voltage-current characteristic: The static voltage-current curve of an alkaline electrolyzer of 3 MW is shown in Figure 12. Based on Equation (4) and Figure 12, the equivalent resistive load according to the current is provided in Figure 13. It has to be noted that the equivalent resistive load decreases according to the current since electrolyzers are high current/low voltage loads as shown in Figure 9. To emphasize better the resistance decrease for a current range included between 2000 and 20,000 A, an additional curve has been inserted in Figure 13. This simple model has been used to study the power quality of an alkaline electrolyzer [10] and to develop controllers applied to a synchronous [14] and classic buck converter [17], full-bridge multi resonant converter [18], soft-switching full-bridge converter [15,16], and three-level interleaved buck converter [23]. However, only in [10,23], details are provided to determine the value of the resistive load. In [10], the authors have expressed the equivalent resistance of the electrolyzer based on the static voltage-current characteristic: The static voltage-current curve of an alkaline electrolyzer of 3 MW is shown in Figure 12. Based on Equation (4) and Figure 12, the equivalent resistive load according to the current is provided in Figure 13. It has to be noted that the equivalent resistive load decreases according to the current since electrolyzers are high current/low voltage loads as shown in Figure 9. To emphasize better the resistance decrease for a current range included between 2000 and 20,000 A, an additional curve has been inserted in Figure 13.
By comparison, in [23], the authors have expressed the value of the resistive load based on a dynamic electrical circuit model (presented in this section). The value of the resistance is determined as follows: where -R tot : sum of the resistances in the equivalent model, taking into account activation losses both at the anode and cathode and ohmic losses. The total resistance value is equal to 0.441 Ω. Based on Equation (5), the resistance value has been computed according to the PEMEL voltage range included between 5 and 8 V (rated stack power of 400 W), as shown in Figure 14. The higher the PEMEL voltage, the lower the resistive load. These results are similar to those introduced in [10].  By comparison, in [23], the authors have expressed the value of the resistive load based on a dynamic electrical circuit model (presented in this section). The value of the resistance is determined as follows: where -: sum of the resistances in the equivalent model, taking into account activation losses both at the anode and cathode and ohmic losses. The total resistance value is equal to 0.441 Ω.
-: reversible voltage of the PEM EL, which is equal to 4.38 V. -: PEMEL voltage (V), depending on the PEMEL current if the latter is controlled instead of the voltage.
Based on Equation (5), the resistance value has been computed according to the PEMEL voltage range included between 5 and 8 V (rated stack power of 400 W), as shown in Figure 14. The higher the PEMEL voltage, the lower the resistive load. These results are similar to those introduced in [10].   By comparison, in [23], the authors have expressed the value of the resistive load based on a dynamic electrical circuit model (presented in this section). The value of the resistance is determined as follows: where -: sum of the resistances in the equivalent model, taking into account activation losses both at the anode and cathode and ohmic losses. The total resistance value is equal to 0.441 Ω.
-: reversible voltage of the PEM EL, which is equal to 4.38 V. -: PEMEL voltage (V), depending on the PEMEL current if the latter is controlled instead of the voltage.
Based on Equation (5), the resistance value has been computed according to the PEMEL voltage range included between 5 and 8 V (rated stack power of 400 W), as shown in Figure 14. The higher the PEMEL voltage, the lower the resistive load. These results are similar to those introduced in [10].  To summarize, based on these two previous works, the calculation of the resistance is based on the static behavior of the PEMEL.  To summarize, based on these two previous works, the calculation of the resistance is based on the static behavior of the PEMEL.

Static Model
The static model of PEMELs is based on its current-voltage characteristic [19][20][21][22]. The parameters of the I-V characteristic curve are determined by using the measurements of a PEMEL single cell.
An example of an I-V characteristic curve for a PEMEL single cell is shown in Figure 15. This curve has been obtained for the following operating conditions: gas pressures of 1 bar and an ambient temperature of 20 • C. The specifications of the PEMEL used for the experiments are provided in Table 1. C 2020, 6, x FOR PEER REVIEW 12 of 21

Static Model
The static model of PEMELs is based on its current-voltage characteristic [19][20][21][22]. The parameters of the I-V characteristic curve are determined by using the measurements of a PEMEL single cell.
An example of an I-V characteristic curve for a PEMEL single cell is shown in Figure 15. This curve has been obtained for the following operating conditions: gas pressures of 1 bar and an ambient temperature of 20 °C. The specifications of the PEMEL used for the experiments are provided in Table  1.   Figure 15, the PEMEL has a nonlinear behavior. It can be modeled by using either a linear or a nonlinear model. First of all, a simple static modeling approach can be developed based on the following expression [20]:    Figure 15, the PEMEL has a nonlinear behavior. It can be modeled by using either a linear or a nonlinear model. First of all, a simple static modeling approach can be developed based on the following expression [20]:  The reversible voltage V int (electrochemical) is important for electrolysis process and hydrogen production, and can be calculated as follows [3]: where: -∆G is the Gibbs energy (238 kJ mol −1 for T = 20 • C) allowing splitting the deionized water into hydrogen and oxygen.
If the water is liquid, ∆G can be determined according to the temperature ( • C): -z is the number of electrons exchanged during the reaction. For H 2 , z = 2.
-F is the Faraday's constant (96 485 C mol −1 ). From Equation (6), an equivalent static electrical circuit for a single PEMEL cell can be deduced as shown in Figure 16. It is composed of an electromotive force representing the reversible voltage (V int ) in series-connected with an equivalent resistance (R 1 ) [21].
The reversible voltage Vint (electrochemical) is important for electrolysis process and hydrogen production, and can be calculated as follows [3]: where: -∆ is the Gibbs energy (238 kJ.mol −1 for T = 20 °C) allowing splitting the deionized water into hydrogen and oxygen.
If the water is liquid, ∆ can be determined according to the temperature (°C): -z is the number of electrons exchanged during the reaction. For H2, z = 2.
-F is the Faraday's constant (96 485 C. mol ). From Equation (6), an equivalent static electrical circuit for a single PEMEL cell can be deduced as shown in Figure 16. It is composed of an electromotive force representing the reversible voltage (Vint) in series-connected with an equivalent resistance (R1) [21].
The electrical power Pcell for a single cell is given by: Hence, the total electrical power (Pel) is obtained by multiplying the power (Pcell) by the number of cells (N): From Figure 16, the energy efficiency of the PEMEL is expressed as the ratio between the electrochemical hydrogen energy (PH2) and electrical power (Pel) as follows: From Equation (11), the higher the cell voltage, the lower the energy efficiency. Besides, the cell voltage efficiency is defined as the ratio between the thermoneutral cell voltage VTN (given in Equation (13)) and the cell voltage, vcell [3]: where -∆ is the change of enthalpy (285.84 kJ. mol −1 ). The electrical power Pcell for a single cell is given by: Hence, the total electrical power (P el ) is obtained by multiplying the power (P cell ) by the number of cells (N): From Figure 16, the energy efficiency of the PEMEL is expressed as the ratio between the electrochemical hydrogen energy (P H2 ) and electrical power (P el ) as follows: From Equation (11), the higher the cell voltage, the lower the energy efficiency. Besides, the cell voltage efficiency is defined as the ratio between the thermoneutral cell voltage V TN (given in Equation (13)) and the cell voltage, v cell [3]: C 2020, 6, 29 13 of 20 where -∆H is the change of enthalpy (285.84 kJ mol −1 ). Finally, the static model of the PEMEL can be modeled by an empirical current-voltage expression. The different parameter values can be found also by using a least-squares regression algorithm. The I-V relationship of the PEMEL is written as follows [41]: where -r is the ohmic resistance of electrolyte (Ω).
-s, t are the coefficients for overvoltage on electrodes It is well known that the overvoltages appearing in Equation (14) are dependent on the temperature. Therefore, the previous equation can be modified to take into consideration the effects of the temperature on the electrode and electrolyte overvoltage. The temperature-dependent I-V model can be expressed as follows [41]: -r i is the parameters for ohmic resistance of electrolyte (i = 1...2).
-s i , t i are the parameters for overvoltage on electrodes (i = 1...3).
-T is the temperature of electrolyte ( • C). The parameters r i , s i , t i are provided in [41].

Dynamic Model
The previously presented models are static and do not take into consideration the dynamics of the PEMEL in case of operating conditions change. These dynamics can be investigated in performing experimental tests by supplying the PEMEL with dynamic current profiles [24]. The realized experimental test bench to analyze the dynamic behavior of PEMELs is shown in Figure 17.  The DC power supply is controlled by using a virtual control panel with a laptop. The specifications of the used PEMEL are reported in Table 1. A test has been carried out with a DC current (from 0 to C 2020, 6, 29 14 of 20 10 A) to supply the PEMEL. The PEMEL stack response is provided in Figure 18. A high time scale for this test has been chosen to see the dynamics of the PEMEL and the steady-state operation.  The equivalent electrical circuit used to model the dynamic behavior of the PEMEL is shown in Figure 19. It consists of an electromotive force representing the reversible voltage (Vint), seriesconnected with a resistor (Rint) modeling the membrane, and an RC branch to emulate the losses in the anode and the dynamics [24]. In the literature, this dynamic model has been reported in other related works [9,42,43]. In [42], the authors have proposed this model based on previous works reported for fuel cells. Experimental results are provided to show the dynamics of a single-cell PEM electrolyzer, but the authors did not determine the parameters of the equivalent circuit. By comparison, in [43], the same model is presented based on [42], but only the static model is considered to be coupled with an interleaved DC-DC buck converter for control purposes. Finally, in [9], this model is presented for an alkaline electrolyzer to investigate the effects of rectifiers on the specific energy consumption and gas quality during dynamic operation. From Figure 18, the PEMEL stack voltage is equal to the open-circuit voltage called Nernst voltage or reversible voltage (V int ) at i el = 0 A. This voltage is approximately equal to 4.33 V. When the DC current is supplied to the PEMEL, an immediate rise in PEMEL stack voltage can be noticed, followed by a slow final rise before reaching its steady-state value. The immediate rise in stack voltage can be modeled as a resistor, representing the ohmic losses in the membrane. Besides, the slow final rise in stack voltage represents the speed of the chemical reactions both at the anode and the cathode. In this test, this final rise lasts 13 s and can change according to the input energy supplying the PEMEL. Indeed, as highlighted in [24], the duration of the final rise to reach the steady-state operation due to the move of the electrons into the anode and cathode lengthens especially as the input energy increases. However, when the stack voltage gets closer to the limit operating voltage (for this PEMEL, the limit voltage is equal to 8 V), the dynamics are faster, and the stack voltage in steady-state operation may drop [24]. To summarize, this slow rise represents the activation overvoltage both at the anode and cathode. Since this dynamic behavior is close to a first-order system, it can be modeled as RC branches. On one hand, the resistors result in activation losses both at the anode and cathode, while the capacitors represent the charge separation into the anode and cathode. Hence, the dynamic behavior of the PEMEL to the final duration required by the charge layers to move when a DC current is supplied can be replicated.
On the other hand, in [24], it has been emphasized that the reaction speed at the anode is much slower than the reaction speed at the cathode. For this reason, the activation phenomena in PEMEL are mainly dominated by the anode reaction.
The equivalent electrical circuit used to model the dynamic behavior of the PEMEL is shown in Figure 19. It consists of an electromotive force representing the reversible voltage (V int ), series-connected with a resistor (R int ) modeling the membrane, and an RC branch to emulate the losses in the anode and the dynamics [24]. In the literature, this dynamic model has been reported in other related works [9,42,43]. In [42], the authors have proposed this model based on previous works reported for fuel cells. Experimental results are provided to show the dynamics of a single-cell PEM electrolyzer, but the authors did not determine the parameters of the equivalent circuit. By comparison, in [43], the same model is presented based on [42], but only the static model is considered to be coupled with an interleaved DC-DC buck converter for control purposes. Finally, in [9], this model is presented for an alkaline electrolyzer to investigate the effects of rectifiers on the specific energy consumption and gas quality during dynamic operation. From Figure 19, the PEMEL stack voltage can be expressed as follows: The dynamic activation overvoltage at the anode can be written as: The couple R1C1 represents the time constant in seconds of the dynamics. The equivalent resistor R1 at the anode can be determined based on the activation overvoltage vACT and the current flowing through the electrolyzer iel. Hence, the constant time at the anode is expressed as follows: By using the static model identification (see Section 3.2), the reversible voltage (Vint) and the membrane resistor (Rint) can be determined. After that, to identify the parameters of the RC branch for the anode reaction, a dynamic model identification has to be considered by investigating the transient operation as shown in Figure 18. A least-squares regression algorithm has been used to determine the time constant of the transient operation, then the parameters of the activation resistance and double-layer capacitor [24]. It has to be noted that the equivalent double-layer capacitor determined for this study is equal roughly to 37 F, while the values reported for PEM fuel cells and alkaline electrolyzers are of the order of a few Farad [44] and milli Farad [9], respectively. These results are high in terms of capacitance values and are assimilated to supercapacitors.

Comparison between Static and Dynamic Model
A comparison between the static model ( Figure 16) and the dynamic model ( Figure 19) providing the voltage response of the PEMEL according to a current profile is given in Figure 20 [24]. It can be noted that the voltage response obtained by the dynamic model allows replicating accurately the real dynamic behavior of a PEMEL. Based on previous work [24], a dynamic model features a maximum error of about 4%, whereas a static model presents a maximum error of about 15%. Unfortunately, the maximum error for the static model occurs just after the step solicitation; it is a typical situation depicted by an abrupt variation of the available power when a RES is used.
The equivalent model developed in Figure 19 has been used for the first time to design the controller of a suitable DC-DC converter to supply a PEMEL, namely a stacked interleaved buck converter (SIBC) as shown in Figure 21 [25]. Compared to a classic interleaved buck converter, the SIBC includes an additional capacitor Cs between the first and the second phase. The capacitor Cs From Figure 19, the PEMEL stack voltage can be expressed as follows: The dynamic activation overvoltage at the anode can be written as: The couple R 1 C 1 represents the time constant in seconds of the dynamics. The equivalent resistor R 1 at the anode can be determined based on the activation overvoltage v ACT and the current flowing through the electrolyzer i el . Hence, the constant time τ a at the anode is expressed as follows: By using the static model identification (see Section 3.2), the reversible voltage (V int ) and the membrane resistor (R int ) can be determined. After that, to identify the parameters of the RC branch for the anode reaction, a dynamic model identification has to be considered by investigating the transient operation as shown in Figure 18. A least-squares regression algorithm has been used to determine the time constant of the transient operation, then the parameters of the activation resistance and double-layer capacitor [24]. It has to be noted that the equivalent double-layer capacitor determined for this study is equal roughly to 37 F, while the values reported for PEM fuel cells and alkaline electrolyzers are of the order of a few Farad [44] and milli Farad [9], respectively. These results are high in terms of capacitance values and are assimilated to supercapacitors.

Comparison between Static and Dynamic Model
A comparison between the static model ( Figure 16) and the dynamic model ( Figure 19) providing the voltage response of the PEMEL according to a current profile is given in Figure 20 [24]. It can be noted that the voltage response obtained by the dynamic model allows replicating accurately the real dynamic behavior of a PEMEL. Based on previous work [24], a dynamic model features a maximum C 2020, 6, 29 16 of 20 error of about 4%, whereas a static model presents a maximum error of about 15%. Unfortunately, the maximum error for the static model occurs just after the step solicitation; it is a typical situation depicted by an abrupt variation of the available power when a RES is used.  Based on the static and dynamic models, the bode diagrams of the transfer function (to control the stack voltage of the PEMEL) of the SIBC (Figure 22), and accordingly, the step response of the system (Figure 23) changes by considering either a static or a dynamic load. Besides, it can be seen in Figure 24 that the step response of the SIBC with a dynamic model is slower due to the anode reaction. As a result, the design and the tuning of the controller require an accurate model, such as the dynamic model presented in Section 3.3. To control PEM electrolyzers efficiently, overshoot must be avoided that may damage the electrolyzer and stability must be ensured without oscillation as highlighted in [25]. Besides, since PEM electrolyzers present slow dynamics, the rapidity is not the most important criterium to design controllers.
The different step responses demonstrate that adopting a simplified model such as the static one could lead to a wrong design of the controller. Indeed, the bode diagram of the converter with the static model shows a gain higher compared to the dynamic model. This implies the response with overshoot requiring a correction that would worsen the dynamic. Contrarily, the real response of the system corresponds to the dynamic model in which no overshoot occurs. It means that the dynamic of the PEMEL enforces the stability of the system; for this reason, it is crucial to take it into account. Besides, this issue is strictly tied to the use of RES since for stationary grid-supplied applications no input power variations occur and the voltage or current supplying the electrolyzer can be varied according to the quantity of hydrogen to be produced.  The equivalent model developed in Figure 19 has been used for the first time to design the controller of a suitable DC-DC converter to supply a PEMEL, namely a stacked interleaved buck converter (SIBC) as shown in Figure 21 [25]. Compared to a classic interleaved buck converter, the SIBC includes an additional capacitor C s between the first and the second phase. The capacitor C s allows blocking the DC component of the current flowing through the second phase. As a result, only the AC component of the current flows through the second phase. Besides, since a couple of power switches (i.e., S 1 and S 4 , and S 2 and S 3 ) are controlled in an opposite way, only the DC current supplies the PEMEL. A low current ripple is one of the most important features required for the DC-DC converter to supply the electrolyzer. Indeed, as demonstrated in the literature, the current ripple leads up to the decrease of the energy efficiency of the electrolyzer [45,46]. converter to supply the electrolyzer. Indeed, as demonstrated in the literature, the current ripple leads up to the decrease of the energy efficiency of the electrolyzer [45,46]. Based on the static and dynamic models, the bode diagrams of the transfer function (to control the stack voltage of the PEMEL) of the SIBC (Figure 22), and accordingly, the step response of the system ( Figure 23) changes by considering either a static or a dynamic load. Besides, it can be seen in Figure 24 that the step response of the SIBC with a dynamic model is slower due to the anode reaction. As a result, the design and the tuning of the controller require an accurate model, such as the dynamic model presented in Section 3.3. To control PEM electrolyzers efficiently, overshoot must be avoided that may damage the electrolyzer and stability must be ensured without oscillation as highlighted in [25]. Besides, since PEM electrolyzers present slow dynamics, the rapidity is not the most important criterium to design controllers.
The different step responses demonstrate that adopting a simplified model such as the static one could lead to a wrong design of the controller. Indeed, the bode diagram of the converter with the static model shows a gain higher compared to the dynamic model. This implies the response with overshoot requiring a correction that would worsen the dynamic. Contrarily, the real response of the system corresponds to the dynamic model in which no overshoot occurs. It means that the dynamic of the PEMEL enforces the stability of the system; for this reason, it is crucial to take it into account. Besides, this issue is strictly tied to the use of RES since for stationary grid-supplied applications no input power variations occur and the voltage or current supplying the electrolyzer can be varied according to the quantity of hydrogen to be produced.  Based on the static and dynamic models, the bode diagrams of the transfer function (to control the stack voltage of the PEMEL) of the SIBC (Figure 22), and accordingly, the step response of the system (Figure 23) changes by considering either a static or a dynamic load. Besides, it can be seen in Figure 24 that the step response of the SIBC with a dynamic model is slower due to the anode reaction. As a result, the design and the tuning of the controller require an accurate model, such as the dynamic model presented in Section 3.3. To control PEM electrolyzers efficiently, overshoot must be avoided that may damage the electrolyzer and stability must be ensured without oscillation as highlighted in [25]. Besides, since PEM electrolyzers present slow dynamics, the rapidity is not the most important criterium to design controllers.
The different step responses demonstrate that adopting a simplified model such as the static one could lead to a wrong design of the controller. Indeed, the bode diagram of the converter with the static model shows a gain higher compared to the dynamic model. This implies the response with overshoot requiring a correction that would worsen the dynamic. Contrarily, the real response of the system corresponds to the dynamic model in which no overshoot occurs. It means that the dynamic of the PEMEL enforces the stability of the system; for this reason, it is crucial to take it into account. Besides, this issue is strictly tied to the use of RES since for stationary grid-supplied applications no input power variations occur and the voltage or current supplying the electrolyzer can be varied according to the quantity of hydrogen to be produced.

Figure 23.
Step response of the SIBC either by considering a static or dynamic model. Figure 23. Step response of the SIBC either by considering a static or dynamic model.

Figure 23.
Step response of the SIBC either by considering a static or dynamic model. Figure 24. Slow response of the system with the dynamic model. Figure 24. Slow response of the system with the dynamic model.

Conclusions
The goal of this article was to underline the importance of using power converters, particularly DC-DC converters, to interface a high DC bus voltage to a low DC voltage required for PEM electrolyzers. With renewable energy sources as supply, green hydrogen can be produced and blended to natural gas or employed in other clean applications. Besides, the modeling of the PEM electrolyzer is mandatory to develop controllers such as energy efficiency or hydrogen flow rate-based controls. In this article, three types of modeling have been investigated (i.e., resistive, static, and dynamic load) showing that the use of RES implies different design constraints compared to stationary grid-supplied converters for electrolyzers. It has been demonstrated that the dynamics play an important role to understand the behavior of the PEM electrolyzer in case of dynamic operating conditions. Indeed, the PEM electrolyzers can be coupled with renewable energy sources and are submitted to dynamic operations due to the weather conditions change. As a result, if these dynamics were taken into account in designing the controller, the performance of the system could be enhanced greatly making hydrogen more attractive also in carbon-based fuels.