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 181 records · Page 10

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↗

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↗

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↗

Leak-tight crack repair for 304L stainless steel using friction surfacing

Here, friction surfacing (FS) was investigated as a method for sealing 50 μm wide thru-cracks in 3-mm-thick 304L stainless steel (304L) plate. Friction surfacing was performed using two sizes of 304L consumable rod, with diameters of 9.52 mm and 12.7 mm, on two substrate conditions (clean and oxidized). Friction surfacing was able to deposit a 600 μm thick coating and repair cracks to a depth of 100 μm - 200 μm below the original surface, when using the 12.7-mm-diameter consumable rod for both substrate conditions. Helium leak rates of 10 -10 atm-cc/s were achieved on crack repairs, designating them as leak-tight by the ANSI N14.5 standard. Optical and scanning electron microscopy were used to investigate the coating microstructure and bond interfaces. In this study, the 12.7-mm-diameter rod performed better with disruption of the oxide layer on the oxidized substrate and providing a more homogenous coating. Mechanical properties of the coated samples were evaluated by performing micro-indentation, tensile, bending, and adhesion testing. Higher hardness in the coatings was observed due to the fine equiaxed grains resulting from dynamic recrystallization. Complete consolidation of plastically deformed material on the substrate and a strong diffusion bond across the interface was observed. This was reflected in tensile and bending properties of the FS coatings being comparable to those of uncoated specimen. No delamination of the FS coating was observed during adhesion testing. The results demonstrate that friction surfacing is a viable option to seal cracks in stainless steels.

36 MATERIALS SCIENCE↗

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗

Interfaces enhanced plasma irradiation resistance in CrMoTaWV/W multilayer films through blocking He diffusion

The performance of plasma-facing materials (PFMs) is one of the key factors that significantly impact the stability of operation in fusion reactors. Herein, a new CrMoTaWV/W (high entropy alloy (HEA)/W) multilayer structure is designed as PFM to investigate its resistance to He plasma irradiation. It was observed that the introduction of the interfaces effectively absorbed plenty of He atoms, preventing them from diffusing into the material and delaying the formation of fuzz incubation zone, therefore, enhancing the resistance to plasma irradiation. The thickness transformed to fuzz in the HEA/W multilayer films was observed to be about two-thirds of those in the CrMoTaWV (HEA) film. Additionally, the fuzz growth rates in HEA/W multilayer films are lower than the average growth rate of bulk W and HEA films combined. These findings highlight a promising new avenue for the exploration of high-performance PFMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unusual electrical conductivity driven by localized stoichiometry modification at vertical epitaxial interfaces

Precise control of lattice mismatch accommodation and interfacial cation diffusion is critical to modulate correlated functionalities in epitaxial heterostructures, particularly when the interface composition is positioned near a compositional phase transition boundary. Here we select La 1-x Sr x MnO 3 (LSMO) as a prototype phase transition oxide and establish vertical epitaxial interfaces with NiO for exploring the strong interplay between strain accommodation, stoichiometry modification, and localized electron transport across the interface. It is found that localized stoichiometry modification overcomes the dead layer problem plaguing LSMO and leads to strongly directional conductivity, as manifested by a more than three orders of magnitude difference between out-of-plane and in-plane conductivity. Comprehensive atomic scale structural characterization and transport measurements reveal that this emerging behavior is created by a compositional change produced by preferential cation diffusion that pushes the LSMO phase transition from insulating into metallic within an ultrathin interface region. This study explores the nature of unusual electrical conductivity at vertical epitaxial interfaces and highlights the beneficial role of controllable cation diffusion that enables emerging functionalities for a broad range of potential applications such as memristors, spintronic devices, and novel nanoelectronic devices using strongly correlated materials.

36 MATERIALS SCIENCE↗

Creep and failure at metal-oxide interfaces

Small-scale bicrystal creep experiments were performed on contacts formed via in situ high-temperature diffusion bonding of metal-oxide interfaces including Ag-ZrO 2 , Pd-ZrO 2 , Pt-ZrO 2 , and Ag-high entropy oxide. This work characterizes deformation and failure at metal-oxide interfaces during mechanical loading. Interfacial sliding can be activated easily, while tensile interfacial creep was not observed at any condition of stress or temperature measured. Plastic strain, instead, localizes within the metal under tensile loading. A variety of mechanisms for plastic strain occur in the metal including lattice dislocation-mediated plasticity, twinning, low-angle grain boundary formation, and low-angle grain boundary creep. Surface and low-angle grain boundary diffusion occur under conditions where no metal-oxide tensile creep is observed, highlighting the significant differences in their interfacial mechanical response. High-temperature interfacial failure occurs when the mean curvature at the contact neck is approximately zero and the applied stresses comparable to brittle fracture stresses. The brittle fracture stresses were measured to be σ ƒ = 180 ± 90 MPa at the Ag-ZrO 2 interface at 225 °C, σ ƒ = 460 ± 160 MPa at the Pd-ZrO 2 interface at 680 °C, and σ ƒ = 640 ± 440 MPa at the Pt-ZrO 2 interface at 1010 °C.

Creep↗

High-Throughput In-Line Deposition of Silicon Oxide for Polycrystalline Silicon Passivating Contacts

Polycrystalline silicon passivating contacts rely on an ultrathin (1–2 nm) silicon oxide layer to minimize recombination at the wafer/oxide interface and regulate dopant diffusion. Traditionally formed by thermal or chemical oxidation, this oxide is herein replaced by silicon oxide deposited via aerosol impact-driven assembly (AIDA), enabling high wafer-per-hour throughput and precise thickness control. In this study, we show that AIDA coatings conformally cover planar or textured substrates and achieve a SiO x /poly-Si(n) structure with an implied open-circuit voltage (iV oc = 726 mV) and contact saturation current density (J 0 = 8.8 fA/cm 2 ). Furthermore, annealing AIDA SiO x films at elevated temperatures desorbs hydroxyl groups while the stoichiometry transitions toward SiO 2 , improving passivation quality. Together, these results highlight AIDA’s potential for scalable, high-throughput manufacturing of advanced passivating contacts, offering a cost-effective alternative to conventional low-pressure chemical vapor deposition and plasma-enhanced chemical vapor deposition-based silicon and oxide processes.

TOPcon↗

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗

Quantitatively controlled electrophoretic deposition of nanocrystal films from non-aqueous suspensions

This study presents a novel method to correlate the mass and charge transfer kinetics during the electrophoretic deposition of nanocrystal films by using a purpose-built double quartz crystal microbalance combined with simultaneous current-measurement. Our data support a multistep process for film formation: generation of charged nanocrystal flux, charge transfer at the electrode, and polarization of neutral nanocrystals near the electrode surface. The polarized particles are then subject to dielectrophoretic forces that reduce diffusion away from the interface, generating a sufficiently high neutral particle concentration at the interface to form a film. The correlation of mass and charge transfer enables quantification of the nanocrystal charge, the fraction of charged nanocrystals, and the initial sticking coefficient of the particles. These quantities permit calculation of the film thickness, providing a theoretical basis for using concentration and voltage as process parameters to grow films of targeted thicknesses.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

On the high-temperature stability of the Al 8 Cu 3 Ce intermetallic in an additively manufactured Al-Cu-Ce-Zr alloy

High-temperature resistant eutectic Al alloys are crucial materials for lightweight and energy efficient design in the automotive and aviation industries. Additive manufacturing offers a pathway to refine eutectic microstructures and develop novel alloys with superior high-temperature strength. High-volume fraction intermetallic Al-Cu-Ce alloys have been developed to deliver high-temperature strength in combination with reduced hot-tearing susceptibility. Zr is added to provide additional strengthening via nanoscale Al 3 Zr precipitation, and to stabilize and avoid coarsening of the Al 8 Cu 3 Ce phase. However, the detailed interaction between Zr and Al 8 Cu 3 Ce remains unexplored. In this work, we show with synchrotron X-ray diffraction that laser powder bed fusion fabricated Al-Cu-Ce and Al-Cu-Ce-Zr alloys contain predominantly the Al 8 Cu 3 Ce intermetallic in the as-fabricated condition. Heat treatment of the Al-Cu-Ce alloy results in the Al 8 Cu 3 Ce → Al 8 Cu 4 Ce phase transformation. In the Al-Cu-Ce-Zr alloy, minor fractions of (Al,Cu,Si) 4 Ce and Al 2 Cu-θ are found in the as-fabricated condition, while Al 8 Cu 3 Ce remains stable during heat treatment. Atom probe microscopy quantifies intermetallic stoichiometries and reveals how Zr is enriched at the Al-matrix/Al 8 Cu 3 Ce interface acting as a diffusion barrier against solute exchange. Calibrated thermodynamic modeling underpins this as a kinetic effect. Here, a qualitative microstructural model summarizes, how Zr stabilizes Al 8 Cu 3 Ce against phase transformations and coarsening.

36 MATERIALS SCIENCE↗

Strain-Driven Surface Reactivity in Magnesium-Ion Battery Cathodes

Successful deployment of a Mg-ion battery requires cathodes that can achieve reversible Mg intercalation and high energy density. Recent theoretical and experimental studies indicated that the overall transport is likely limited by sluggish Mg transport at the cathode-electrolyte interface and not Mg diffusion through bulk. Here, in this work, we investigated the surface electrochemical activity of Mg ions by using a spinel-structured manganese oxide thin-film model system and in situ X-ray scattering. In combination with post-mortem microscopy analysis, we found that magnesium insertion was more favorable than subsequent extraction near the surface of the Mg x Mn 2 O 4 film, resulting in overmagnesiation, and eventually amorphization of the surface. This structural irreversibility and high overpotential required for Mg extraction could explain significant voltage hysteresis and Mg surface enrichment previously observed in bulk cathodes. Density functional theory calculations suggested that the tendency for the Mg surface enrichment could be associated with Mg diffusion kinetics, which varies with the strain state evolved due to constrained film volume change during Mg insertion and extraction. Particularly, out-of-plane Mg migration was predicted to be favorable in the tensile strain rather than in the compressive case.

25 ENERGY STORAGE↗

Accumulation of Antifreeze Proteins on Ice Is Determined by Adsorption

Antifreeze proteins (AFPs) enable the survival of organisms in freezing environments by binding to ice crystals and suppressing their growth. The rate at which AFPs accumulate on ice is determined by an interplay between their diffusion from the aqueous solution and their adsorption to the ice crystal. To interrogate the relative importance of these two processes, we combined fluorescence experiments with reaction–diffusion modeling. We found that as AFPs accumulate on ice, their concentration in the aqueous solution does not develop a gradient but remains equal to its bulk concentration throughout our experiments. Furthermore, these findings led us to conclude that AFP accumulation on ice crystals, which are smaller than 100 μm in radius, is not limited by the diffusion of AFPs, but by the kinetics of AFP adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficiency of bulk perovskite-sensitized upconversion: Illuminating matters

Photon upconversion via triplet–triplet annihilation could allow for the existing efficiency limit of single junction solar cells to be surpassed. Indeed, efficient upconversion at subsolar fluences has been realized in bulk perovskite-sensitized systems. Many questions have remained unanswered, in particular, regarding their behavior under photovoltaic operating conditions. Here, we investigate the impact of repeated and continuous illumination on bilayer perovskite/rubrene upconversion devices. We find that variations of the underlying perovskite carrier recombination dynamics greatly impact the upconversion process. Trap filling and triplet sensitization are in direct competition: more saturated trap states in the perovskite and, thus, longer underlying perovskite photoluminescence lifetimes allow for an increased number of carriers to diffuse to the perovskite/rubrene interface and undergo charge extraction to the triplet state of rubrene. As a result, the upconversion efficiency is greatly influenced by the underlying trap density: the upconverted photoluminescence intensity increases by two orders of magnitude under continuous illumination for 4 h. This shows that the upconversion efficiency is difficult to define for this system. Importantly, these results indicate that perovskite-sensitized upconversion devices exhibit peak performance under continuous illumination, which is a requirement for their successful integration into photovoltaics to help overcome the Shockley–Queisser limit in single junction solar cells.

VanOrman, Zachary A. (ORCID:000000027773012X)↗

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗

Diffuse energy transport and coda-wave interferometry for resonant transmission between reverberant structures

Approximate analytic and numerical theories are developed with which to model, and compare with laboratory measurements, the diffuse ultrasonic field transmitted from one elastic body to another through a resonant coupling. Particular focus is on the sensitivity of the transmitted field to perturbations in the natural frequency of the coupling, and on the manifestation of such perturbations as time delays in the second body. Immediate application is to measurements of slow dynamic elastic nonlinearity. It is found that the developed approximate theories do a good job of predicting the time dependence of mean square transmitted signal and signal spectra. They also predict the time delays, both their erratic character and their means. Analytic dependence of these time shifts on the fractional frequency perturbation df/f is derived.

42 ENGINEERING↗