Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion model”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Experimental and modeling study of the effect of confinement on the thermal decomposition of organic materials

Simultaneous Thermal Analysis is widely used to study thermal decomposition of a variety of organic materials, and the choice of experimental conditions is very important to obtain reliable mass loss and heat flow curves. Here in this work we investigate the effect of confinement on the thermal decomposition of organic materials. As an example, we explore high-density polyethylene (HDPE) and one of its pyrolysis products, namely eicosane (C 20 H 42 ), through experiments and modeling in open crucibles and crucibles closed with a pierced lid. A new model is developed for evaporation and sublimation from open and pinhole thermal-analysis pans. It considers the kinetic resistance of evaporation at the liquid-vapor interface, inhibition of evaporation by buildup of vapor inside the pan, thermal expansion of vapor inside the pan, and diffusion of vapor from the pan. The model is validated using simultaneous heat-flow and mass-loss measurements of n-eicosane evaporation for various pinhole sizes, down to 50 µm. The importance of product evaporation inhibition for measuring polymer decomposition was demonstrated using HDPE. This study sheds light on the effect of confinement on the mass loss rate of organic materials.

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

Toward a high-fidelity tritium transport modeling for retention and permeation experiments

We report Tritium Migration Analysis Program version 8 (TMAP8), the latest version of TMAP, was developed within the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE). Created at Idaho National Laboratory (INL), MOOSE is an open-source, dimension-agnostic, fully coupled, and fully implicit multiphysics platform featuring massively parallel computation capabilities. Using TMAP8, tritium transport in a divertor monoblock was analyzed to elucidate the effects of pulsed operation (up to fifty 1,600 s plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. With its built-in Message Passing Interface capability, TMAP8 can, in under 2 h, simulate tritium transport in three different layered materials (i.e., tungsten, copper, and copper-chromium-zirconium alloy) in 2D geometry, using a single device/computer with 10 cores. The MOOSE-based TMAP8 code can leverage other MOOSE tools developed under the Nuclear Energy Advanced Modeling and Simulation program to perform tritium and thermal transport in complex geometries and multiphysics environments. And via its massively parallel computation, MOOSE will enable the fusion pilot plant designers to conduct high-fidelity multiphysics modeling for the design of the divertor and blanket systems as well as for the safety analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Spontaneous Polarization Suppression of Exciton–Exciton Annihilation in Rhombohedral-Stacked Bilayer Molybdenum Disulfide

Rapid exciton–exciton annihilation (EEA) in two-dimensional semiconductors limits access to high-density excitonic regimes essential for efficient optoelectronic operation under strong excitation. Here, we show that EEA is strongly suppressed in rhombohedral (3R)-stacked MoS 2 bilayers relative to nonpolar 2H bilayers, and that this suppression is quantitatively consistent with repulsive dipole–dipole interactions between layer-polarized excitons induced by spontaneous polarization. Using ultrafast pump–probe spectroscopy, we measure an EEA rate of γ EEA = (5.03 ± 0.99) × 10 –3 cm 2 s –1 in 3R bilayers, which is approximately 18.2-fold smaller than that in monolayers and 2.9-fold smaller than that in nonpolar 2H bilayers. Despite the higher exciton diffusivity recently reported for 3R relative to 2H bilayers, the reduced EEA rate in 3R indicates a rate-limited regime governed by the close-encounter annihilation probability rather than diffusion. A rate-limited annihilation model incorporating a dipole–dipole repulsive potential captures the observed ratio γ EEA,3R /γ EEA,2H ≈ 0.35 for an exciton–exciton encounter distance of ∼1.3 nm, consistent with the bilayer exciton Bohr radius. These results indicate that spontaneous polarization in 3R-stacked bilayers can suppress nonlinear excitonic losses and provide a route toward high-density excitonics.

dipolar excitons↗
DOE OSTI2025

Spin-on deposition of amorphous zeolitic imidazolate framework films for lithography applications

Amorphous zeolitic imidazolate framework (aZIF) films have been recently introduced as resists for electron beam and extreme ultraviolet lithography. aZIFs are also being considered for separation applications, including thin film membranes. However, the reported methods for aZIF deposition are currently based on highly empirical trial-and-error approaches that hinder control of film composition, thickness and uniformity as well as scale-up and transferability to different coating geometries. This work presents a method for depositing aZIF films with controllable thickness using dilute precursors mixed immediately before encountering the substrate. Importantly, the method is amenable to quantitative analysis by computational fluid dynamics to extract intrinsic deposition rates and limiting reactant transport diffusivities, enabling predictive physics-based modeling of the deposition process. This allows the deposition method to be adapted for spin coating on silicon wafers to prepare high-quality aZIF films with consistently controlled thickness. Using this approach, high-resolution resist performance and wafer-scale use for beyond extreme-ultraviolet lithography of aZIF films is demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Room temperature buffer gas beam of metastable state titanium atoms

We produce beams of neutral titanium (Ti) atoms in their metastable 3d3(F4)4sa5F5 state by laser ablation into He, N2, and Ar buffer gases. The high temperatures associated with the ablation process populate the a5F5 level without the need for optical pumping. Remarkably, we observe that Ti atoms in the a5F5 state survive thousands of collisions with He and Ar buffer gas atoms without being quenched to lower-energy states. We study the yield of Ti atoms when ablated into buffer gases of varying species and pressure, quantify quenching rates and diffusion cross sections based on simple models, and provide insight into optimal design parameters for an ablation cell. Using a 3.3 cm ablation cell with interchangeable exit apertures, we produce metastable atom beams and quantify their brilliance and velocity distributions as functions of buffer gas pressure.

Schrott, Jackson↗
DOE OSTI2025

Multipoint Correlations in Poisson Media

Multipoint correlations in randomly heterogeneous composite media are determined by the probability that a set of points belong to specific phases. They determine a wide range of macroscopic transport properties such as conductivity, dielectric constant, diffusion coefficient, and transmittance. The Poisson model—a random tesselation of space by hyperplanes—provides realistic descriptions of heterogeneous media in, e.g., radiation-transport applications; yet, until now, it has lacked closed-form expressions for its multipoint correlations. We resolve this problem by presenting an exact solution for the multipoint correlations in the Poisson model. By comparing it to Monte Carlo simulations of four-point correlations in three dimensions, we demonstrate the accuracy of our solution. In conclusion, we visualize the multipoint correlations and discuss their features.

Amorphous materials↗
DOE OSTI2024

A Model to Account for the Effects of Load Ratio and Hydrogen Pressure on the Fatigue Crack Growth Behavior of Pressure Vessel Steels

A phenomenological model for estimating the effects of load ratio R and hydrogen pressure P H2 on the hydrogen-assisted fatigue crack growth rate (HA–FCGR) behavior in the transient and steady-state regimes of pressure vessel steels is described. The “transient regime” is identified with crack growth within a severely embrittled zone of intense plasticity at the crack tip. The “steady-state” behavior is associated with the crack growing into a region of comparatively lower hydrogen concentration located further away from the crack tip. The model treats the effects of R and P H2 as being functionally separable. In the transient regime, the effects of the hydrogen pressure on the HA–FCGR behavior were negligible but were significant in the steady-state regime. The hydrogen concentration in the steady-state region is modeled as being dependent on the kinetics of lattice diffusion, which is sensitive to pressure. Experimental HA–FCGR data from the literature were used to validate the model. The new model was shown to be valid over a wide range of conditions that ranged between –1 ≤ R ≤ 0.8 and 0.02 ≤ P H2 ≤ 103 MPa for pressure vessel steels.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Generative AI in Supply Chain Management: Applications, Challenges, and Future Directions

Supply chain management (SCM) is undergoing rapid transformation due to increasing global complexity, demand volatility, and operational disruptions. Generative Artificial Intelligence (GenAI) has emerged as a powerful paradigm capable of synthesizing data, simulating operational scenarios, and enabling adaptive decision-making across supply chain networks. This paper presents a survey of GenAI’s role in SCM, focusing on its applications in predictive analytics, autonomous logistics, and fraud detection. Unlike traditional AI systems that rely primarily on predictive analytics, GenAI models, including large language models, generative adversarial networks, and diffusion-based architectures, enable the creation of synthetic supply chain scenarios and autonomous optimization strategies. This survey provides (1) a taxonomy of GenAI techniques for supply chain applications, (2) a comparative analysis of generative AI approaches with traditional machine learning, reinforcement learning, and blockchain-based methods, and (3) a discussion of key challenges such as data privacy, interpretability, and integration with legacy enterprise systems. Furthermore, we outline open research problems and propose directions for future research toward autonomous, resilient, and sustainable AI-driven supply chains.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2024

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Analysis of Polarized Dust Emission Using Data from the First Flight of SPIDER

Using data from the first flight of Spider and from the Planck High Frequency Instrument, we probe the properties of polarized emission from interstellar dust in the Spider observing region. Component-separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses of diffuse Galactic dust emission spanning the full Spider region demonstrate (i) a spectral energy distribution that is broadly consistent with a modified-blackbody (MBB) model with a spectral index of β d = 1.45 ± 0.05 (1.47 ± 0.06) for E (B)-mode polarization, slightly lower than that reported by Planck for the full sky; (ii) an angular power spectrum broadly consistent with a power law; and (iii) no significant detection of line-of-sight polarization decorrelation. Tests of several modeling uncertainties find only a modest impact (~10% in σ r ) on Spider's sensitivity to the cosmological tensor-to-scalar ratio. The size of the Spider region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature T d = 19.6 K, an analysis of two independent subregions of that field results in inferred values of β d = 1.52 ± 0.06 and β d = 1.09 ± 0.09, which are inconsistent at the 3.9σ level. Furthermore, a joint analysis of Spider and Planck 217 and 353 GHz data within one subregion is inconsistent with a simple MBB at more than 3σ, assuming a common morphology of polarized dust emission over the full range of frequencies. This evidence of variation may inform the component-separation approaches of future cosmic microwave background polarization experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulative Prediction of Solar Illuminance and Application of the Du-Sharples Model in Estimating Adapted Daylighting Metrics for an Urban Environment

The practice of daylighting in indoor spaces can significantly reduce electricity consumption and carbon emissions, improve human productivity, and enhance mood and cognitive perception. This work discussed the recent developments in daylighting science and practice, computed the periodic variations in average diurnal daylight levels for each month, quantified in terms of global horizontal illuminance and diffuse horizontal illuminance, for Kolkata, India, a city with tropical wet and dry climate, with two empirical luminous efficacy models of estimating solar illuminance, and assessed daylighting metrics with the Du-Sharples model. A program was formulated that could compute and generate daylight data with monthly-hourly solar irradiation data and the Du-Sharples model was utilized to predict dirt-corrected daylighting metrics for three glazing transmittance values and five elemental carbon deposition levels on glazing material. The highest monthly average global horizontal illuminance is recorded in April (64.05 klx for Littlefair model and 66.82 klx for Muneer-Kinghorn model) and the highest monthly average diffuse horizontal illuminance is recorded in July (33.23 klx for Littlefair model and 30.63 klx for Muneer-Kinghorn model). Further, the computed yearly average global and diffuse horizontal illuminance levels agree well with a previous study that applied the Perez model. Yearly average horizontal work surface illuminance level remained >1.5 klx for window-towall area ratio >30 %. The approach adopted in this work and the temporal variation charts of computed exterior daylight level data may assist building service engineers, architects, and indoor lighting practitioners in making informed policy decisions at different stages of building planning.

Engineering↗
DOE OSTI2025

Detailed characterization of runaway electron driven whistler waves in low-density DIII-D discharges

RE-driven whistler waves during quiescent DIII-D shots have been investigated further. The waves are confirmed to be mostly perpendicularly propagating and are observed for the first time with frequencies up to 700 MHz. Phase-spectral analysis has been used to infer their toroidal mode numbers, n, which are expected to scale with the wavenumber, k, of the mode. Though we derive a theoretical scaling of k ≈ 4n, the measured mode numbers are found to exhibit a very weak dependence on k. In addition, increases in synchrotron emission have been found to consistently lag whistler wave bursts by roughly 3–5 ms, suggesting the waves are causing pitch-angle scattering, since the emitted synchrotron radiation is a strong function of the REs' perpendicular energy. The stronger the wave bursts, the greater the subsequent increase in synchrotron emission. A predator-prey model is used to describe these nonlinear wave-particle interactions, from which the wave damping rates and the loss parameter can be inferred. The damping rates are found to be of the order of (1.6 ± 0.8) × 10 4 /s, and the unitless loss parameter is found to be approximately 2, suggesting that the loss mechanism is diffusive. These observations will serve to validate models of RE-driven waves in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2024

Hadronization of heavy quarks

Heavy-flavor hadrons produced in ultrarelativistic heavy-ion collisions are a sensitive probe for studying hadronization mechanisms of the quark-gluon-plasma. In this paper, we survey how different transport models for the simulation of heavy-quark diffusion through a quark-gluon plasma in heavy-ion collisions implement hadronization and how this affects final state observables. Utilizing the same input charm-quark distribution in all models at the hadronization transition, we find that the transverse-momentum dependence of the nuclear modification factor of various charm hadron species has significant sensitivity to the hadronization scheme. In addition, the charm-hadron elliptic flow exhibits a nontrivial dependence on the elliptic flow of the hadronizing partonic medium.

Atomic, Nuclear and Particle Physics↗
DOE OSTI2024

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗