SEARCH · Search NASA

Results for “Microstructure simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

360 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗

DOE OSTI2025

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗

DOE OSTI2024

Comprehensive new insights on the potential use of SiC as plasma-facing materials in future fusion reactors

Abstract The performance of silicon carbide as an alternative plasma facing material (PFM) was studied at various irradiation conditions relevant to ion energies and fluxes of a fusion reactor. This analysis involves detailed modeling of subsurface plasma/material interactions, sputtered particle transport above the surface and redeposition, and related changes in material composition and microstructure induced by steady-state and Edge Localized Mode ion fluxes. Transition of a crystalline SiC surface to semi-crystalline and amorphous phases was analyzed based on advanced modeling of DIII-D tokamak experiments where SiC was irradiated in single- and multiple- L-mode and H-mode discharges. This analysis shows that displacement damage, particle deposition/redeposition, and D accumulation on the SiC divertor surface can lead to significant microstructural changes that result in enhanced sputtering erosion in comparison with the original crystalline material. However, the resulting total net erosion rate for a full-coverage, advanced tokamak, SiC coated divertor may well be acceptably low. Moreover, the C sputtering yield from the evolved SiC surface can be seven times lower than from a pure graphite surface; this would imply significantly reduced tritium co-deposition rates in a D-T tokamak reactor, compared with a pure carbon surface. It was also determined that chemical sputtering of both C and Si should not result in any noticeable effect on the net erosion, for attached plasma regimes. Our results thus show encouraging results overall for use of SiC as a PFM in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Interpretation of wide-angle x-ray diffraction patterns from mesophase pitch-based carbon fibers – a simulation and experimental study

Carbon fiber is a critical material in a wide range of industries, where it is highly valued for its high specific strength/stiffness, excellent wear resistance, efficient electrical and thermal transport properties, chemical resistance, and low coefficient of thermal expansion. The properties of a specific carbon fiber are closely tied to its structural characteristics at all length scales. Here, in this work, we applied wide-angle x-ray diffraction to a set of heat-treated mesophase pitch-based carbon fibers, with the goal of elucidating the crystalline structures as a function of fiber orientation. To assist with analysis and interpretation of the experimental data, we employed diffraction pattern simulations using the scalar and vector forms of the Debye scattering equation to determine the influence of basal plane orientation, crystalline ordering (turbostratic-graphitic), and basal plane asymmetry on the diffraction patterns. The results presented here suggest that growth of the transverse crystallites in mesophase pitch-based carbon fiber is fixed until graphitization temperatures are reached. The work completed here provides a framework for the analysis of carbon fiber and other oriented carbon-based materials via diffraction.

36 MATERIALS SCIENCE↗

DOE OSTI2025

An integrated approach to examine fuel-cladding chemical interaction in HT9/U-10Zr metallic fast reactor fuels: Coupling machine learning with electron microscopy and local mechanical properties analysis

The metallic U-Zr nuclear fuel alloy has garnered renewed interest as a promising candidate for next-generation sodium-cooled fast reactors. Recent studies and technology assessments have identified several areas requiring improvements, enhanced knowledge, and reliable data to strengthen the U-Zr fuel design basis for qualification and commercial applications. One of the most challenging phenomena impacting this fuel system’s performance is fuel-cladding chemical interaction (FCCI). This work aimed to harvest FCCI data by examining selected HT9/U-10Zr (wt. %) fuel samples of prototypic full-length fuel pins through an integrated approach. This approach integrated scanning electron microscopy (SEM) microstructure characterization with localized mechanical properties examination to deepen understanding of FCCI phenomenon in HT9/U-10Zr fuel system. Particularly, this study focused on MFF fuel pins irradiated at Fast Flux Test Facility (FFTF), which aimed to qualify metallic fuel as a driver fuel for FFTF and to assess its viability for larger-scale fast reactors. Electron microscopy provided high confidence in detecting and distinguishing the different FCCI layers, while small-scale mechanical testing (SSMT) probed the mechanical properties of these layers. SEM examination of a MFF-2 pin 192167, with a time averaged inner cladding temperature (TICT) slightly over 500°C, revealed minimal cladding-side FCCI (cladding wastage). In contrast, significantly thicker cladding wastage comprising two distinct sublayers was observed in samples from the thermally hot MFF-3 pin 193045 and MFF-5 pin 195011 where the TICT ranged from 610-635°C. SSMT indicated complete embrittlement in the sublayer adjacent to the fuel and a tendency toward embrittlement in the other sublayer. Additionally, a new machine learning method was developed, validated, and used to quantify cladding wastage thickness. The machine learning method reliably predicted the wastage thickness across various fuel pins and sample cross-sections. Furthermore, the available cladding wastage data from HT9/U-10Zr fuel system demonstrated a strong temperature dependency. However, the dataset remains small, and ongoing research activities are essential to further understand the FCCI phenomenon and develop a reliable FCCI model for enhanced fuel performance simulation under various conditions.

36 - MATERIALS SCIENCE↗

DOE OSTI2024

Tailoring Piezoelectricity of 3D Printing PVDF-MoS 2 Nanocomposite via In Situ Induced Shear Stress

3D printing of unique structures with tunable properties offers significant advantages in the fabrication of complex and customized electronic devices. Here, this study introduces a process-microstructure–property-guided manufacturing route to fabricate PVDF-2D MoS 2 piezoelectric nanocomposites with tunable piezoelectric properties without having a postprocess. We control PVDF’s microstructure through direct ink writing (DIW) 3D printing while tuning PVDF-MoS 2 interfacial strain by controlling rheology and 3D printing parameters, such as nozzle size and printing speed. Our approach demonstrates tunable piezoelectricity in PVDF-MoS 2 , achieving a 15-fold increase in the piezoelectric coefficient (d 33 ) at a printing-induced shear stress of 6685 Pa. This enhancement arises from the electrostatic interactions between PVDF and MoS 2 and the filler distribution and alignment caused by the in situ shear stress in 3D printing, as confirmed by XPS and Raman mapping analyses. Our findings advance the understanding of piezoelectric mechanisms in PVDF-based nanocomposites, laying the foundation for 3D printing of piezoelectric sensors in wearable device applications with enhanced performance and customization capabilities.

2D MoS2↗

DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗

DOE OSTI2025

Experiment to Investigate the On-set of Transient Fission Gas Release Under Simulated LOCA Conditions

Continuing efforts to better understand the phenomenon of transient fission gas release in irradiated commercial light water reactor fuel undergoing a simulated loss of coolant accident (LOCA) will be discussed in this work. During heating associated with LOCA conditions there is an additional release of fission gas beyond what is typically released during steady state operation. The exact timing of this release could impact the cladding balloon and rupture as the LOCA progresses in an actual accident. Recent experiments performed on the upgraded Severe Accident Test Station at Oak Ridge National Laboratory have expanded the current assessment of transient fission gas release. This system evaluates fission gas release by sweeping fission gas out of the LOCA test furnace and into cryogenically cooled cold traps where radioactive noble gas fission products are measured by gamma spectrometry. Additionally, the sweep gas can be monitored by a gas mass spectrometer to provide additional data on noble gas release and changes in the sweep gas atmosphere during heating. Several different conditions have been evaluated that explored the response of irradiated fuel to different thermal ramp conditions and different overpressure conditions. This further confirmed the impact of an overpressure in suppressing transient fission gas release. Further experiments are planned to evaluate the on-set temperature of transient fission gas release under a known overpressure. This test will progressively step through different temperature conditions while maintaining hydrostatic pressure on the fuel that is representative of high burnup fuel conditions. The available results from this type of testing will be discussed. Post-test microscopy will be compared to other transient fission gas release and standard LOCA tests to further understand the microstructural source of transient fission gas release.

Harp, Jason [ORNL] (ORCID:0000000253458440)↗

DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale molecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

DOE OSTI2024

Roadmap on data-centric materials science

Science is and always has been based on data, but the terms ‘data-centric’ and the ‘4th paradigm’ of materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of artificial intelligence and its subset machine learning, has become pivotal in addressing all these challenges. This Roadmap on Data-Centric Materials Science explores fundamental concepts and methodologies, illustrating diverse applications in electronic-structure theory, soft matter theory, microstructure research, and experimental techniques like photoemission, atom probe tomography, and electron microscopy. While the roadmap delves into specific areas within the broad interdisciplinary field of materials science, the provided examples elucidate key concepts applicable to a wider range of topics. The discussed instances offer insights into addressing the multifaceted challenges encountered in contemporary materials research.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗

DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗

DOE OSTI2025

Residual strain and microstructure evolution in 316L stainless steel fabricated by hybrid additive/subtractive manufacturing

The extreme thermal gradients in additive manufacturing (AM) result in residual strains. Subsequent machining operations can relieve the AM strains but can introduce additional strains. AM and machining on the same platform, also called interleaving, during hybrid manufacturing has been shown to improve productivity. However, the impacts of these operations on the resulting residual strains are not well understood. In this investigation, neutron diffraction was used to measure residual strains in hybrid manufactured 316L stainless steel where blown-powder laser directed energy deposition (DED), and machining operations were interleaved in 1-inch intervals. We observe that the dominant texture developed during AM, resulting from different spot sizes and deposition parameters, significantly affects residual strain evolution. A dominant <011> texture along the build direction does not show any strain localization during interleaved deposition and machining whereas a dominant <001> texture shows strong tensile strain localization consistent with the pauses for machining at the end of AM step. Using coolant during machining can further alter the residual strain evolution. Our work shows that a combination of AM parameters and machining strategy can be used to influence the nature and magnitude of residual strains in components fabricated using convergent manufacturing technologies.

Hybrid manufacturing↗

DOE OSTI2024

Radiation damage study of POCO ZXF-5Q graphite for neutrino production targets using 4.5 MeV helium ions

To address the challenges of increased beam power and target survivability associated with next-generation particle production beam lines, high dose, high-energy proton beam conditions are simulated using irradiation from low-energy ion beams. A low-energy ion irradiation study of POCO ZXF-5Q graphite under conditions similar to those of the NuMI NT-02 neutrino production target at the Fermi National Accelerator Laboratory is reported. Helium ion irradiation was performed at 100 ∘ C to a maximum damage level of 0.9 displacements per atom (DPA). Irradiation induced hardening, swelling of the irradiated region, inter-plane lattice expansion, and intraplane lattice contraction with increasing ion fluence was observed using micromechanical (nanoindentation, atomic force microscopy) and electron microscopy (high-resolution imaging, selected area diffraction) characterization. Similar changes were also observed in post irradiation examination of the NT-02 target indicating that ion irradiation can be a valuable tool for estimating radiation damage in proton beam targets. Caution must be exercised though, because the hardening, lattice alteration, and swelling occur to different magnitudes for a given damage level. The observed hardening and embrittlement were greater for ion irradiated graphite. For He ion irradiated samples the lattice spacing changes were smaller at low damage levels (78% less expansion and 71% less contraction at 0.1 DPA) and larger at high damage levels (38% more expansion and 5% more contraction at 0.9 DPA) relative to that observed in the NT-02 target. The magnitude of swelling was 8.5× greater under ion irradiation which is influenced by the differing damage gradients and inclusion of implanted He ions in the region of interest.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Irradiation-induced formation of G-phase precipitates and M 2 X carbides in self-ion irradiated HT-9

Ferritic-martensitic steels with high chromium content are a promising material group for advanced nuclear systems due to their high temperature strength and good irradiation tolerance. HT-9 is an optimized and often-studied alloy in this group, but additional studies are required on its radiation response under extreme conditions to be experienced in various types of nuclear reactors, especially with respect to phase stability under irradiation. Self-ion irradiation of HT-9 by 5 MeV Fe ions was used to simulate neutron-induced behavior reaching peak doses of 100 and 300 dpa at temperatures ranging from 450 to 550 °C. M 23 C 6 carbides that existed prior to irradiation were found to remain stable under all examined irradiation conditions. As irradiation progressed at 450 and 500 °C, however, formation of spherical-like G-phase precipitates and needle-like M 2 X carbides was observed. G-phase precipitates were found to be enriched in Ni, Si, and Mn, and show no interface segregation, whereas needle-like M 2 X carbides were rich in Cr and Mo and clearly displayed interface segregation of Ni and Si. M 2 X carbide formation is believed to be assisted by vacancies, while G-phase precipitation is thought to be assisted by interstitials. Finally, this difference in defect-mediated formation leads to a difference in distribution with depth. M 2 X carbides are distributed over shallower depths than that of G-phase precipitates, consistent with defect imbalance predictions that consider the influence of the injected interstitial effect.

atom probe tomography↗

DOE OSTI2025

Kinetics and Thermodynamics of Sr Permeation in CeO 2 -Based Barrier Layers for Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) convert steam to hydrogen efficiently at high temperatures. However, during operation, the diffusion of cations or impurities through the cells due to electrode degradation can cause unwanted secondary phases to form, which may degrade device performance. Here, in this study, we use atomistic and mesoscale simulations coupled with experimental analysis to study the diffusion of Sr through the Gd-doped CeO 2 (GDC) barrier layer used to protect the yttria-stabilized zirconia (YSZ) electrolyte in SOECs. From our atomistic calculations, we find Sr diffusion to be negligibly slow in bulk GDC; however, surface diffusion is much more favorable. Subsequent mesoscale simulations show that Sr diffusion is activated when the porosity of GDC exceeds ∼10% and significantly exceeds diffusion in bulk and grain boundary regions. We also find that SrO-based species can accumulate at GDC surfaces; however, SrO aggregation and coarsening will be limited by the large lattice mismatch between GDC and SrO. Energy-dispersive X-ray spectroscopy (EDS) and electron diffraction confirm that Sr can accumulate within GDC pores and form disperse Sr-containing secondary phases. Altogether, Sr diffusion in dense GDC is unlikely to give rise to thick SrO layers, which would severely limit device performance. The formation of Sr-containing secondary phases can largely be avoided by restricting the porosity of the GDC layer as much as possible.

36 MATERIALS SCIENCE↗