Computational Analysis of Solid Oxide Fuel Cell with Functionally Graded Electrodes

Ceramic high-temperature fuel cells and electrolysers are efficient energy conversion systems for electrical power generation and hydrogen production. Their core components are constituted by a stack of electroactive Solid Oxide Cells (SOCs) in which the electrochemical reactions take place. Thanks to their flexibility, the same stack can be alternatively operated in both fuel cell and electrolysis modes. However, the insufficient durability of SOCs still constitutes a major limitation for the technology. The present study addresses this issue and aims to bring some new insights on the effect of the Solid Oxide Fuel Cell (SOFC) versus Solid Oxide Electrolysis Cell (SOEC) operating modes on degradation of a typical Ni-YSZ//YSZ//CGO//LSCF-CGO cell. The electrochemical degradations are generally attributed to several underlying phenomena such as electrode microstructural evolution, material chemical decomposition or electroactive sites poisoning by contaminants. Among them, it is generally considered that Ni agglomeration in the Ni-YSZ cermet and Lanthanum Strontium Cobalt Ferrite (LSCF) material destabilization are two prevalent mechanisms involved in the cell performance deterioration. Therefore, these two mechanisms have been specifically investigated by a coupled approach of long-term testing in both SOFC and SOEC modes (1000 ≤ t (h) ≤9000) and post-test characterizations. The experimental results have then been analyzed in the frame of an in-house multi-scale model with the purpose to interpret them and to quantify the effect of material ageing on cell performances. The extent of Ni agglomeration has been characterized by three-dimensional electrode reconstructions obtained by X-ray nano-holotomography at European Synchrotron Radiation Facility (ESRF) (on the new Nano-Imaging beamline ID16A-NI). The new set-up and protocol enable the reconstructions of valuable 3D volumes (Fig. 1) with a large field of view (~50 µm) along with a high spatial resolution (~50 nm). The electrode morphological properties, which have been measured on the 3D volumes, have revealed a substantial Ni coarsening over time at 850°C and 750°C, whereas no Ni depletion was detected at the electrolyte interface. The increase of the Ni particle size is found to induce a decrease in both (i) the density of Triple Phase Boundary (TPBs) lines and (ii) the interfacial surface area between Ni and gas. Moreover, it was found that the Ni/YSZ interfacial surface area does not evolve during the experiments. This statement indicates that the ceramic backbone in the cermet prevents a massive Ni agglomeration at the SOFC/SOEC operating temperature. The compilation of all experimental data have allowed fitting the parameters of a physically-based law for Ni coarsening that was introduced in the modelling framework. The simulations have revealed that Ni agglomeration explains around 20-25% of the electrochemical degradation at 850°C after 1000 hrs of operation. However, the electrode microstructural evolution is found not to be affected by the cell polarizations. Therefore, the mechanism cannot explain the higher degradation rates recorded in electrolysis mode compared to the fuel cell ones. To explain the impact of the operating modes (SOFC or SOEC) on the degradation rates, several post-test analyses (i.e. Scanning Electron Microscopy, Transmission Electron Microscopy, X-ray µfluorescence and µdiffraction techniques) have been employed to investigate the phase reactivity in the region of the CGO barrier layer. The characterizations have revealed that Sr diffusion across the barrier layer and formation of SrZrO3 secondary phase occur mainly during electrolysis operation, whereas the process is very limited in fuel cell mode (Fig. 2). As a consequence, LSCF destabilization is found not to be involved in the degradation of cell performances during fuel cell operation while it could explain the highest degradation rates recorded in electrolysis mode. The post-test analyses have also revealed a diffusion and an accumulation of Co in the region of the barrier layer which is concomitant with the formation of SrZrO3. The formation of these Co-rich segregates in contact with SrZrO3 grains have been identified as cobalt-ferrite type compound. The in-house multi-scale model has been used to interpret the role of the cell operating mode on the LSCF destabilization mechanism. The cell polarization curves and the local quantities within the O2 electrode have been computed in both fuel cell and electrolysis modes. The simulations have shown that the electrolysis operation leads to a strong depletion of oxygen vacancies in the LSCF material. It has been proposed that the depletion in oxygen vacancies under electrolysis polarization could drive the Sr release from the structure, and in turn, could explain the experimental results. Based on this proposition, a possible mechanism for the LSCF destabilization and SrZrO3 formation has been detailed. Figure 1

Solid oxide fuel cells (SOFCs) attract attentions as high-efficiency energy-conversion devices. In order to establish high-performance SOFCs, the slow electrode kinetics, which is often a main cause of the energy loss, should be improved. In present SOFCs, mixed ionic-electronic conducting (MIEC) oxides, such as (La,Sr)CoO3 -δ and (La,Sr)(Co,Fe)O3 -δ , are typically used as the cathode materials. In such an MIEC cathode, it is generally believed that electrochemical oxygen reduction reaction takes place mainly at the oxide surfaces (double phase boundary, DPB) and the contribution of the reaction at the electrode/electrolyte/gas interface (triple phase boundary, TPB) is insignificant. Therefore, researches on the kinetics of SOFC cathodic reaction have mostly concerned to investigate the surface reaction on the electrode oxides. In our previous works, the reaction distributions in MIEC SOFC electrodes were experimentally evaluated by using operando micro X-ray absorption spectroscopy (operando μ-XAS) [1, 2]. For more quantitative and precise investigation of the reaction distribution, we proposed a patterned thin film electrode, which is a kind of the columnar electrode simplifying the microstructures of a porous electrode [2]. However, the patterned thin film electrode in our previous work did not have triple phase boundaries, and it was thus impossible to discuss the contribution of TPB reaction to the total reaction. In this work, a patterned thin film electrode having triple phase boundaries was fabricated. (La,Sr)CoO3- d (LSC) was chosen as the model MIEC oxide for SOFC cathodes. The reaction distribution in this new type of the patterned thin film electrode was experimentally determined by using operando μ-XAS, and the results were compared with ones obtained with the old type of the patterned thin film electrode without TPB. It was found the effective reaction area tended to shrink by the introduction of TPB to the patterned thin film electrode. This fact might suggest that the reaction at TPBs in SOFC cathodes is not negligible even when MIEC oxides are used as electrode materials.

A micro scale model of a solid oxide fuel cell (SOFC) involving the mass transfer together with the electrochemical reaction, the electron and ion transports through respective cylindrical shaped electron-and ion-conducting particles inside the electrodes was mathematically developed. The predicted cell performance was showed according to the operating and design condition. The effects of micro-scale electrode geometry on the cell performance were also taken into account. This present study reveals the working mechanisms of SOFC at the micro-scale level, while demonstrating the use of micro-scale relations to enhance the SOFC performance. The accuracy of the presented model was validated by comparing to already existing experimental results from the available literatures.

The functional layer of a flat‐tubular solid oxide fuel cell (SOFC) is examined using a three‐dimensional microscale electrode model. SOFC electrodes essentially include two types of layers: a structural layer and a functional layer. The structural layers, which are the anode support layer and the cathode current collector layer, are composed of large particles with a high porosity that facilitates gas diffusion. The functional layers consist of small particles with a low porosity that increases the triple phase boundary (TPB) reaction area and reduces the activation overpotential. In the model, the particle diameter and functional layer thickness are adjusted and analyzed. The effects of the two parameters on the performance of the functional layer are monitored in the contexts of several multilateral approaches. Most reactions occurred near the electrode–electrolyte interface; however, an electrode design that included additional TPB areas improved the electrode performance. The role of the functional layer in a flat‐tubular SOFC is examined as a function of the functional layer particle size and thickness. The performance of a cell could be enhanced by preparing a functional layer using particles of optimal size and thickness, and by operating the device under conditions optimized for these parameters.

In order to develop a predictive model of real cell performance, firm relationships and assumptions need to be established for the definition of the physical and microstructure parameters for solid oxide fuel cells (SOFCs). This study explores the correlations of microstructure parameters from a microscale level, together with mass transfer and electrochemical reactions inside the electrodes, providing a novel approach to predict SOFC performance numerically. Based on the physical connections and interactions of microstructure parameters, two submodel correlations (i.e., porosity–tortuosity and porosity–particle size ratio) are proposed. Three experiments from literature are selected to facilitate the validation of the numerical results with experimental data. In addition, a sensitivity analysis is performed to investigate the impact of the adopted submodel correlations to the SOFC performance predictions. Normally, the microstructural inputs in the numerical model need to be measured by experiments in order to test the real cell performance. By adopting the two submodel correlations, the simulation can be performed without obtaining relatively hard-to-measure microstructural parameters such as tortuosity and particle size, yet still accurately mimicking a real-world well-structured SOFC operation. By accurately and rationally predicting the microstructural parameters, this study can eventually help to aid the experimental and optimization study of SOFC.

In the present study, numerical simulation of a solid oxide fuel cell (SOFC) button cell with a multi-step oxygen reduction reaction (ORR) mechanism in a composite cathode was developed for better understanding of processes within the SOFC. The ORR mechanism involves both surface charge transfer through triple phase boundary (3PB) and bulk charge transfer via a two phase boundary (2PB) pathway. The developed multi-physics simulation considered the charge conservation and species transportation through porous media. The developed simulation was first applied to a dense thin film study only involving surface adsorption and the 2PB pathway. By calibrating the simulation with measurements of dense thin film cathodes, the kinetic rates of rate-determining steps in ORR model were decided. Finally, the simulations of full SOFC button cell were calibrated with measured cell performance (polarization curves and impedance behavior) at three different air/fuel utilization conditions. It was found that it is beneficial to determine the kinetic rates of rate-determining steps in a multi-step charge transfer model before calibrating the simulations for full cell system. With the determined rates, the various performance of different cells can be attributed to the microstructure versus material variations. Therefore, the present study not only helps better understand the details of processes in SOFC but also provides guidance for electrode design and manufacturing.

Solid Oxide Fuel Cells (SOFCs) are energy conversion devices known for their significantly higher power density compared to other fuel cell types. However, their high operating temperatures pose challenges related to thermal stability. To address this, research is focusing on Low-Temperature SOFCs (LT-SOFCs), which function at lower temperatures and exhibit enhanced electrochemical performance. While various electrode materials are utilized in SOFCs, platinum (Pt) stands out for its excellent electronic conductivity and catalytic activity. Unfortunately, at the operating temperatures of SOFCs, Pt tends to agglomerate, leading to a rapid reduction in the triple phase boundary (TPB) and a subsequent decline in electrochemical reactions. In this study, LT-SOFCs were fabricated with a Praseodymium Oxide (PrOx) capping layer applied to a porous Pt cathode using sputtering, with various thicknesses achieved by adjusting the deposition time. The electrochemical performance of the LT-SOFCs was measured at 500oC. Additionally, the degradation behavior of the LT-SOFCs was assessed by applying a constant voltage of 0.5 V for 48 hours. Scanning Electron Microscopy (SEM) analysis was also conducted on the PrOx capping layer thin films under the same operating conditions.

The nickel(Ni)-patterned electrode enables researchers to quantify the triple-phase boundary (TPB) length and Ni surface area, exclude the influence of bulk gas diffusion and clearly separate the active regions of the chemical/electrochemical reactions. This tool has been studied for identifying the reaction mechanism in SOFC but rarely in SOEC. In this study, the Ni-patterned electrodes with the stripe width of 100 μm were tested in both the solid oxide fuel cell (SOFC) and solid oxide electrolysis cell (SOEC) modes at the atmosphere of H2O/H2 and further H2O/CO2/H2/CO. The experimental test shows the stability of the Ni-patterned electrode in the H2O/H2 atmosphere was much poorer than that in the CO/CO2 atmosphere due to the nickel reacting with steam. Thus, the Ni-patterned electrode needed to operate at a temperature of ≤700oC with the inlet H2/H2O molar ratio of >7 and the current-applied time shouldn’t be over 5 h to keep the accuracy of the TPB length and Ni surface area. The effects of the temperature, partial pressure of H2O and H2 are investigated. The activation energy of the Ni-patterned electrodes was 0.386 eV. The electrochemical performance had a positive correlation with both the partial pressure of H2 and H2O, especially showed higher sensitivity to the partial pressure of H2O. Further, a possible mechanism of the H2O/H2 electrochemical conversion was proposed according to the experimental data and existing literature, which contains two-step charge-transfer reaction: H(Ni)+O2-(YSZ)OH-(YSZ)+(Ni)+e- and H(Ni)+OH-(YSZ)H2O(YSZ)+(Ni)+e-. Based on the experiments, an analytical calculation was performed to investigate the rate-determining steps in the Ni-patterned electrodes by bridging the connection between the kinetic parameters and operating conditions. The results indicate the rate-limiting steps may be different for the SOFC and SOEC modes. In SOFC mode, the H2 electrochemical oxidation could depend more on the charge-transfer reaction H(Ni)+O2-(YSZ)→OH-(YSZ)+(Ni)+e-, however in SOEC mode, the rate-determining step of the H2O electrochemical reduction should be H2O(YSZ)+(Ni)+e-→OH-(YSZ)+H(Ni). In this way, a reversible mechanism could be applied for reversible SOFC (RSOFC) to describe the conversion between hydrogen and steam by different charge-transfer-reaction domination in the SOFC and SOEC mode. Figure 1

YSZ ( yttria-stabilized zirconia)-modified LSM composite cathodes were fabricated by formation of YSZ film on triple phase boundary (TPB) of LSM/YSZ/gas. The YSZ coating film greatly enlarged electrochemical reaction sites from the increase of additional TPB. The composite cathode was formed on thin YSZ electrolyte (about 30 Um thickness) supported on an anode and then I-V characterization and AC impedance analyses were performed at temperature between . As results of the impedance analysis on the cell at with humidified hydrogen as the fuel and air as the oxidant, R1 around the frequency of 1000 Hz represents the anode Polarization. R2 around the frequency of 100Hz indicates the cathode polarization, and R3 below the frequency of 10 Hz is the resistance of gas phase diffusion through the anode. The cell with the composite cathode produced power density of at air and oxygen atmosphere, respectively. The I-V curve could be divided into two parts showing distinctive behavior. At low current density region (part I) the performance decreased steeply and at high current density region (part II) the performance decreased gradually. At the part I the performance decrease was especially resulted from the large cathode polarization, while at the part H the performance decrease related to the electrolyte polarization.

Implementation of nano-catalyst materials into solid oxide fuel cell (SOFC) electrodes to improve performance and stability has been widely studied. Addition of the nano-catalysts into an electrode structure serves to enhance the electrochemical performance of the SOFC by increasing the Triple Phase Boundary (TPB) area, improving redox stabilization, and modifying reaction kinetics of hydrocarbon gases that cause anode degradation due to carbon deposition.Typical Ni-based cermet anodes suffer from anode deactivation due to carbon build up under hydrocarbon flows. Carbon builds up and covers the TPB area which results in poor electrochemical performance. Larger carbon deposits can block pores within the anode microstructure which leads to gas diffusion issues. Carbon buildup also causes mechanical stresses to the electrode due to volumetric changes which can lead to fracture and complete failure of the cell.This work studied nano-catalyst decoration of the Ni-based cermet anodes with catalysts that promote internal reforming to protect against coking. The addition of active metal components and ceramic reforming promoters were investigated. Multi-component systems with several catalysts were also examined. Uniform incorporation of nano-catalyst into the anode microstructure was achieved through a patented liquid phase surfactant assisted process (using various catechol surfactants). Deposition loading densities and distribution of nanoparticles was controlled by altering the surfactant and catalyst solution concentrations.This work defined initial morphology and studied the coarsening/annealing kinetics of nano-catalysts used for enhanced solid oxide fuel cell (SOFC) operation. Nano-catalyst depositions were performed on flat single crystal surfaces using a bio-inspired surface modification technique. 3D profiles of the nanoparticles at various coarsening stages were obtained using atomic force microscopy (AFM). X-ray photoelectron spectroscopy (XPS) was used to track the changes in composition. Figure 1 shows AFM imaging of single component cerium oxide and cobalt oxide depositions on flat single crystal yttrium stabilized zirconia (YSZ) before and after 24 h of annealing in reducing conditions.Impregnated SOFCs were tested for 200+ h in constant load while operating under high methane fuel streams. Fuel compositions were selected to have high hydrocarbon content to accelerate degradation. Deposition techniques were optimized for multi-component internal reforming catalysts to obtain sustained performance enhancement. Current-voltage-power (I-V-P) measurements and electrochemical impedance spectroscopy (EIS) results were used for evaluation. Post-mortem microstructure and chemistry characterizations were used in analysis.Impregnated SOFCs demonstrated greater than 50% sustained improvement in harsh environments in long-term tests where the anode was subjected to 40% CH4 for 200+ h. Acknowledgements: Work supported by a subcontract under the GE Aerospace’s US DOE- ARPA-E project (DE-AR0001344) entitled “FueL CelL Embedded ENgine (FLyCLEEN)”. Material characterization and imaging work was made possible with the support of the West Virginia University Shared Research Facilities. Figure 1

Solid Oxide Fuel Cells (SOFCs) are widely considered to be an ideal power source across a range of industries in the near future due to their high efficiency, high power density, and fuel flexibility. However, for widespread usage of SOFCs to be feasible, certain disadvantages, such as their high operating temperature, electrochemical losses, and high cost of manufacturing, must be solved. In particular, the field of materials research occupies a unique position of being able to quickly and simultaneously solve many of these issues through the development of material compositions and structures that do not currently exist in the industry. Especially relevant to solid oxide fuel cells in particular, is the impact that the electrodes’ geometries and compositions have on the operating performance. The electrodes’ ability to quickly move reactants to the interface of the electrolyte and remove products to the bulk flow comes into play as the fuel cell operates at higher current densities and concentration losses become prominent. Therefore, the microstructures present throughout the electrodes’ depth plays a vital role in improving overall cell performance. However, current fuel cell research has not yet effectively captured and defined the connections of the microstructure parameters and measured their impact on cell behavior. Additionally, homogeneous electrode design and implementation have dominated the research and commercial space thus far. The technique of functional material grading, which has begun to see wide research and commercial usage for such problems as acoustic impedance matching, thermal control, and mechanical design will be leveraged for the enhancement of SOFC electrodes. Functionally graded electrodes have been used in solid oxide fuel cell research in recent years in an effort to improve the cell performance by altering the microstructure including porosity, particle size, and composition of electronic and ionic conductors near the triple phase boundary region. Normally the functional grading process adds additional fabrication steps making it less desirable as an optimization process and easy manufacturing method. However, by using an additive manufacturing process it eliminates the need for multiple discrete graded layers to be deposited to obtain a graded electrode functional layer. This study seeks to explore the correlations of geometry and structure parameters from the meso-scale through to the micro-scale level, together with mass transfer, ionic and electronic transport, and gas-surface electrochemical reactions inside the electrodes to guide future additive manufacturing of SOFCs. Outlined is the development of an effective medium model for simulating the performance of fuel cell microstructures as a function of specific operating ranges in temperature, pressure, and in fuels used. Smooth and continuous linear and nonlinear functional gradation profiles of porosity and material composition are studied within the framework of an effective medium boundary value problem where the electrochemical relations, such as the Butler-Volmer equation, percolating conduction paths, and gas diffusion are solved using the finite difference method. The form of these relations for homogeneous electrodes was analyzed by Costamagna, et. al [1], who also defined outstanding problems with regards to optimizing SOFC electrodes [2]. Globally optimal linear and nonlinear gradation profiles and cell parameters are derived as a function of particle sizes, ratio of conductivities, and desired operating conditions and scored based on their reduction of cell overpotentials.[1] Costamagna, P., P. Costa, and V. Antonucci, Micro-modelling of solid oxide fuel cell electrodes. Electrochimica Acta, 1998. 43(3-4), pp.375-394.[2] Costamagna, P., P. Costa, and E. Arato, Some more considerations on the optimization of cermet solid oxide fuel cell electrodes. Electrochimica Acta, 1998. 43(8), pp.967-972. Figure 1

This paper describes the fabrication and investigation of morphologically stable model electrode structures with well-defined and sharp platinum/yttria-stabilized zirconia (YSZ) interfaces to study geometric effects at triple phase boundaries (TPBs). A nanosphere patterning technique using monodispersed silica nanoparticles, which are applied to the YSZ surface by the Langmuir–Blodgett method, is employed to deposit nonporous platinum electrodes containing close-packed arrays of circular openings 300–400 nm in diameter through which the underlying YSZ surface is exposed to the gas phase. These nanostructured dense Pt array cathodes exhibited better structural integrity and thermal stability at the solid oxide fuel cell (SOFC) operating temperature of when compared to porous sputtered Pt electrodes. More importantly, electrochemical studies on geometrically well-defined Pt/YSZ sharp interfaces demonstrated that the cathode impedance and cell performance both scale almost linearly with the aerial density of TPB length. These controlled experiments also demonstrated that when normalized with respect to TPB length, the performance of different cells with different TBP densities agree well each other, indicating that TPB length governs cell performance especially in the activation polarization regime, as expected. Cells with a higher TPB density achieved better fuel cell performance in terms of higher power density and lower electrode impedance.

Gadolinium doped ceria (GDC) gains increasing attention as a promising material for both solid oxide fuel cells (SOFCs) and solid oxide electrolysis cells (SOECs). The GDC-based electrodes demonstrate positive features, e.g. enhanced electrochemical performance in low operation temperatures, resistivity to sulfur poisoning and carbon deposition, etc. As GDC is a mixed electronic-ionic conductor in reducing atmosphere, electrochemical reactions can be supported by the double phase boundary (DPB) at the GDC surface. Several studies demonstrated that the DPB reaction may dominate the electrochemical performance over triple phase boundary (TPB) reaction in Ni-GDC electrodes, where it can account for as much as 95 % of the total reaction depending on the cermet microstructure [1]. Recent studies even show that pure ceria electrodes can achieve substantial power density [2,3]. In this context, the microstructure design principle for GDC-based electrodes should focus on maximizing GDC surface area. Particularly, GDC electrodes with infiltrated GDC nanoparticles or electro-catalytic nanoparticles may achieve competitive performance compared with the conventional Ni-YSZ electrodes. Although many strategies have been proposed to improve performance of GDC-based electrodes, it is difficult to quantitatively compare all the reported results. Furthermore, the correlation between microstructure and performance of pure GDC electrode are still scarce.Several different types of GDC-based electrodes are fabricated and tested in the unified characterization procedure in this study. Four basic designs are considered, i.e. Ni-GDC, GDC-perovskite composites, GDC electrodes with deposited Ni nanoparticles and pure GDC electrodes. Several strategies to modify the GDC surface area are discussed. In addition, the substantial enhancements in the DPB density and electrochemical performance were achieved by precipitating Ni-nanoparticles directly on sub-micron GDC powders at 2wt.% loading. The performance of electrodes was evaluated by electrochemical impedance spectra measurements with various humidity and temperature conditions. The results showed clear dependance between GDC surface area and electrode performance. At the same time, it should be noted that the electrode active thickness depends on the electrode microstructure. In particular, it is in the range of 5 µm for the pure nano-GDC electrode. When the GDC electrode is thicker than the active thickness, it results in large Ohmic resistance due to the electronic conduction in the porous electrode.[1] A. Nenning, M. Holzmann, J. Fleig, A.K. Opitz, Mater. Adv. 2 (2021) 5422–5431.[2] M. Ouyang, A. Bertei, S.J. Cooper, Y. Wu, P. Boldrin, X. Liu, M. Kishimoto, H. Wang, M. Naylor Marlow, J. Chen, X. Chen, Y. Xia, B. Wu, N.P. Brandon, J. Energy Chemistry. 56 (2021) 98–112.[3] W. Jung, K.L. Gu, Y. Choi, S.M. Haile, Energy Env. Sci. 7 (2014) 1685–1692.

Solid oxide fuel cells (SOFCs) are getting more importance with their promising future, due to their high-energy conversion efficiency, low-emissions and flexibility of usable fuel type. The performance and the lifetime of SOFCs are keenly dependent on electrode microstructure. In order to recognize the microstructural evolution and its degradation kinetics in SOFC cermet cathodes during long exposure time (up to 20 000 h) under realistic operating conditions ( T= 850 °C, J= 190-250 ), investigations on porous LSM/YSZ cathodes were conducted. The 3D-tomography technique (FIB/SEM) offers extensive data about the microstructures of various cathode aged during different operating times (2 500 h, 15 000 h and 20 000 h), which allows an exact quantification of particle size distribution, phase-connectivity, tortuosity factor and Triple Phase Boundary Length (TPBL). With the increasing operating time no significant 3D microstructural changes in the cathode were noticed in the obtained data. However, additional qualitative X-ray fluorescence measurement, indicate a clear presence of chrome contamination on aged cathodes, which may be the main degradation mechanism in the SOFC cathode. Keywords: chrome poisoning, long operating time, particle size distribution, SOFC, triple phase boundary, SEM, XRF

Integrated gasification fuel cell (IGFC) systems combining coal gasification and solid oxide fuel cells (SOFC) are promising for highly efficient and environmentally friendly utilization of coal for power production. Most IGFC system analyses performed to-date have used nondimensional thermodynamic SOFC models that do not resolve the intrinsic constraints of SOFC operation. In this work a quasi-two-dimensional (2D) finite volume model for planar SOFC is developed and verified using literature data. Special attention is paid to making the model capable of supporting recent SOFC technology improvements, including the use of anode-supported configurations, metallic interconnects, and reduced polarization losses. Activation polarization parameters previously used for high temperature electrolyte-supported SOFC result in cell performance that is much poorer than that observed for modern intermediate temperature anode-supported configurations; thus, a sensitivity analysis was conducted to identify appropriate parameters for modern SOFC modeling. Model results are shown for SOFC operation on humidified H2 and CH4 containing syngas, under coflow and counterflow configurations; detailed internal profiles of species mole fractions, temperature, current density, and electrochemical performance are obtained. The effects of performance, fuel composition, and flow configuration of SOFC performance and thermal profiles are evaluated, and the implications of these results for system design and analysis are discussed. The model can be implemented not only as a stand-alone SOFC analysis tool, but also a subroutine that can communicate and cooperate with chemical flow sheet software seamlessly for convenient IGFC system analysis.