Reversible and Irreversible Potentials and an Inaccuracy in Popular Models in the Fuel Cell Literature

: Modeling is an integral part of fuel cell design and development. This paper identiﬁes a long-standing inaccuracy in the fuel cell modeling literature. Speciﬁcally, it discusses an inexact insertion, in popular models, of cell/stack current into Nernst’s equation in the derivation of output (load) voltage. The origin of the inaccuracy is traced to the nature of reversible and irreversible potentials (equilibrium and non-equilibrium states) in the cell. The signiﬁcance of the inaccuracy is explained in the context of the electrochemistry and thermodynamics of the fuel cell.


Introduction
Modeling the operation of fuel cells is an indispensable part of fuel cell research. In many papers in the literature involving the modeling of solid oxide fuel cells (SOFCs) and proton exchange membrane fuel cells (PEMFCs), an inaccuracy is found in the cell/stack voltage expression (see, for example, references ). We argue that the problem arises because of an inexact intermediate step in the derivation of the output voltage. The origin of this inaccuracy can be traced to the substitution in the Nernst equation of activity (or concentration or partial pressure) with a function of the load current. An analysis of the flaw is presented in the remainder of this paper.

Background: The Nernst Equation
The Nernst equation, which is the cornerstone of fuel cell thermodynamics, provides an expression for the reversible thermodynamic potential, also known as the equilibrium voltage or the open-circuit electromotive force (EMF), of the fuel cell [31]: where E 0 is the reference (standard) EMF at unit activity and atmospheric pressure; i and j are the numbers of reactant and product species; a represents the activity; c i is the stoichiometric coefficient of species i; R is the universal gas constant; F is Faraday's constant; n is the number of electrons transferred for each molecule of the fuel participating in the reaction; and T is the temperature. For a hydrogen-oxygen fuel cell (e.g., solid oxide fuel cell or proton exchange membrane fuel cell), hydrogen and oxygen are the reactants, and the product is water (or steam). The reference EMF, E 0 , depends on the temperature, T: Assuming perfect gas properties for hydrogen at the anode channel, we have where p anode is the partial pressure of hydrogen in the anode channel; n anode H 2 = n anode H 2 (t) is the number of moles of hydrogen in the anode channel; and V anode is the anode volume. Taking the derivative with respect to time, we can express the rate of change of partial pressure in terms of the molar flow rate: with q anode H 2 (t) defined as q anode Splitting q anode H 2 (t) into its three components, namely, the input flow rate, the flow rate that takes part in the reaction, and the output flow rate, we have and thus (omitting the t for simplicity) The process of oxygen inflow-reaction-outflow at the cathode channel can be described by In SOFCs, water vapor (steam) forms at the anode: In PEMFCs, water forms at the cathode. Of course, when liquid water is produced in a fuel cell (as in PEMFCs), we have Assuming that "the molar flow of any gas through the valve is proportional to its partial pressure inside the channel" [1], we have for an SOFC q out and q out at the anode and q out at the cathode, with K H 2 , K H 2 O , and K O 2 being constants. From electrochemistry we know that where I FC = I FC (t) is the cell current. Upon the introduction of a constant [1,13], Equation (20) yields and Now, inserting Equations (17) and (22) into Equation (13) and taking the Laplace transform (we use the notation that the transform of x(t) is X (s)), we get [12] which, upon the substitution of p H 2 (0) = 0 and becomes Similarly, for oxygen, Equations (14), (19) and (23) yield which, with the substitution of p O 2 (0) = 0 and gives Finally, for steam, we have from Equations (15), (18), and (22) where (Equations (26), (29) and (30) were used in Figure 2 of reference [1].) Now, applying the inverse Laplace transform to Equation (26) allows us go from the s-domain back to the time domain: The steady-state partial pressure, then, is obtained from the above equation as (A simpler derivation of Equation (33) is given in Appendix A.1) The steady-state partial pressures of oxygen and steam are obtained similarly: and Now, plugging the above three Equations (33)-(35), as well as the definitions of u and r HO (Equations (A4) and (A5), respectively, from Appendix A.2) into the Nernst EMF expression (Equation (6) of Section 2), we obtain, after some algebra (see Appendix A.2), our familiar Equation (8).

Analysis of the Inaccuracy
The terminal (load) voltage is generally obtained by subtracting the following types of losses (or "irreversibilities") from E Nernst : • activation loss; • concentration loss; • ohmic loss; and • losses due to fuel crossover and internal current.
Equation (8), or one of its many (implicit or explicit) variants, has been the mainstay of a line of fuel cell modeling research for about two decades now (examples abound; a few representative articles are listed in the References section). Not all of the models in references  show the explicit form of Equation (8), but they all use some form of this equation, by expressing-directly or indirectly-the Nernst voltage as a function of, among other variables, current.
The problem with Equation (8) is that it mixes equilibrium and non-equilibrium expressions. The Nernst voltage (EMF) is the reversible thermodynamic potential that applies only to the equilibrium condition of the cell; the equilibrium is lost when current is drawn from the cell. In other words, the Nernst voltage is, by definition, the open-circuit EMF and cannot therefore be expressed in terms of the cell (load) current or current density.
That the use of I FC in Equation (8) is questionable can also be seen from the fact that setting I FC = 0 causes the right side to be mathematically undefined. The effect of Equation (8) on the output voltage can be seen in the SOFC polarization curve (Figure 1) obtained from the following relationship: where V is the output (load) voltage; E Nernst is given by Equation (8); and r is the ohmic resistance of the cell (cell parameter values, taken from references [1,6,12], are given in Table 1). The upper curve (red) in Figure 1 shows the Nernst voltage, E Nernst , computed from Equation (8), while the lower one plots the output voltage, V. Equation (36), like Equation (12) of reference [1] and Equation (6) of reference [6], considers only the ohmic loss out of the four types of losses mentioned earlier.   Figure 1 has the horizontal (current) axis starting at 1 A, not zero, because in Equation (8), the Nernst voltage is undefined at I FC = 0. Figure 1 shows an increase in the Nernst voltage with an increasing load current. At first sight, it may not be immediately clear how the model represented by Equation (8) produces the V-values in this figure (it may not be impossible for the polarization characteristics in Figure 1 to be obtained from actual measurements of a physical fuel cell stack under specific operating conditions; much depends on fuel flow, fuel utilization, thermal effects, diffusion, back-diffusion, load current and voltage requirements, and not all complex interactions are fully understood). It can, however, be argued that Equation (8) produces this polarization curve by holding both u and r HO constant to mimic a constant fuel utilization ratio and a constant hydrogen-oxygen ratio of input flow rates. We need to account for the changes in the V-values without having to make E Nernst a function of current.
For a head-to-head comparison of the polarization produced by Equation (8) with the "correct" polarization, we need the "correct" Nernst voltage (that is, the Nernst voltage at open circuit), which, unfortunately, cannot be obtained from Equation (8). ( (6)

Discussion
If a non-zero current is to be considered, the proper equation to use is the Butler-Volmer equation, not the Nernst equation. The current-voltage relationship at a non-zero current (away from equilibrium) is perhaps best investigated using reaction kinetic expressions with microscopic reversibility. When modeling specific polarization behaviors, it is important to distinguish between the "reversible" potential and the "irreversibilities".
The difference between reversible and irreversible potentials and its connection to the Nernst equation are important in the study of solid oxide electrolysis cells (e.g., [32]), too.
In reference [33], the partial pressures of hydrogen and oxygen of a proton exchange membrane fuel cell are expressed (in Equations (2) and (3) of that paper [33]) in terms of current, but before these partial pressures are plugged into the Nernst equation Equation (1) of that paper [33]) to obtain the reversible potential, the current, correctly, is set to zero. Reference [33] also reported using non-zero values of current in the Nernst equation.
Reference [12] used Equation (8) to generate both training and test data, and therefore, none of the conclusions in that paper are invalidated by the inaccuracy discussed here.

Conclusions
An inaccuracy in the fuel cell modeling literature involving the nature of the reversible thermodynamic potential was identified in this paper. An analysis of the inaccuracy was provided and its significance explained. To summarize, there is nothing wrong with the Nernst equation that provides the reversible (equilibrium) voltage of a fuel cell; what is wrong is the mixing up of equilibrium and non-equilibrium conditions by inserting cell (load) current into the Nernst voltage expression. It can be argued that expressing the Nernst voltage as a function of current is a "quick-and-dirty" trick that makes modeling easy for us. However, by taking such a shortcut, we compromise scientific rigor for the sake of practical convenience.