Search NASA⌕ Search

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.

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 OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2025

Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data

Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.

MPPD↗
DOE OSTI2024

A chemo-mechanical model for describing sorption hysteresis in a glassy polyurethane

Abstract Hysteretic sorption and desorption of water is observed from 0 to 95% relative humidity and 298–333 K on a glassy polyurethane foam. It is postulated that sorption-induced swelling of the glassy polyurethane increases the concentration of accessible hydrogen-bonding adsorption sites for water. The accessibility of sites is kinetically controlled due to the restricted thermal motions of chains in the glassy polymer, causing a difference in accessible site concentrations during sorption and desorption. This discrepancy leads to hysteresis in the sorbed concentrations of water. A coupled chemo-mechanical model relating volumetric strain, adsorption site concentration, and sorbed water concentration is employed to describe water sorption hysteresis in the glassy polyurethane. This model not only describes the final mass uptake for each relative humidity step, but also captures the dynamics of water uptake, which exhibit diffusion and relaxation rate-controlled regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Amorphization and siliconization of silicon carbide as a first wall material

The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dark Photons from Perturbative Decay of a Misaligned Higgs Field

We reconsider the production of dark photons $A'$ as dark matter, from the perturbative decay of a dark Higgs field $h$, that is stochastically misaligned from the minimum of its potential during inflation. This is a simple and predictive framework for generating the $A'$ relic abundance. It is constrained by structure formation, since the $A'$ are initially boosted, and inflationary isocurvature fluctuations, which require small quartic couplings $λh^4$. We identify $A'$ masses between 100 eV and 1 GeV and gauge couplings $g\sim 10^{-15}-10^{-10}$ that are consistent in this scenario, and which become more tightly constrained if a generic level of kinetic mixing is present. The favored parameter region could be tested through future CMB or Lyman-$α$ observations, and, in the presence of kinetic mixing, by direct detection experiments or diffuse soft gamma-ray searches.

Cline, James M. [McGill U., Montreal (main)]↗
DOE OSTI2025

Durability of PGM catalyst MEAs of polymer electrolyte membrane fuel cells for heavy-duty vehicles

Polymer electrolyte membrane fuel cells (PEMFCs) are promising power sources for heavy-duty vehicles (HDVs) owing to cleanliness and efficiency. However, the degradation of membrane electrode assemblies (MEAs) under HDV conditions remains a huge challenge. Here, this work investigated MEA durability under HDV conditions using a US Department of Energy standard accelerated stress test for 180,000 cycles (equivalent to 1 million miles of HDV operation). Effects of catalyst Pt content on MEA durability were examined using homemade 30% Pt/C (H-Pt/C) and commercial 46% Pt/C (C-Pt/C) catalysts. Both MEAs experienced H 2 /air and H 2 /O 2 performance loss over cycles. Analysis with scanning transmission electron microscopy, X-ray diffraction, inductively coupled plasma mass spectrometry, and mercury intrusion porosimetry revealed severe degradation of Pt nanoparticles (NPs), support structures, and the catalyst layer. Two degradation stages for NPs were proposed: Ostwald ripening dominated the initial 60,000 cycles, followed by combined Ostwald ripening and particle migration. Measurements with ion chromatography, high-frequency resistance, and oxygen-diffusion resistance revealed degradation of membrane and ionomer, respectively.

25 ENERGY STORAGE↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Solution-phase sample-averaged single-particle spectroscopy of quantum emitters with femtosecond resolution

Here, the development of many optical quantum technologies depends on the availability of solid-state single quantum emitters with near-perfect optical coherence. However, a standing issue that limits systematic improvement is the significant sample heterogeneity and lack of mechanistic understanding of microscopic energy flow at the single emitter level and ultrafast timescales. Here we develop solution-phase single-particle pump-probe spectroscopy with photon correlation detection that captures sample-averaged dynamics in single molecules and/or defect states with unprecedented clarity at femtosecond resolution. We apply this technique to single quantum emitters in two-dimensional hexagonal boron nitride, which suffers from significant heterogeneity and low quantum efficiency. From millisecond to nanosecond timescales, the translation diffusion, metastable-state-related bunching shoulders, rotational dynamics, and antibunching features are disentangled by their distinct photon-correlation timescales, which collectively quantify the normalized two-photon emission quantum yield. Leveraging its femtosecond resolution, spectral selectivity and ultralow noise (two orders of magnitude improvement over solid-state methods), we visualize electron-phonon coupling in the time domain at the single defect level, and discover the acceleration of polaronic formation driven by multi-electron excitation. Corroborated with results from a theoretical polaron model, we show how this translates to sample-averaged photon fidelity characterization of cascaded emission efficiency and optical decoherence time. Our work provides a framework for ultrafast spectroscopy in single emitters, molecules, or defects prone to photoluminescence intermittency and heterogeneity, opening new avenues of extreme-scale characterization and synthetic improvements for quantum information applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Sulfide‐Based Solid Electrolyte With High Humid Air Tolerance for Long Lifespan All‐Solid‐State Sodium Batteries

Abstract Sulfide‐based superionic conductors present great promise to achieve high energy density and safety for all‐solid‐state sodium batteries (ASSSBs). However, the poor electrolyte/electrode interface compatibility and humid air stability seriously hinder their deployment in ASSSBs. Herein, a series of high‐performance Na 3‐□ Sb 1‐4x (SnWCaTi) x S 4 sulfide‐based solid electrolytes (SSEs) are reported by coupling the vacancy effect with configurational entropy, which displays an excellent interface stability against sodium metal and an extraordinary tolerance toward the moist atmosphere, even for water. The optimized electrolyte effectively inhibits the detrimental mixed ion‐electron conducting interphase formation, achieving the ultra‐stable operation of Na–Na symmetric cell up to 1000 h. Furthermore, the Na + diffusion kinetics is obviously enhanced by increasing the Na sites local anisotropy and Na vacancies. Eventually, the assembled TiS 2 //Na 5 Sn ASSSBs deliver a remarkable reversible capacity of 211.6 mAh g −1 at 0.5C with a long‐term cycling performance of 450 cycles at room temperature. More importantly, it achieves a steady running up to 100 cycles at 1C even if this electrolyte is placed in the air with a dew temperature of 13.8 °C for 30 min, the highest values in the state‐of‐the‐art sulfide‐based ASSSBs. The well‐designed SSEs open a new avenue for realizing the advanced and powerful ASSSBs.

Guo, Yayu↗
DOE OSTI2025

Measurement of Zero-Frequency Fluctuations Generated by Coupling between Alfvén Modes in the JET Tokamak

We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with theoretical and numerical predictions of zonal modes pumped by Alfvén modes, and is correlated with an increase in the deep core ion temperature, temperature gradient, confinement factor H 89,P , and a reduction in the main ion heat diffusivity. Despite the energetic particle transport induced by the Alfvén eigenmodes, the generation of a zero-frequency fluctuation that can suppress the turbulence leads to an overall improvement of confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overcoming limitations in propanedehydrogenation by codesigning catalyst-membrane systems

Propylene production through propane dehydrogenation (PDH) is endothermic, and high temperatures required to achieve acceptable propane conversions lead to low selectivity and severe carbon-induced deactivation of conventional catalysts. We developed a catalyst-membrane system that removes the hydrogen by-product and can thus achieve propane conversions that exceed equilibrium limits. In this codesigned system, a silica/alumina (SiO 2 /Al 2 O 3 ) hollow-fiber hydrogen membrane was packed with a selective platinum-tin (Pt 1 Sn 1 /SiO 2 ) PDH catalyst on the tube side with hydrogen diffusing from the tube to the shell side. We demonstrate that the catalyst-membrane system can achieve propane conversions >140% of the nominal equilibrium conversion with a propylene selectivity >98% without deactivation of the system components. Here, we also show that by introducing oxygen on the shell side of the catalyst-membrane system, we can couple the endothermic PDH reaction on the tube side with exothermic hydrogen oxidation on the shell side. This coupling results in higher rates of hydrogen transport, leading to further enhancements in the propane conversion as well as desired thermoneutral system operation.

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

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗