Polymer electrolyte membrane (PEM) fuel cells convert the chemical energy of hydrogen and oxygen directly into electrical energy. Waste heat and water are the reaction by-products, making PEM fuel cells a promising zero-emission power source for transportation and stationary co-generation applications. In this study. a mathematical model of a PEM fuel cell stack is formulated. The distributions of the pressure and mass flow rate for the fuel and oxidant streams in the stack are determined with a hydraulic network analysis. Using these distributions as operating conditions, the performance of each cell in the stack is determined with a mathematical, single cell model that has been developed previously. The stack model has been applied to PEM fuel cell stacks with two common stack configurations: the U and Z stack design. The former is designed such that the reactant streams enter and exit the stack on the same end, while the latter has reactant streams entering and exiting on opposite ends. The stack analysed consists of 50 individual active cells with fully humidified H-2 or reformate as fuel and humidified O-2 or air as the oxidant. It is found that the average voltage of the cells in the stack is lower than the voltage of the cell operating individually, and this difference in the cell performance is significantly larger for reformate/air reactants when compared to the H-2/O-2 reactants. It is observed that the performance degradation for cells operating within a stack results from the unequal distribution of reactant mass flow among the cells in the stack. It is shown that strategies for performance improvement rely on obtaining a uniform reactant distribution within the stack, and include increasing stack manifold size, decreasing the number of gas flow channels per bipolar plate, and judicially varying the resistance to mass flow in the gas flow channels from cell to cell. Copyright (C) 2004 John Wiley Sons, Ltd.