In this work, a three-dimensional, steady-state, multi-agglomerate model of cathode catalyst layer in polymer electrolyte membrane (PEM) fuel cells has been developed to assess the activation polarization and the current densities in the cathode catalyst layer. A finite element technique is used for the numerical solution to the model developed. The cathode activation overpotentials, and the membrane and solid phase current densities are calculated for different operating conditions. Three different configurations of agglomerate arrangements are considered, an in-line and two staggered arrangements. All the three arrangements are simulated for typical operating conditions inside the PEM fuel cell in order to investigate the oxygen transport process through the cathode catalyst layer, and its impact on the activation polarization. A comprehensive validation with the well-established two-dimensional "axi-symmetric model" has been performed to validate the three-dimensional numerical model results. Present results show a lowest activation overpotential when the agglomerate arrangement is in-line. Formore realistic scenarios, staggered arrangements, the activation overpotentials are higher due to the slower oxygen transport and lesser passage or void region available around the individual agglomerate. The present study elucidates that the cathode overpotential reduction is possible through the changing of agglomerate arrangements. Hence, the approaches to cathode overpotential reduction through the optimization of agglomerate arrangement will be helpful for the next generation fuel cell design. (c) 2007 Elsevier B.V. All rights reserved.