Search NASASearch

SEARCH · Search NASA

Results for “AI-driven simulation”

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.

439 records · Page 13

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES

A review of United States energy-only generator interconnection service policy and considerations for reform

Grid interconnection has emerged as a significant obstacle to the development of new electricity resources. There is growing interest in energy-only interconnection, which is an interconnection service option meant to allow the interconnection of new generators without ensuring their energy deliverability during all hours through the transmission system to customers. This approach potentially avoids upfront congestion-related transmission upgrades but could increase curtailment risk. Interest in energy-only interconnection is shaped by incomplete understanding of how interconnection policy functions in different jurisdictions, a knowledge gap that makes it difficult to determine how energy-only interconnection might be better used or re-designed. In this paper, we provide a regulatory review of energy-only interconnection in U.S. interconnection policy and practice, identifying jurisdictions in which rules are close to -or farther from-the theoretical concept of energy-only interconnection service. We find substantial jurisdictional differences in how energy-only interconnection is implemented, driven by differences in resource adequacy frameworks, real-time transmission operations, and state-level procurement practices. U.S. regulators have preferred local jurisdictional flexibility over federal prescription of interconnection study methods and procedures, which also contributes to differences among regions. Such findings raise fundamental questions about whether competition policies in electricity markets should extend beyond spot energy markets and into more prescriptive guidelines around interconnection rules and market entry. This paper sheds light on tensions that energy-only interconnection raises in allowing generators to access the transmission system on an as available basis and discusses how controlling thresholds for congestion-related network upgrades may be a barrier to electricity market entry.

Gorman, Will

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID

Aragonite saturation horizon variability along North Pacific seamounts and implications for deep-sea coral reefs

The 2014 discovery of living deep-sea coral reefs along the Northwest Hawaiian Islands (NWHI) and lower Emperor Seamount Chain (ESC), despite the North Pacific’s shallow aragonite saturation horizon (ASH) and high CaCO3 dissolution rates, underscores the need to understand the local seawater chemistry that allows these reefs to persist. We investigated seawater carbonate chemistry using discrete samples along the NWHI and ESC from two cruises ~1 year apart (08/26/21 – 09/26/21, 09/09/22 – 10/24/22). Across the two cruises, ASH depth difference ranged from 15 to 77 m. Since the Pacific ASH shoals by 1–2 m yr?¹, this long-term trend cannot explain the magnitude of ASH change observed between cruises. Temperature-salinity plots similarly indicate no major shifts in intermediate water masses that could account for these changes. Instead, ASH depth variability was primarily governed by localized biogeochemical processes, namely changes in intermediate water respiration and CaCO3 dissolution. Indicators for dissolution (TA*) and respiration (AOU) suggest changes in ASH depth were driven by changes in dissolution at the northern- and southern-most sites, whereas respiration exerted stronger control at central sites. Combining 2021 and 2022 data with data from 2014 – 2019 (excluding 2018) revealed high interannual ASH variability, by as much as >200 m at one site. Deep-sea coral reefs across the NWHI and ESC currently reside close to the ASH depth and likely experience seasonal and interannual shifts between under- and supersaturation. As ocean acidification progresses however, persistent exposure to undersaturated seawater could further threaten these critical deep-sea ecosystems.

coral reefs

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)

Catalytic disproportionation on carbon superstructures enables long-life, high-loading Li–S batteries

Electrocatalysis has been widely explored as an effective strategy to accelerate polysulfide (PS) conversion and suppress the shuttle effect in lithium–sulfur (Li–S) batteries. However, the underlying mechanisms remain elusive, and electrocatalytic reactions are inactive during cell resting. In this work, we reveal and quantitatively analyze a previously unrecognized sulfur reduction route (SRR) driven by catalytic disproportionation at the carbon cathode surface—fundamentally distinct from conventional electrocatalysis. Unlike conventional stepwise pathways, this SRR enables high-order polysulfides (Sₓ²⁻, x = 5–8) to directly convert into S₈ and Li₂S₂, bypassing low-order intermediates. This sulfur-reduction shortcut is systematically elucidated through high-performance liquid chromatography, revealing the intrinsic catalytic contribution of carbon frameworks and the dynamic evolution of PS species. We demonstrate that carbon superstructures (CSS-0.5), assembled from nanosheet subunits with abundant N/O functionalities and interconnected charge-migration channels, synergistically promote this catalytic process. Benefiting from these features, CSS-0.5 delivers superior electrochemical performance under practical conditions, enabling high sulfur loading (6.0 mg cm⁻²) pouch cells with 80.5% capacity retention over 210 cycles. This study provides the first quantitative evidence of electrocatalytic disproportionation in Li–S batteries, offering mechanistic insights and design principles for advanced sulfur cathodes.

25 ENERGY STORAGE

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS

Microstructure-dependent particulate populations in ODS steels revealed by single-particle ICP-TOF-MS

In this work, we evaluated the ability of single-particle inductively coupled plasma time-of-flight mass spectrometry (SP-ICP-TOF-MS) to characterize particulates generated from 300 series stainless steels and oxide dispersion–strengthened (ODS) steels. Particulates were produced using two different rotary tools (cutting disk and grinding bit) to assess the influence of particle generation on effluent particulates. Across all materials and preparation techniques, SP-ICP-TOF-MS detected a substantial population of particles containing the expected constituent elements, with average compositions closely matching the nominal alloy values. A secondary population of Y-bearing particles was evident in ODS steels; their formation is strongly favored in particulates generated using the grinding bit. Specifically, in PM2000, Y-bearing particles were predominantly colocated with Al, whereas in 14YWT, Y-bearing particles were less frequently Fe-free but showed stronger colocation with Ti. These trends are consistent with the known microstructural differences between Y–Al–O oxide particles in PM2000 and finer Y–Ti–O nanoclusters in 14YWT. Together, these results demonstrate that SP-ICP-TOF-MS can resolve distinct particulate populations in complex, multicomponent steels at the femtogram mass scale, enabling particle-resolved insight into microstructure-driven phase separation and preparation-dependent particle generation.

Brubaker, Zach E. [Oak Ridge National Laboratory (

Valorization of consolidated bioprocessing residues for bioplastics

This study demonstrates an organic solvent-free processing strategy to valorize consolidated bioprocessing (CBP) residues, from switchgrass and poplar biomass, into functional poly(butylene succinate) (PBS)-based biocomposites using high-shear homogenization (HSH). HSH transformed the switchgrass and poplar CBP residues (CBP-R) into fine, uniformly distributed particles and microfibers. The composites of PBS with homogenized switchgrass residues (H-CBP-R-SG) or homogenized poplar residues (H-CBP-R-P) at a 70/30 weight ratio exhibited improved processability and mechanical integrity, with the Young's modulus for the PBS/H-CBP-R-SG and PBS/H-CBP-R-P nearly doubling to 0.66 ± 0.07 GPa and 0.65 ± 0.04 GPa, respectively, compared to neat PBS (0.36 ± 0.02 GPa). Dynamic Mechanical Analysis (DMA) reveals a significant suppression of the tan δ peak magnitude, indicating that HSH-mediated physical activation facilitates stress transfer in composites typical of covalent chemical grafting systems. While the transition to a stiffness-dominated profile reduces ductility, the resulting composites exhibit the dimensional stability and resistance to thermal warping required for high-fidelity FDM 3D printing and injection molding. Beyond material performance, comprehensive techno-economic analysis (TEA) and life cycle assessment (LCA) confirmed that diverting CBP residues into composite products can improve the economic viability of the biorefinery without substantially increasing biorefinery global warming potential (GWP). At a 30 wt% blend ratio, incorporating residuals into PBS yielded a minimum selling price for the composite of $\$4.07$ per kg compared to the conventional bioplastic price of $\$5.00$ per kg. This approach aligns with circular bioeconomy principles by converting waste streams into value-added products. Furthermore, this innovative strategy addresses key challenges in bioplastic development, including cost, compatibility, and performance, while simultaneously advancing waste minimization strategies for sustainable manufacturing systems.

09 BIOMASS FUELS

Socio-Techno-Economic Feasibility of Deep Geothermal with EGS for Residential District Heating in the Northeastern United States

In this paper, we present socio-techno-economic simulation results for enhanced geothermal system (EGS) district heating systems in the northeastern United States. We applied the geospatial Distributed Geothermal Market Demand (dGeo) tool coupled with GEOPHIRES to evaluate the feasibility of geothermal deep direct-use, utilizing EGS reservoirs for providing residential heating with district heating systems. The results were combined with JEDI results to assess social impacts. Simulation results assuming doublet EGS reservoirs with target production temperature of 80 degrees C and newly constructed district heating systems indicate a wide range in system levelized cost of heat (LCOH) values ranging from ~$15/MMBtu to over $1,000/MMBtu, with most the attractive regions urban areas with high total thermal demand and high thermal demand density. We estimate about 60 GWth in potential installed capacity with LCOH values under $50/MMBtu and 100 GWth in potential capacity with LCOH values under $100/MMBtu, respectively.

15 GEOTHERMAL ENERGY

Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

active matter

Terahertz‐Nanoscale Visualization of the Microscopic Spin‐Charge Architecture of Colossal Magnetoresistive Switching

Resolving sub-10 nm spin switching and the associated terahertz (THz) electrodynamics during the colossal magnetoresistance (CMR) transition is a definitive frontier in reaching the fundamental spatial, temporal, and energy-dissipation limits of spin-electronics. Yet, simultaneous control of high magnetic field, cryogenic environment, and nanometer resolution has remained an elusive benchmark for THz nanoscopy, leaving the local THz dynamics of these transitions largely unexplored. Here, we overcome these limitations by utilizing a custom-built cryogenic magneto-THz scattering-type scanning near-field optical microscopy (cm-THz-sSNOM) to resolve the near-field THz spectroscopic evolution of the magnetic field-driven CMR transition in a manganite single crystal. Our measurements provide a nanoscale visualization of the THz conductivity, capturing the moment that magnetic-field-induced spin switching triggers the transition from an antiferromagnetic insulator to a ferromagnetic metal. An ellipsoidal near-field model reveals a multi-scale transition initiated by 1–2 nm isolated spin-flip sites at low magnetic fields, which coalesce into ∼15 nm conducting regions as the threshold field is approached. These results provide an nano-THz view of CMR switching, establishing an analysis framework for mapping spin–charge–lattice–orbit–coupled dynamics at spatial scales that transcend the nominal sSNOM resolution.

Haeuser, Samuel [Iowa State University, Ames, IA (

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 twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 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.

Sitaraman, Hariswaran

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE

Decorated Honeycomb Lattice and Successive Magnetic Transitions of the Delafossite Derivative K 4 Ni 5 Te 3 O 16

Here, we report the synthesis and characterization of the delafossite derivative K 4 Ni 5 Te 3 O 16 , which hosts a decorated honeycomb Ni 2+ lattice. Single-crystal X-ray diffraction, electron microscopy, and neutron diffraction establish a noncentrosymmetric Pmm2 crystal structure that combines regular honeycomb units with elongated motifs formed by Ni 2+ trimers along extended edges. Magnetic susceptibility reveals strong antiferromagnetic interactions, while heat capacity measurements identify two successive magnetic transitions at 30 and 12 K. Neutron diffraction shows that these transitions correspond to a progression from partial to full long-range magnetic order on distinct Ni sites, reflecting competition among the various exchange interactions. The resulting magnetic point group, mm2.1′, inherits the lattice polarity and permits high-order magnetoelectric coupling. Consistent with this symmetry, a quadratic magnetic-field-induced polarization is experimentally observed below the magnetic ordering temperature, likely arising from exchange-driven magnetostriction coupled to the polar lattice. These results establish K 4 Ni 5 Te 3 O 16 as an interesting platform for engineering B-site ordering in delafossite derivatives to realize complex quantum magnetism and functional responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Temperature as a control knob on spin-orbit coupling

Spin-orbit coupling (SOC) governs many physical phenomena. Through ab initio molecular dynamics simulations, Lu and Sun demonstrated that structural disorder at elevated temperatures substantially reduces the effective SOC contribution that stabilizes band inversion, driving a topological-to-normal-insulator transition in Bi2Se3. Their work identifies temperature as a meaningful control parameter for SOC-mediated topology and other properties.

Liang, Liangbo [ORNL] (ORCID:0000000311990049)