Abstract:Solid oxide fuel cell (SOFC) and electrolyzer cell (SOEC) are the main focus for researchers in recent times to diminish greenhouse gas emissions and meet cleaner energy demand, as these devices have highest energy conversion rates and efficiencies [1]. The major hinderance in commercializing SOEC is due to its faster cell degradation under higher current density, compared against its SOFC counterpart. The delamination of oxygen electrode (OE) from the electrolyte (EL) is one of the major causes for the cell degradation [2]. We employed a multi-physics numerical model to understand the electrochemical performance of a SOEC half-cell more precisely at a microscale. In this study, widely used mixed ionic and electronic conductor (MIEC) La 1-x Sr x Co 1-y Fe y O 3-d (LSCF) and Gadolinium-doped Ceria (GDC) composite is used as the oxygen electrode and, pure ionic conductor 8 mol% Y 2 O 3 doped ZrO 2 (YSZ) is used as the electrolyte. A GDC buffer layer is tucked in between to prevent undesirable side reactions between OE and EL during long term polarization operation. Until now, a homogenous OE structure has commonly been used with effective physicochemical properties to predict the performance of SOFC/SOEC [3,4]. Only a few researchers have worked on heterogenous microstructure model [5] and simulated the electrochemical performance of the cell. Synthetic microstructure of LSCF-GDC porous electrode as shown in Figure 1 are constructed with an opensource software Dream 3D with particle size ranging from 0.45um - 1.5um. We developed in-house Matlab code for labelling domains, interfaces (2PB and 3PB) and digitally removing the unconnected phases to reduce numerical instabilities while solving the mathematical model in COMSOL Multiphysics 6.0. Mesh in these domains is constructed using Iso2Mesh toolbox [6] and translated as mphtxt file which is later imported into COMSOL. This practice reduces the computational time and memory. We employed Nernst-Planck theory for solving the transport of species like vacancy, hole and electron within LSCF inside OE. Multi-step charge transfer at two phase boundary (2PB) and triple phase boundary (3PB) interfaces for oxygen evolution reaction are formulated by Butler-Volmer expressions [4]. A physics-based impedance model under DC bias is developed to separate the electrochemical contributions from 2PB and 3PB transport pathways. The voltage-current curve during long-term polarization test is coupled with the model for validation and approximation of process parameters. Furthermore, we will simulate the crack propagation at OE-buffer layer interface in various directions and orientations to understand the failure mechanism on the synthetically constructed microstructure. References: A. Laguna-Bercero, “Recent advances in high temperature electrolysis using solid oxide fuel cells: a review”, J. Power Sources , vol. 203, no. 4, 2012 N. Rashkeev and M. V. Glazoff, “Atomic-scale mechanisms of oxygen electrode delamination in solid oxide electrolyzer cells”, Int. J. Hydrogen Energy , vol. 37, no. 2, pp. 1280–1291, 2012 Pakalapati, K. Gerdes, H. Finklea, M. Gong, X. Liu and I. Celik, “Micro scale dynamic modeling of LSM/YSZ composite cathodes”, Solid State Ionics, vol. 258, pp. 45-60, 2014 Cook, J. Wrubel, Z. Ma, K. Huang, X. Jin, “Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells”, Journal of The Electrochemical Society, vol. 168, no. 11, pp. 114510, 2021 Su, Z. Zhong and Z. Jiao, “A novel multi-physics coupled heterogenous single-cell numerical model for solid oxide fuel cell based on 3D microstructure reconstructions”, Energy & Environmental Science, vol. 15 pp. 2410 – 2424, 2022 Fang and D. Boas, “Tetrahedral mesh generation from volumetric binary and grayscale images”, Proceedings of IEEE International Symposium on Biomedical Imaging , 2009, pp. 1142–1145 Figure 1

Novel 3-D resistor network simulation method for mixed ionic and electronic conducting electrodes

Multiphysics Modeling for Solid Oxide Electrolyzer Cell with Heterogenous Synthetic Microstructure

Using resistor network models to predict the transport properties of solid-state battery composites

Reconstruction and optimization of LSCF cathode microstructure based on Kinetic Monte Carlo method and Lattice Boltzmann method

Solid Oxide Fuel Cell Porous Cathode Micro Modeling and Performance Prediction

The CO2 electrolysing mechanism in single-phase mixed-conducting cathode of solid oxide cell

Thermo-Electro-Chemo-Mechanical Coupled Modeling of Solid Oxide Fuel Cell with LSCF-GDC Composite Cathode

Oxygen Transport Kinetics of Mixed Ionic-Electronic Conductors by Coupling Focused Ion Beam Tomography and Electrochemical Impedance Spectroscopy

Modeling of cooperative defect transport and thermal mismatch in a planar solid oxide fuel cell

A General Approach for Electrochemical Impedance Spectroscopy Simulation using Transient Mechanistic SOFC Model

A fast regression model for the interpretation of electrochemical impedance spectra of Intermediate Temperature Solid Oxide Fuel Cells

Direct Two-Dimensional Electrochemical Impedance Spectra Simulation for Solid Oxide Fuel Cell

Coupled Electrochemical-Thermo-Stress Analysis for Methanol-Fueled Solid Oxide Fuel Cells

Simulation of Two-Dimensional Electrochemical Impedance Spectra of Solid Oxide Fuel Cells Using Transient Physical Models

Electrode Performance and Analysis of Reversible Solid Oxide Fuel Cells with Proton Conducting Electrolyte of BaCe0.5Zr0.3Y0.2O3−δ

A Novel Equivalent Circuit for GDC-based Solid Oxide Cells Considering the Variations of Cell Resistances under Load

Optimization of interconnect flow channels width in a planar solid oxide fuel cell

Assessing the Versatility and Robustness of Pore Network Modeling to Simulate Redox Flow Battery Electrode Performance

Theoretical Modeling of Air Electrode Operating in SOFC Mode and SOEC Mode: the Effects of Microstructure and Thickness

Simulation of Electrochemical Impedance Spectra of Solid Oxide Fuel Cells Using Transient Physical Models

Simulation of Eis Spectra and Polarization Curves Based on Ni/Ysz Patterned Anode Elementary Reaction Models