- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in Power Electronics

Volume 2013 (2013), Article ID 352765, 9 pages

http://dx.doi.org/10.1155/2013/352765

## Modeling, Control, and Simulation of a Solar Hydrogen/Fuel Cell Hybrid Energy System for Grid-Connected Applications

Power Electronic Group (PEG), ENIS, BP W, 3038 Sfax, Tunisia

Received 10 January 2013; Revised 18 February 2013; Accepted 18 February 2013

Academic Editor: Jose Pomilio

Copyright © 2013 Tourkia Lajnef et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Different energy sources and converters need to be integrated with each other for extended usage of alternative energy, in order to meet sustained load demands during various weather conditions. The objective of this paper is to associate photovoltaic generators, fuel cells, and electrolysers. Here, to sustain the power demand and solve the energy storage problem, electrical energy can be stored in the form of hydrogen. By using an electrolyser, hydrogen can be generated and stored for future use. The hydrogen produced by the electrolyser using PV power is used in the FC system and acts as an energy buffer. Thus, the effects of reduction and even the absence of the available power from the PV system can be easily tackled. Modeling and simulations are performed using MATLAB/Simulink and SimPowerSystems packages and results are presented to verify the effectiveness of the proposed system.

#### 1. Introduction

At present, most of energy demand in the world relies on fossil fuels such as petroleum, coal, and natural gas that are being exhausted very fast. One of the major severe problems of global warming is one of these fuels combustion products, carbon dioxide; these are resulting in great danger for life on our planet [1].

Fossil fuels can have as an alternative some renewable energy sources like solar, wind, biomass, and so; among them on the photovoltaic (PV) generator which converts the solar radiation into electricity, largely used in low power applications. The photovoltaic generator is chosen for its positive points including being carbon free and inexhaustible; moreover, it does not cause noise for it is without moving parts and with size-independent electric conversion efficiency [2].

Nevertheless, the power generated by a PV system is influenced by weather conditions; for example, at night or in cloudy periods, it would not generate any power or application. In addition, it is difficult to store the power generated by a PV system for future use. The best method to overcome this problem is to integrate the PV generator with other power sources such as an electrolyser, hydrogen storage tank, FC system, or battery due to their good features such as high efficiency response, modular production, and fuel flexibility [3, 4]. Its coordination with a PV system could be successful for both grid-connected and stand-alone power applications. Thanks to the rapid response capability of the fuel cell power system, the photovoltaic fuel cell hybrid system can be able to overcome the inconvenience of the intermittent power generation. Furthermore, unlike a secondary battery, the FC does not only store energy but also produce electricity for unlimited time to support the PV power generator. Hence, the coordination between the FC power system and the photovoltaic generator becomes necessary in order to smooth out the PV power fluctuations.

This paper focuses on developing a simulation model to design and size the hybrid system for a variety of loading and meteorological conditions. This simulation model is performed using Matlab and SimPowerSystems and results are presented to verify the effectiveness of the proposed system.

#### 2. Modeling

##### 2.1. A Dynamic Model of PV Generator

PV arrays are built up with combined series/parallel combinations of PV solar cells, which allow extracting the characteristic parameters of the one-diode equivalent model for a single solar cell.

The PV cell output voltage is a function of the photocurrent that is mainly determined by load current depending on the solar irradiation level during the operation [2]: The parameters used in the mathematical model of the PV system are shown in the abbreviations section.

##### 2.2. A Dynamic Model of PEMFC

The PEMFC is an electrochemical device which allows the electric energy conversion of the chemical energy contained in a reaction between a fuel, the hydrogen, and an oxidizer, the oxygen. The temperature effects have been taken into account in the typical range of low temperature PEM (60–100°C) and a thermal behaviour submodel has been introduced.

A bias voltage is applied across the electrochemical cell in order to induce electrochemical reactions at both electrodes. Water is introduced at the anode and dissociated into oxygen, protons, and electrons. The protons are driven by an electric field through the PEM to the cathode where they combine with the electrons arriving from the external circuit to form hydrogen gas. The assessment of these two half reactions produces water, heat, and electricity as Figure 1 shows.

Many models have been proposed to simulate the fuel cell in the literature [4, 5]. This model is built by utilizing the relationship between the output voltage and potential pressure of hydrogen, oxygen, and water. Figure 2 shows the detailed PEMFC model, which is then embedded into the SimPowerSystems of MATLAB-controlled voltage source and integrated into the overall system.

The proportional relationship of the molar flow of gas through a valve with its partial pressure can be expressed as

For hydrogen molar flow, the derivative of the partial pressure can be calculated using the perfect gas equation as follows [5]
The relationship between the hydrogen flow and the FC system current can be written as [5]
Using (4) and (3), the hydrogen partial pressure can be rewritten in the (*s*) domain as [5]
where
Typically, the fuel cell output voltage is obtained from the sum of three effects [5, 6], the Nernst potential, the cathode and anode activation overvoltage, and the ohmic overvoltage due to internal resistance:
where
Now, the Nernst voltage in terms of gas molarities may be expressed as [5]

The FC system consumes hydrogen obtained from the on-board high-pressure hydrogen tanks according to the power demand. A feedback control strategy is used to control hydrogen flow rate according to the output power of the FC system. This feedback control is achieved where FC current how the output is taken back into the input while converting the hydrogen into molar form [5, 7].

The quality of hydrogen found from the hydrogen tank is given by:

##### 2.3. An Electrolyser Model

Electrolysis of water is the dissociation of water molecules into hydrogen and oxygen gas. The electrochemical reaction of water electrolysis is given by

The rate of hydrogen reacting is directly proportional to the electrical current in the equivalent electrolysis circuit [8, 9], given by

The relation between the real hydrogen flow rate and the theoretical one is defined as the Faraday’s efficiency. In general, it is assumed to be more than 99%. The Faraday efficiency is expressed by [9, 10]

Figure 3 shows the simulation model of the electrolyser implemented in MATLAB/Simulink.

##### 2.4. A Hydrogen Storage System

The amount of hydrogen required by the PEMFC is sent directly from the electrolyzer system based on the relationship between the output power and the hydrogen requirement of the PEMFC system. The remaining amount of hydrogen is sent to the storage tank [11]. The parameters used in the hydrogen storage system are listed in the abbreviations section.

In this study, the dynamic of the storage is obtained as follows [10, 11]:

Neither the compression dynamics nor the compression energy requirements are accounted for in our calculations. All auxiliary power requirements such as pumps, valves were ignored in the dynamic model. The Simulink version of the hydrogen storage model is depicted in Figure 4.

##### 2.5. An Averaged Model of the Converters

The simulation of power electronic systems behavior using semiconductor refined models gives an accurate, but unaffordable results. The structure of the switch cell is given in Figure 5; it has two basic switching cells (P-cell and N-cell).

Each cell consists of one controlled switch and one diode. The two active switches are directly controlled by external control signals. The diodes are indirectly controlled by the state of the controlled switches and the circuit conditions. The load is presented by an inductor and a voltage source (). The DC loop inductance is modeled by an inductor . Depending on the sign of the load current , only two devices (one controlled device and one diode) operate simultaneously [12–15].

The proposed averaged model of the PWM-switch is presented in Figure 5(b). This model contains a controlled voltage source () and a controlled current source (). The PWM-switch is the only nonlinear element which is supposed to be responsible for the nonlinear behavior of the converter, Considering as the switching period of the controlled switches and () the duty ratio which is the ratio of the on-time value () of the upper controlled switch () and the switching period .

In Figure 5(b), the current source () and the voltage source () are given by where and are the time averaged values of the instantaneous terminal waveforms and , respectively, over one cycle .

In this paper we have chosen MOSFET’s devices for the controlled switch and the PIN devices for diodes. We notice that the load current is considered constant and equal to the averaged value of the real current in the load over the switching period . During devices turn-on and turn-off phases, the different voltage magnitudes are given by

The averaged values of the voltage across the device of the voltage source and current are obtained by integrating the voltage and the current evolutions over one cycle .

##### 2.6. The PV System Model and Integration to the Overall Hybrid Model

As shown in Figure 6, the majority of the PVFC system comprises a solar-cell model, a PEM fuel cell generator, and a water electrolyser. The parameters of the PV systems are listed in Table 1. The PEMFC system parameters are given in Table 2.

DC-DC converters are widely used in PV generating systems as an interface between the PV generator and the water electrolyser. An optimized controller for the DC-DC must find the optimum duty cycle which leads the PV generator as close as possible to its MPPT and ensure that the working point at the water electrolyser is a safety point. Since the fuel cell stack operates at a low DC voltage range (102 V in this paper), the DC-DC converter must boost the DC voltage and invert it to the AC grid frequency (230 V/50 Hz here) for grid-connected operations. To keep the DC buses fixed at 400 V, we chose to use a hysteresis regulator.

Hence, the fuel cell control problem is translated into an output current control requirement, to be realized by the DC/DC converter, in order to ensure optimal operation for a given fuel flow rate. Under these conditions, the cell output power is directly related to its fuel consumption at the selected optimum operating point of the *V*-*I* characteristic.

Operating the fuel cell at different output power levels requires suitable variation of its input flow rate, to be realized by the overall control system of the cell. The power demand requirement of the fuel cell is translated into a current demand input by dividing with stack output voltage:

##### 2.7. Control Strategy of the Hybrid Power System

The proposed control strategy is based on a power decoupling strategy in the frequency domain of the power source. This energy-management strategy fulfills the fast energy demands if the load respects the integrity of each source. The inner control loop is tuned to drive the PEMFC current. can be considered constant during current regulation. With this assumption, a classical proportional-integral (PI) controller is proposed (Figure 7).

The implementation structure of the loop regulation with hysteresis is given in Figure 8 scheme. This cascaded control structure shows that it could satisfy the different requirements of the hybrid system both from the point of view of the feed load and from the point of view of the component limits.

#### 3. Results and Discussion

In this section we present simulation results for the coupling between the PV/FC and the PEM electrolyser through the DC-DC converter controller. The PV generator is a 2880 W power plant at 1000 (w/m^{2}) solar radiation. From Figures 9 and 10, we observe that the current and the available power from the PV generator decreases because of radiation variation.

Initially, the total power generated by the PV generator is sent to the electrolyser through a DC-DC converter to generate hydrogen. The hydrogen produced by the electrolyser causes the pressure of the storage tank to vary as shown in Figure 11. These results are obtained for two values of solar radiation 1000 (w/m^{2}) and 800 (w/m^{2}). We notice that the final value of pressure is used as initial condition in the following simulation.

After this, the total power generated by the PV system is sent to grid via a DC/AC converter. At 1000 (w/m^{2}) solar radiation, the PV generator used in this study is capable of delivering 2880 W while the FC delivers 200 W, but at 800 (w/m^{2}) solar radiation (at 1 s) the power produced by the PV system tends to 2300 W. So the generated power by the PV is less than the demand; power will be supplied from the FC system. The power produced by the FC system tends to 500 W at 1 s of time simulation.

The internal voltage of the FC system decreases when the FC output power increases. This relation between power and the voltage of the FC system authenticates the reliability of the FC model. The transient response of the FC system voltage to the load changes varies according to the amount of power supplied by the FC system as shown in Figure 12. The power produced by the FC system is given in Figure 13.

The amount of hydrogen moles consumed by the FC system is proportional to the power drawn from the FC system. The hydrogen flow to the FC system per second is depicted in Figure 14.

Figure 15 shows the hydrogen storage tank pressure variation corresponding to the amount of hydrogen extracted from the storage tank. It is evident that the hydrogen storage tank pressure decreases with time as more and more hydrogen is extracted from the storage tank. The pressure variation of storage hydrogen is as illustrated in Figure 16.

#### 4. Conclusion

In this paper, a PV/FC generator and PEM electrolyser have been described for a PV/FC system intended for grid-connected operations. Special attention has been paid to the modeling of temperature dependence, concentration over potential, and limiting current in the PEM electrolyser model. Then, the power conditioning system, including the DC/DC and DC/AC converters, is presented and typical waveforms are shown from its simulation in MATLAB/Simulink.

#### Abbreviations

: | Ideality or completion factor |

: | PV cell reverse saturation current [A] |

: | PV cell output current [A] |

: | Short-circuit cell current (representing insolation level [A]) |

: | Boltzmann’s constant [j/°K] |

: | Voltage factor |

: | The number of parallel strings |

: | The number of series cells per string |

: | Electron charge [] |

: | Series resistance of PV cell [] |

: | PV cell temperature [°K] |

: | PV cell voltage corresponding to maximum power [] |

: | Open-circuit voltage [] |

: | Terminal voltage for PV cell [V] |

: | Molar mass of hydrogen [kg kmol^{−1}] |

: | Hydrogen moles per second delivered to the storage tank [kmol/s] |

: | Pressure of tank [pascal] |

: | Initial pressure of the storage tank [pascal] |

: | Universal (Rydberg) gas constant [j/(kmol K)] |

: | Operating temperature [°K] |

: | Volume of the tank [m^{3}] |

: | Compressibility factor as a function of pressure. |

#### References

- K.-S. Ro and S. Rahman, “Two-loop controller for maximizing performance of a grid-connected photovoltaic-fuel cell hybrid power plant,”
*IEEE Transactions on Energy Conversion*, vol. 13, no. 3, pp. 276–281, 1998. View at Publisher · View at Google Scholar · View at Scopus - Y. Sukamongkol, S. Chungpaibulpatana, and W. Ongsakul, “A simulation model for predicting the performance of a solar photovoltaic system with alternating current loads,”
*Renewable Energy*, vol. 27, no. 2, pp. 237–258, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. J. Hwang, D. Y. Wang, N. C. Shih, D. Y. Lai, and C. K. Chen, “Development of fuel-cell-powered electric bicycle,”
*Journal of Power Sources*, vol. 133, no. 2, pp. 223–228, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. J. Hwang, W. R. Chang, F. B. Weng, A. Su, and C. K. Chen, “Development of a small vehicular PEM fuel cell system,”
*International Journal of Hydrogen Energy*, vol. 33, no. 14, pp. 3801–3807, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Y. El-Sharkh, A. Rahman, M. S. Alam, P. C. Byrne, A. A. Sakla, and T. Thomas, “A dynamic model for a stand-alone PEM fuel cell power plant for residential applications,”
*Journal of Power Sources*, vol. 138, no. 1-2, pp. 199–204, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. Padullés, G. W. Ault, and J. R. McDonald, “Integrated SOFC plant dynamic model for power systems simulation,”
*Journal of Power Sources*, vol. 86, no. 1, pp. 495–500, 2000. View at Publisher · View at Google Scholar · View at Scopus - J. Hamelin, K. Agbossou, A. Laperrière, F. Laurencelle, and T. K. Bose, “Dynamic behavior of a PEM fuel cell stack for stationary applications,”
*International Journal of Hydrogen Energy*, vol. 26, no. 6, pp. 625–629, 2001. View at Publisher · View at Google Scholar · View at Scopus - M. J. Khan and M. T. Iqbal, “Dynamic modeling and simulation of a small wind-fuel cell hybrid energy system,”
*Renewable Energy*, vol. 30, no. 3, pp. 421–439, 2005. View at Publisher · View at Google Scholar · View at Scopus - O. Ulleberg,
*Stand-alone power systems for the future: optimal design, operation and control of solar-hydrogen energy systems [Ph.D. thesis]*, Norwegian University of Science and Technology, 1998. - K. Sapru, N. T. Stetson, S. R. Ovshinsky et al., “Development of a small scale hydrogen production-storage system for hydrogen applications,” in
*Proceedings of the 32nd Intersociety Energy Conversion Engineering Conference*, pp. 1947–1952, August 1997. View at Scopus - H. Görgün, “Dynamic modelling of a proton exchange membrane (PEM) electrolyzer,”
*International Journal of Hydrogen Energy*, vol. 31, no. 1, pp. 29–38, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. J. N. Spruijt, D. M. O'Sullivan, and J. B. Klaassens, “PWM-switch modeling of DC-DC converters,”
*IEEE Transactions on Power Electronics*, vol. 10, no. 6, pp. 659–664, 1995. View at Publisher · View at Google Scholar - J. Chen, R. W. Erickson, and D. Maksimovic', “Averaged switch modeling of boundary conduction mode DC-to-DC converters,” in
*Proceedings of the 27th Annual Conference of the IEEE Industrial Electronics Society (IECON '01)*, vol. 2, pp. 844–849, Denver, Colo, USA, November 2001. View at Publisher · View at Google Scholar - A. Ammous, K. Ammous, M. Ayedi, Y. Ounajjar, and F. Sellami, “An advanced PWM-switch model including semiconductor device nonlinearities,”
*IEEE Transactions on Power Electronics*, vol. 18, no. 5, pp. 1230–1237, 2003. View at Publisher · View at Google Scholar · View at Scopus - S. Abis, K. Ammous, H. Morel, and A. Ammous, “Advanced averaged model of PWM-switch operating in CCM and DCM conduction modes,”
*International Review of Electrical Engineering*, vol. 2, no. 4, pp. 544–556, 2007.