Exchange current, Butler-Volmer model

Charge transport is modeled by Ohm s law (Equation (3.10)) and the charge conservation equation (Equation (3.68)), while the current density distribution at the electrode/electrolyte interface is modeled through the Butler-Volmer equation (Equation (3.102)). It should be noted that, contrarily to Section 3.7, Equation (3.102) is here derived from Equation (3.37) rather than Equation (3.39), because the former allows for a better agreement between experimental and simulated results. Equations (3.40)-(3.42) are used to model, the exchange current density, the activation overpotential, and the ideal potential drop at the electrode/electrolyte interface, respectively. Heat transfer is modeled through Equation (3.6), and the appropriate heat terms for each domain. 

Additional parameters specified in the numerical model include the electrode exchange current densities and several gap electrical contact resistances. These quantities were determined empirically by comparing FLUENT predictions with stack performance data. The FLUENT model uses the electrode exchange current densities to quantify the magnitude of the activation overpotentials via a Butler-Volmer equation [1], A radiation heat transfer boundary condition was applied around the periphery of the model to simulate the thermal conditions of our experimental stack, situated in a high-temperature electrically heated radiant furnace. The edges ofthe numerical model are treated as a small surface in a large enclosure with an effective emissivity of 1.0, subjected to a radiant temperature of 1 103 K, equal to the gas-inlet temperatures. 

Eor very high exchange current densities (i.e., rapid reactions), a linearized form of Eq. 27 can often be used. For very slow reaction kinetics, either the anodic or cathodic term dominates the kinetics, and so the other term is often ignored, yielded what is known as a Tafel equation for the kinetics. Often, more complicated expressions than that of Eq. 27 are used. For example, if the elementary reaction steps are known, one can write down the individual steps and derive the concentration dependence of the exchange current density and the kinetic equation. Other examples include accounting for surface species adsorption or additional internal or external mass transfer to the reaction site [9]. All of these additional issues are beyond the scope of this chapter, and often an empirically based Butler-Volmer equation is used for modeling the charge transfer in porous electrodes. 

Electrode Kinetic and Mass Transfer for Fuel Cell Reactions For the reaction occurring inside a porous three-dimensional catalyst layer, a thin-film flooded agglomerate model has been developed [149, 150] to describe the potential-current behavior as a function of reaction kinetics and reactant diffusion. For simplicity, if the kinetic parameters, such as flie exchange current density and diffusion limiting current density, can be defined as apparent parameters, the corresponding Butler-Volmer and mass diffusion relationships can be obtained [134]. For an H2/air (O2) fuel cell, considering bofli the electrode kinetic and the mass transfer, the i-rj relationships of the fuel cell electrode reactions within flie catalyst layer can be expressed as Equations 1.130 and 1.131, respectively, based on Equation 1.122. The i-rj relationship of the catalyzed cathode reaction wifliin the catalyst layer is... 

Regarding catalyst aging. Darling and Meyers have proposed a mechanistic model, based on empirical parameters, of the Pt oxidation/dissolution in a PEMFC cathode, largely cited by experimentalists in subsequent papers. By using classical Butler-Volmer equations written in terms of the CL electrode potential and empirical parameters (e.g. symmetry factors, zero-exchange current and... 

The rate of anodic hydrogen oxidation is proportional to the current density of hydrogen oxidation (iH2)- As discussed in Section 11.3, this current density according to the electrochemical model follows a Butler-Volmer-type equation (Eqs. (10a) and (10b)) with a concentration-dependent exchange current... 

---