Next Article in Journal
Information Entropy of Tight-Binding Random Networks with Losses and Gain: Scaling and Universality
Next Article in Special Issue
Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion
Previous Article in Journal
Desalination Processes’ Efficiency and Future Roadmap
Previous Article in Special Issue
A LiBr-H2O Absorption Refrigerator Incorporating a Thermally Activated Solution Pumping Mechanism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Analysis of a Proton Exchange Membrane Fuel Cell Based Syngas

1
Department of Physics, Jimei University, Xiamen 361021, China
2
Department of Building and Real Estate, The Hong Kong Polytechnic University, Hong Kong, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(1), 85; https://doi.org/10.3390/e21010085
Submission received: 31 October 2018 / Revised: 24 December 2018 / Accepted: 5 January 2019 / Published: 18 January 2019
(This article belongs to the Special Issue Advances in Applied Thermodynamics III)

Abstract

:
External chemical reactors for steam reforming and water gas shift reactions are needed for a proton exchange membrane (PEM) fuel cell system using syngas fuel. For the preheating of syngas and stable steam reforming reaction at 600 °C, residual hydrogen from a fuel cell and a certain amount of additional syngas are burned. The combustion temperature is calculated and the molar ratio of the syngas into burner and steam reformer is determined. Based on thermodynamics and electrochemistry, the electric power density and energy conversion efficiency of a PEM fuel cell based syngas are expressed. The effects of the temperature, the hydrogen utilization factor at the anode, and the molar ratio of the syngas into burner and steam reformer on the performance of a PEM fuel cell are discussed. To achieve the maximum power density or efficiency, the key parameters are determined. This manuscript presents the detailed operating process of a PEM fuel cell, the allocation of the syngas for combustion and electric generation, and the feasibility of a PEM fuel cell using syngas.

1. Introduction

The hydrogen fuel cell is a promising technology for electric vehicle applications [1]. Compared with battery-based electric vehicles, hydrogen fuel cell-based electric vehicles can easily deliver a much higher range, and the fueling of hydrogen only takes a few minutes. However, it is still challenging to store hydrogen in vehicles with small volume and mass. The currently used hydrogen cylinder is bulky while the metal hydride is relatively heavy. One option to solve the hydrogen storage problem is to supply hydrogen onboard by the steam reforming of ethanol [2,3] or diesel [4]. This option has been demonstrated to be feasible and is practically used. However, hydrogen production from nonrenewable resources results in pollutant emission from the life cycle point of view. Biomass is an important renewable resource and is carbon neutral. The use of biomass-derived fuel in fuel cell vehicles can achieve green transportation with minimal pollutant emission. When biomass-derived fuel is used in a proton exchange membrane (PEM) fuel cell, the syngas production [5], steam reforming, and water gas shift reactions are most necessary to provide pure hydrogen. Although a recent study showed that it is possible to directly use biomass in a fuel cell without converting biomass into hydrogen rich gases [6], the power density is too low for practical applications. Two conversion paths, anaerobic digestion and gasification, are available to convert biomass into hydrogen-rich gases, which contains hydrogen, methane, carbon monoxide, carbon dioxide, and other gases [7]. As methane cannot be electrochemically oxidized in the PEM fuel cell and the carbon monoxide can easily poison the noble metal catalyst in the fuel cell, the further steam reforming and the water gas shift are necessary to convert the residual methane and carbon monoxide to hydrogen, respectively. As the steam reforming reaction is endothermic, the addition of CaO into glycerol steam reforming is proposed [8], since the conversion of CaO into CaCO3 can provide enough heat for the endothermic reaction. However, the new product CaCO3 may cover the catalyst of the steam reforming reaction and thus limits the long-term durability. The high-temperature exhaust heat from a Stirling engine can be used for steam reforming [9], but the operation of the PEM fuel cell is dependent on the engine.
It is better if the steam reforming and the combustion of residual hydrogen from a PEM fuel cell are integrated. The combustion of residual fuel from a fuel cell and directly fed fuel of a burner can generate enough heat for the reforming [10]. A membrane reactor of reforming is employed for the improvement of the system efficiency [11]. The effect of the molar ratio of steam and hydrocarbon on the steam reforming product [12] and system efficiency [13] is present. The molar ratio of the fuel into the burner and steam reformer has almost not been discussed, and it is considered as a parameter in this paper which will affect the performance of the system.
To improve the performance of a PEM fuel cell, various strategies have been proposed. A gas diffusion layer with micropores of the PEM fuel cell is suggested [14] because the low thermal conductivity of it can accelerate the transportation of water while the operating temperature of the fuel cell may be unstable. The addition of silicates to the proton exchange membrane will improve the thermal conductivity, water uptake, and proton conductivity of the membrane [15]. The work in Reference [16] exhibits high water uptake and proton conductivity of the activated carbon/Nafion hybrid composite. The ohmic polarization loss is influenced by the proton conductivity and gas humidity [17]. If there is analytic relation between the proton conductivity and composition of the membrane, the performance of a PEM fuel cell can be optimized. The concentration polarization loss mainly occurs in the cathode [18]. The flow field of the bipolar plate is designed to reduce the loss [19]. The mathematical expressions of activation, ohm, and concentration overpotentials in engineering are used in this study.

2. The Auxiliary Systems outside a PEM Fuel Cell

The schematic diagram of the system is shown in Figure 1, where the syngas and water are preheated from 25 °C to 600 °C in heat exchanger 1 (HE1); the steam reforming reaction of methane and water vapor in the gases is taken place in heat exchanger 2 (HE2); and the carbon monoxide is eliminated after passing the high-temperature water gas shift (HTS), the low-temperature water gas shift (LTS), and the preferential oxidation reactions (PROX). The air flowing into the cathode of a fuel cell is preheated by the gases leaving the fuel cell in heat exchanger 3 (HE3). The residual hydrogen leaving the anode of a fuel cell and a certain amount of additional syngas are burned in the after burner (AB); the high-temperature combustion product is used to supply enough heat for the preheating and steam reforming reactions in HE1 and HE2, respectively.

2.1. The Heat Needed in HE1 and HE2

The water gas shift and preferential oxidation reactions are exothermal; the heat about them will not be discussed. Assuming the heat released from the residual air and water leaving the cathode is enough for the preheating of air flowing into the fuel cell, the process of heat transfer in HE3 will also not be discussed.
If the mole flow rate of syngas into HE1 is n (mol·s−1) and the amount of water in the flow is twice of that needed for the steam reforming and water gas shift reactions, the steam reforming and water gas shift reactions are, respectively,
CH 4 + H 2 O CO + 3 H 2    Δ h > 0
and
CO + H 2 O CO 2 + H 2    Δ h < 0
where Δ h is the enthalpy change of gases if one mole C H 4 is consumed in Equation (1) or one mole CO is consumed in Equation (2). The mole rate of water added in the syngas is
n H 2 O = 2 n ( 2 x CH 4 + x CO ) n x H 2 O ,
where x k is the mole fraction of k in the syngas; the specific components of the syngas are shown in Table 1 [20]. The heat needed for the preheating of the syngas and water is
Q ˙ 1 = n k ( x k 298.15 873.15 C k , m d t ) + n H 2 O 298.15 873.15 C H 2 O , m d t + 2 n ( 2 x CH 4 + x CO ) L H 2 O , m ,
where C k , m is the molar heat capacity of k, the value of which is shown in Table 2 [21,22]; t is temperature; and L H 2 O , m is the latent heat of one mole water. The heat needed for the steam reforming reaction is
Q ˙ 2 = n x CH 4 ( 3 h H 2 + h CO h CH 4 h H 2 O ) ,
where h k is the enthalpy of k per mole at 600 °C and 1atm and is calculated as
h k = h k 0 + 298.15 873.15 C k , m d t ,
where h k 0 is the enthalpy of k per mole at 25 °C and 1 atm, the value of which is shown in Table 2.

2.2. The Molar Ratio of Syngas into the Burner and Steam Reformer

The relative low temperature of a PEM fuel cell will reduce the cost of materials and will present a short start time and a fast transient response [23]. In this paper, the working temperature of the PEM fuel cell is under 100 °C and the water can be separated from the gases which will enter into the anode of the fuel cell. If CO is reacted fully in the water gas shift reactions, the gases flowing into the fuel cell is composed of H 2 , CO 2 , and N 2 by combining Equations (1) and (2).
The total electrochemical reaction in a fuel cell is
H 2 + 0.5 O 2 H 2 O + Electricity .    Δ h < 0
If the hydrogen utilization factor at the anode is u H 2 , the mole flowing rate of residual hydrogen is determined, and the total enthalpy of residual H 2 , CO 2 , and N 2 leaving the fuel cell is
H ˙ c e l l = n ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) h H 2 ( T ) + n ( x CO 2 + x CH 4 + x CO ) h CO 2 ( T ) + n x N 2 h N 2 ( T ) ,
where h k ( T ) is the enthalpy of k per mole at temperature T and 1 atm.
The combustion of residual hydrogen may not be enough for the preheating and steam reforming of the syngas, so extra syngas is added in the burner. The combustion equations in the burners are, respectively,
H 2 + 0.5 O 2 H 2 O ,
CH 4 + 2 O 2 CO 2 + 2 H 2 O ,
and
CO + 0.5 O 2 CO 2 .
Assuming the molar ratio of the syngas into the burner and steam reformer is x and the air in the burner is just enough for the combustion, the enthalpy of gases flowing into the burner unit time is
H ˙ i n = H ˙ c e l l + n x k x k h k 0 + 1 2 n [ ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) + x x H 2 + x x CO + 4 x x CH 4 ] ( h O 2 0 + 78 21 h N 2 0 )
Combining Equations (8)–(12), the component of the combustion product is determined. If the combustion temperature is T C , one can obtain the enthalpy of gases leaving the burner as
H ˙ o u t = n [ ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) + x x H 2 + x x H 2 O + 2 x x CH 4 ] h H 2 O ( T C ) + n ( 1 + x ) ( x CO 2 + x CH 4 + x CO ) h CO 2 ( T C ) + n ( 1 + x ) x N 2 h N 2 ( T C ) + 1 2 n [ ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) + x x H 2 + x x CO + 4 x x CH 4 ] 78 21 h N 2 ( T C ) .
Assuming the burner is adiabatic, the enthalpy of gases flowing into the burner is equal to that leaving the burner:
H ˙ i n = H ˙ o u t .
Combining Equations (8) and (12)–(14), x is a function of u H 2 , T , and T C . T C is not random; it should be larger than 600 °C, and the heat released by the combustion product must be enough for the preheating and steam reforming reaction in HE1 and HE2 as
n [ ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) + x x H 2 + x x H 2 O + 2 x x CH 4 ] ( T C T C C H 2 O , m d t + L H 2 O , m ) + n ( 1 + x ) ( x CO 2 + x CH 4 + x CO ) T C T C C CO 2 , m d t + n ( 1 + x ) x N 2 T C T C C N 2 , m d t + 1 2 n [ ( x H 2 + 4 x CH 4 + x CO ) ( 1 u H 2 ) + x x H 2 + x x CO + 4 x x CH 4 ] 78 21 T C T C C N 2 , m d t Q ˙ 1 + Q ˙ 2 ,
where T C is the temperature of the combustion product leaving HE1. If the rate of heat released from the combustion product is as close to that needed for the preheating and steam reforming reaction as possible, Equation (15) is another relationship between x , u H 2 , T , and T C .
Combining Equations (3)–(6), (8), and (12)–(15), x and T C are both independent of n , while they are functions of u H 2 , T , and T C . If T and T C are given, one can obtain the curves of x and T C varying with u H 2 as shown in Figure 2 and Figure 3, respectively. If the hydrogen utilization factor is smaller, there are more residual hydrogen leaving the fuel cell, the combustion temperature is higher, and the amount of syngas in the burner is less as shown in Figure 2 and Figure 3, so the cost of the steam reforming reaction is lower. However, the hydrogen utilization factor at the anode of a fuel cell will affect the electric power of the fuel cell [24]; there is a specific hydrogen utilization factor under the maximum electric power. The energy conversion efficiency is dependent on both the amount of syngas in the burner and the electric power need to be introduced and optimized. Figure 2 and Figure 3 also show that if the working temperature of the PEM fuel cell is higher, the temperature of the combustion product is higher and the amount of syngas needed for the burner is less.

3. The Electric Power of a PEM Fuel Cell Based Syngas

As defined above, the molar rate of the hydrogen used in Equation (7) is
n ( x H 2 + 4 x CH 4 + x CO ) u H 2 .
According to thermodynamics, the available energy released by the reaction is
n ( x H 2 + 4 x CH 4 + x CO ) u H 2 Δ g ( T ) ,
where Δ g ( T ) is the Gibbs function change of the gases if one molar hydrogen is consumed in Equation (7), and it can be calculated as
Δ g ( T ) = Δ h ( T ) T Δ s ( T ) ,
where Δ h ( T ) and Δ s ( T ) are, respectively, the enthalpy and entropy changes of the gases if one molar hydrogen is consumed in Equation (7). The enthalpy change of the gases is
Δ h ( T ) = h H 2 O ( T ) h H 2 ( T ) 0.5 h O 2 ( T ) .
The entropy change of the gases is
Δ s ( T ) = s H 2 O 0 + 298.15 T C H 2 O , m t d t ( s H 2 0 + 298.15 T C H 2 , m t d t ) 0.5 ( s O 2 0 + 298.15 T C O 2 , m t d t ) + R ln p H 2 p O 2 0.5 p H 2 O ,
where s k 0 is the entropy of per mole k at 25 °C and 1 atm; R is the universal gas constant; p k is the partial pressure of k at the electrode; the partial pressure of liquid water is 1;
p H 2 = p a ( 1 p s p a ) / [ 1 + x a 2 ( 1 + δ a δ a 1 ) ] ;
and
p O 2 = p c ( 1 p s p c ) / [ 1 + x c 2 ( 1 + δ c δ c 1 ) ] ,
where p a (atm) and p c (atm) are pressures at anode and cathode, respectively [25], p s (atm) is the saturation pressure of water and is dependent on the temperature [26], x a and x c are dry gas molar ratios at anode and cathode, respectively, and δ a and δ c are stoichiometry coefficients.
As shown by the equivalent circuit of a PEM fuel cell in Ref. [27], there is internal resistance of a fuel cell. The heat produced by the internal resistance is
I ( V a c t + V o h m + V c o n ) ,
where I is the intensity of the electric current of the PEM fuel cell and V a c t is the activation overpotential. The general expression of it is the Butler–Vollmer equation [28]. There is another pattern reported by Tafel [29] based on the results of the experiments and is defined as
V a c t = R T n e F ( α a + α c ) α a α c ln i i 0 ,
where n e is the number of electrons transferred through the external circuit if one molecule of hydrogen is reacted; F is Faraday’s constant; α a and α c represent the anode and cathode charge transfer coefficients; i = I / A c is the current density of the PEM fuel cell; A c is the effective surface area of the bipolar plate; i 0 is the exchange current density of the electrodes [30] and is defined as
i 0 = 1.27 × 10 8 exp ( 2.06 p O 2 ) ;
and V o h m and V c o n are ohm and concentration overpotentials [31]. The expressions of them are [27]
V o h m = i δ m e m σ m e m
and
V c o n = i ( β 1 i i L ) β 2 ,
respectively, where δ m e m is the thickness of membrane and σ m e m is the membrane conductivity [27]. It can be defined as
σ m e m = ( 0.005139 μ m e m 0.003260 ) exp [ 1268 ( 1 303 1 T ) ] ,
where μ m e m is the water content and is determined by water vapor activity, i L is the limiting current density, β 1 is dependent on the partial pressure of oxygen and the temperature at the cathode [30], and β 2 is constant.
According to energy conservation, the electric power of a PEM fuel cell is
P e = n ( x H 2 + 4 x CH 4 + x CO ) u H 2 Δ g ( T ) I ( V a c t + V o h m + V c o n ) .
Based on Faraday’s law, the intensity of the electric current in Equation (29) can be presented as
I = n ( x H 2 + 4 x CH 4 + x CO ) u H 2 n e F .
Combining Equations (18)–(30) and Table 3 [25,30,32], the electric power density “ P e * = P e / A c ” is dependent on “ n * = n / A c ”, u H 2 , and T . If T is given, one can obtain the maximum power density varying with n * by optimizing the hydrogen utilization factor, as shown by the solid curve in Figure 4. The hydrogen utilization factor under the maximum power density is shown by the solid curve in Figure 5. There is a certain electric current density under the maximum power density of the PEM fuel cell [24]. This means that there is an optimal product of n * and u H 2 from Equation (30). When n * is smaller than 0.17 (mol·s−1·m−2), the product of n * and u H 2 is always smaller than the optimum value, even though u H 2 = 1 in Figure 5 and the electric power density of the PEM fuel cell is a monotonous increasing function of the molar flow rate of the syngas in Figure 4. When n * 0.17 (mol·s−1·m−2), the hydrogen utilization factor is decreasing with the increase of the molar flow rate of syngas to reach the optimal product of them in Figure 5, and the electric power density is constant in Figure 4.
There is a linear relationship between x and u H 2 in Figure 2, so one can also obtain x under the maximum power density, as shown by the solid curve in Figure 6. The graphs of Figure 5 and Figure 6 are similar. When there is no residual hydrogen from the fuel cell, the molar ratio of the syngas into the burner and steam reformer is 0.4. As the amount of residual hydrogen increases, the molar ratio x is decreasing. When the hydrogen utilization factor is about 0.61, there is no need to add syngas into the burner. As discussed in Section 2.2, the total energy conversion efficiency of the hybrid system of a PEM fuel cell and steam reformer is dependent on the hydrogen utilization factor at the anode of the fuel cell, so the efficiency of the hybrid system will be further optimized.

4. The Total Energy Conversion Efficiency

The energy conversion efficiency of the hybrid system of a PEM fuel cell and steam reformer is
η = P e n ( 1 + x ) ( x H 2 q LHV ( H 2 ) + x CO q LHV ( CO ) + x CH 4 q LHV ( CH 4 ) ) = ( x H 2 + 4 x CH 4 + x CO ) u H 2 [ Δ g ( T ) n e F ( V a c t + V o h m + V c o n ) ] ( 1 + x ) ( x H 2 q LHV ( H 2 ) + x CO q LHV ( CO ) + x CH 4 q LHV ( CH 4 ) ) ,
where q LHV ( k ) is the lower heating value of per molar k. Combining Equations (18)–(31), η is dependent on n * , u H 2 , x, and T , while there is a linear relationship between u H 2 and x. If temperature T is given, one can obtain the maximum efficiency varying with the molar flowing rate of the syngas by optimizing the hydrogen utilization factor in the fuel cell, as shown by the dash curve in Figure 7. The electric power density, hydrogen utilization factor, and x under the maximum efficiency are, respectively, shown by the dash curves in Figure 4, Figure 5 and Figure 6. If the values of u H 2 and x derived in Section 3 are substituted into Equation (31), one can obtain the efficiency under the maximum power density as shown by the solid curve in Figure 7. The efficiency of the system is decreasing with the increase of the molar flowing rate of syngas.
When n * is smaller than 0.065 (mol·s−1·m−2), the working conditions of the maximum efficiency and the maximum power density are the same: The hydrogen should be reacted totally in the fuel cell, and the molar ratio of the syngas into the burner and steam reformer is 0.4. When 0.065 n * 0.14 (mol·s−1·m−2), there should be residual hydrogen leaving the fuel cell and x is smaller than 0.4 for the maximum efficiency of the system. When n * is larger than 0.14 (mol·s−1·m−2), the hydrogen utilization factor is constant and about 0.61; there is no need to add syngas into the burner.

5. Conclusions

A model of the PEM fuel cell based syngas is established, and the residual hydrogen leaving the fuel cell and extra syngas are burned in a burner to supply enough high-temperature heat for the steam reforming reaction. There is a linear relationship between the hydrogen utilization factor in the fuel cell and the molar ratio of the syngas into the burner and steam reformer; if the hydrogen utilization factor is higher, the molar ratio is larger, and vice versa. Based on the thermodynamics and electrochemistry, the expressions of the electric power and energy conversion efficiency of the system are derived. For the maximum power density or maximum efficiency of a PEM fuel cell based syngas, the optimal hydrogen utilization factor in the fuel cell and the molar ratio of the syngas into the burner and steam reformer are, respectively, determined. In the future, the waste heat from a PEM fuel cell and the water gas shift reactions can be considered to be utilized to improve the energy conversion efficiency of the PEM fuel cell.

Author Contributions

Conceptualization, X.Z., Q.L., H.L., X.C., and S.S.; software, X.Z.; Writing—review and editing, M.N. and X.Z.

Funding

This work was funded by the National Natural Science Foundation, grant number 11605073 and 51601068; Natural Science Foundation of Fujian Province, grant number 2017J05009; and Fujian Province for Distinguished Young Scholars, grant number B17164.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A c Effective surface area of the bipolar plate
α a Charge transfer coefficient at anode
α c Charge transfer coefficient at cathode
β 1 Parameter in the expression of overpotential
β 2 Constant in the expression of overpotential
C k , m Molar heat capacity of k
Δ h Enthalpy change of gases
Δ h ( T ) Enthalpy change of gases
Δ g ( T ) Gibbs function change of gases
Δ s ( T ) Entropy change of gases
δ a Stoichiometry coefficient
δ c Stoichiometry coefficient
δ m e m Membrane thickness
F Faraday constant
η Efficiency of the hybrid system
h k Molar enthalpy of k at 600 °C
h k 0 Molar enthalpy of k at 25 °C
h k ( T ) Molar enthalpy of k at temperature T
H ˙ c e l l Enthalpy of gases leaving the fuel cell
H ˙ i n Enthalpy of gases into burner
H ˙ o u t Enthalpy of gases leaving the burner
I Electric current of a PEM fuel cell
i 0 Exchange current density of the electrodes
i Current density of a PEM fuel cell
i L Limiting current density of a PEM fuel cell
L H 2 O , m Latent heat of one mole water
μ m e m Water content
n Mole flow rate of syngas into HE1
n H 2 O Mole rate of water added into syngas
n * Mole flow rate of syngas into HE1 unit area
n e Number of electrons
P e * Electric power density of a PEM fuel cell
P e Electric power of a PEM fuel cell
p c Pressure at the cathode
p a Pressure at the anode
p s Saturation pressure of water
p k Partial pressure of k at the electrode
Q ˙ 2 The heat needed in HE2 unit time
Q ˙ 1 The heat needed in HE1 unit time
q L H V ( k ) Lower heating value of k per molar
R Universal gas constant
s k 0 Molar entropy of k at 25°C
σ m e m Membrane conductivity
T Operating temperature of a PEM fuel cell
T C Combustion temperature
t Temperature
T C Temperature of gases leaving HE1
u H 2 Hydrogen utilization factor in a fuel cell
V a c t Activation overpotential
V o h m Ohm overpotential
V c o n Concentration overpotential
x c Dry gas molar ratio at cathode
x a Dry gas molar ratio at anode
x The molar ratio
x k Molar fraction of k in syngas

References

  1. Brey, J.J.; Carazo, A.F.; Brey, R. Exploring the marketability of fuel cell electric vehicles in terms of infrastructure and hydrogen costs in Spain. Renew. Sustain. Energy Rev. 2018, 82, 2893–2899. [Google Scholar] [CrossRef]
  2. Ham, S.W.; Jo, S.Y.; Dong, H.W.; Jeong, J.W. A simplified PEM fuel cell model for building cogeneration applications. Energy Build. 2015, 107, 213–225. [Google Scholar] [CrossRef]
  3. Izurieta, E.M.; Adrover, M.E.; Pedernera, M.N.; López, E. Ethanol processor design for hydrogen production. Kinetic analysis and process integration. Ind. Eng. Chem. Res. 2018, 57, 13615–13626. [Google Scholar] [CrossRef]
  4. Maximini, M.; Engelhardt, P.; Brenner, M.; Beckmann, F.; Moritz, O. Fast start-up of a diesel fuel processor for PEM fuel cells. Int. J. Hydrog. Energy 2014, 39, 18154–18163. [Google Scholar] [CrossRef]
  5. Peng, W.X.; Wang, L.S.; Mirzaee, M.; Ahmadi, H.; Esfahani, M.J.; Fremaux, S. Hydrogen and syngas production by catalytic biomass gasification. Energy Convers. Manag. 2017, 135, 270–273. [Google Scholar] [CrossRef]
  6. Zhao, X.; Liu, W.; Deng, Y.; Zhu, J.Y. Low-temperature microbial and direct conversion of lignocellulosic biomass to electricity: Advances and Challenges. Renew. Sustain. Energy Rev. 2017, 71, 268–282. [Google Scholar] [CrossRef]
  7. Guan, T.; Chutichai, B.; Alvfors, P.; Arpornwichanop, A. Biomass-fuelled PEM FC systems: Evaluation of two conversion paths relevant for different raw materials. Energy Convers. Manag. 2015, 106, 1183–1191. [Google Scholar] [CrossRef]
  8. Wess, R.; Nores-Pondal, F.; Laborde, M.; Giunta, P. Single stage H2 production, purification, and heat supply by means of sorption-enhanced steam reforming of glycerol. A thermodynamic analysis. Chem. Eng. Sci. 2015, 134, 86–95. [Google Scholar] [CrossRef]
  9. Shin’ya, O.; Itaru, T.; Shunsuke, K.; Akira, H.; Seizi, S. Exergy analysis of the woody biomass Stirling engine and PEM-FC combined system with exhaust heat reforming. Int. J. Hydrog. Energy 2008, 33, 2289–2299. [Google Scholar]
  10. Purnima, P.; Jayanti, S. A high-efficiency, auto-thermal system for on-board hydrogen production for low temperature PEM fuel cells using dual reforming of ethanol. Int. J. Hydrog. Energy 2016, 41, 13800–13810. [Google Scholar] [CrossRef]
  11. Salemme, L.; Menna, L.; Simeone, M.; Volpicelli, G. Energy efficiency of membrane-based fuel processors-PEM fuel cell systems. Int. J. Hydrog. Energy 2010, 35, 3712–3720. [Google Scholar] [CrossRef]
  12. Comas, J.; Dieuzeide, M.; Baronetti, G.; Laborde, M.; Amadeo, N. Methane steam reforming and ethanol steam reforming using a Ni(II)-Al(III) catalyst prepared from lamellar double hydroxides. Chem. Eng. J. 2006, 118, 11–15. [Google Scholar] [CrossRef]
  13. Ersoz, A.; Olgun, H.; Ozdogan, S. Reforming options for hydrogen production from fossil fuels for PEM fuel cells. J. Power Sources 2006, 154, 67–73. [Google Scholar] [CrossRef]
  14. Ito, H.; Heo, Y.; Ishida, M.; Nakano, A.; Someya, S.; Munakata, T. Application of a self-supporting micro-porous layer to gas diffusion layers of proton exchange membrane fuel cells. J. Power Sources 2017, 342, 393–404. [Google Scholar] [CrossRef]
  15. Mishra, A.K.; Bose, S.; Kuila, T.; Kim, N.H.; Lee, J.H. Silicate-based polymer-nanocomposite membranes for polymer electrolyte membrane fuel cells. Prog. Polym. Sci. 2012, 37, 842–869. [Google Scholar] [CrossRef]
  16. Chien, H.C.; Tsai, L.D.; Lai, C.M.; Lin, J.N.; Zhu, C.Y.; Chang, F.C. Characteristics of high-water-uptake activated carbon/Nafion hybrid membranes for proton exchange membrane fuel cells. J. Power Sources 2013, 226, 87–93. [Google Scholar] [CrossRef]
  17. Kim, H.; Kim, K. Numerical study on the effects of gas humidity on proton-exchange membrane fuel cell performance. Int. J. Hydrog. Energy 2016, 41, 11776–11783. [Google Scholar] [CrossRef]
  18. Liu, J.X.; Guo, H.; Ye, F.; Ma, C.F. Two-dimensional analytical model of a proton exchange membrane fuel cell. Energy 2017, 119, 299–308. [Google Scholar]
  19. Kahraman, H.; Orhan, M.F. Flow field bipolar plates in a proton exchange membrane fuel cell: Analysis & modeling. Energy Convers. Manag. 2017, 133, 363–384. [Google Scholar]
  20. Zhang, X.; Wang, Y.; Liu, T.; Chen, J. Theoretical basis and performance optimization analysis of a solid oxide fuel cell-gas turbine hybrid system with fuel reforming. Energy Convers. Manag. 2014, 86, 1102–1109. [Google Scholar] [CrossRef]
  21. Shin, Y.; Park, W.; Chang, J.; Park, J. Evaluation of the high temperature electrolysis of steam to produce hydrogen. Int. J. Hydrog. Energy 2007, 32, 1486–1491. [Google Scholar] [CrossRef]
  22. Zhang, X.; Li, G.; Li, J.; Feng, Z. Numerical study on electric characteristics of solid oxide fuel cells. Energy Convers. Manag. 2007, 48, 977–989. [Google Scholar] [CrossRef]
  23. Strahl, S.; Costa-Castelló, R. Temperature control of open-cathode PEM fuel cells. IFAC-PapersOnLine 2017, 50-1, 11088–11093. [Google Scholar]
  24. Zhang, X.; Ni, M.; He, W.; Dong, F. Theoretical analysis and optimum integration strategy of the PEM fuel cell and internal combustion engine hybrid system for vehicle applications. Int. J. Energy Res. 2015, 39, 1664–1672. [Google Scholar] [CrossRef]
  25. Rowe, A.; Li, X. Mathematical modeling of proton exchange membrane fuel cells. J. Power Sources 2001, 102, 82–96. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, Y.; Mawardi, A.; Pitchumani, R. Numerical studies on an air-breathing proton exchange membrane (PEM) fuel cell stack. J. Power Sources 2007, 173, 264–276. [Google Scholar] [CrossRef]
  27. Zhang, X.; Guo, J.; Chen, J. The parametric optimum analysis of a proton exchange membrane (PEM) fuel cell and its load matching. Energy 2010, 35, 5294–5299. [Google Scholar] [CrossRef]
  28. Chan, S.H.; Low, C.F.; Ding, O.L. Energy and exergy analysis of simple solid-oxide fuel-cell power systems. J. Power Sources 2002, 103, 188–200. [Google Scholar] [CrossRef]
  29. Larminie, J.; Dicks, A. Fuel Cell Systems Explained; Wiley: New York, NY, USA, 2003. [Google Scholar]
  30. Ay, M.; Midilli, A.; Dincer, I. Exergetic performance analysis of a PEM fuel cell. Int. J. Energy Res. 2006, 30, 307–321. [Google Scholar] [CrossRef]
  31. Pukrushpan, J.; Stefanopoulou, A.; Peng, H. Modeling and control for PEM fuel cell stack system. In Proceedings of the American Control Conference, Anchorage, AK, USA, 8–10 May 2002; IEEE: Piscataway, NJ, USA, 2002. [Google Scholar]
  32. Wang, Y.; Chen, K.; Mishler, J.; Cho, S.; Adroher, X. A review of polymer electrolyte membrane fuel cells: Technology, applications, and needs on fundamental research. Appl. Energy 2011, 88, 981–1007. [Google Scholar] [CrossRef]
Figure 1. The schematic diagram of a proton exchange membrane (PEM) fuel cell based syngas, where abbreviations of HE1, HE2, and HE3 are heat exchangers, SR is a steam reformer, HTS and LTS are water gas shift reactions, PROX is a preferential oxidizer, and AB is an auxiliary burner.
Figure 1. The schematic diagram of a proton exchange membrane (PEM) fuel cell based syngas, where abbreviations of HE1, HE2, and HE3 are heat exchangers, SR is a steam reformer, HTS and LTS are water gas shift reactions, PROX is a preferential oxidizer, and AB is an auxiliary burner.
Entropy 21 00085 g001
Figure 2. Curves of x varying with u H 2 , where T C = 30 °C.
Figure 2. Curves of x varying with u H 2 , where T C = 30 °C.
Entropy 21 00085 g002
Figure 3. Curves of T C varying with u H 2 , where T C = 30 °C.
Figure 3. Curves of T C varying with u H 2 , where T C = 30 °C.
Entropy 21 00085 g003
Figure 4. The power density varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are, respectively, the maximum power density and the power density under the maximum efficiency.
Figure 4. The power density varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are, respectively, the maximum power density and the power density under the maximum efficiency.
Entropy 21 00085 g004
Figure 5. The hydrogen utilization factor varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are, respectively, the hydrogen utilization factor under the maximum power density and the maximum efficiency.
Figure 5. The hydrogen utilization factor varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are, respectively, the hydrogen utilization factor under the maximum power density and the maximum efficiency.
Entropy 21 00085 g005
Figure 6. Curves of x varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are x under the maximum power density and the maximum efficiency, respectively.
Figure 6. Curves of x varying with n * , where T = 70 °C and T C = 30 °C. The solid and dash curves are x under the maximum power density and the maximum efficiency, respectively.
Entropy 21 00085 g006
Figure 7. The efficiency varying with n * , where T = 70 °C and T C = 30 °C. The dash and solid curves are, respectively, the maximum efficiency and the efficiency under the maximum power density.
Figure 7. The efficiency varying with n * , where T = 70 °C and T C = 30 °C. The dash and solid curves are, respectively, the maximum efficiency and the efficiency under the maximum power density.
Entropy 21 00085 g007
Table 1. The composition of syngas [20].
Table 1. The composition of syngas [20].
Component k H 2 CH 4 CO CO 2 H 2 O N 2
Mole fraction of k in syngas: xk0.130.010.160.050.360.29
Table 2. Thermodynamic parameters of the chemical components [21,22], where (g) and (l) refer to gas and liquid phases, respectively.
Table 2. Thermodynamic parameters of the chemical components [21,22], where (g) and (l) refer to gas and liquid phases, respectively.
Component k h k 0
(J·mol−1)
s k 0
(J·mol−1·K−1)
L H 2 O , m
(J·mol−1)
Molar Heat Capacity C k , m
(J·mol−1·K−1)
N2029.12
O20205.13825.8911 + 0.0129874t − 0.0000038644t2
CH4−75,00014.1555 + 0.0755466t − 0.0000180032t2
CO2−393,80026.0167 + 0.0435259t − 0.0000148422t2
CO−110,50026.8742 + 0.006971t − 0.0000008206t2
H20130.69529.0856 − 0.0008373t + 0.0000020138t2
H2O (g)−241,80030 + 0.01071t + 33000/t2
H2O (l)−285,80069.94040,70075.44
Table 3. Parameters used in the system.
Table 3. Parameters used in the system.
ParameterValue
Number of electrons, n e 2
Faraday constant, F (C mol−1) 96485
Universal gas constant, R (J·mol·K−1)8.314
Pressure at the anode, p a (atm)3 [13]
Pressure at the cathode, p c (atm)5 [13]
Anode stoichiometry, δ a 1.5 [20]
Cathode stoichiometry, δ c 3 [20]
Dry gas molar ratio at anode, x a ( x N 2 + x CO 2 + x CH 4 + x CO ) / ( x H 2 + x CO + 4 x CH 4 ) [13]
Dry gas molar ratio at cathode, x c 3.762 (air)
Charge transfer coefficient at the anode, α a 0.5 [20]
Charge transfer coefficient at the cathode, α c 1 [20]
Membrane thickness, δ m e m (cm)0.018 [13]
μ m e m 14 [21]
Constant, β 2 2 [20]
Limiting current density, i L (A cm−2)2 [20]
T = 70 °C: p s (atm); β 1 0.3071; 0.2048
q LHV ( k ) (kJ mol−1): k = H2; CO; CH4241.9; 283.2; 803.7

Share and Cite

MDPI and ACS Style

Zhang, X.; Lin, Q.; Liu, H.; Chen, X.; Su, S.; Ni, M. Performance Analysis of a Proton Exchange Membrane Fuel Cell Based Syngas. Entropy 2019, 21, 85. https://doi.org/10.3390/e21010085

AMA Style

Zhang X, Lin Q, Liu H, Chen X, Su S, Ni M. Performance Analysis of a Proton Exchange Membrane Fuel Cell Based Syngas. Entropy. 2019; 21(1):85. https://doi.org/10.3390/e21010085

Chicago/Turabian Style

Zhang, Xiuqin, Qiubao Lin, Huiying Liu, Xiaowei Chen, Sunqing Su, and Meng Ni. 2019. "Performance Analysis of a Proton Exchange Membrane Fuel Cell Based Syngas" Entropy 21, no. 1: 85. https://doi.org/10.3390/e21010085

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop