A modeling study of porous composite microstructures for solid oxide fuel cell anodes

High operating temperature affords solid oxide fuel cells (SOFCs) fuel flexibility, favorable kinetics, enhanced transport properties and high quality heat for secondary thermodynamic cycles.1 Elevated temperatures, however, present challenges to ensuring cell durability.2 Coarsening of nickel in the Ni-YSZ anode is one such challenge that can lead to significant performance degradation. The fundamental driving force Δμ for coarsening phenomena is given by Eq. 1.3 The expression, which assumes isotropic interfacial energy, indicates that the chemical potential μ of a curved interface differs from that of a flat interface μ 0 by the factor VmγH, where H is mean curvature, γ interfacial energy and Vm molar volume. Δμ = μ - μ 0 = VmγH (1) Eq. 1 makes evident that variations in interface curvature, i.e., interface shape, directly correspond to variations in chemical potential. Furthermore, it follows from Eq. 1 that microstructural coarsening is fundamentally related to variations in interface shape. The link between interface shape and coarsening propensity is well known and methods have been developed to directly measure interfacial curvature distributions for complex morphologies in 3-D.4,5 In this work, these techniques are applied to 3-D images obtained by x-ray nanotomography to study coarsening in microstructures obtained from Ni-YSZ solid oxide fuel cell anodes aged for varying lengths of time. This work builds on a previous characterization of the same samples.6 The curvature measurements collected suggest that a combination of coarsening and morphological instability, which leads to channel pinch-off, severely diminish the connectivity of the anode’s Ni network and, in addition, lead to a significant loss of both nominal and effective triple-phase boundary (TPB) density. Both the loss of connectivity and TPB density coincide with cell performance degradation that was observed in the aging experiment7 from which the anode samples were taken. References Singhal, S. C.; Kendall, K. High Temperature Solid Oxide Fuel Cells, Elsevier, 2003.Yokokawa, H.; Tu, H.; Iwanschitz, B.; Mai, A. Fundamental mechanisms limiting solid oxide fuel cell durability. J. Power Sources 2008, 400-412.Voorhees, P. W. The Theory of Ostwald Ripening. J. Stat. Phys. 1985, 231-252.Glicksmann, M. E.; Voorhees, P. W. Ostwald Ripening and Relaxation in Dendritic Structures. Metall. Trans. A, 1984, 995-1001.Mendoza, R.; Savin, I.; Thornton, K.; Voorhees, P. W. Topological complexity and the dynamics of coarsening. Nat. Mater., 2004 385-388.Nelson, G. J.; Grew, K. N.; Izzo, J. R. Jr.; Lombardo, J. J.; Harris, W. M.; Faes, A.; Hessler-Wyser, A.; Van herle, J.; Wang, S.; Chu, Y. S.; Virkar, A. V.; Chiu, W. K. S. Three-dimensional microstructural changes in the Ni-YSZ solid oxide fuel cell anode during operation. Acta Mater., 2012 3491-3500. 7. Faes, A.; Hessler-Wyser, A.; Presvytes, D.; Vayenas, C. G.; Van herle, J. Nickel-Zirconia Anode Degradation and Triple Phase Boundary Quantification from Microstructural Analysis. Fuel Cells, 2009 841-851. Acknowledgements Financial support from the National Science Foundation (Award CBET-1134052) and an Energy Frontier Research Center on Science Based Nano-Structure Design and Synthesis of Heterogeneous Functional Materials for Energy Systems (HeteroFoaM Center) funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences (Award DE-SC0001061) is gratefully acknowledged. X-ray nanotomography was performed with Dr. Steve Wang at the Advanced Photon Source supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under contract No. DE-AC02-06CH11357. Ni-YSZ samples were provided by MER Dr. J. Van herle (Ecole Polytechnique Federale de Lausanne, Switzerland).

Solid oxide fuel cells are among the most thermodynamically efficient energy conversion devices. They are critical for the future if we are to meet rising energy demands and increasing selectivity over our future energy mix. Traditionally, solid oxide fuel cell performance has been assessed using bulk parameters, and the microstructure-bulk properties relationship is yet to be well understood. Here, we use time-lapse 3D imaging to probe 3D electrode structure and move towards intelligently designing structures with specific performance attributes. Using 3D imaging allows for characterisation and quantification of parameters that are not well described in 2D, such as triple-phase boundary lengths, tortuosity and phase percolation. In this work, the parameters for an optimised porous nickel-scandia stabilized zirconia anode (Ni-ScSZ) were calculated using an electrochemical model to maximise triple-phase boundary length and density, gas permeation and electronic and ionic transport. Ni-ScSZ scaffolds were then produced using screen-printing and impregnation and characterized with FIB-SEM to obtain 3D microstructural parameters such as actual percolation of the phases, tortuosity and connectivity of the phases, and the length and density of triple-phase boundaries (Figure 1). The experimental microstructural parameters were then compared to the modelled optimised microstructure, allowing for an examination of the accuracy of the manufacturing parameters. Symmetrical cells were set up to compare the anode performance with the performance predicted by the electrochemical model. Early iterations of impregnated scaffolds show triple-phase boundary densities increased by approximately an order of magnitude and allowed for a critical assessment of the electrochemical model. The degraded electrode morphology was also characterized in 3D to link its changes during cycling to performance degradation. This approach enables us to move towards designing optimised electrodes with specific performance attributes. Figure 1

Design of optimal microstructures for infiltrated solid oxide fuel cell (SOFC) electrodes is a complicated process because of the multitude of the electrochemical and physical phenomena taking place in the electrodes in different temperatures, current densities and reactant flow rates. In this study, a stochastic geometric modeling method is used to create a range of digitally realized infiltrated SOFC electrode microstructures to extract their geometry-related electrochemical and physical properties. Triple Phase Boundary (TPB), active surface density of particles along with the gas transport factor is evaluated in those realized models to adapt for various infiltration strategies. Recently, additive manufacturing or freeze type casting methods enable researchers to investigate the performance of directional electrodes to get the maximum electrochemical reaction sites, gas diffusivity and ionic conductivity simultaneously. A series of directional backbones with different amount of virtually deposited electrocatalyst particles are characterized in the first step. The database of microstructural parameters (inputs) and effective geometric properties (outputs) is used to train a range neural network. A microstructure property hull is created using the best neural network model to discover the range of effective properties, their relative behaviour and optimum microstructure. The characteristics of models is shown that there is not any contradiction between the high level of TPB and contact surface density of particles, but the highest amount of gas diffusivity can be found in the microstructures with lower level of reaction sites. Also increasing the contact surface density has a negative effect on gas transport but the high level of TPB density is feasible in the full range of microstructures. In the other hand, TPB density and gas diffusion into the models are inversely related, although there are a limited number of microstructures with high level of reaction sites and acceptable gas diffusivity. Finally, using a simple optimization process, the microstructures with the highest level of reaction sites and gas transport factor are identified which have the backbone porosity of about 50%, and extremely higher gain growth rate normal to the electrolyte. Additive manufacturing and 3D printing methods will enhance researchers in the future to create the real directional electrodes on the base of these proposed models.

Microstructural stability of supported metal catalysts: A phase field approach

3D microstructure-performance relationships in Ni-YSZ anodes for electrolyte-supported cells are investigated in terms of the correlation between the triple phase boundary (TPB) length and polarization resistance (Rpol). Three different Ni-YSZ anodes of varying microstructure are subjected to eight reduction-oxidation (redox) cycles at 950 °C. In general the TPB lengths correlate with anode performance. However, the quantitative results also show that there is no simplistic relationship between TPB and Rpol. The degradation mechanism strongly depends on the initial microstructure. Finer microstructures exhibit lower degradation rates of TPB and Rpol. In fine microstructures, TPB loss is found to be due to Ni coarsening, while in coarse microstructures reduction of active TPB results mainly from loss of YSZ percolation. The latter is attributed to weak bottlenecks associated with lower sintering activity of the coarse YSZ. The coarse anode suffers from complete loss of YSZ connectivity and associated drop of TPBactive by 93%. Surprisingly, this severe microstructure degradation did not lead to electrochemical failure. Mechanistic scenarios are discussed for different anode microstructures. These scenarios are based on a model for coupled charge transfer and transport, which allows using TPB and effective properties as input. The mechanistic scenarios describe the microstructure influence on current distributions, which explains the observed complex relationship between TPB lengths and anode performances. The observed loss of YSZ percolation in the coarse anode is not detrimental because the electrochemical activity is concentrated in a narrow active layer. The anode performance can be predicted reliably if the volume-averaged properties (TPBactive, effective ionic conductivity) are corrected for the so-called short-range effect, which is particularly important in cases with a narrow active layer.

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

Polarization resistance of solid oxide fuel cell electrodes consists of resistances due to mass transfer and chemical reaction. Therefore, higher ion conductivity, electron conductivity, and gas diffusivity and larger triple-phase boundary (TPB) density lead to lower polarization resistance. To increase effective conductivities and diffusivity, tortuosity factors in each mass-transfer path should be reduced. In this study, magnetic alignment method is developed to align the mass-transfer paths in the electrode microstructure. Since Ni as an electron-conducting material in anode has high magnetic permeability, the particles are polarized and aligned by a magnetic field. Therefore, the tortuosity factors can be reduced. In addition, this Ni alignment is expected to reduce the tortuosity factors of yttria-stabilized zirconia (YSZ) and pore phases as well. However, there is a risk for a decrease in TPB density. The effect of magnetic field on resultant microstructure and polarization resistance of Ni-YSZ anode are investigated by impedance measurement of anode symmetrical cells and by reconstruction of three-dimensional microstructure. The tortuosity factors of Ni and YSZ decreased and TPB density decreased by the magnetic field. As a result, the polarization resistance of anode was kept nearly unchanged. The effect of reducing Ni volume fraction is also investigated.

Geometrical modeling of microstructure of solid oxide fuel cell composite electrodes

Effect of powder morphology on the microstructural characteristics of La0.6Sr0.4Co0.2Fe0.8O3 cathode: A Kinetic Monte Carlo investigation

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

3D imaging and quantification of interfaces in SOFC anodes

Numerical study on TPB density and percolation properties of microstructure reconstruction of nickel/yttria stabilized zirconia cermet anode based on discrete element method

Quantitative estimation of triple phase boundaries in solid oxide fuel cell electrodes via artificial neural network

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.

Evaluation of electrochemical reaction mechanisms of La0.6Sr0.4CoO3-δ-Gd0.1Ce0.9O2-δ composite cathodes by 3D numerical simulation