Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source Code”

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 451 records · Page 25

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of fast particle experiments in the first MAST Upgrade experimental campaigns

Abstract MAST-U is equipped with on-axis and off-axis neutral beam injectors (NBI), and these external sources of super-Alfvénic deuterium fast-ions provide opportunities for studying a wide range of phenomena relevant to the physics of alpha-particles in burning plasmas. The MeV range D-D fusion product ions are also produced but are not confined. Simulations with the ASCOT code show that up to 20% of fast ions produced by NBI can be lost due to charge exchange (CX) with edge neutrals. Dedicated experiments employing low field side (LFS) gas fuelling show a significant drop in the measured neutron fluxes resulting from beam-plasma reactions, providing additional evidence of CX-induced fast-ion losses, similar to the ASCOT findings. Clear evidence of fast-ion redistribution and loss due to sawteeth (ST), fishbones (FB), long-lived modes (LLM), Toroidal Alfvén Eigenmodes (TAE), Edge Localised Modes (ELM) and neoclassical tearing modes (NTM) has been found in measurements with a Neutron Camera (NCU), a scintillator-based Fast-Ion Loss Detector (FILD), a Solid-State Neutral Particle Analyser (SSNPA) and a Fast-Ion Deuterium- α (FIDA) spectrometer. Unprecedented FILD measurements in the range of 1–2 MHz indicate that fast-ion losses can be also induced by the beam ion cyclotron resonance interaction with compressional or global Alfvén eigenmodes (CAEs or GAEs). These results show the wide variety of scenarios and the unique conditions in which fast ions can be studied in MAST-U, under conditions that are relevant for future devices like STEP or ITER.

Physics↗

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

Abstract The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA’s versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

59 BASIC BIOLOGICAL SCIENCES↗

Advancing Neutrino Simulation Modeling with MARLEY: Insights from the NNSA-MSIIP Internship

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. However, there is no experimental data to validate the current event generator MARLEY (Model of Argon Reaction Low Energy Yields) that can model supernova neutrinos. To validate the model, I developed a new muon capture feature within the MARLEY simulation framework by coding key functions in C++, Python, and ROOT. I generated and analyzed one million simulated muon capture events and comparing the results to experimental data. To improve the model, I worked on optimizing the model’s parameters to improve its precision using statistical methods.

Wong, Baker [Fermilab]↗

Demonstration of Optimal Benchmark Selection Website and Validation of the q c Coverage Metric Using HEU-SOL-THERM-013-003 Experiment

In the work documented in this interim report, the experiment selection toolkit web site was demonstrated and q C coverage metric methodology was validated for IEU-MET-FAST-002-001, MIX-COMP-THERM 004-004, and HEU-SOL-THERM-013-003 experiments. 𝑞 𝐶 is an information-theoretic measure based on mutual information that quantifies the ability of candidate benchmark experiments to reduce the bias and uncertainty of a target criticality safety application. The metric and an accompanying open-source Python toolkit with a web-based interface were tested against a benchmark set of 425 experiments drawn from the International Criticality Safety Benchmark Evaluation Project Handbook. The interface is hosted at https://edim.covdef.com. It accepts sensitivity data files produced by the TSUNAMI-IP module of the SCALE code system and supports both (i) deterministic analysis using the ENDF/B-VII.0 covariance library and (ii) stochastic analysis based on user-supplied keff samples. Demonstrations on representative applications across a range of material composition, spectrum, and form show that q C -guided benchmark selection achieves greater uncertainty reduction with fewer experiments and yields more stable posterior bias and uncertainty estimates than traditional similarity coefficient ( c k )–based selection, while also capturing valuable low-ck experiments that one-to-one metrics overlook.

Abdel-khalik, Hany S. [Indiana Univ.-Purdue Univ. ↗

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗

HALLUFIELD: DETECTING LLM HALLUCINATIONS VIA FIELD-THEORETIC MODELING

A research-focused Python package that implements our hallucination-detection method for large language models(LLMs). The code computes stability signals from LLM predictions across various hyperparameters sweeps and combines free-energy/entropy–style metrics to flag likely hallucinations, with tunable thresholds for batch scoring. The repo includes evaluation scripts, config files, and example notebooks to reproduce benchmark results and ablations; it depends on standard open-source libraries (PyTorch, Hugging Face) and runs on CPU/GPU. The repository examples only use public models/datasets only.

Bhattarai, Manish [Los Alamos National Labs]↗

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗

PlantCV v4: Image analysis software for high‐throughput plant phenotyping

PlantCV is an open-source Python project aimed at developing tools to address a range of image-based, plant phenotyping questions. PlantCV has been used for more than 10 years to automate trait collection from image data, and the newest release, PlantCV version 4, continues to lower the barrier to entry for users without substantial coding experience through extensive example use-case tutorials and simplified installation. In addition to usability, we document added functionality since the release of PlantCV v2, including support for more image types such as fluorescence, thermal, and hyperspectral data. Finally, we describe the development of a new subpackage focused on morphological trait measurements like leaf angle, and demonstrate its utility as compared to more manual methods of data collection.

Schuhl, Haley [Donald Danforth Plant Science Cente↗

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Outflowing photoionized plasma in Circinus X-1 using the high-resolution X-ray spectrometer Resolve onboard XRISM and the radiative transfer code cloudy

High-resolution X-ray spectroscopy is a key to understanding the mass inflow and outflow of compact objects. Spectral lines carry information about the ionization, density, and velocity structures through their intensity ratios and profiles. They are formed in non-local thermodynamic equilibrium conditions under the intense radiation field from the compact objects; thus radiative transfer (RT) calculation is a requisite for proper interpretations. We present such a study for a low-mass X-ray binary, Circinus X-1, from which the P Cygni profile was discovered using the X-ray grating spectrometer onboard Chandra. We observed the source using the X-ray microcalorimeter onboard XRISM at an orbital phase of 0.93–0.97 and revealed many spectral features unidentified before: the higher series transitions ($n{\rightarrow }1$; $n>2$) of highly ionized (H- and He-like) S, Ca, Ar, and Fe in emission and absorption, the Fe K$\alpha$ and K$\beta$ inner-shell excitation absorption of mildly ionized (O- to Li-like) Fe, and resolved fine-structure level transitions in the Fe Ly$\alpha$ and He$\alpha$ complexes. They blend with each other at different velocity shifts on top of apparently variable continuum emission that changed its flux by an order of magnitude within a 70 ks telescope time. Despite such complexity in the observed spectra, most of them can be explained by a simple model consisting of the photoionized plasma outflowing at $\sim$300 km s$^{-1}$ and the variable blocking material in the line of sight of the incident continuum emission from the accretion disk. We demonstrate this with the aid of the RT code cloudy for the line ratio diagnostics and spectral fitting. We further constrain the physical parameters of the outflow and argue that the outflow is launched close to the outer edge of the accretion disk and can be driven radiatively by being assisted by the line force calculated using the RT simulation.

X-rays: binaries↗

Analysis of a Prototype Multi-Detector Fast-Neutron Radiography Panel

A multi-detector fast neutron radiography panel was built using the previous work on scalable neutron radiography using the IDEAS ROSSPAD readout module. A new aluminum housing was built to accommodate a large number of detectors tiled together. Additional changes to startup and processing code were made to operate the detector as one cohesive unit. Spatial resolution of the full panel using Cs-137 gammas was reported to be 0.42 line pairs per centimeter at 90% MTF and 2.09 line pairs per centimeter at 10% MTF. Here, three neutron radiographs generated using a Cf-252 fission neutron source were used to determine the spatial resolution of the panel for neutrons. The experiments had 90% MTF values of 0.24, 0.3, and 0.27 line pairs per centimeter and 10% MTF values of 1.30, 1.46, and 1.40 line pairs per centimeter. An example neutron radiograph was also used to prove that the radiography panel can perform true neutron radiography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

CHARACTERIZATION OF FATIGUE BEHAVIORS OF NOTCHED 316L DED AM SPECIMENS

ASME Codification of Additive Manufacturing • Integration of AM into ASME Codes and Standards • The ASME goal is to have AM requirements in ASME Code Cases preceding the 2025 Edition. • The ASME Special Committee on AM has drafted criteria for two Code Cases for Additive Manufacturing. • AM Construction of Pressure Equipment using the Direct Energy Deposition Process with Wire Feedstock. • Includes Gas Metal Arc Welding. • Time-independent material properties. • Status - Criteria endorsed by AM Committee. • AM Construction of Pressure Equipment using the Powder Bed Fusion AM Process. • Includes Laser and Electron Beam Energy Sources. • Austenitic and Nonferrous materials. • Time-independent material properties. • Status – Approval ballot circulating to the AM Committee.

Krentz, Timothy M.↗

Demonstration and Automation of Reflected Target Optical Measurement for Heliostats

Accurate optical surfaces are a primary driver of concentrated solar power plant performance. Errors in pointing and tracking mirrors, the canting of individual mirror facets, and the surface slope of the mirror itself can be caused by errors during assembly, transportation, wind loading, gravity, and many other sources. The tools that exist to measure these error sources today largely rely on fringe deflectometry (SOFAST, QDec, others), or photogrammetry with targets attached to the mirror surface. Since 2022, NREL has been developing a measurement method called the Reflected Target Non-intrusive Assessment (ReTNA) system. This system differs from most established methods in that we perform deflectometry with a pattern of coded targets, identified in space with photogrammetry. Reflected target systems have several advantages over traditional fringe deflectometry systems. Firstly, they can be operated in bright or ambient lighting, a challenge for fringe systems that use a projector and screen. Reflected target systems also can use a much lighter and less expensive target than projector-based systems. Lastly, 2D slope measurement can be solved from a single image, which leads to several advantages for accommodating faster measurements and smaller sized targets. These advantages make ReTNA particularly well-suited for applications where there are space or lighting constraints, like performing heliostat quality assurance on an assembly line. It's also useful when a lightweight, flexible system is needed, like for heliostat developers to quickly measure a new heliostat design at different orientations, to observe gravitational effects on the mirror surface shape. In the last year, significant improvements were made to this tool to make it more useful for these applications. These improvements were focused around validation of the ReTNA measurement system, and automation of the setup and measurement process. First, we present an improved ReTNA layout, for use on the heliostat assembly line. Next, we detail the various changes to the ReTNA software and computer vision methods to automate data collection in this new setup, and lessons learned from this process. The goal with this new setup is to perform a full heliostat surface characterization without removing the mirror from the assembly line. Lastly, we share results from several ReTNA validation studies undertaken over the last year. These include repeated ReTNA measurement on demonstration mirror facets, comparisons with other optical measurement tools, and some studies aimed at quantifying the uncertainty of ReTNA measurement under various constraints (mirror-target spacing, camera resolution, etc.). These results are compared with 2024 HelioCon performance targets, and our planned next steps for the ReTNA measurement system are presented.

CSP↗

Radiological Releases from Novel Fuel Forms in Advanced Reactors During Severe Accidents for Consequence Analyses

Various advanced reactor developers are exploring the potential for reductions in the size of physical security forces and emergency planning zones. These reductions are based on robust fuel forms and inherently safe reactor designs. However, such reductions in physical protection measures could increase the risk of sabotage. To assess the possibility of reducing these measures, sabotage-induced radiological consequence analyses were carried out. These analyses considered accident scenarios that were beyond design basis accidents and overly conservative (Shah, 2025a; Shah, 2025b; Shah and Hartanto, 2026), yielding very large release fractions. These fractions, which can be used to evaluate physical protection and emergency planning requirements, have been crudely determined and applied as demonstrations for a sodium-cooled fast reactor (SFR) (Shah and Hartanto, 2025a), a high-temperature gas-cooled reactor (HTGR) (Shah and Hartanto, 2025b), a heat pipe–cooled reactor (HPR) (Shah and Hartanto, 2025c), and a molten salt–cooled reactor (MSR) (Shah et al., 2026). A Sandia National Laboratories (SNL) team used MELCOR—a fully integrated severe accident analysis code—to demonstrate the code’s capability to analyze advanced (i.e., not light water–cooled) reactors (including a fluoride salt–cooled high-temperature reactor [FHR]) and calculate radiological releases to the environment during severe accidents (Wagner et al., 2022a, 2022b, 2022c, 2023a, and 2023b). Although the analyses were carried out to demonstrate MELCOR’s growing capability, the release source terms were estimated for advanced reactors, providing valuable insights into the accident progression and radiological releases. These findings from prior SNL studies, including estimated source terms and related sensitivity studies, were leveraged to derive source terms for postulated sabotage-induced accidents. Insights from these sensitivity studies informed the scaling of SNL’s estimated source terms for the defined accident scenarios. The derived release fractions for the severe accident scenarios for the respective reactor designs can be used to perform more nuanced dose consequence analyses to evaluate the reactors’ physical protection and emergency planning zone requirements. These analyses are in accordance with the risk-informed, performance-based approach proposed under 10 CFR Part 53. This study builds on the prior source term analyses and associated sensitivity studies by SNL to derive time-dependent and design-informed release fractions. Section 2 describes the diverse advanced reactor designs analyzed by the SNL team. Section 3 discusses the severe accident analyses, the release fractions calculated, and the limitations and assumptions of the demonstration project. Section 4 presents the release percentages derived for the hypothetical sabotage-induced severe accidents at the advanced reactors. Section 5 summarizes the study’s findings and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗