SEARCH · Search NASA

Results for “generative 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Modified Sand’s Time Incorporating Li-Ion Transport Across the SEI: Basis for Understanding Li Dendrite Formation and Li-Metal Battery Electrolyte Selection

Abstract Understanding the initiation of lithium dendrites remains elusive, largely due to the intricate role of the solid electrolyte interphase (SEI) which forms on the Li surface during electrodeposition. Many studies have utilized the classical Sand’s equation to estimate the onset time when lithium dendrites begin to form. The Sand’s equation provides the time when the cation (Li+) concentration at the electrode-electrolyte interface approaches zero under diffusion-limited conditions in galvanostatic Li electrodeposition. However, recent experimental studies have revealed that the observed lithium dendrite onset time deviates considerably from the Sand’s time. Here, we show that this deviation from classical theory is likely due to the transport of Li+ ions through the SEI - a transport limitation that is much more dominant in controlling dendrite formation. We develop a ‘modified’ Sand's equation, incorporating the SEI layer and the diffusional transport across it to predict Li dendrite onset times. To validate this approach, we conducted Li electrodeposition experiments at various current densities using two distinct organic electrolytes. Analysis of the results demonstrates that the modified Sand's equation provides a more accurate prediction of dendrite onset times, highlighting the importance of incorporating SEI into transport models of Li plating in next-generation rechargeable Li-metal batteries.

Ma, Yuanman (ORCID:0000000200444811)↗

DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗

DOE OSTI2026

Time-integrated Southern-sky Neutrino Source Searches with 10 yr of IceCube Starting-track Events at Energies Down to 1 TeV

In the IceCube Neutrino Observatory, a signal of astrophysical neutrinos is obscured by backgrounds from atmospheric neutrinos and muons produced in cosmic-ray interactions. IceCube event selections used to isolate the astrophysical neutrino signal often focus on the morphology of the light patterns recorded by the detector. The analyses presented here use the new IceCube Enhanced Starting Track Event Selection (ESTES), which identifies events likely generated by muon–neutrino interactions within the detector geometry, focusing on neutrino energies of 1–500 TeV with a median angular resolution of 1.4°. Selecting for starting-track events filters out not only the atmospheric-muon background but also the atmospheric-neutrino background in the southern sky. This improves IceCube’s muon–neutrino sensitivity to southern-sky neutrino sources, especially for Galactic sources that are not expected to produce a substantial flux of neutrinos above 100 TeV. In this work, the ESTES sample was applied for the first time to search for astrophysical sources of neutrinos, including a search for diffuse neutrino emission from the Galactic plane. No significant excesses were identified from any of the analyses; however, constraining limits are set on the hadronic emission from TeV gamma-ray Galactic plane objects and models of the diffuse Galactic plane neutrino flux.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2025

Dynamical Ion Correlations in Polymer and Membrane Electrolytes: Recent Advances and Open Questions

Ion transport in polymer and membrane electrolytes is often described using uncorrelated models like the Nernst–Einstein (NE) approximation, which link conductivity, transference number, and salt diffusivity to individual ion mobilities. However, growing evidence shows that dynamical ion correlations─arising from ion pairing, polymer-mediated coupling, confinement, and electrostatic interactions─can substantially alter macroscopic transport properties. The Onsager transport formalism offers a rigorous framework to quantify these effects via cross-correlated ionic displacements. This Perspective highlights how such correlations manifest across dry and hydrated systems, including salt-doped polymers, block copolymers, hydrogels, and ion-exchange membranes. Here, we examine key physical origins, the influence of salt chemistry and polymer structure, and open questions in measuring and controlling correlated motion. Embracing correlation-aware design enables new pathways to optimize conductivity, transference number, and selectivity─shifting the focus from individual ion mobility to engineered collective dynamics for next-generation electrolytes in energy storage and separations.

Charge transport↗

DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Li1−xNiO2 Many-body DMC Benchmark Dataset

The dataset contains all numerical data generated in support of the manuscript “Many‑body Benchmark of Electronic Charge and Spin Densities for Li1–xNiO2” (Journal of Chemical Theory and Computation, DOI: 10.1021/acs.jctc.5c02097, URL: https://pubs.acs.org/doi/10.1021/acs.jctc.5c02097). The materials included in this repository are: 1. Data files used to produce all figures and tables in the main manuscript and supporting information. 2. Benchmark density‑functional theory (DFT) datasets used for the charge‑ and spin‑density analyses. 3. Reference many‑body diffusion Monte Carlo (DMC) calculations and associated input/output files.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗

DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗

DOE OSTI2026

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗

DOE OSTI2025

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Insights Into Seismicity Associated With Flexibly Operating Enhanced Geothermal System From Real‐Time Distributed Acoustic Sensing

Enhanced Geothermal Systems (EGS) have the capacity to broaden the accessible resource pool for geothermal power generation. Traditionally viewed as a “baseload” resource, their flexible operation might also enable dispatchable load‐following generation and long‐term energy storage, aligning them with the evolving landscape of decarbonized electricity systems. However, increasing permeability and extracting energy during EGS operations can induce microseismic events; for many prior EGS efforts, some associated seismicity has been observed. While energetically beneficial, the flexibility of EGS operations prompts our inquiry into whether new types of operations will yield previously unseen seismicity patterns. We demonstrate the use of distributed acoustic sensing (DAS) with real‐time edge computing to monitor seismicity during a pilot test of a cyclically operated EGS facility at the Blue Mountain geothermal field. Our focus lies in uncovering seismicity insights from the real‐time microseismic catalog, particularly during load‐following dispatchability tests simulating flexible EGS operation. Here, we find that variations in pore pressure consistently correlate with seismicity, and that controlling pressure cycles during flexible operations appears to constrain microseismic activity during subsequent cycles. The spatio‐temporal evolution of microseismic clouds recorded during cyclic injection cycles fits diffusive models over our available observation period. Additionally, seismicity elevation lags behind pore pressure increases, likely due to pressure diffusion to the fracture system boundary. Through real‐time monitoring, we offer novel insights into seismicity associated with flexibly operating EGS. Our findings suggest that leveraging DAS and edge computing can inform EGS operations and help mitigate induced seismicity.

Chamarczuk, Michal [Rice Univ., Houston, TX (Unite↗

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗

DOE OSTI2026

Hydrogen uptake in graphite matrix at high temperature

Tritium management is a critical challenge for the next generation of nuclear reactors, such as Fluoride Salt Cooled High Temperature Reactors (FHRs) and High Temperature Gas-cooled Reactors (HTGRs), due to the higher production rate (up to 10,000 times) than conventional Light Water Reactors (LWRs). Graphitic materials employed as moderator, reflector, and fuel pebbles offer a potential pathway for tritium recovery by serving as a sink for tritium. Prediction of uptake capacity under reactor relevant conditions remains a challenge due to a lack of low partial pressure data and significant inter-grade variability of graphite. This study addresses these gaps by providing a comprehensive characterization of hydrogen (as a tritium surrogate) uptake and release behavior in the A3-3 graphite matrix (GM) used in fuel pebbles. Uptake measurements are performed at reactor relevant temperatures of 600- 800 °C, 1-200 Torr hydrogen pressure, and 15-120 min equilibration time, followed by thermal desorption spectroscopy up to 1100 °C. Uptake experiments at different equilibration times demonstrate the role of kinetics in hydrogen uptake, which can be modeled as a diffusion-with-trapping process. In the thermodynamics limit, the Sips adsorption model is shown to capture the uptake in A3-3 GM well. Our campaign provides a set of new results for hydrogen uptake in A3-3, including limiting uptake capacity at 600 °C, apparent diffusion coefficient at 600 °C, and the first estimates of the FHR/HTGR relevant (600 °C, 20 Pa partial pressure) equilibrium uptake capacity and time to saturation. Desorption data highlights a new site for hydrogen uptake, not observed in nuclear graphite, which we attribute to the non-graphitized binder. Using the Kissinger method, we estimate activation energy for release from the desorption peaks, confirming the activation energy for release from the basal planes and providing the first estimate for the activation energy of release from the binder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗