Development of a 1D Finite-Volume Model for the Simulation of Solid Oxide Fuel Cells

This work presents the development and validation of a 1D finite-volume model for the simulation of planar solid oxide cells (SOCs), developed for integration in more complex systems and process simulations. The model allows to investigate the temperature, composition, and current density profiles along the channel. In this work, the Fick’s equations typically used to calculate the concentration overpotential due to H2 and H2O diffusion in the electrode are improved compared to 1D SOC models available in the literature. In particular, the approximate analytical solution of the dusty gas model (DGM) equations allows for a better definition of H2 and H2O mixture diffusion coefficients, which are relevant, for instance, in the case of solid oxide fuel cells (SOFCs) fed with reformate gas mixtures. Differently from other 1D models available in the literature, the model developed is validated using experimental SOFC polarization curves covering a wide range of operating conditions in terms of molar fraction of H2 (21–93%) and H2O (7–50%) in the fuel, temperature (550–750 °C), and fuel utilization factor (exceeding 90%), demonstrating that 1D SOC models retain a good description of the physical processes occurring within the cell. While this work focuses on a co-flow SOFC configuration, the model can simulate a counter-flow configuration and electrolysis operation without modifying the model equations.