Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

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.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling the aqueous interface of amorphous TiO 2 using deep potential molecular dynamics

Amorphous titanium dioxide (a-TiO 2 ) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO 2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO 2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO 2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO 2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO 2 , leading to an ~10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti 2 –O b H) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti–O w H) due to fast Ti–O w H 2 → Ti–O w H proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO 2 in electrochemical environments. Moreover, the procedure of generating the a-TiO 2 -interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides.

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

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Effect of lithium diffusion into Ga 2 O 3 thin films

The integration of lithium based compounds (e.g., Li:NiO/Ga 2 O 3 , LiGa 5 O 8 /Ga 2 O 3 ) in Ga 2 O 3 based pn-heterojunctions raises concerns about interface stability, i.e., Li diffusion effects on Ga 2 O 3 properties. In this work the ex-situ diffusion of Li is investigated in three different Ga 2 O 3 epilayers systems [(001) κ-Ga 2 O 3 and (−201) β-Ga2O3 heteroepitaxy on (001) α-Al 2 O 3 , and (010) β-Ga 2 O 3 homoepitaxy] at relevant temperatures for the synthesis / processing of Li-based epilayers. It is here experimentally demonstrated and quantified the Li diffusion in all the investigated Ga 2 O 3 epilayers systems and Li bulk (D Li,bulk ) and 2D defects (D Li,2D ) diffusion coefficients are provided. In the case of the (010) β-Ga 2 O 3 homoepitaxial layer (nominally free of structural defects), hybrid functional theory calculations foresee a diffusion mechanism mediated by Ga vacancies (VGa). Moreover, in the (010) β-Ga 2 O 3 homo-layer a significant effect on its functional properties (e.g., additional Raman vibrational modes, induced conductivity in an otherwise insulating sample) upon the Li-diffusion process is experimentally highlighted and tentatively related to the passivation of acceptor defects (i.e., formation of V Ga -nLi complexes).

Defects↗
DOE OSTI2020

Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock

The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Giant interfacial spin Hall angle from Rashba-Edelstein effect revealed by the spin Hall Hanle process

The Rashba-Edelstein effect (REE), which generates interfacial spin polarization and subsequent spin current, is a compelling spin-charge conversion mechanism for spintronic applications, since it’s not limited by the elemental spin-orbit couplings. In this work we demonstrate REE at Pt/ferroelectric interfaces by showing a positive correlation between polarization and effective spin Hall angle in the recently elucidated spin Hall Hanle effects (SHHE), in which a Larmor precession of spin polarization in a diffusion process from the interface manifest as magnetoresistance and Hall effect. We show that REE leads to a large enhancement of the effective spin Hall angle of ferroelectric interface Pt/h-LuFeO 3 compared with that of Pt/Al 2 O 3 , without obvious differences in the spin relaxation time. Modeling using SHHEs involving REE as an additional source of interfacial polarization suggests that REE can lead to an interfacial spin Hall angle (0.4 ± 0.1) in Pt/h-LuFeO 3 that is one order of magnitude larger than the bulk value of Pt. Our results demonstrate that a ferroelectric interface can produce large spin-charge conversion and that SHHEs are a sensitive tool for characterizing interfacial spin-transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Irradiation-enhanced Interactions at UO 2 /Al 2 O 3 /Al Interfaces

Combustion synthesis is used to prepare thin UO 2 films on aluminum alloy substrates. This simple preparation method involves electrospraying uranyl nitrate + acetylacetone + 2-methoxyethanol solution on the substrate, followed by a short annealing at 350 or 550 °C. The irradiation of films with a 40 Ar 2+ ion beam (energy of 1.7 MeV and fluences of 7.7 × 10 16 and 1.3 × 10 17 ions/cm 2 ) is conducted to investigate irradiation-induced restructuring processes. High-resolution transmission electron microscopy (TEM) and X-ray photoelectron spectroscopy (XPS) investigations show that the annealing temperature significantly influences the crystallinity and stability of materials during irradiation. A small amount of Mg in the alloy substrate diffuses into the amorphous Al 2 O 3 interfacial layer between the film and the substrate. Local thermal spikes from the incoming ions facilitate the irradiation-induced mixing of immiscible Al 2 O 3 and UO 2 for the materials prepared at 350 °C. This mass diffusion produces relatively large cavities at the interface. Selective diffusion of a more significant amount of Mg for the materials prepared at 550 °C suppresses the mixing of the Al 2 O 3 interlayer with the film but forms Mg y U 1–y O 2±x solid solutions during irradiation. Local thermal heating triggers the precipitation of a discontinuous crystalline MgO layer close to the film surface. As a result, the enhanced and selective diffusion of Mg into the film makes the materials prepared at 550 °C more robust and mechanically stable than those prepared at 350 °C.

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

Critical length scales for chemical heterogeneity at Cu/Nb 3D interfaces by atom probe tomography

Cu/Nb nanocomposites containing sharp, two-dimensional (2D) interfaces have outstanding strength but limited deformability. In contrast, Cu/Nb with three dimensional (3D) biphase interfaces exhibiting crystallographic, topological, and chemical variations in all spatial dimensions overcomes this limitation by simultaneously enhancing material strength and deformability. Here, while structural characterization of 3D interfaces has been performed to understand their mechanical behavior, three dimensional chemical characterization of such interfaces is lacking. In this work we quantify the local chemistry of 3D interfaces in Cu/Nb nanocomposites using atom probe tomography (APT). Our analysis demonstrates chemical heterogeneities along all spatial dimensions in 3D interfaces, establishes the length scale of such features, and quantifies the morphology of 3D interfaces. 3D interface heterogeneities form by surface diffusion during physical vapor deposition (PVD), suggesting that deposition parameters can be used to control interface structure and provide unique ways to explore processing-structure-property relationships in interface-dominated nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2024

A dynamically bare metal interface enables reversible magnesium electrodeposition at 50 mAh cm −2

Understanding and facilitating pure magnesium nucleation/growth electrodeposition behavior with ultrahigh Coulombic efficiency is complicated by the phenomenon of solid electrolyte interphase (SEI) formation in state-of-the-art, halogen-free magnesium electrolytes. Defining the electrolyte properties necessary to achieve ideal electrodeposition/stripping (E/S) thus remains elusive. Here, in this work, we reveal for the first time, rapid magnesium electrodeposition behavior that forms densely aligned, micron-sized thin platelets by establishing a dynamic bare magnesium/electrolyte interface during high-rate net plating. This effectively “SEI-free” interface allows facile magnesium diffusion and migration in stripping with near-unity E/S efficiency under demanding conditions over long-term cycling. The intrinsic electrolyte stability of the salt/solvent at the molecular level is the key to forming such an interface. The efficacy of the dynamic bare interface and an electrodeposited, free-standing magnesium anode is demonstrated in a high-areal-capacity full cell. These findings provide new design principles and fundamental understanding of interfacial chemistry in multivalent metal batteries.

Mg full cell↗
NASA NTRS1988

The surface morphology of crystals melting under solutions of different densities

Examples of solids melting beneath liquids are described for cases where the bulk liquid volume is stabilized against convection by a positive vertical temperature gradient, either with, or without local density inversion at the melting interface. The examples include ice melting beneath brine or methanol solutions and tin or lead melting under molten Sn-20 wt pct Pb or Pb-20 wt pct Sn, respectively. Without density inversion the melting is slow, purely diffusion controlled and the interfaces are smooth; with convection assisted melting the rate increases by some two orders of magnitude and the interfaces develop a rough profile - in the case of ice both irregular and quasi-steady state features are observed. The observations are discussed in terms of prevailing temperature and concentration gradients.

Fang, Dacheng↗
DOE OSTI2023

Understanding hydrogen electrocatalysis by probing the hydrogen-bond network of water at the electrified Pt–solution interface

Rational construction of the electrode–solution interface where electrochemical processes occur is of paramount importance in electrochemistry. Efforts to gain better control and understanding of the interface have been hindered by lack of probing methods. Here, in this study, we show that the hydrogen evolution and oxidation reactions (HER/HOR) catalysed by platinum in base can be promoted by introduction of N-methylimidazoles at the platinum–water interface. In situ spectroscopic characterization together with simulations indicate that the N-methylimidazoles facilitate diffusion of hydroxides across the interface by holding the second layer of water close to platinum surfaces, thereby promoting the HER/HOR. We thus propose that the HER/HOR kinetics of platinum in acid and base is governed by diffusion of protons and hydroxides, respectively, through the hydrogen-bond network of interfacial water by the Grotthuss mechanism. Moreover, we demonstrate a 40% performance improvement of an anion exchange membrane electrolyser by adding 1,2-dimethylimidazole into the alkali fed into its platinum cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Understanding Process–Structure Relationships during Lamination of Halide Perovskite Interfaces

Fabrication of halide perovskite (HP) solar cells typically involves the sequential deposition of multiple layers to create a device stack, which is limited by the thermal and chemical incompatibility of top contact layers with the underlying HP semiconductor. One emerging strategy to overcome these restrictions on material selection and processing conditions is lamination, where two half-stacks are independently processed and then diffusion bonded to complete the device. Lamination reduces the processing constraints on the top side of the solar cell to allow new device designs, expanded use of deposition methods, and self-encapsulation of devices. While laminated perovskite solar cells with high efficiencies and novel interlayer combinations have been demonstrated, there is a limited understanding of how the lamination process parameters affect the diffusion-bond quality and material properties of the resulting HP layer. In this study, we systematically vary temperature, pressure, and time during lamination and quantify the resulting impacts on bonded area, grain domain size, and photoluminescence. A design of experiments is performed, and statistical analysis of the experimental results is used to quantitatively evaluate the resulting process–structure–property relationships. The lamination temperature is found to be the key parameter controlling these properties. Furthermore, a temperature of 150 °C enables successful bonding over 95% of the substrate area and also results in increases in apparent grain domain size and photoluminescence intensity. Based on these insights, the lamination temperature of functional perovskite solar cell devices is varied, demonstrating the importance of the resulting bond quality on device performance metrics.

14 SOLAR ENERGY↗
DOE OSTI2022

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li 7 La 3 Zr 2 O 12 solid electrolytes

Batteries based on solid-state electrolytes, including Li 7 La 3 Zr 2 O 12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. Here these examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

25 ENERGY STORAGE↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Optical processing for semiconductor device fabrication

A new technique for semiconductor device processing is described that uses optical energy to produce local heating/melting in the vicinity of a preselected interface of the device. This process, called optical processing, invokes assistance of photons to enhance interface reactions such as diffusion and melting, as compared to the use of thermal heating alone. Optical processing is performed in a 'cold wall' furnace, and requires considerably lower energies than furnace or rapid thermal annealing. This technique can produce some device structures with unique properties that cannot be produced by conventional thermal processing. Some applications of optical processing involving semiconductor-metal interfaces are described.

Sopori, Bhushan L.↗
DOE OSTI2022

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗