# On the Performance of Hybrid PV/Unitized Regenerative Fuel Cell System in the Tropics

Authors: Salwan Dihrab; Tamer Khatib; Kamaruzzaman Sopian; Habeeb Al-Ani; Saleem H. Zaidi
Journal: International Journal of Photoenergy
 (2012)

---

## Abstract

Solar hydrogen system is a unique power system that can meet the power requirements for future energy demands. Such a system uses the hydrogen as the energy carrier, which produces energy through the electrolyzer with assistance of the power from the PV during the sunny hours, and then uses stored hydrogen to produce energy through the fuel cell after sunset or on cloudy days. The current study has used premanufactured unitized regenerative fuel cells in which the electrolyzer and the fuel cell function within one cell at different modes. The system components were modeled and the one-day real operational and simulated data has been presented and compared. The measured results showed the ability of the system to meet the proposed load, and the total efficiency was about 4.5%.


---

## Body

## 1. Introduction

Photovoltaic (PV) power generation, which directly converts solar radiation into electricity, contains a lot of significant advantages, such as being inexhaustible and pollution-free, silent and with no rotating parts, and its size-independent electricity conversion efficiency. A positive environmental effect of photovoltaic is replacing electricity generated in a more pollutant way, or providing electricity where none was available before. While increasing the penetration of solar photovoltaic devices, various antipollution apparatuses can be operated by solar PV power, for example, water purification by electrochemical processing or stopping desert expansion by photovoltaic water pumping with tree implantation [1].The hybrid PV/hydrogen systems have units producing hydrogen called electrolyzers, which are used to generate hydrogen with the assistance of the PV as an energy source. This hydrogen will be preserved in suitable storage medium before being converted to DC power through a unit called fuel cell. Fuel cells are electrochemical devices that convert chemical energy, typically from hydrogen, directly into electrical energy. The introduction of fuel cell systems into the power generation market will not only supply clean renewable energy to millions of users, but it will help to reduce the dependence on oil [2].Such a system includes a source of power (PV modules), a hydrogen generator (electrolyzer), a storage medium to store the hydrogen and hydrogen utilization units (fuel cell) which are able to make a continuous, and self-dependant or uninterruptable source of energy. The PV size should be designed properly, in order to have excess power over the load requirement. This excess power will be used to generate hydrogen through the generation units, and it will then utilize it after sunset or at a time of low solar irradiation. The load in this system will be powered in the short-term by the PV modules, and the hydrogen system will act as a long-term supplier. An illustration of short- and long-term supplying is shown in Figure1.Figure 1
Long- and short-term power supplying for solar hydrogen systems [3].The modeling of a PV-powered hydrogen generation system is needed to solve the electrical equations for the PV and the DC/DC converters (if they exist) and the electrochemical equation for the electrolyzers. Fischer [4] compared the experimental results of a small system consisting of a solar cell generator, water electrolysis, and power conditioning with calculated results of a system simulation model, while Vidueira et al. [5] examined the performance of the solar hydrogen system for producing hydrogen via an electrolysis generator, in order to satisfy the hydrogen consumption of the first two fuel cell buses in Madrid. Some other researchers used various software to simulate the system performance, such as Park et al. [6, 7], who used (PSCAD/EMTDC) software to simulate a hydrogen production system performance.Galli and Stefanoni [8] investigated and tested some commercial solar-hydrogen technologies, preliminary performance results, and the control of the system for solar-powered regenerative fuel cell system in Italy, while Hedström et al. [9] did so in Sweden, Abaoud and Steeb [10] in the German-Saudi bilateral program HYSOLAR, Voss et al. [11] in Germany, Lutfi and Veziroǧlu [12] in Pakistan, Chaparro et al. [13] in Spain, and Little et al. [14] in UK. Torres et al. [15] also simulated a PV-hydrogen-fuel-cell hybrid system in Mexico.The main objective of this paper is to present the performance of a hybrid PV/unitized regenerative fuel cell system under Malaysian weather conditions. A PV/unitized regenerative fuel cell system is built in the laboratories of the solar energy research institute (SERI)/Universiti Kebangsaan Malaysia. The results presented in this paper are based on a real operation conditions. The comparison between the operational and simulated data is also presented in this paper.


## 2. System Description

The proposed system mainly consists of two major subsystems, which are the PV modules (Kyosera, KC85T) located outdoor and the hydrogen system located inside the lab. The main components of the system used in this work are 3 same oriented PV modules (the method used by Markvart [16] is based on the maximum operating current, and the lowest solar irradiation was used to predict how many modules needed to be used) with 87 W maximum power and 0.65 m2 active area, unitized regenerative fuel cell (URFC), gases storages and water tanks, humidifiers, water traps, DC/DC converter (WD150-400), maximum power point tracker (MPPT) (STECA PR 1010) to ensure maximum PV power output, and bus bar to transfer the current either from the PV to the load and electrolyzer in the sunny hours or from the fuel cell to the load during the sunset. The system starts in the morning, when the PV current is high enough to power the load, while the excess current sends power to URFC stack, which works in electrolysis mode, receiving water from small water pump. The generated gases are then sent to the storage tanks. When the PV current fails to power the load, the system will then switch to fuel cell mode after purging the lines with nitrogen. The humidifiers used for fuel cell mode were PERMA PURE MH Series humidifiers which are shell and tube moisture exchangers that allow the transfer of water vapor between a liquid water supply and a flowing gas stream. Water is absorbed into the walls of the Nafion tube and transferred to the dry gas stream. This transfer is driven by the difference in partial pressures of water vapor on opposing sides. Deionized water was filled into the tank which is equipped with heater (2 kW) to heat up the water to the desired temperature through the temperature controller. Nitrogen for purging was used when switching from electrolyzer mode to fuel cell mode and at system shutdown time at which the nitrogen supplies to both hydrogen and oxygen piping line using solenoid valve. The solenoid valves will open and allow nitrogen to flow through the fuel and oxidant. A schematic of the hydrogen system is shown in Figure 2.Figure 2
The hydrogen system prototype.The laboratory was equipped with a fuel cell test station capable of measuring the stack temperature, hydrogen and air flow rate measurement, and the stack operating current and voltage. The subsystem contains humidification bottles for the fuel cell mode to humidify the reactant gases. The humidification can be bypassed if a dry gas feed into the stack is desired. The gas flow rates are controlled manually by a set of valves for each side of the stacks. Moreover, a data acquisition system was incorporated to generate polarization curves and time history of the stacks performance, the temperature at various locations in the system, and the solar irradiation.


## 3. System Modeling

### 3.1. PV Module

The model given by Duffie and Beckman [17] was used to model the PV modules’ current and voltage output. At fixed temperatures and solar radiation, the current-voltage characteristic of a PV module can be modeled using
(1)I=IL-ID=IL-Io(e(V+IRs)/a-1)-V+IRsRsh.
The power is given by
(2)PPV=IV.In most cases, the shunt current can be ignored because the shunt resistance is so high that the term goes to zero, particularly for monocrystalline solar cells, so (1) becomes
(3)I=IL-ID=IL-Io(e(V+IRs)/a-1)
and in terms of voltage becomes
(4)V=a·ln(IL-I+IoIo)-IRs.
The values of IL, Io, Voc, and the parameter a are functions of the cell temperature. The model presented by Duffie and Beckman [17] was used to find these parameters and consequently the operating voltage at any load current.


### 3.2. Electrolyzer Mode

The electrode kinetics of an electrolyzer modeled using the semiempirical current-voltage relationships is presented by Ulleberg [18]. The basic form of the I-V curve is given by (5), and that curve varies for a given temperature:
(5)Uel=Uo-RTnFln(PH2·P020.5)+rAI+slog(tAI+1).
Ulleberg and Mørner [19] used the following equations to count the temperature dependence of the ohmic resistance (r) and the overvoltage coefficients (s and t):
(6)t=t1+t2T+t3T2,r=r1+r2T,s=s1+sT+s3T22.
The electrolyzer used in this study had been tested under different operating temperatures, and with and without ohmic losses to predict the ohmic losses (as well as the rest of the losses). The results are then plotted, and, through curve fitting, the parameters can be found. Table 2 listed the overpotential parameters for the electrolyzer. Once the eight parameters for the electrolyzer mode were found, (5) can be used to describe the I-U for the stack. Table 1 listed the electrolyzer mode parameters.Table 1
Electrolyzer mode parameters.
Parameter
Value
Unit
r
1
9.5 × 10−5
Ω·m2
r
2
−4.7 × 10−7
Ω·m2·°C−1
s
1
0.1689
V
s
2
0.0008
V·°C−1
s
3
−1.33 × 10−5
V·°C−2
t
1
0.33
m2·A−1
t
2
5.48
m2·A−1·°C
t
3
242
m2·A−1·°C2Table 2
Fuel cell mode parameters.
T,oC
Oxygen
Air
E
rev
. V
U
o V
b, mV/dec
i
o, mA/cm2
E
rev
. V
U
o V
b, mV/dec
i
o, mA/cm2
30
1.18
1.028
86.6
0.0016
1.17
0.972
91.1
0.0012
70
1.17
0.999
71.8
0.0035
1.16
0.962
79
0.0028


### 3.3. Fuel Cell Mode

The current-voltage characteristics of a PEM fuel cell used in this study are an empirical equation used by Ulleberg [20] that takes overvoltages due to activation into account. However, the concentration losses were neglected since the operating current density was lower than 200 mA/cm2. The following equation was used to predict the fuel cell mode voltage as a function of the current density:
(7)Ufc=Uo-blogi-ri.
The open-circuit voltage can be measured or calculated from
(8)Uo=Erev+blogio.
Table 2 listed the real operation results for the URFC working on the fuel cell mode at two different cell temperatures and 1 bar pressure fully humidified hydrogen and oxidants.In order to simulate the performance of the system, the electric and thermoelectric equations of each part of the system were solved using MATLAB code. The input parameters were the solar irradiation and location parameters (whether data), proposed load profile, regenerative fuel cell specifications, and the hydrogen and oxygen tanks specifications. The DC/DC controller efficiency was set to 95%.


## 3.1. PV Module

The model given by Duffie and Beckman [17] was used to model the PV modules’ current and voltage output. At fixed temperatures and solar radiation, the current-voltage characteristic of a PV module can be modeled using
(1)I=IL-ID=IL-Io(e(V+IRs)/a-1)-V+IRsRsh.
The power is given by
(2)PPV=IV.In most cases, the shunt current can be ignored because the shunt resistance is so high that the term goes to zero, particularly for monocrystalline solar cells, so (1) becomes
(3)I=IL-ID=IL-Io(e(V+IRs)/a-1)
and in terms of voltage becomes
(4)V=a·ln(IL-I+IoIo)-IRs.
The values of IL, Io, Voc, and the parameter a are functions of the cell temperature. The model presented by Duffie and Beckman [17] was used to find these parameters and consequently the operating voltage at any load current.


## 3.2. Electrolyzer Mode

The electrode kinetics of an electrolyzer modeled using the semiempirical current-voltage relationships is presented by Ulleberg [18]. The basic form of the I-V curve is given by (5), and that curve varies for a given temperature:
(5)Uel=Uo-RTnFln(PH2·P020.5)+rAI+slog(tAI+1).
Ulleberg and Mørner [19] used the following equations to count the temperature dependence of the ohmic resistance (r) and the overvoltage coefficients (s and t):
(6)t=t1+t2T+t3T2,r=r1+r2T,s=s1+sT+s3T22.
The electrolyzer used in this study had been tested under different operating temperatures, and with and without ohmic losses to predict the ohmic losses (as well as the rest of the losses). The results are then plotted, and, through curve fitting, the parameters can be found. Table 2 listed the overpotential parameters for the electrolyzer. Once the eight parameters for the electrolyzer mode were found, (5) can be used to describe the I-U for the stack. Table 1 listed the electrolyzer mode parameters.Table 1
Electrolyzer mode parameters.
Parameter
Value
Unit
r
1
9.5 × 10−5
Ω·m2
r
2
−4.7 × 10−7
Ω·m2·°C−1
s
1
0.1689
V
s
2
0.0008
V·°C−1
s
3
−1.33 × 10−5
V·°C−2
t
1
0.33
m2·A−1
t
2
5.48
m2·A−1·°C
t
3
242
m2·A−1·°C2Table 2
Fuel cell mode parameters.
T,oC
Oxygen
Air
E
rev
. V
U
o V
b, mV/dec
i
o, mA/cm2
E
rev
. V
U
o V
b, mV/dec
i
o, mA/cm2
30
1.18
1.028
86.6
0.0016
1.17
0.972
91.1
0.0012
70
1.17
0.999
71.8
0.0035
1.16
0.962
79
0.0028


## 3.3. Fuel Cell Mode

The current-voltage characteristics of a PEM fuel cell used in this study are an empirical equation used by Ulleberg [20] that takes overvoltages due to activation into account. However, the concentration losses were neglected since the operating current density was lower than 200 mA/cm2. The following equation was used to predict the fuel cell mode voltage as a function of the current density:
(7)Ufc=Uo-blogi-ri.
The open-circuit voltage can be measured or calculated from
(8)Uo=Erev+blogio.
Table 2 listed the real operation results for the URFC working on the fuel cell mode at two different cell temperatures and 1 bar pressure fully humidified hydrogen and oxidants.In order to simulate the performance of the system, the electric and thermoelectric equations of each part of the system were solved using MATLAB code. The input parameters were the solar irradiation and location parameters (whether data), proposed load profile, regenerative fuel cell specifications, and the hydrogen and oxygen tanks specifications. The DC/DC controller efficiency was set to 95%.


## 4. Results and Discussion

Among the testing days, a relatively low solar radiation day was selected, on which the total daily radiation was about 4.1 kWh (refer to Figure3). The system first turned on the electrolyzer mode ELM to generate hydrogen. Figure 4 shows the measured and simulated PV current output and the current given to the load. Its clear to see that the current from the PV during the sunny hours was much higher than the load current, and therefore the excess current goes to the stack to generate hydrogen. The simulated and measured data was almost with relatively low error, especially at a low operating current, and slightly increased when the current increased due to nonconsiderable ohmic loses for the simulation data.Figure 3
Solar irradiation for the selected day (November 19, 2009).Figure 4
PV current output and load current for the selected day (November 19, 2009).The electrolyzer current keeps increasing as the PV current increases and then decreases till reaching to zero at about 4 PM. At that time, the stack switched to fuel cell mode FCM after purging the stack with inert nitrogen. The fuel cell mode current is higher than the load current and the load voltage is 12 V DC, while the fuel cell mode voltage is less than 6 V. Because of this the stack should generate a higher current to meet the power required for the load after DC/DC controls the voltage to the desired load amount. Figure5 shows the stack current for the electrolyzer and fuel cell mode.Figure 5
Fuel cell and electrolyzer mode measured current.The Fuel cell and electrolyzer mode voltage are shown in Figure6. The voltage of the stack during the fuel cell mode is around 6 V (7 cells stack). The measured electrolyzer voltage was slightly higher than the simulated values because the PV output voltage was slightly higher than the electrolyzer required voltage. The bus bar current available at the system is either from the PV current or from the stack through the fuel cell mode. Figure 7 shows the bus bar current distribution for the selected day.Figure 6
Fuel cell mode voltage and electrolyzer mode voltage.Figure 7
Measured fuel cell and electrolyzer mode bus available current.The variation of the stored hydrogen inside the hydrogen tank is shown in Figure8. The initial pressure value was set at 1 bar. The simulation results showed acceptable predictable values comparing to the real data, with a RMS error of 0.07 bars. This was mainly because the gases were assumed to behave as an ideal gas and the real amount of the produced and consumed amount of hydrogen was not the same as the simulated values.Figure 8
Simulated and measured hydrogen tank pressure.The summary of the full measured results for the selected day is shown in Figure9 and all the results are listed in Table 3.Table 3
Summary of the selected day measured and simulation results.
System
Simulated
Measured
Solar input
4.08
4.08
kWh/m2
PV output
0.83
0.78
kWh
PV  efficiencya
10.3
9.7
%
User load
0.37
0.37
kWh
PV to loadb
28.85
31.43
%
Load powered by PVc
65.26
66.94
%
FCM power
0.13
0.12
kWh
FCM operating hours
7
7
hr
Total H2consumedd
0.25
0.27
kWh
FCM energy efficiencye
50.1
45.35
%
ELM power
0.46
0.51
kWh
ELM operating hours
9
9
hr
Total H2produced
0.39
0.35
kWh
ELM energy efficiencyf
84.48
68.98
%
Total system Efficiencyg
4.55
4.55
%
Round trip efficiencyh
42.31
31.28
%
aPV output to the solar input.
b(Load − Fuel cell mode power)/PV output.
c(Load − Fuel cell mode power)/load.
dBased on HHV for hydrogen.
eFuel cell mode power to the hydrogen consumed.
fHydrogen produced to the electrolyzer mode power.
gLoad to the solar input.
hFuel cell mode efficiency and electrolyzer mode  efficiency.Figure 9
Summary of measured data for the selected day.


## 5. Conclusion

The following can be concluded from the current study.(1)
The load will be powered by the PV array at the solar time, and the extra power will go to the electrolyzer to produce hydrogen. After sunset, or at the time of low irradiation, this stored hydrogen will be used to power the load through the fuel cell.(2)
Comparing the simulated results with measured data shows the ability of the cell modeling to predict the cell performance with an acceptable level of error. The difference between the measured and simulated data increased as the operating current increases.(3)
The performance of the stack works on electrolyzer or fuel cell mode was slightly lower when compared with the discrete electrolyzer and fuel cell, which is an expected result since the membrane electrode assembly catalysts were chosen based on compromise between fuel cell and electrolyzer mode performance.



---

*Source: 942784-2012-12-10.xml*