Search NASA⌕ Search

SEARCH · Search NASA

Results for “generative 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 235 records · Page 13

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗

Dust Survival in Galactic Winds

This repository contains three-dimensional volumetric data from an Eulerian hydrodynamical simulation (conducted on a uniform Cartesian grid) generated by the Cholla hydrodynamics code. The datasets contain snapshots (full-grid, projections, and slices) in the HDF5 format of a multi-phase medium in which a hot, diffuse, dust-free background wind accelerates a cool, dense cloud of gas and dust. This scenario is intended to represent a supernova-driven galactic outflow, in which hot supernova winds are thought to accelerate cool interstellar medium material out of the galactic disk into the surrounding circumgalactic medium. There are three separate datasets for simulations corresponding to three cloud evolutionary scenarios: long-term cloud survival (surv), marginal cloud survival (disr), and cloud destruction (dest). Projection and slice images of the simulations are also included in this repository.

79 ASTRONOMY AND ASTROPHYSICS↗

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗

Systematic correction of the density functional theory spectra via a quantum Monte Carlo approach

Numerical outputs and driver scripts supporting auxiliary-boson corrected diffusion Monte Carlo (ABCDMC) benchmarks on second-row neutral atoms, cations, and dications, plus a C2 molecule active-space study. Includes NIST reference energies, ABCDMC timestep extrapolation summaries, orbital generation inputs, singles-only CASCI driver scripts, an O-atom basis-set convergence study, and the PyQMC boson code snapshot used for the calculations.

36 MATERIALS SCIENCE↗

Decoupling of Nitrogen and Oxygen Impurities in Nitrogen Doped SRF Cavities

The performance of superconducting radiofrequency (SRF) cavities is critical to enabling the next generation of efficient high-energy particle accelerators. Recent developments have focused on altering the surface impurity profile through \textit{in-situ} baking, furnace baking, and doping to introduce and diffuse beneficial impurities such as nitrogen, oxygen, and carbon. However, the precise role and properties of each impurity are not well understood. In this work, we attempt to disentangle the role of nitrogen and oxygen impurities through time-of-flight secondary ion mass spectrometry of niobium samples baked at temperatures varying from 75-800 $^\circ$C with and without nitrogen injection. From these results, we developed treatments recipe that decouple the effects of oxygen and nitrogen in doping treatments. Understanding how these impurities and their underlying mechanisms drive further optimization in the tailoring of impurity profiles for high-performance SRF cavities.

43 PARTICLE ACCELERATORS↗

Decoupling of Nitrogen and Oxygen Impurities in Nitrogen Doped SRF Cavities

The performance of superconducting radiofrequency (SRF) cavities is critical for enabling the next generation of efficient high-energy particle accelerators. Recent developments have focused on altering the surface impurity profile through \textit{in-situ} baking, furnace baking, and doping to introduce and diffuse beneficial impurities such as nitrogen, oxygen, and carbon. However, the precise role and properties of each impurity are not well understood. In this work, we attempt to disentangle the role of nitrogen and oxygen impurities through time-of-flight secondary ion mass spectrometry using niobium cavity cutouts baked at temperatures varying from 120-800$^\circ$C with and without nitrogen injection. Understanding the underlying mechanisms behind these impurities will drive further optimization in the tailoring of nitrogen and oxygen profiles for high-performance SRF cavities.

43 PARTICLE ACCELERATORS↗

Tuning Cation (Dis)Order in Cr‐Based Li‐Excess Oxide Cathode Materials to Improve Li+ Transport Properties

Abstract Li‐excess disordered rocksalts (DRXs) hold promise as next‐generation cathodes for Li‐ion batteries due to their high capacity and energy density, along with the potential to eliminate the need for Co and Ni. However, due to their low Li + diffusivity, DRXs need to be pulverized into nanoparticles to achieve high performance. Herein, a new strategy for overcoming this limitation is demonstrated, involving the design of as‐synthesized partially disordered oxides with a structure that lies between ordered layered and fully disordered, exhibiting a varying degree of local (dis)order. This unique structure activates new Li + diffusion channels, improving percolation and transport properties. This strategy allows a large content of Li + to be accessed in material with large, micron‐sized particles through highly reversible Cr 3+/6+ and O redox, yielding a first discharge capacity of 286 mAh g −1 (881 Wh kg −1 ). The Li + percolation network is further improved by substituting Ti with a mixture of multiple metals, which appears to locally decrease the migration barrier through lattice distortion. Proper tuning of the chemical composition, especially the content of metals with empty d orbitals, is established as a crucial factor for controlling the degree of disorder and mitigating voltage fade and hysteresis growth upon cycling.

Moździerz, Maciej↗

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗

SCALE 6.3 Validation: Nuclear Criticality Safety

The KENO codes have been used extensively around the world since their initial development in the 1960s. Therefore, a large number and wide range of validation reports have been generated for KENO by a variety of organizations. This section provides a synopsis of recent validation reports generated at ORNL, along with a brief list of some older validation reports generated in Oak Ridge, Tennessee. Some of the oldest reports are not attributed directly to ORNL, as the Y-12 Plant (currently the Y-12 National Security Complex), the K-25 Gaseous Diffusion Plant (now the East Tennessee Technology Park), and the X-10 Site (ORNL) were operated as a single entity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗

Polyethylene Upcycling to Diacids Using Acid-Only Activation

There is a pressing need to develop plastic recycling technologies. Chemical upcycling converts low-value waste plastics into higher-value products. Here, polyethylene (PE) upcycling is accomplished by using acid-only activation of PE (PEAA), thus eliminating costly and toxic organic solvents. Acid mixtures of chlorosulfuric acid (CSA) and sulfuric acid (SA) were used to effectively sulfonate PE, allowing for subsequent facile depolymerization. Molecular deconstruction of sulfonated PE using H 2 O 2 in the presence of an Fe(III) catalyst achieved molar yields of C 2 –C 4 diacids exceeding 40% (based on estimates of acids generated from sulfonates). The sulfonation step using PEAA is hypothesized to follow a shrinking core mechanism (SCM), in which PE particles quickly reach a saturated state in the outer sulfonated layer as CSA diffuses radially inward. This proposed mechanism is supported by the fact that, regardless of the extent of sulfonation, similar molar yields of products and ratios of liberated carbon to sulfur are observed. Modeling and experimentation show that deviations from the SCM occur by lowering the Damköhler number (Da) or increasing CSA solubility in PE. In both cases, CSA diffusion is enhanced, which favors kinetically limited sulfonation.

36 MATERIALS SCIENCE↗

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗

Refining fast calorimeter simulations with a Schrödinger Bridge

Machine learning-based simulations, especially calorimeter simulations, are promising tools for approximating the precision of classical high energy physics simulations with a fraction of the generation time. Nearly all methods proposed so far learn neural networks that map a random variable with a known probability density, like a Gaussian, to realistic-looking events. In many cases, physics events are not close to Gaussian and so these neural networks have to learn a highly complex function. We study an alternative approach: Schrödinger bridge Quality Improvement via Refinement of Existing Lightweight Simulations (SQuIRELS). SQuIRELS leverages the power of diffusion-based neural networks and Schrödinger bridges to map between samples where the probability density is not known explicitly. We apply SQuIRELS to the task of refining a classical fast simulation to approximate a full classical simulation. On simulated calorimeter events, we find that SQuIRELS is able to reproduce highly non-trivial features of the full simulation with a fraction of the generation time.

Calorimeter methods↗

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↗

A Straightforward Model for Quantifying Local pH Gradients Governing the Oxygen Evolution Reaction

The production and consumption of protons by an electrocatalyst will, under certain conditions, generate localized microenvironments with properties distinct from those of the bulk solution. These local properties are particularly impactful for reactions involving proton-coupled electron transfer, where the generation of locally basic or acidic environments may significantly influence the energy efficiency and reaction selectivity of the electrocatalyst. Whereas local pH environments have been observed and characterized in reductive half-reactions, including the CO 2 reduction and hydrogen evolution reactions, the incompatibility of conventional techniques and materials has limited studies in oxidative half-reactions, including the oxygen evolution reaction (OER), which provides the reducing equivalents for solar-to-fuels electrolysis. With the straightforward parameters bulk pH, buffer composition and pK a , and mass transport, we develop a model for describing local pH as a function of current density regardless of the microscopic details of the mechanism. Using an acid-stable PbO x OER catalyst, we observe the formation and dissipation of pH gradients during the OER and validate the model with voltammetric and potentiometric studies. Here, the model predicts how local acidic environments can develop over a narrow OER current density window, thus providing further motivation for the development of OER catalysts that are stable to acid, even when operating in basic aqueous conditions. More generally, the model is not restricted to the OER and is useful for determining the onset of local pH gradients for other electrocatalytic reactions that involve the consumption or generation of protons in energy conversion reactions.

Anions↗

Domestically Manufactured Fuel Cells for Heavy Duty Applications

This project, funded by the U.S. Department of Energy under Award DE-EE0009248, advances the domestic development of fuel-cell stack technologies for heavy-duty transportation applications. The work directly supports the Million-Mile Fuel Cell Truck (M2FCT) initiative by improving performance, durability, and manufacturability—key requirements for next-generation hydrogen fuel-cell systems.

08 HYDROGEN↗