Reduced dimensional computational models of polymer electrolyte membrane fuel cell stacks

  • Authors:
  • Paul Chang;Gwang-Soo Kim;Keith Promislow;Brian Wetton

  • Affiliations:
  • Mathematics Department, 121-1984 Mathematics Road, UBC, Vancouver, BC, Canada V6T 1Z2;Ballard Power Systems, 9000 Glenlyon Parkway Burnaby, BC Canada V5J 5J8;Department of Mathematics, MSU East Lansing, MI 48824, USA;Mathematics Department, 121-1984 Mathematics Road, UBC, Vancouver, BC, Canada V6T 1Z2

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.45

Visualization

Abstract

A model of steady state operation of polymer electrolyte membrane fuel cell (PEMFC) stacks with straight gas channels is presented. The model is based on a decoupling of transport in the down-channel direction from transport in the cross-channel plane. Further, cross-channel transport is approximated heuristically using one-dimensional processes. The model takes into account the consumption of reactants down the channel, the effect of membrane hydration on its conductivity, water crossover through the membrane, the electrochemistry of the oxygen reduction reaction, thermal transport within the membrane electrode assembly (MEA) and bipolar plates to the coolant, heat due to reaction and condensation and membrane resistance, electrical interaction between unit cells due to in-plane currents in the bipolar plates, and thermal coupling of unit cells through shared bipolar plates. The model corresponds to the typical operation with counter-flowing reactant gas streams. The model is a nonstandard system of non-smooth boundary value differential algebraic equations (DAEs) with strong, nonlocal coupling. A discretization of the system and a successful iterative strategy are described. Some preliminary analysis of the system and iterative strategy is given, using simple, constant coefficient, linear versions of the key components of the model. Representative computational results, validation against existing experimental data and a numerical convergence study are shown.