Research Papers

Reaction Kinetics at the Triple-Phase Boundary in PEM Fuel Cells

[+] Author and Article Information
P. Berg

Faculty of Science, UOIT, 2000 Simcoe Street N, Oshawa, ON, L1H 7K4, Canadapeter.berg@uoit.ca

A. Novruzi

Department of Mathematics and Statistics, University of Ottawa, 585 King Edward Avenue, Ottawa, ON, K1N 6N5, Canadanovruzi@uottawa.ca

O. Volkov

Department of Mathematics and Statistics, University of Ottawa, 585 King Edward Avenue, Ottawa, ON, K1N 6N5, Canadaovolkov@uottawa.ca

J. Fuel Cell Sci. Technol 5(2), 021007 (Apr 11, 2008) (10 pages) doi:10.1115/1.2821599 History: Received July 03, 2007; Revised July 10, 2007; Published April 11, 2008

A mathematical model is developed to describe the reaction dynamics in the vicinity of the triple-phase boundary (TPB), which is an important part of the pore scale structure of the catalyst layer in proton exchange membrane fuel cells. The model incorporates coupled diffusion, migration, and reaction phenomena of the chemical components in an undersaturated air pore and ionomer. One challenging feature of the work is the description of the TPB by a system of nonlinear partial differential equations (PDEs), coupling bulk, and surface-diffusion phenomena, which offers an approach to study the rarely investigated proton surface diffusion along the air pore surface. A numerical technique is implemented, taking into account the particular form of the domain, in order to solve the nonlinear PDE system efficiently. Several numerical results are discussed, including a sensitivity analysis with respect to the physical reference case and geometric parameters. The results indicate that surface diffusion might play a major role for the reaction kinetics, but only if the air pore is void of liquid water. In contrast, the formation of liquid water in the gas pores will turn surface diffusion into bulk diffusion, with the latter resembling the Grotthus mechanism.

Copyright © 2008 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

Schematic picture of a C∕Pt particle in the CL and geometric modeling of the triple-phase contact point

Grahic Jump Location
Figure 2

Diffusion fluxes considered in this model

Grahic Jump Location
Figure 3

Surface diffusion

Grahic Jump Location
Figure 4

Reference case solution. Ω1 and Ω2 each occupy a right angle sector. Solutions of oxygen, vapor, water and protons are given in a dimensionless form, the electric potential in Volt.

Grahic Jump Location
Figure 5

Dimensionless electrochemical reaction term S¯c (see Eq. 17) along the TPB Γ01∪Γ02 for different current densities i0ref. The coordinate −1 corresponds to the boundary Γ1, and the coordinate +1 corresponds to the boundary Γ2. Triangles and circles correspond to high and low current densities, respectively.

Grahic Jump Location
Figure 6

Dimensionless reaction term c¯Oc¯+exp(αcφ¯) along the ionomer-C∕Pt interface. Comparison between the (oxygen-only) model of O’Hayre (5) and the results calculated with the full model presented in this article for different values of the kinetic reaction parameter k0.

Grahic Jump Location
Figure 7

Current as a function of different reference values

Grahic Jump Location
Figure 8

Current as a function of the ionomer fraction angle α

Grahic Jump Location
Figure 9

Water sorption isotherm on the air pore/ionomer interface Γ12 and in the air pore Ω1

Grahic Jump Location
Figure 10

Flux balance at the origin




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In