SEARCH · Search NASA

Results for “Diffusion Couple”

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.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Selective Recovery of Critical Minerals from Simulated Electronic Wastes Via Reaction‐Diffusion Coupling

Abstract Atom‐ and energy‐efficient chemical separations are urgently needed to meet the surging demand for critical materials that has strained supply chains and threatened environmental damage. In this study, we used reaction‐diffusion coupling to separate iron, neodymium, and dysprosium ions from model feedstocks of permanent magnets, which are typically found in electronic wastes. Feedstock solutions were placed in contact with a hydrogel loaded with potassium hydroxide and/or dibutyl phosphate, resulting in complex precipitation patterns as the various metal ions diffused into the reaction medium. Specifically, we observed the precipitation of up to 40 mM of iron from the feedstock, followed by the enrichment of 73 % dysprosium, and the extraction of >95 % neodymium product at a further distance from the solution‐gel interface. We designed a series of experiments and simulations to determine the relevant ion diffusivities, D Nd =5.4×10 −10 and D Dy =5.1×10 −10 m 2 /s, and precipitation rates, k Nd =1.0×10 −5 and k Dy =5.0×10 −3 m 9 mol −3 s −1 , which enabled a numerical model to be established for predicting the distribution of products in the reaction medium. Our proof‐of‐concept study validates reaction‐diffusion coupling as an effective and versatile approach for critical materials separations, without relying on ligands, membranes, resins, or other specialty chemicals.

Wang, Qingpu [Physical and Computational Sciences ↗

DOE OSTI2025

A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction

Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.

Wang, Ziyan [Energy Geosciences Division Lawrence ↗

DOE OSTI2025

Diffusion and phase formation in the γ-uranium-technetium system

Phase formation in the U-Tc binary system at 800 °C was investigated using a diffusion couple experiment. Scanning electron microscopy (SEM) and energy dispersive X-ray spectroscopy (EDS) identified four novel potential intermetallic phases - U 7 Tc 3 , U 13 Tc 12 , U 3 Tc 5 , and UTc 4 . Diffusion coefficients were calculated for the intermetallic phases using the Boltzmann-Matano method and were respectively found to be – 120, 38.2, 15.6, and 1.51 × 10 −18 m 2 /s. Tc also exhibits a solid-solution phase with high penetration into the U with a diffusion coefficient of ∼ 10 −14 m 2 /s. Furthermore, these findings expand the number of known U-Tc phases and provide the first diffusion coefficients for the U-Tc system, and contribute valuable data to the broader field of actinide metallurgy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Chemical interaction and compatibility of uranium mononitride and alumina forming austenitic stainless steel

Uranium mononitride (UN) and alumina forming austenitic (AFA) stainless steel are a potential fuel-cladding combination for the lead-cooled fast reactor (LFR). Chemical compatibility between UN and AFA steel needs to be verified before implementation in a nuclear reactor. Diffusion couple experiments at 823 K and 1023 K were conducted for nonirradiated UN samples in contact with as-cast (no thermally grown Al 2 O 3 ) and preoxidized (with thermally grown Al 2 O 3 ) AFA for 500 and 1000 h in an inert environment. Preoxidized AFA exhibited little to no interaction with all UN samples tested at both 823 K and 1023 K, displaying the stability and capability of the Al 2 O 3 layer to prevent chemical interaction and inter-diffusion with UN. Chemical interaction occurs between UN and as-cast AFA. At 1023 K, an aluminum and nitrogen rich phase (likely AlN) formed along the interface of as-cast AFA and UN samples. At 823 K the AlN phase was not prominently observed due to the reduced diffusivity of aluminum through AFA. The aluminum and nitrogen-enriched phase was also observed in a high temperature pressure-assisted test sample of UN and as-cast AFA thermally treated at 1373 K. Finally, in UN samples doped with a low weight percent of UO 2 (< 3 wt%), AlN was not detected along the interface at either temperatures, and an Al 2 O 3 layer likely formed along the interface and prevented further chemical interaction between UN and as-cast AFA.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2026

Thermochemical Stability of Ca2Yb8(SiO4)6O2 Apatite in Presence of Molten Calcium-Magnesium-Aluminosilicate (CMAS)

Thermochemical stability of ytterbium silicon oxyapatite Ca 2 Yb 8 (SiO 4 ) 6 O 2 (CYbS) in the presence of molten calcium-magnesium aluminosilicate (CMAS) has been investigated at elevated temperatures for consideration as a thermal and environmental barrier coating (T/EBC) material. CYbS apatite powder was synthesized from the constituent oxides via a solid-state reaction method. Hot-pressed apatite substrates were exposed to molten CMAS at 1200, 1300, and 1400 °C for 1, 10, and 50 h. Development of phases in the interaction region of the heat-treated specimens was monitored using scanning electron microscopy, transmission electron microscopy, high-angle annular dark-field imaging, selected area electron diffraction, and energy dispersive x-ray spectroscopy. Monoclinic cyclosilicate Ca 3 Yb 2 (Si 3 O 9 ) 2 formed from interaction of CYbS apatite with CaO in the CMAS melt at the apatite-CMAS reaction front and continued to nucleate and grow within the residual CMAS in diffusion couples annealed for 1 to 50 h at 1200 °C as well as in those heat treated at 1300 °C for 1 h. Residual CMAS was depleted of Ca when cyclosilicate was present. Dendritic wollastonite (CaSiO3) was observed within the residual CMAS in couples annealed at 1200 and 1300 °C. Ingress of molten CMAS, because of its exponential decrease in viscosity, occurred through open pores and along the grain boundaries of the apatite substrates without any detectable chemical reaction at 1300 and 1400 °C. Results of this study indicate that Ca 2 Yb 8 (SiO 4 ) 6 O 2 apatite has the potential to mitigate the CMAS corrosion up to about 1200 °C but not at higher temperatures.

Narottam P Bansal↗

DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗

DOE OSTI2026

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗

DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

DOE OSTI2026

Interfacial Void Formation and Self-Healing in Oxide Scales on Al-containing High-Entropy Alloy

The exceptional high-temperature oxidation resistance of Al-containing high-entropy alloys (HEAs) is often attributed to the formation of a protective α-Al 2 O 3 scale. However, the dynamic, atomic-scale mechanisms governing the stability of this scale—including interfacial void formation and the often-postulated but rarely visualized “self-healing” capacity—remain poorly understood. Herein, we reveal the complex evolution of the triple-layer oxide scale on an Al 10 CoCrFeNi HEA through combined electron microscopy and diffraction study. We show that interfacial voids are an inherent consequence of the scaling process, originating from two distinct mechanisms: the Kirkendall effect at the interface between the γ-Al 2 O 3 /α-Al 2 O 3 and alloy driven by cationic diffusion imbalance and volumetric contraction due to phase transformations at the spinel/Cr 2 O 3 interface. Crucially, we provide microstructural evidence consistent with an intrinsic self-healing response. This process is driven by coupled inward diffusion of oxygen and outward diffusion of metal cations, leading to the in-situ formation of transient θ-Al 2 O 3 and spinel phases that partially fill and seal the voids. Here, these results provide atomic-scale insights into the phase evolution, defect formation, and self-repair of oxide scales in HEAs—highlighting pathways to enhance their oxidation resistance in extreme environments.

High-entropy alloy↗

DOE OSTI2026

Binder-Free Graphite Anodes for Next-Generation High-Performance Lithium-Ion Batteries

High-energy density anodes are crucial for next-generation lithium-ion batteries (LIBs) particularly for electric vehicle (EV) applications. Sluggish lithium-diffusion kinetics coupled with conventional anode fabrication processes containing polymeric binders hinder fast-charging capabilities and high-energy density of graphite. Herein, we introduce a binder-free graphite anode fabrication strategy using the electrospinning technique that contains ~2.41% carbon nanotubes (CNTs). Our strategy relies on the formation of an interconnecting conductive CNT network coupled with an ultrathin N-doped carbon coating on graphite particles from sacrificial binders. This combination enhances both structural integrity and electrical conductivity and, in turn, improves fast-charging capabilities and high energy density of LIBs. The binder-free graphite anode achieves ~335.0 mAh g–1 capacity at C/3 rate over 400 cycles with capacity retention of >95% and average Coulombic efficiencies >99.95%. These promising results suggest that the binder-free anode fabrication with a multifunctional design approach could elevate the energy-density limits of the graphite anodes, solving high-energy density requirements of EVs, and potentially provides a path forward for the development of economically feasible energy storage systems for various applications.

Ozcan, Muca [ORNL] (ORCID:0000000320020474)↗

DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗

DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A Transferable Force Field for Predicting Adsorption and Diffusion of Water in Cationic Zeolites with Coupled Cluster Accuracy

We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.

Adsorption↗

DOE OSTI2026

FROMP-from-Monomer

This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.

Chua, Lauren [Massachusetts Inst. of Technology (M↗

DOE OSTI2025

Asymmetric rotations slow down diffusion under confinement

Translation and rotation are the two most fundamental forms of diffusion, yet their coupling mechanism is not clear, especially under confinement. Here, we provided evidence of the coupling between rotation and translation using a substituted benzene molecule as an example. A counterintuitive behavior was observed where the movement of the smaller molecule with an asymmetric shape was unexpectedly slower than the larger one with a symmetric shape in confined channels of zeolite. We showed that this diffusion behavior was caused by the presence of the specific and selective interaction of the asymmetric guest with the pores, which increased the local restricted residence time, thus inhibiting the translation under confinement, as further confirmed by dynamic breakthrough curves, uptake measurements, quasi-elastic neutron scattering, and 2 H solid-state NMR techniques. Our work correlated asymmetric rotation and diffusion under a confined environment, which enriched our understanding of the coupling between rotation and translation and could shed light on a fundamental understanding of the diffusion process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Chemo-Mechanics of α-V 2 O 5 During Lithiation and Implications for Rechargeable Battery Cathodes

Chemo-mechanical degradation of layered oxide electrodes is strongly influenced by crystallographic anisotropy, local stress evolution, and ion insertion, yet the intrinsic mechanical response of layered materials remains incompletely understood. Indeed, most prior studies have focused on polycrystalline materials but single crystals enable direct observation of coupling between anisotropic ion diffusion and mechanical response. This study aims to determine how crystallographic anisotropy and lithiation affect deformation, fracture, and mechanical properties in single-crystal V 2 O 5 , and compares this behavior with polycrystalline counterparts. Polycrystalline V 2 O 5 thin films and single-crystal α-V 2 O 5 were studied using nanoindentation, scanning electron microscopy, focused ion beam cross-sectioning, and Raman spectroscopy. Single crystals were tested in pristine and chemically lithiated states, including experiments in which crystals were first plastically deformed via nanoindentation and subsequently lithiated. Polycrystalline films exhibited significantly higher hardness and elastic modulus than single crystals. Single crystals indented normal to the exposed (001) basal plane exhibited pronounced anisotropic deformation, including directional slip, crystallographically-guided cracking, anisotropic crack propagation, interlayer separation, and shear localization. Lithiation caused substantial softening, reduced hardness and modulus, and suppressed displacement bursts during nanoindentation, while previously indented regions showed crack formation and growth upon lithiation. Mechanical behavior of α-V 2 O 5 is strongly governed by crystallographic anisotropy and further altered by lithiation, with pre-existing deformation serving as a strong driver of fracture during ion insertion. These findings illuminate the coupling among ion insertion, deformation, and fracture in layered oxides and provide a basis for understanding and mitigating mechanical failure in electrochemical energy-storage materials.

Anisotropy↗