Search NASA⌕ Search

SEARCH · Search NASA

Results for “O CODES”

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 415 records · Page 23

Status of SAS4A/SASSYS-1 Software Development and Application (FY2024)

SAS4A/SASSYS-1 is a simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. With its origin as SAS1A in the late 1960s, the SAS series of codes has been under continuous use and development for over fifty years. It has been identified as a critical element of safety analysis capabilities for the U.S. Department of Energy and is utilized within industry to perform the transient safety analyses required to support the licensing of Liquid Metal-cooled Fast Reactors (LMFRs). This report summarizes the code development and update activities carried out during FY2024. In FY2024, programmatic activities focused on key improvements to software useability, such as enhanced user interfaces for reactivity feedback modeling, improvements in stability/useability of the Code Manual, and improvements to the acceptance testing infrastructure, including automation of acceptance testing and generation of the Acceptance Testing Report. To support end user applications, an open training was held, a semi-public forum was maintained, and a practical benchmarking and validation matrix was developed which allowed limitations of existing testing capabilities to be assessed. The existing fuel models were also enhanced with improved modeling capabilities and testing for the oxide and annular fuel models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗

Generation of meter-scale hydrogen plasmas and efficient, pump-depletion-limited wakefield excitation using 10 GeV electron bunches

High repetition rates and efficient energy transfer to the accelerating beam are important for a future linear collider based on the beam-driven plasma wakefield acceleration scheme (PWFA-LC). This paper reports the first results from the Plasma Wakefield Acceleration Collaboration (E300) that are beginning to address both of these issues using the recently commissioned FACET-II facility at SLAC national accelerator laboratory. Here we have generated meter-scale hydrogen plasmas using time-structured 10 GeV electron bunches from FACET-II, which hold the promise of dramatically increasing the repetition rate of PWFA by rapidly replenishing the gas between each shot compared to the hitherto used lithium plasmas that operate at 1–10 Hz. Furthermore, we have excited wakes in such plasmas that are suitable for high gradient particle acceleration with high drive-bunch to wake energy transfer efficiency- a first step in achieving a high overall energy transfer efficiency. We have done this by using time-structured electron drive bunches that typically have one or more ultra-high current (> 30 kA) femtosecond spike(s) superimposed on a longer (~0.4 ps) lower current (< 10 kA) bunch structure. The first spike effectively field-ionizes the gas and produces a meter-scale (30–160 cm) plasma, whereas the subsequent beam charge creates a wake. The length and amplitude of the wake depends on the longitudinal current profile of the bunch and plasma density. We find that the onset of pump depletion, when some of the drive beam electrons are nearly fully depleted of their energy, occurs for hydrogen pressure ≥ 1.5Torr. We also show that some electrons in the rear of the bunch can gain several GeV energies from the wake. These results are reproduced by particle-in-cell simulations using the QPAD code. At a pressure of ~2Torr, simulation results and experimental data show that the beam transfers about 60% of its energy to the wake.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multigrid reduction in time with Richardson extrapolation

The advent of exascale computing will leave many users with access to more computational resources than they can simultaneously use, e.g., billion-way parallelism. In particular, this is true for time-dependent simulations that limit parallelism to the spatial domain. One method to add parallelism in time to existing simulation codes and thus take advantage of ever larger compute resources is Multigrid Reduction in Time (MGRIT). The goal is to achieve a smaller time-to-solution through parallelism in time. In this paper, MGRIT is enhanced with Richardson extrapolation in a cost-efficient way to produce a parallel-in-time method with improved accuracy. Overall, this leads to a large improvement in the accuracy per computational cost of MGRIT.

97 MATHEMATICS AND COMPUTING↗

Exascale models of stellar explosions: Quintessential multi-physics simulation

The ExaStar project aims to deliver an efficient, versatile, and portable software ecosystem for multi-physics astrophysics simulations run on exascale machines. The code suite is a component-based multi-physics toolkit, built on the capabilities of current simulation codes (in particular Flash-X and Castro), and based on the massively parallel adaptive mesh refinement framework AMReX. It includes modules for hydrodynamics, advanced radiation transport, thermonuclear kinetics, and nuclear microphysics. The code will reach exascale efficiency by building upon current multi- and many-core packages integrated into an orchestration system that uses a combination of configuration tools, code translators, and a domain-specific asynchronous runtime to manage performance across a range of platform architectures. The target science includes multi-physics simulations of astrophysical explosions (such as supernovae and neutron star mergers) to understand the cosmic origin of the elements and the fundamental physics of matter and neutrinos under extreme conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic state switching in FeGa microstructures

This work demonstrates that magnetoelectric composite heterostructures can be designed at the length scale of 10µms that can be switched from a magnetized state to a vortex state, effectively switching the magnetization off, using electric field induced strain. This was accomplished using thin film magnetoelectric heterostructures of Fe 81.4 Ga 18.6 on a single crystal (011) [Pb(Mg 1/3 Nb 2/3 )O 3 ] 0.68 -[PbTiO 3 ] 0.32 (PMN-32PT) ferroelectric substrate. Here, the heterostructures were tripped from a multi-domain magnetized state to a flux closure vortex state using voltage induced strain in a piezoelectric substrate. FeGa heterostructures were deposited on a Si-substrate for superconducting quantum interference device magnetometry characterization of the magnetic properties. The magnetoelectric coupling of a FeGa continuous film on PMN-32PT was characterized using a magneto optical Kerr effect magnetometer with bi-axial strain gauges, and magnetic multi-domain heterostructures were imaged using x-ray magnetic circular dichroism—photoemission electron microscopy during the transition to the vortex state. The domain structures were modelled using MuMax 3 , a micromagnetics code, and compared with observations. The results provide considerable insight into designing magnetoelectric heterostructures that can be switched from an ‘on’ state to an ‘off’ state using electric field induced strain.

36 MATERIALS SCIENCE↗

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CROCUS Optical All Precipitation Gauge Data at Argonne National Laboratory Prairie Site

The APG (Optical Scientific Inc. All-Precipitation Gauge 815-DS) dataset contains one-minute measurements of precipitation rate, precipitation accumulation, air temperature, and present weather detection, both in 4680 format and decoded. Data were collected at the Argonne Testbed for Multiscale Observational Science (ATMOS), a 20-acre prairie site at Argonne National Laboratory in Lemont, Illinois. The data is presented as daily NetCDF (.nc) files, each containing approximately 24 hours of observations. Files follow the naming convention of: the project (CROCUS), location (atmos), instrument name (apg), data level (raw, a1), and date (year, month, day). The NetCDF format can be accessed using common scientific software such as Python using xarray, netCDF4 or act-doe.

54 ENVIRONMENTAL SCIENCES↗

Vertical control of DIII-D discharges with strong negative triangularity

Through predictive modeling validated by a series of experiments on DIII-D, the vertical stability of low β diverted plasmas with strong negative triangularity (NT) ($\delta\sim-0.6$) is assessed. As a result of their unique magnetic geometry, NT plasmas feature larger Shafranov shifts and more elongated inner flux surfaces than positive triangularity counterparts, typically leading to enhanced vertical instability growth rates. However, coupling with the non-conformal vessel DIII-D wall reduces these growth rates to controllable values, providing a path forward for stabilizing strongly NT plasma with a diverted geometry. These discharges are used to validate GSdesign (part of the TokSys code suite) stability calculations in strong NT plasmas on DIII-D, with errors of no more than ~20% observed between modeled and experimentally measured growth rates. Additions of diagnostic noise and power supply tuning to the TokSys model are needed to accurately capture the time dependence of DIII-D NT discharges, assisting with the design of control schemes specific to the DIII-D poloidal field coils. Lastly, implications of these results on a future NT reactor are briefly described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of suprathermal ions on neutron yield in the pre-DT phase of ITER operation

An assessment of neutron production during the pre-fusion power-operation (PFPO) phase has been carried out for a representative set of plasma scenarios predicted by the ITER Research Plan. A range of heating systems, namely neutral beam injection (NBI) (hydrogen), electron cyclotron resonance heating (ECRH), and ion cyclotron resonance heating (ICRH) are planned to be used for PFPO studies in helium, hydrogen, and mixed hydrogen–helium plasmas. Fast ions (protons and 3He) originating from NBI and ICRH systems can increase neutron production in PFPO plasmas by directly interacting with intrinsic Be impurities or through secondary processes, as also evidenced at JET. The generation of fast ions in ITER PFPO scenarios has been modelled using the ASTRA-NBI and TORIC-SSFPQL codes. A significant impact of the synergy between hydrogen NBI and hydrogen-minority ICRH on neutron production in helium plasmas is reported. In addition, the stability of the toroidicity-induced Alfvén eigenmodes (TAE) is analyzed for PFPO plasmas with a high pressure of suprathermal ions and a weak reversed shear. The possible impact of sawtooth oscillations and TAEs on neutron production is discussed, based on a linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Carbon dioxide, water vapor and methane soil efflux (soil respiration) in a Pinus palustris restoration site in Georgetown, SC

This dataset contains processed data from a combination of survey flux chambers and long-term automated flux chambers. Biweekly soil flux measurements were conducted from June 2023 through December 2025 at a longleaf pine restoration site in Georgetown, SC. Processed, QAQC’d data can be found in the file: 1_DATA_ESS_DOE_HR_RS_HB3_QAQC_Survey_Data_20260223.csv. Raw and working data files (.json, & .81x format) from LI-COR equipment are included for reference and can be accessed using SoilFluxPro software. CSV metadata files describe the raw data and modifications made using SoilFluxPro v5 and Matlab R2024b, as well as formatting and units for processed CSVs. Matlab code is included for reading in the processed CSVs. This research was performed as part of the project: “Improving models of stand and watershed carbon and water fluxes with more accurate representations of soil-plant-water dynamics in southern pine ecosystems”, which examines in part the effects hydraulic redistribution on soil efflux of carbon dioxide, water vapor and methane, as well as soil moisture and temperature in a southern pine ecosystem with sandy soils and high water table.

CARBON DIOXIDE FLUX↗

Development Of Thermodynamic and Kinetic Simulation Tools and Testing Procedures for Enhanced Durability of Concrete Containing Industrial By-Products

This project developed screening tools that enable evaluation of alternative cementitious binders that create concretes to significantly reduce energy and emissions while remaining cost competitive on both initial and long-term costs. The team began with the viewpoint that acceptance of new cementitious binder products has a substantially greater chance of successful implementation when capital investment is not excessive and the end product has customers that have experience and resources to use this. As such, a binder system that is based on portland cement with blended using industrial by-products (alternative cementitious materials) has the potential for dramatic and meaningful impact. The team has focused on developing implementable solutions in specifications and current practice. This however requires three main factors: 1) ability to screen byproducts and alternative materials for success, 2) ability to ‘treat’ materials chemically to enhance kinetics, and 3) ability to provide predictions of performance of both binders and concrete from first principles. The project developed/refined a state of the art and scientifically based screening test for SCM called the pozzolanic reactivity test. The team developed kinetic models to simulate these materials as well as experimental approaches to alter selected reactions. Simulation tools were developed that enable the performance of concrete to be predicted based on the chemistry and reactivity of the cement and alternative SCM. Specifically, this project: • enhanced the kinetic reactivity models for use in multi-scale computational programs that use thermodynamics to predict reaction products. 2 • developed scaling models to extend thermodynamic modeling to link these models with pore structure. This enables strength, transport property, and coupled transport prediction. • developed tools to predict performance of cementitious materials using the pozzolanic reactivity test and chemical composition. The predicted properties are consistent with AASHTO R101 and the CEB-FIP model code and can be measured using associated test procedures. • demonstrated mechanical and fracture based modeling tools that thermodynamic predictions and inputs to predict concrete service life. These results have been used to demonstrate the value of enabling specifications to include ASTM C 595 cement as well ASTM C150 cement. In addition, these products are being used to expedite the evaluation of alternative SCM to aide in determining which materials have potential value and what ‘compositions’ of blended cements merit further investment.

42 ENGINEERING↗

Understanding Binary Systems—a Comparison between COSMIC and MESA

Abstract We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion temperature clamping in Wendelstein 7-X electron cyclotron heated plasmas

The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

General Relativistic Neutrino-driven Turbulence in One-dimensional Core-collapse Supernovae

Abstract Convection and turbulence in core-collapse supernovae (CCSNe) are inherently three-dimensional (3D) in nature. However, 3D simulations of CCSNe are computationally demanding. Thus, it is valuable to modify simulations in spherical symmetry to incorporate 3D effects using some parametric model. In this paper, we report on the formulation and implementation of general relativistic neutrino-driven turbulent convection in the spherically symmetric core-collapse supernova code GR1D . This is based upon the recently proposed method of Simulated Turbulence in Reduced Dimensionality (STIR) in Newtonian simulations from Couch et al. (2020). When the parameters of this model are calibrated to 3D simulations, we find that our general relativistic formulation of STIR requires larger turbulent eddies to achieve a shock evolution similar to the original STIR model. We also find that general relativity may alter the correspondence between progenitor mass and successful versus failed explosions.

Astronomy & Astrophysics↗

Data for Myers-Pigg et al. (2026), "Short-term coastal forest responses to a hurricane-scale freshwater and saltwater flooding experiment"

Coastal upland forests are exposed to intensifying precipitation regimes and sea level rise, increasing tree mortality and transforming these coastal forests into wetland ecosystems. Despite these well-known risks, the differing degrees to which hydrological, biogeochemical, and biological components of upland forests respond to novel salinity exposure is relatively unknown. The Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experiment decouples two distinct disturbances associated with hydrological extremes: (1) flooding from heavy precipitation and (2) exposure to saline conditions from storm surge. This dataset includes data reported in Myers-Pigg et al. (2025), which analyzed data from the first TEMPEST flooding treatment in 2022. This includes: - Colored dissolved organic matter in porewaters - Soil temperature and oxygen - Groundwater temperature and chemistry - Dissolved organic carbon concentrations in porewaters - Soil-to-atmosphere CH4 and CO2 fluxes - Soil temperature, water content, and electrical conductivity - Root-influenced CH4 and CO2 flux - Tree sap flow velocity - The R analytical code and documentation about the computational environmental in which it was run (the "sessionInfo.txt" file) All data files are plain-text comma separated value (CSV) and no special software is required to read them.

54 ENVIRONMENTAL SCIENCES↗

LOCOMOTIVES - Comprehensive Impact and Cost Assessment Framework of Carbon Lowering Approaches for the US Rail Freight System

The goal of this project is to develop a tool to aid railroads and other stakeholders assess and approach the decarbonization of freight rail operations by identifying new, viable low-carbon energy storage and conversion systems for future locomotive systems and how they should be deployed on the existing US freight rail network. In the first quarter, the project focused on collecting data, establishing a simulation workflow, and engaging industry through the creation of the Industry Advisory Board (IAB). In the second quarter, the project focused on selecting fuel pathways and powertrain technologies, setting performance targets, conducting a techno-economic analyses, and developing the simulation framework that would serve as the backbone of the future toolhead. The third quarter involved developing an industry-oriented interactive dashboard powered by a five-step sequential framework, as well as holding industry advisory board meetings as per the initial technology-to-market plan. In the remaining project quarters, the NUFRIEND dashboard were fine-tuned with the help of IAB member feedback and in-depth scenario analyses were conducted to support the techno-economic analysis of energy sources. Additionally, dashboard documentation, project insights, and open-source code on GitHub were prepared and released. Throughout the project, the team completed testing and analysis of all model components, integrated all initial test scenarios, and conducted stakeholder engagement. Lower-carbon drop-in fuels can be deployed as admixtures and are considered uniform across the network at a desired penetration rate, while hydrogen and battery-electric technology deployment poses a more complex problem as they require significant investments to be made in the siting of refueling/charging facilities and the replacement of locomotive fleets. Thus, strategies for locating and sizing refueling/charging facilities on a railroad’s network to meet their energy demands were developed to inform deployment decisions. To address this challenge, the Northwestern University Freight Rail Infrastructure & Energy Network Decarbonization (NUFRIEND) framework presents a five-step sequential framework to select O-D paths, locate facilities, reroute flows, size facilities, and evaluate the deployment for alternative energy sources that require locomotive powertrains to be converted and new refueling infrastructure to be deployed. The NUFRIEND Framework is an industry-oriented tool for simulating the deployment of new energy technologies across the US freight rail network. The framework provides a comprehensive network-level optimization and scenario simulation tool for decarbonizing the freight rail sector, addressing the uncertainties surrounding technological developments by supporting sensitivity analyses for different operational and technological parameters through a transparent and flexible input module. It offers practical alternatives to diesel locomotives and can be applied for any railroad considering the specific network structure and freight demand, outputting evaluation metrics for the associated emissions and costs relative to diesel operations. A number of relevant simulation scenarios were run and analyzed for key insights on the value of different alternative technologies for freight rail decarbonization. The project developments and findings have been presented at numerous conferences and events.

08 HYDROGEN↗