Document Type : Full Lenght Research Article
Authors
Faculty of Engineering, The University of Guilan
Abstract
Keywords
Main Subjects
A B S T R A C T In the present study, the operating conditions of an anode-supported high temperature direct internal reforming (DIR) planar solid oxide fuel cell (SOFC) are numerically analyzed. SOFCs are the energy conversion devices that produce electricity and heat directly from a gaseous fuel by the electrochemical combination of that fuel with an oxidant. A planar SOFC consists of an interconnect structure and a three-layer region often referred to as the PEN (Positive electrode/Electrolyte/Negative electrode) composed of two ceramic electrodes, anode and cathode, separated by a dense ceramic electrolyte. A 1D mathematical model based on the species and the energy balances is conducted to investigate the concentration of the gas components in the channels and the cell temperature distribution in both solid structures and gas channels. In addition, the electrochemical model of the cell is applied to determine its performance parameters. A synthesis gas with 10% pre-reformed methane as fuel stream is fed to the cell and the water-gas shift, the methane-steam reforming and the electrochemical reactions are considered as mass and heat sources/sinks. The effects of the air ratio and the fuel utilization factor on the temperature field and the performance of the cell are examined through parametric analysis to determine the optimum operating conditions of the cell. The investigations are performed to ensure no cracking will occur due to cell crucial temperature difference. The results indicate that increasing the fuel utilization factor from 0.45 to 0.85 leads to increase the value of the maximum temperature difference across the PEN layer from 101K to 145K while it decreases from 177K to 124K by increasing the air ratio from 6.5 to 12.5. In order to achieve the maximum power density by considering the maximum allowable temperature difference across the cell, the value of the fuel utilization factor and the air ratio are obtained 0.85 and 9, respectively. © 2016 Published by Semnan University Press. All rights reserved. |
History: Received 8 February 2014 Received in revised form 24 April 2014 Accepted 27 September 2015
Keywords: Numerical Analysis Fuel cell DIR-PSOFC
|
The Numerical analysis of an anode-supported high temperature DIR-PSOFC operating conditions with considering the maximum allowable temperature difference
A.Etemadi*, Saba Ghorbani, Mehdi Masoumpour, Mostafa Dadkhah, Kazem Atashkari
Faculty of Engineering, The University of Guilan- Guilan, Iran
|
Journal of Heat and Mass Transfer Research
Journal homepage: http://jhmtr.journals.semnan.ac.ir |
|
|
1. Introduction
Today, concerns of the environmental and energy crisis, which are due to the rising energy demands, especially in developing countries, result in more attention to the alternative energy conversion methods, e.g. the wind turbines, geothermal energy, photovoltaic, fuel cells and etc. the Fuel cells, as a promising technology to produce the clean power based on renewable energy resources, have attracted researchers to explore the various aspects of its technology in order to develop it as a feasible and reliable power generation method. The Fuel cells are electrochemical devices that produce electricity via electrochemical reactions between fuel and air flows without occurring any combustion and related exergy destructions. Among different types of fuel cells, Solid Oxide Fuel Cell, due to its high operating temperature and capability to work with a wide variety of fuels, appears to be the most appropriate one for power generation purposes, whether using in distributed power generation systems or combining with conventional power plant cycles to improve their efficiencies. From the geometric viewpoint, there are two main types of SOFCs: tubular and planar which are presented in Fig.1. It is necessary to do theoretical investigations where sealing problems and cracking of the solid structure which are due to the high temperature gradients and consequent thermal stresses particularly in power demand fluctuations are probable. Thus, detailed numerical simulations are necessary in order to illuminate the effect of the physical or chemical phenomena that occurring in the gas channels of the Planar SOFCs (PSOFC).
Fig.1. Schematic of tubular and planar SOFCs.
|
In 1994, Achenbach [1] proposed a three dimensional dynamic modeling of a planar SOFC stack with taking into account the radiation heat transfer mechanism. Moreover, anode gas recycling has also been considered and results indicated the higher plant efficiency and more uniform current density distribution. Ferguson et al [2] investigated the effect of gas channels flow configuration on the planar SOFC performance through a three dimensional modeling. A quasi-two and three dimensional simulation was done by Iwata et al [3] on the mass, the heat and charge transfer in a planar SOFC. Results show an increasing the cell temperature which is due to the fact that increasing the cell power can be controlled by increasing of the operating pressure. Yakabe et al [4] made a computational analysis on momentum and heat transfer of a planar SOFC unit by using syngas as fuel in two cases of external and internal reforming of methane. Their results demonstrated that the internal reforming reaction causes a drop in temperature at the inlet which produces a tensile stress in electrolyte structure. Aguiar et al [5] made a steady state study on IIR-SOFC modeling by using methane as fuel to find the relation between different reforming catalysts and generated power. The results show that the use of typical metal-based (e.g. Ni) steam reforming catalysts leads to full methane consumption and undesirable local cooling in the first one third of the reformer length. Also Aguiar et al [6] presented the steady state modeling of a planar high temperature anode supported DIR-SOFC by applying syngas as fuel in both co-flow and counter-flow operating conditions. They reported an output voltage of 0.66V, a power density of 0.33Wcm−2, and a fuel efficiency of 47% in case of using a fuel mixture with 10% pre-reformed methane and 75% fuel utilization, inlet fuel and air temperatures of 1023 K, average current density of 0.5Acm−2, and an air ratio of 8.5. In addition, Aguiar et al [7] extended their previous research to a dynamic modeling in order to determine the response of open-loop system against changes in the operating conditions and design appropriate closed-loop system proportional to changes in power demands. Colpan et al [8] analyzed the performance of a planar high temperature anode supported DIR-SOFC thermodynamically by using syngas as fuel. By considering fuel channel gas recirculation, they found that the recirculation ratio has no significant effect power output and electrical efficiency of the cellin case of low current density conditions. Meng Ni et al [9] compared the performance of methane fed planar SOFC with two different types of electrolyte: a) using proton conducting electrolyte and b) using oxygen ion conducting electrolyte, and proved that the latter one has higher performance regardless of previous results, which is due to lower ohm over potential. Jeon et al [10] accomplished a two-dimensional computational modeling of a planar high temperature anode supported SOFC and solved the conservation of ions, mass, momentum and energy equations. As an important point, their results demonstrate that the contribution of potential gain due to temperature rising was considerably high. Iwai et al [11] simulated a planar high temperature anode supported DIR-SOFC numerically to investigate the cell temperature distribution and the local current density. The gas channels are considered as porous media and the analysis was done for the different pre-reforming rates, considering the current density is produced by the carbon monoxide electrochemical reaction. Results illustrate that the current density produced due to hydrogen oxidation is quietly higher than twice of which produced by carbon monoxide.
In the current study, the gaseous components and the temperature distribution of one cell of a planar anode supported high temperature SOFC is
Table .1. Reactions considered in the mass and energy balances and rate of them
|
Reactions |
Rate of reactions |
Steam reforming reaction |
||
Water-gas shift reaction |
||
Hydrogen oxidation reaction |
|
|
Oxygen reduction reaction |
|
|
Overall cell reaction |
numerically investigated. The governing equations including mass fractions for the fuel and air flow in gas channels and also energy balances in the domain are solved numerically along with the appropriate boundary conditions by using a numerical procedure through a code using MATLAB in order to determine the mass fractions and the temperature distribution in the gas channels and solid structures. The main advantage of the present work is that the thermal conductivity of the fuel and air components are function of their temperature in gas channels and are not considered to be constant. In order to prevent the emergence of severe thermal gradient and to find out optimal operating condition, the effects of the SOFC operating conditions on the temperature field during the cell are investigated.
2. Mathematical model
In this study, 1D steady-state solid oxide fuel cell has been developed. A complete solid oxide fuel cell is not only consisting of a single cell, but also the several cells, which are connected both in series and/or in parallel. The developed model includes the mass, thermal and electrochemical analysis, as well as the kinetic models of hydrocarbon reactions. Also the pressure drop along the cell channel is assumed to be zero [6].
For the purpose of modeling, it is considered that the SOFC is composed of the fuel and air channels and of the electrochemical region represented by anode, electrolyte, and cathode. The air which has been comprised of and is considered as cathode channel reactant and 10% pre-reformed methane with steam to carbon 2 is considered as fuel for anode channel. Two types of chemical reactions (methane reforming and water - gas shift reaction) occur at the anode while the cells are supplied with the natural gas because the cell operating temperature was sufficiently high [12]. Table 1 presents the cell reactions and the rate of them.
2.1. Mass balance
Both air and fuel flows are assumed to be steady and laminar. The Molar fluxes in the gas channels are considered to be convective in flow direction and the gas flow in the same direction. Mass balances are derived in plug-flow form along the longitudinal direction, because of high Peclet numbers [6]. By Applying the Faraday’s law, the rates of the electrochemical reactions are related to the electric current [1] which can be expressed by:
(1)
|
Steam reforming and water-gas shift reactions occur on the anode surface. The reaction rates for steam reforming and water-gas shift reactions are calculated from Achenbach [1] and Aguiar [6]
(2) |
|
(3) |
where , and is the equilibrium constant for shift reaction which is expressed by
(4) |
The mass balance equation for air channel can be written as
(5) |
And for the fuel channel can be written as
, |
(6) |
Where is the stoichiometry coefficient of component i in reaction k. And the boundary conditions are
(7) |
2.2. Energy balance
In a fuel cell, the current and temperature distributions are strongly coupled. Therefore, knowledge of the predominant heat transfer processes is required. Heat transfer processes include heat release from electrochemical reactions, electrical resistances, convective heat transfer between cell components and fuel and air gas streams, heat conduction through cell components and radiation between solid parts.
Energy balance has to be established for the fuel and air flows and the solid structure separately. It is assumed that all reaction enthalpies are released at the electrode/electrolyte interfaces or at three phase boundaries (TPB). Therefore, the reaction enthalpies appear in the source terms of the solid structure balances.
The energy balances for the gas flows include the convective heat transfer from channel walls to the gas and heat capacity of the gases with respect to their composition. Thermal flux in the PEN and interconnector structures is conductive. And due to the high temperatures, radiation heat transfer between PEN and interconnect structures is taken into account. It is assumed that the cell is in the center of a stack, so there is no heat transfer through the external walls. Thus the energy balances are calculated through Eqs. (8)–(11).
Air channel
|
(8) |
Fuel channel
|
(9) |
PEN structure
|
(10) |
Interconnector
|
(11) |
And the boundary conditions are written as
(12) |
|
(13) |
Heat transfer coefficients between gas channels and solid structures were calculated by using a constant Nusselt number of 3.09 [6] and the channel hydraulic diameter. Therefore, these coefficients are expressed by
(14) |
||
(15) |
||
Where and are the thermal conductivity of air and fuel flow, respectively.
2.3. Electrochemical model
The reversible output voltage of the SOFC is expressed by the Nernst equation [1] as a function of the Gibbs free energy change during electrochemical reaction. This is called open circuit voltage of the cell:
(16) |
where is the Gibbs free energy change from reactants to products ( ), F=96485.3 is the Faraday constant ( ), R=8.314 is the universal gas contact ( ), T is the electrochemical reaction temperature ( ) and , and are the reacting gas pressures ( ) of hydrogen, oxygen and vapor, respectively. Is expressed by:
(17) |
|
(18) |
Where , and are the Gibbs free energy ( ) of vapor, hydrogen and oxygen respectively. is the corresponding enthalpy ( ) and is the corresponding entropy ( ). The gas enthalpy and entropy is calculated by referring to the JANAF Table [13].
Reversible voltage is the maximum potential that can be produced by a cell under specific operating conditions. Actually, the output voltage of the cell, due to over potentials, is lower than the reversible voltage and can be calculated from Eq. (19):
(19)
|
Where is ohmic loss and , , and are the anode and cathode concentration and activation overpotential losses, respectively.
2.4. Overpotentials
Ohmic, concentration and activation overpotentials are the three types of losses that occur during the operation of the cell. Definitely, the ohmic loss is defined as the resistance of the cell solid structures against the conduction of ions and electrons which is expressed by:
(20) |
Where j is the current density and is the internal resistance of the cell, which contains electronic and ionic resistances. is written as
(21) |
Where , and are thickness of anode, cathode and electrolyte, respectively. Also and are electronic conductivities of the anode and cathode, and is ionic conductivity of the electrolyte [6].
Concentration overpotential is the other loss in the cell performance, which is produced due to the difference between species partial pressures in the channel flow and three phase boundary (TPB). In order to occurring electrochemical reaction that occurs on the TPB, the gaseous species must diffuse through porous electrodes which called Knudsen diffusion. In addition, each species should diffuse through the other species via bulk diffusion. Due to these facts, the values of the local concentration species on the TPB differ with their values in the bulk fuel flow and, consequently, the species partial pressures in the Nernst Equation should be rewritten for TPB which expressed in Eq. (22).
(22)
|
The difference between Eqs. (16) and (22) produces the concentration overpotential
(23) |
Where, the first and the second terms on the right hand side are anodic and cathodic concentration overpotentials, respectively.
Partial pressure of water, hydrogen and oxygen at the three- phase boundary can be calculated by
(24) |
|
(25) |
|
(26) |
Where and are average effective diffusivity coefficients at anode and cathode, respectively and is written as
(27) |
is the material porosity and is tortuosity factor. is expressed by
(28) |
Bulk diffusivity, is caused by fluid–fluid intermolecular collisions as considered by the kinetic theory of gases. Knudsen diffusivity is the result of fluid-wall collisions, when the pore size is so small that it is less than the mean path of the fluid.
Activation overpotential is extra potential necessary to defeat the energy barrier and transfer an electron between electronic the ionic conductors. This overpotential reflects the kinetics involved with the electrochemical reactions.
The activation overpotential is expressed by the Butler-Volmer equation
(29) |
Where is the transfer coefficient (usually =0.5) [9], n is the number of electrons participating in electrochemical reaction and is the exchange current density, expressed as
(30) |
Where and values are presented in the Table 2.
3. Result and discussion
The anode supported design is the most favorable one for SOFC [14], therefore this study is carried out on an anode support solid oxide fuel cell with internal reformer. In order to verify the results, the obtained results are compared with the data from Aguiar[6] which is represented in table 4. Gas composition and temperature profile of each layer show s an acceptable agreement with Aguiar’s results. The comparison between Aguiar’s results and the present work with constant thermal conductivity is presented in Fig.2. Range of parameters involved in this modeling and default values of geometry and operating condition of fuel cell are given in table 3.
The simulated results of current-voltage characteristic are shown in Fig.3. As it can be seen, the cathode activation overpotential has a significant magnitude in comparison with other overpotentials. Anode activation overpotential and ohmic loss are the next reasons causing the voltage drop.
Fig.4 shows the magnitude of overpotential in fuel cell. The vertical axis shows overpotential and horizontal axis shows three operating temperatures: 973, 1023 and 1073. At 973K, activation overpotential accounted for 72.5% and concentration overpotential accounted for only 3% of the total of losses as mentioned in the legend. Also, these results show that the amount of over potential decreases by increasing the work temperature. At 973K, overpotential contribution will be increased to 9.3%.
Table2. Data for exchange current density equation
|
k |
E |
Cathode |
||
Anode |
Fig.2. Comparison of Aguiar’s result and present work for (a) Mole fraction of fuel components and (b) PEN temperature |
Fig.3. Cell voltage and power density as a function of current density at1023K for a totally reformed fuel mixture. |
Table3. the operating and geometrical parameters used in the mathematical model |
|||
Geometrical parameter |
Value |
Operating parameter |
Value |
Cell length (m) |
0.4 |
Air and fuel input temperature (˚K) |
1024 |
Cell width (m) |
0.1 |
Air and fuel input pressure (bar) |
1 |
Fuel and air channel height (mm) |
1 |
Average current density (A/cm2) |
0.5 |
Anode thickness (µm) |
500 |
air ratio |
8.5 |
cathode thickness (µm) |
50 |
Fuel utilization |
75% |
electrolyte thickness (µm) |
50 |
steam / carbon |
2. |
interconnect thickness (µm) |
500 |
pre reforming rate |
10% |
3.1. Gas composition and temperature distribution
Fig.5 illustrates temperature curves of each layer along fuel cell length. The Cell temperature falls steadily in the initial part of a fuel cell as a result of the sharply endothermic reforming reaction, as expected. In Fig.6, distributions of reforming and shift reaction rates are shown. It can be seen that by reducing the influence of reforming reaction, the temperature will raise over the cell.
Fig.7 shows the mole fraction distribution in the fuel channel. The mole fraction of water and methane at the beginning of the cell are prominent because the steam reforming reaction is dominant. As a result of reforming, the concentration of hydrogen and monoxide carbon will rise gradually and the concentration of methane and water will be reduced. To put it more simply, 80 percent of the input methane will be consumed in 22% of the elementary input cell and hydrogen molar fraction will be more than doubled. But due to the electrochemical reaction, the hydrogen concentration near the outlet reduces and the concentration of water will be increased.
3.2. Effect of air ratio
The Air ratio is one of the most effective parameters for controlling the cell temperature and preventing the occurrence of a temperature gradient. The Inlet air flow rate is adjusted to keep maximum temperature difference across the fuel cell lower than 150K [15]. Values of the maximum PEN temperature difference in the air ratio between 6.5 and 12.5 are presented in Table 5.
Figs. 8 and 9 show the effect of this parameter on the temperature and local current density distribution across the cell. Increasing of the air ratio will decrease the temperature difference throughout the channel and current density distribution becomes more uniform. But it should
Fig. 4. The Magnitude of overpotential in the fuel cell at three operating temperatures: 973, 1023 and 1073. |
Table4. The comparison results obtained with Aguiar [6] |
||
Parameter |
Aguiar [2] |
Present |
Output voltage (V) |
0.66 |
0.65 |
PEN temperature (˚K) min max |
980 1134 |
985 1138 |
Fuel Channel temperature (˚K) min max |
970 1127 |
981 1135 |
Interconnect temperature (˚K) min max |
976 1120 |
983 1131 |
Local current density (A/m2) min max |
2385 9385 |
2751 8517 |
be considered that at high pressure, increase in the work of the compressor and the cycle efficiency will eventually come down. Hence it is necessary to choose it carefully.
Fig.5. Fuel and air channels, PEN structure, and interconnect temperature along the cell length for the conditions in Table 2. |
Fig.6. The distributions of thereforming and shift reaction rates |
Table 5. PEN temperature difference for SOFC model of Figs 2 and 3
Air ratio |
temperature difference(˚K) |
6.5 |
177.8 |
8.5 |
154.1 |
10.5 |
127.2 |
12.5 |
110.1 |
Fig.8. The Effect of air ratio on the PEN temperature curves. |
Fig.9. The Effect of air ratio on the local current density. |
3.3. Effect of fuel utilization factor
Figs 10 and 11 illustrate the effect of changing fuel utilization. Fuel utilization factor is the ratio of fuel consumed to produce electricity to the total amount of inlet fuel. Fig 8 presents the PEN temperature profile along the cell. It indicates that temperature will rise gradually with increase of the utilization factor. As a result of increasing temperature, the activation loss (that causes major losses in the fuel cell) will be reduced. Thus, the operating voltage and subsequently power density will be increased.
In general, by increasing this parameter power density will be increased, but the excessive increase of fuel utilization factor causes reduction in the amount of hydrogen in the terminal of cell and consequently the concentration over potential will be increased. Also in high fuel utilization, the temperature difference is higher due to unbalanced current density distribution in the cell.
4. Conclusions
In this study, an anode-supported high-temperature DIR-SOFC fed by 10% pre-reformed methane is studied. This model considers the mass transport and heat transfer, over potentials and ohmic loss which the heat transfer process is based on the exothermic and endothermic reactions occurring in fuel channel. The model can effectively predict the distribution of each layer temperature curve, current density and gas molar fraction along the channels.
A parametric analysis was performed to investigate the effects of some operational parameters on temperature difference across the cell. The investigations indicate that the temperature difference grows from 101K to 145K by increasing the fuel utilization factor from 0.45 to 0.85 and deduces from 177K to 124K by increasing the air ratio from 6.5 to 12.5. Because of the fact that the cell power density can dwindled by increasing the anode activation loss which itselfis temperature dependent, the maximum power density by considering the maximum allowable cell temperature difference is considered as objective and the values of the fuel utilization factor and air ratio are obtained 0.85 and 9, respectively.
Fig.10. The Effect of fuel utilization factor on the PEN temperature curves. |
Fig.11. The Effect of fuel utilization factor on the operating voltage and open-circuit potential. |
Nomenclature |
||
specific heat capacity |
(kJ kg−1 K−1) |
|
hydraulic diameter of a rectangular cross section gas channel |
(m) |
|
bulk diffusivity coefficient |
(m2 s-1) |
|
effective diffusivity coefficient |
(m2 s-1) |
|
molecular diffusivity coefficient |
(m2 s-1) |
|
Knudsen diffusivity coefficient |
(m2 s-1) |
|
Potential |
(V) |
|
open-circuit potential |
(V) |
|
activation energy of the exchange current density |
(kJ mol−1) |
|
Faraday’s constant |
(C mol−1) |
|
air channel height |
(m) |
|
fuel channel height |
(m) |
|
heat transfer coefficient for cathode |
(kJm−2 s−1 K−1) |
|
heat transfer coefficient for anode |
(kJm−2 s−1 K−1) |
|
local current density |
(Am−2) |
|
average current density |
(Am−2) |
|
exchange current density |
(Am−2) |
|
thermal conductivity |
(kJm−1 s−1 K−1) |
|
pre-exponential factor of the exchange current density |
(Am−2) |
|
number of electrons participating in the electrochemical reaction |
|
|
Nusselt number |
|
|
pressure |
(Pa) |
|
partial pressure of component i |
(Pa) |
|
partial pressure of component i at three phase boundary |
(Pa) |
|
rate of reaction i |
(molm−2 s−1) |
|
gas constant |
(kJ mol−1 K−1) |
|
total cell resistance, including both ionic and electronic resistances |
(Ωm2) |
|
temperature |
(K) |
|
air velocity |
(ms−1) |
|
fuel velocity |
(ms−1) |
|
system width |
(m) |
|
transfer coefficient |
|
|
enthalpy change of the electrochemical reaction |
(kJ mol−1) |
|
electrode porosity |
|
|
Stefan–Boltzmann constant |
(Wm−2 K−4) |
|
electronic conductivity of the anode |
(Ω−1 m−1) |
|
electronic conductivity of the cathode |
(Ω−1 m−1) |
|
ionic conductivity of the electrolyte |
(Ω−1 m−1) |
|
Density |
(kgm−3) |
|
thickness |
(m) |
|
Tortuosity |
|
References
[1] Achenbach E. Three-dimensional and time-dependent simulation of a planar solid oxide fuelcell stack. J Power Sources 1994;49:333–48.
[2] Ferguson JR, Fiard JM, Herbin R. Three-dimensional numerical simulation for variousgeometries of solid oxide fuel cells. J Power Sources 1996;58:109–22.
[3] Iwata M, Hikosaka T, Morita M, Iwanari T, Ito K, Onda K, et al. Performance analysis of planar-type unit SOFC considering current and temperature distributions. Solid State Ionics. 2000;132:297-308.
[4] Yakabe H, Ogiwara T, Hishinuma M, Yasuda I. 3-D model calculation for planar SOFC. J Power Sources 2001;102:144–54.
[5] Aguiar.P, Chadwick.D, Kershenbaum.L, Modeling of an indirect internal reforming solid oxide fuel cell, Chem. Eng. Sci. 57 (2002)1665–1677.
[6] Aguiar P, Adjiman CS, Brandon NP. Anode-supported high temperature direct internal reforming solid oxide fuel cell. I: modelbased steady-state performance. J Power Sources 2004;138:120–36.
[7] Aguiar P, Adjiman CS, Brandon NP. Anode-supported high temperature direct internal reforming solid oxide fuel cell: II. Modelbased dynamic performance and control. JPower Sources 2005;147: 136–47.
[8] Colpan CO, Dincer I, Hamdullahpur F. Thermodynamic modeling of direct internalreforming solid oxide fuel cells operating with syngas. J Hydrogen Energy 2007;32,787-795.
[9] Meng Ni, Dennis Y.C. Leung, Michael K.H. Leung. Modeling of methane fed solid oxide fuel cells: Comparison between proton conducting electrolyte and oxygen ion conducting electrolyte J Power Sources 183 (2008) 133–142
[10] Dong HyupJeon A comprehensive CFD model of anode-supported solid oxide fuel cells .ElectrochimicaActa 54 (2009) 2727–2736
[11] Iwai H, Yamamoto Y, Saito M, Yoshida H. Numerical simulation of hightemperaturedirect-internal-reforming planar solid oxide fuel cell. Energy 2011:1–10.
[12] Singhal SC, Kendall K, editors. High temperature solid oxide fuel cells: fundamentals, design and applications. Oxford: Elsevier; 2003. p. 12–7.
[13] M.W. Chase, NIST-JANAF thermochemical tables, fourth edition, J. Phys. Chem. Ref. Data, Monograph 9 (1998) 1–1951.
[14] Meng Ni, Dennis YC Leung, Michael KH Leung. Electrochemical modeling and parametric study of methane fed solid oxide fuel cells. J Energy Conversion and Management;2009.50.268- 278.
[15] Bove, R. and S. Ubertini, Modeling Solid Oxide Fuel Cells: Methods, Procedures and Techniques. 2008: Springer.