Carbon dioxide reduction in solid oxide electrolyzer cells utilizing nickel bimetallic alloys infiltrated into Gd0.1Ce0.9O1.95 (GDC10) scaffolds
Not Available
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Overview -------------------------- This repository contains datasets from the manuscript **"Enhanced Generalizability to Deep-Learning Quantification of 3D Microstructural Characteristics through Microstructurally Aware Augmentation of Scarce Data"** (*William F. Kent, Rochan Bajpai, Rachel C. Kurchin, William K. Epting, Harry W. Abernathy, Paul A. Salvador. Submitted 2026*). The methods are also described in the dissertation **Data Intensive Analysis of Solid Oxide Cell Microstructures** (*Doctoral dissertation, Carnegie Mellon University, 2025*). The datasets here are trained convolutional neural network (CNN) models for predicting key microstructural properties of solid oxide cell (SOC) electrodes from low-res, 2-channel 3D images, as well as some helpful code. The parameters for input images are provided in the paper. Sample data is provided in the file `Combined_anode_aug_dual_1k_examples` - that particular data was used to train `anode_all_aug.pth` and will work most accurately with that model. Please familiarize yourself with all caveats on accuracy and applicability, as detailed in the associated paper. Usage -------------------------- The basic usage is as follows, assuming `model_fn` is the path to the .pth file, and `X` is 2-channel input image(s) of the proper dimensions (either one image of shape `[2,12,24,24]`, or a batch of N input images of shape `[N,2,12,24,24]`): from CNN_inferencer import load_model_for_inference model = load_model_for_inference(model_fn) y_predicted = model(X) The model object automatically handles input scaling and output de-scaling based on the way the models were trained - in other words, pass in a 2-channel micro-CT image, and it will output microstructural property values in real units. ## Other model object attributes Note that model has useful attributes other than its forward pass model(X). * `model.output_descaler` - returns the output descaler object. Model does the de-scaling when generating inferences, but you may want to re-use this de-scaler on other values to e.g. compare predictions to ground truth from already-scaled training data. * `model.prop_names` - Gives the property names of the predicted y values, in order. Only exists if there's an output scaler as part of the model object, which there will be in the models provided here. ## Usage with sample data Here is a short script to use with the included sample data. from CNN_inferencer import display_predictions, load_model_for_inference, calculate_mape, parity_plot import h5py import numpy as np model_fn = 'anode_all_aug.pth' data_fn = 'Combined_anode_aug_dual_1k_examples.h5' N_samples = 200 figure_outdir = '.' model = load_model_for_inference(model_fn) with h5py.File(data_fn,'r') as f: XX = f['X'] #These are the 2-channel 3D images yy = f['y'] #These are the ground-truth microstructural properties, but they have been scaled for training - need to de-scale below N = XX.shape[0] #How many images total in the input data file #Run inferences on N_samples random samples from XX. #Run in a batch, much more efficient than one at a time. ii = np.random.choice(N,N_samples,replace=False) ii.sort() y_pred = model(XX[ii]) #Get the original/true (but normalized/scaled) values from the training dataset... #Because they were normalized, they are not in real units yet. So let's also de-scale them using model.output_scaler. y_true = model.output_scaler.transform(yy[ii]) #Let's display actual values for just 5 random ones for i in np.random.choice(N_samples,5,replace=False): display_predictions(y_true[i], y_pred[i], model.prop_names) #Make parity plots for each property (ground truth vs predicted values) #Also label each plot with the mean abs. percent error (MAPE) of the predicted values for i,key in enumerate(model.prop_names): mape = calculate_mape(y_true[:,i], y_pred[:,i]) parity_plot(y_true[:,i], y_pred[:,i], figure_outdir, key, extra_title=f' ({mape:.2f}% MAPE)')
System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.
The ability to recover the oxygen reduction reaction of poisoned metal oxide surfaces, central to many energy related applications, is demonstrated by controlling relative surface acidity.
Polarization modeling based on the measured I-V curves is a useful tool for analyzing SOFC performance and diagnosing performance limitations. Using the polarization model, the overpotential of a cell is usually separated into different polarization contributions, including area specific ohmic polarization, activation polarization, and anodic and cathodic concentration polarizations. In this study, we aim to experimentally validate the accuracy of the polarization model. As a case study, anode-supported cells with LSM-YSZ/YSZ/Ni-YSZ configuration are used. By curve-fitting the I-V curves of a cell tested under different cathode oxygen partial pressures into the polarization model, the contributions of different polarization losses are quantified. To validate the polarization model, independent experiments/studies including electrochemical impedance spectroscopy, cathode kinetics study, ex-situ diffusivity measurement, and tortuosity measurement using 3D reconstructed anode are conducted. The results from polarization modeling and from validation studies are in good agreement, thereby validating the polarization model for SOFC performance analysis.
An individual solid oxide cell (SOC) constructed of a sandwich configuration including in the following order: an oxygen electrode, a solid oxide electrolyte, a fuel electrode, a fuel manifold, and at least one layer of mesh. In one embodiment, the mesh supports a reforming catalyst resulting in a solid oxide fuel cell (SOFC) having a reformer embedded therein. The reformer-modified SOFC functions internally to steam reform or partially oxidize a gaseous hydrocarbon, e.g. methane, to a gaseous reformate of hydrogen and carbon monoxide, which is converted in the SOC to water, carbon dioxide, or a mixture thereof, and an electrical current. In another embodiment, an electrical insulator is disposed between the fuel manifold and the mesh resulting in a solid oxide electrolysis cell (SOEC), which functions to electrolyze water and/or carbon dioxide.
Although the individual effects of airborne Cr and S contaminants on SOFC cathode performance degradation have been extensively studied, the combined effects of Cr and S contaminants remain largely unexplored. Under the real SOFC operating condition where the Cr and S species coexist, their effects may compete, affecting the poisoning behavior. Here, our investigation reveals that the combined Cr and S poisoning behavior of LSCF remains different from those of individual Cr and S poisonings, while the combined poisoning mechanism of LSM is equivalent to the sum of those of individual Cr and S effects. For LSCF electrode, gaseous Cr species are deposited mainly at the LSCF/GDC interface by electrochemical reduction, rather than forming SrCrO 4 on LSCF surfaces (as in the case of Cr-only poisoning), indicating no reaction between Cr vapors and SrO on the LSCF surface. Thermodynamic analysis demonstrates that the SrO on LSCF surface absorbs SO 2 (g) and thereby loses the Cr-gettering effect, allowing Cr vapors to flow through LSCF and to reach the LSCF/GDC interface where the Cr deposition occurs. Unlike LSCF, LSM electrode shows cumulative effects of Cr and S, as Cr accumulation occurs at the triple-phase boundary and S absorption takes place at localized Sr-rich regions.
The chromium (Cr) evaporation behavior of several different types of iron (Fe)-based AFA alloys and benchmark Cr 2 O 3 -forming Fe-based 310 and Ni-based 625 alloys was investigated for 500 h exposures at 800 °C to 900 °C in air with 10% H 2 O. The Cr evaporation rates from alumina-forming austenitic (AFA) alloys were ~5 to 35 times lower than that of the Cr 2 O 3 -forming alloys depending on alloy and temperature. The Cr evaporation behavior was correlated with extensive characterization of the chemistry and microstructure of the oxide scales, which also revealed a degree of quartz tube Si contamination during the test. As a result, long-term oxidation kinetics were also assessed at 800 to 1000 °C for up to 10,000 h in air with 10% H 2 O to provide further guidance for SOFC BOP component alloy selection.
While Cr poisoning of the oxygen reduction reaction (ORR) at SOFC cathodes is widely agreed to involve deactivation of oxygen exchange sites, the degradation mechanism remains ambiguous. Here, we selected an alternative cathode material Pr 0.1 Ce 0.9 O 2− δ , free of Sr segregation, to systematically investigate the effect of Cr-induced degradation in ORR. We expand on our previous studies in which the acidity/basicity of binary additives was found to be a strong indicator of the rate of oxygen surface exchange, by electrochemically investigating the ORR activity of the PCO cathode by impedance spectroscopy. Serial infiltration with acidic Cr-based oxides was found to degrade ORR activity as reflected in a 20-fold increase in area specific resistance (ASR) without corresponding changes in activation energy, with the opposite trend obtained with Ca-based oxides. Detected changes in total capacitance, attributed to changes in surface capacitance, also suggest depressed/enhanced PCO surface redox behavior with Cr/Ca-based infiltrants. Taken together, these results point to a more universal source of ORR poisoning/activation, based on acidity/basicity, rather than physical blocking of active sites. With this improved understanding, one can expect progress to be made in the coming years in optimizing means for protecting catalytic surfaces from degradation and/or improving their performance.
This paper addresses the use of Ce 0.8 Gd 0.2 O 2–δ (GDC) infiltration into the Ni–(Y 2 O 3 ) 0.08 (ZrO 2 ) 0.92 (YSZ) fuel electrode of solid oxide cells (SOCs) for improving their electrochemical performance in fuel cell and electrolysis operation. Although doped ceria infiltration into Ni–YSZ has recently been shown to improve the electrode performance and stability, the mechanisms defining how GDC impacts electrochemical characteristics are not fully delineated. Furthermore, the electrochemical characteristics have not yet been determined over the full range of conditions normally encountered in fuel cell and electrolysis operation. Here we present a study of both symmetric and full cells aimed at understanding the electrochemical mechanisms of GDC-modified Ni–YSZ over a wide range of fuel compositions and temperatures. Single-step GDC infiltration at an appropriate loading substantially reduced the polarization resistance of Ni–YSZ electrodes in electrolyte-supported cells, as measured using electrochemical impedance spectroscopy (EIS) at various temperatures (600–800 °C) in a range of H 2 O–H 2 mixtures (3–90 vol% H 2 O). Fuel-electrode-supported cells had significant concentration polarization due to the thick Ni–YSZ supports. Here, a distribution of relaxation times approach is used to develop a physically-based electrochemical model; the results show that GDC reduces the reaction resistance associated with three-phase boundaries, but also appears to improve oxygen transport in the electrode. Increasing the H 2 O fraction in the H 2 –H 2 O fuel mixture reduced both the three-phase boundary resistance and the gas diffusion resistance for Ni–YSZ; with GDC infiltration, the electrode resistance showed less variation with fuel composition. GDC infiltration improved the performance of fuel-electrode-supported full cells, which yielded a maximum power density of 2.28 W cm –2 in fuel cell mode and an electrolysis current density at 1.3 V of 2.22 A cm –2 , both at 800°C.
Oxygen electrodes are provided, comprising a perovskite compound having Formula (I), Sr(Ti 1-x Fe x-y Co y )O 3-δ wherein 0.90 ≥ x ≥ 0.40 and 0.02 ≥ y ≥ 0.30. Electrochemical devices comprising such oxygen electrodes are also provided, comprising a counter electrode in electrical communication with the oxygen electrode, and a solid oxide electrolyte between the oxygen electrode and the counter electrode. Methods of using such electrochemical devices are also provided, comprising exposing the oxygen electrode to a fluid comprising O 2 under conditions to induce the reaction O 2 +4e − →2O 2− , or to a fluid comprising O 2− under conditions to induce the reaction 2O 2− →O 2 +4e − .
Not Available
Not Available