This paper presents a new computational framework for modeling chemically reacting flow in anode-supported solid-oxide fuel cells (SOFC). Depending on materials and operating conditions, SOFC anodes afford a possibility for internal reforming or catalytic partial oxidation of hydrocarbon fuels. An important new element of the model is the capability to represent elementary heterogeneous chemical kinetics in the form of multistep reaction mechanisms. Porous-media transport in the electrodes is represented with a dusty-gas model. Charge-transfer chemistry is represented in a modified Butler-Volmer setting that is derived from elementary reactions, but assuming a single rate-limiting step. The model is discussed in terms of systems with defined flow channels and planar membrane-electrode assemblies. However, the underlying theory is independent of the particular geometry. Examples are given to illustrate the model.

Monolithic three-way catalysts are applied to reduce the emission of combustion engines. The design of such a catalytic converter is a complex process involving the optimization of different physical and chemical parameters (in the simplest case, e.g., length, cell densities or metal coverage of the catalyst). Numerical simulation can be used as an effective tool for the investigation of the catalytic properties of a catalytic converter and for the prediction of the performance of the catalyst. To attain this goal, a two-dimensional flow-field description is coupled with a detailed surface reaction model (gas-phase reactions can be neglected in three-way catalysts). This surface reaction mechanism (with C₃H₆ taken as representative of unburnt hydrocarbons) was developed using sub-mechanisms recently developed for hydrogen, carbon monoxide and methane oxidation, literature values for C₃H₆ oxidation, and estimates for the remaining unknown reactions. Results of the simulation of a monolithic single channel are used to validate the surface reaction mechanism. The performance of the catalyst was simulated under lean, nearly stoichiometric and rich conditions. For these characteristic conditions, the oxidation of propene and carbon monoxide and the reduction of NO on a typical Pt/Rh coated three-way catalyst were simulated as a function of temperature. The numerically predicted conversion data are compared with experimentally measured data. The simulation further reveals the coupling between chemical reactions and transport processes within the monolithic channel.