In this study a steady-state three-dimensional computational fluid dynamics (CFD) model of a proton exchange membrane fuel cell is developed and presented for a single cell. A complete set of conservation equations of mass, momentum, species, energy transport, and charge is considered with proper account of electrochemical kinetics based on Butler-Volmer equation. The catalyst layer structure is considered to be agglomerate. This model enables us to investigate the flow field, current distribution, and cell voltage over the fuel cell which includes the anode and cathode collector plates, gas channels, catalyst layers, gas diffusion layers, and the membrane. The numerical solution is based on a finite-volume method in a single solution domain. In this investigation a CFD code was used as the core solver for the transport equations, while mathematical models for the main physical and electrochemical phenomena were devised into the solver using user-developed subroutines. Three-dimensional results of the flow structure, species concentrations and current distribution are presented for bipolar plates with square cross section of straight flow channels. A polarization curve is obtained for the fuel cell under consideration. A comparison between the polarization curves obtained from the current study and the corresponding available experimental data is presented and a reasonable agreement is obtained. Such CFD model can be used as a tool in the development and optimization of PEM fuel cells. (C) 2007 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.