The objective of the current work is to develop an electrochemical model based on elementary step charge transfer that occur at the anode and cathode of the HT-PEM fuel cell. In order to develop the mechanistic steps, dissociative adsorption of H_{2} and O_{2} is considered. In addition, to account for performance drop in presence of CO, adsorption of CO is also considered in the elementary step reactions. Assuming the protonation reaction to be rate limiting at the anode three phase interface and OH formation to be rate limiting at the cathode three phase interface, Butler-Volmer expressions are derived for the charge transfer reactions at the respective electrodes.The exchange current density at the anode follows

i=i_{oa }[exp (β_{a}fη_{a}) - exp (-β_{c}fη)],

and at the cathode follows

i=i_{oc }[exp (-β_{c}fηc) - exp ((1+β_{a})fη_{c})].

Here, the exchange current densities at the anode and cathode are respectively given by

i_{0a} =i*x [(K_{H2}p_{H2})^{1/4}]/[1+(K_{CO}p_{CO})+(K_{H2}p_{H2})^{1/2})]

and

ioc = i*x [(p_{H2}/p^{0})^{1/4} (K_{O2}p_{O2})^{3/8}]/[1+(K_{O2}p_{O2})^{1/2}+(p_{H2}/p^{0}) exp (- ΔG*/RT)].

The equilibrium constants appearing the expression for the exchange current densities are evaluated from the adsorption desorption equilibrium of the reactant species (H_{2}, O_{2}, and CO). The adsorption rates are calculated from the sticking coefficients. The activation energy for the desorption reactions is calculated using Unity Bond Index Quadratic Exponential Potential method (UBI-QEP) and the frequency factor for the desorption reactions is calculated from transition state theory.

The electrochemical model is incorporated into a quasi-two dimensional numerical transport model of the cell.

The transport model accounts for the species transport in the flow channels and diffusive transport across the thickness of the electrode (GDL + CL). The diffusive transport in the porous electrode is modeled using Dusty Gas Model (DGM). The resulting governing equations for the species transport are solved using the ODE solver CVode. The electrochemical model, i.e., the Butler-Volmer equations and the cell voltage equation, forms a set of algebraic equations, which are solved using a damped Newton iteration for every call made by the CVode solver.

The whole model is implemented in C++. The model developed can predict the current density and species profiles along the cell length and the thickness of the catalyst layer required to achieve best cell performance. The current density profile predicted by the model along the length is validated by comparing against the experimental measurements reported in the literature. For demonstrating the capability of the model, a comparison between the simulated IV curves and the experimental measurements (Li *et. al, 2003 ECS*) is given in the Figure.