PNNL

Abstract

A two-dimensional numerical model is presented for the efficient computation of the steady-state current density, species concentration, and temperature distributions in planar solid oxide fuel cell stacks. The model reduction techniques, engineering approximations, and numerical procedures used to simulate the stack physics while maintaining adequate computational speed are discussed. The results of the model for benchmark cases with and without on-cell methane reformation are presented with comparisons to results from other research described in the literature. The capabilities, performance, and scalability of the model for the study of large multi-cell stacks are then demonstrated.