Engineering Papers⌕ Search

Engineering topics

Uberuaga, Blas P.

Publications and source records attributed to Uberuaga, Blas P..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Interstitialcy-mediated transport in Hematite

Defect transport mechanisms drive the performance of materials in both energy storage and energy generating systems. Iron oxide systems are of particular interest with use cases in both realms. Previous studies of a-Fe2O3 (hematite) determined that the primary point defects at elevated temperature are O vacancies and Fe interstitials, and these defects determine transport rates at temperatures >850°C. However, in advanced nuclear reactor designs, operating temperatures will range between 400-800°C. In this study, we utilize isotopic tracers and the nanoscale resolution of atom probe tomography to determine diffusivities at lower temperatures relevant to nuclear energy generating systems. Supported by a density functional theory informed rate theory model, we find a change from vacancy to interstitialcy-mediated transport < 800°C in the anion lattice.

Yano, Kayla H.↗
DOE OSTI2024

The role of defect charge, crystal chemistry, and crystal structure on positron lifetimes of vacancies in oxides

Density functional theory based positron lifetime (PL) calculations for cation and oxygen monovacancies in a range of oxides—hematite, magnetite, hercynite, and alumina—have been conducted to compare the impact of defect chemistry and crystal structure on the predicted lifetimes. The role of defect charge state has also been examined. A comparison across the same type of crystalline structure but different composition shows that oxygen vacancies only induce a slight increase in the positron-electron overlap and thus barely modify the PL as compared to the bulk. A much more substantial increase of PL is observed for cation monovacancies, regardless of crystal structure or the elemental nature of the vacancy, which we ascribe to an enhanced localization of charge density around the vacant site. The structural and compositional richness of the oxide leads to longer defect PLs, with defected hercynite exhibiting the longest PLs. The charge state of cation monovacancies modifies only by a small percentage the positron localization, relegating to secondary importance the metal defect's oxidation state in modifying the lifetime of positrons within vacancy traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Directly resolving surface vs. lattice self-diffusion in iron at the nanoscale using in situ atom probe capabilities

Surface self-diffusion studies on metals under elevated reaction conditions are limited, as it is inherently challenging to unambiguously follow atomic transport across highly-reactive surfaces. Here, quantitative and mechanistic insight into thermally induced atomic transport processes in bcc α-iron at the sub-nanometer level was achieved using isotopic tracer techniques coupled with in situ atom probe tomography (APT) capabilities. Specifically, using a reactor directly connected to the APT, needle-shaped specimens fabricated from epitaxial thin films with an embedded 57 Fe tracer layer were annealed in Ar at 500 °C and 350 °C for 1 hour. Furthermore, the tracer was positioned at various depths in the APT specimen by field evaporation, enabling targeted and simultaneous analysis of lattice and surface diffusion. 57 Fe concentration profiles reveal lattice self-diffusion occurs at 500 °C on the order of ~7 – 9 monolayers, while lattice diffusion is not resolvable at 350 °C. Considerable surface transport was, however, observed at both conditions, where atomic transport over the specimen surface led to the formation of a thin (≤1 nm), isotopically-intermixed layer at the surface. Further, the observed isotopic redistributions at 500 °C were convoluted by additional processes occurring in the subsurface, such as atomic intermixing in correlation with lattice diffusion. However, surface diffusion was determined to be the primary transport process at 350 °C and was thereby quantified. Ultimately, these results demonstrate the significance of surface self-diffusion as a short circuit pathway. More broadly, this approach has the potential to provide detailed insight into (self-)diffusion mechanisms across various materials while targeting site-specific reactions under elevated reaction conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Multi–length scale characterization of point defects in thermally oxidized, proton irradiated iron oxides

A key for the success of safe nuclear power generation system is to consider structural materials that are economical, meet mechanical property needs, possess good corrosion resistance, and are radiation tolerant. Nevertheless, fundamental insights that elucidate the details of radiation damage on materials corrosion performance are lacking. This includes the behavior of surface oxides which often regulate corrosion. For example, it is unclear how non-equilibrium point defects, oxide structure, mass transport in oxides, and subsequent oxidation behavior are altered by the radiation. Here, in this work, some of the effects of proton irradiation on the attributes of point defects, iron oxide microstructures, and the physical nature of the oxidation product were correlated with corrosion behavior. Iron oxides, fabricated by thermal oxidation in air at 400°C and 800°C for 1 h, were subjected to 200 keV, 0.03 dpa (displacements per atom) of proton irradiation, and subjected to corrosion reactivity assessment using AC and DC electrochemical methods. Experimental methods that target materials information at different length scales, such as positron annihilation spectroscopy (atomistic), transmission electron microscopy (mesoscopic), and electrochemical methods (macroscopic), were coupled to shed light on the impact of radiation-induced defect modifications and structural changes in oxides on corrosion reactivity which ultimately affects durability in harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2023

Nanoscale mapping of point defect concentrations with $\mathrm{4D-STEM}$

Vacancies are missing atoms in a crystalline material, and occur both at equilibrium (varying with temperature) and out of equilibrium such as when crystalline materials are damaged with radiation or corrosion. While we know of their importance, particularly regarding diffusive mechanisms, it is not straightforward to experimentally measure their concentration or to visualize them directly. Traditionally, measurements such as positron annihilation spectroscopy and to some degree X-ray diffraction can measure average concentrations, but generally lack the ability to visualize or quantify a heterogenous concentration of vacancies that can occur at the level of individual defects and microstructural features. Here, we present a method to map vacancy concentrations and their distribution using local lattice parameter measurements with a high-resolution electron microscope. Our method utilizes a Au thin film as a model to demonstrate the method via four-dimensional scanning transmission electron microscopy (4D-STEM) by correlating the differences between changes in lattice parameter and the volumetric thermal expansion during in situ heating experiments. The vacancy mapping methodology is also applied to non-equilibrium defects accumulated in pure Al via knock-on electron beam irradiation. Our method demonstrates the ability to map point defect concentrations in heterogeneous systems in situ with nanometer spatial resolution. The result is a technique that can provide direct measurements of vacancy concentrations at the level of individual defects in studies of materials in and out of equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interstitial hydrogen enhances the mobility of some grain boundaries in tungsten

Segregation of interstitials at a grain boundary (GB) is known to generally lower its mobility. This phenomenon, called 'solute-drag', has important ramifications on the process of recrystallization and microstructural evolution. Here, we present predictions from molecular dynamics (MD) simulations which demonstrate that interstitial hydrogen in tungsten can in fact increase the mobility of some GBs which exhibit shear coupling. Assuming a disconnection-based mechanism, activation energies and pre-factors for disconnection nucleation are predicted from simulations of shear-coupled motion. In GBs where enhanced mobility is predicted, interstitial H reduces both the activation energy and the pre-factor for disconnection nucleation, thus effectively increasing the mobility. For GBs with diminished mobility, MD predicts that presence of interstitial H reduces the pre-factor and, in some cases, increases the activation energy. The reduction in the activation energy inferred from MD simulations are confirmed by nudged elastic band calculations. Temperature-dependent structural transitions are observed for some GBs, and the effect of interstitial H is found to change with the changes in structure. The effect of interstitial H is predicted to be complex and highly variable, providing some plausible explanations for experimental observations on the recrystallization of tungsten in presence of H-loaded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗