Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code development and validation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SAS4A/SASSYS-1 Validation with FFTF Tests

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. In preparation for the Passive Safety Testing program, a series of individual reactivity feedback tests were performed which were designed to simulate and validate specific features and reactivity feedbacks of the FFTF core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Multiscale, mechanistic modeling of irradiation-enhanced silver diffusion in TRISO particles

Tristructural isotropic (TRISO) particles are under consideration for use in several proposed advanced nuclear reactor concepts. The silicon carbide (SiC) layer in TRISO acts as a barrier to prevent the release of the fission products. However, despite remarkable retention, silver (Ag) release has been observed from intact particles, which requires investigation since the Ag isotope ( 110m Ag) has a long half-life. Previous work focused on developing a multiscale, mechanistic model for Ag diffusion accounting for temperature and microstructure effect and has been successfully validated. In this work, we expand the previous model to account for irradiation-enhanced Ag diffusivity in SiC and improve its accuracy over a wider grain size and temperature ranges relevant for advanced reactor conditions. A temperature, grain size, and flux dependent diffusivity is therefore derived using the mesoscale code MARMOT and implemented in the fuel performance code BISON. The irradiation-enhanced Ag diffusivity in SiC is compared against experimental data and validated using BISON against Ag release measurements from the Advanced Gas Reactor Fuel Development and Qualification Program (AGR-1 and AGR-2). Herein, we quantify the impact of SiC grain size, irradiation, and temperature on Ag release. In agreement with previous studies, we find accounting for SiC grain size improves agreement between BISON predictions and experimental observations for most cases. In conclusion, we also find that accounting for irradiation improves agreement for cases where Ag release was underestimated, but the impact was less significant than accounting for microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Loss of Flow Conditions in a Modern Pool-type SFR and In-Pile Experiment Design

The metallic fuel safety performance under unprotected design-basis transients is a key consideration for the deployment of advanced sodium fast reactors (SFRs). Reliable data are needed to validate advanced safety codes, reduce uncertainty in cladding failure thresholds, and strengthen confidence in licensing approaches. To address this need, this report develops blueprints for a conceptual sodium loss-of-flow (LOF) experiment in the Mk-IIIR loop at the Transient Reactor Test Facility (TREAT), providing the technical foundation for future integral testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

VARI3D & PERSENT: Perturbation and Sensitivity Analysis (Revision 5)

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Study of the heat pipe vapor flow using the MOOSE application Sockeye and the CFD code Nek5000

Heat pipes are efficient heat transfer devices used in various applications, including nuclear microreactors. Heat pipe-cooled microreactors offer advantages in size and cost, which can significantly accelerate their deployment and adoption. Developing accurate models of heat pipes is crucial for a heat pipe-cooled microreactor design and operation. In this work a verification and validation of the heat pipe application Sockeye is presented. Results are compared against RANS simulations, performed with the code Nek5000, and the SPHERE experiment, conducted at INL.

42 ENGINEERING↗
DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Staged Z-Pinch Modeling with Hydra and Chicago Codes (CRADA Final Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC (“Contractor”), as manager and operator of Lawrence Livermore National Laboratory (“LLNL”) and Magneto Inertial Fusion Technologies, Inc. ("MIFTI" or "Participant"), to develop staged Z-pinch (SZP) modeling with HYDRA and CHICAGO codes. The contractor performed numerical simulations of the SZP to reproduce and validate Participant’s MACH2 simulations. Contractor used the state-of-the-art, radiation-magneto-hydrodynamic code HYDRA, as well as the hybrid particle-in-cell code CHICAGO. This work helped assess the physical and numerical uncertainties associated with the simulations and the overall fusion scheme. HYDRA is a well-established and validated code. This work has helped to independently assess the SZP concept to give confidence to the Participant and to the broader fusion community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ultrafast Surface Phosphor Thermometry for Pulsed-power and Hostile Environments

Modern concepts for next generation pulsed power (NGPP) are slated to deliver up to ten times the energy of Z today. An increase of this magnitude is concerning insofar that Z currently exhibits sizable amounts of inner magnetically insulated transmission line (MITL) loss current on the order of 5-10%. Loss phenomenon in these systems are complex and electrode heating and subsequent thermal desorption are a leading cause. Rapid heat-driven thermal desorption of contaminants scales as the square of the current. Therefore, even a modest doubling of drive current would yield an ~ 4X in non-linear surface electrode heating, quickening thermal desorption-based current loss. Exacerbating these physics is a current inability to measure ultra fast heating rates (>20°C/ns), which are paramount to benchmarking and code validation critical to NGPP design – as an empirical approach is not viable. Therefore, Ultrafast Photoluminescent Surface Heating Optical Thermometry (UP-SHOT) was developed as a new diagnostic for measurement of GHz-scale electrode heating. The discovery of UP-SHOT leveraged expertise in Engineering Science, Material Science, Pulsed-Power, and the Center for Integrated Nanotechnologies. This report includes information on: 1) The preparation of zinc oxide (ZnO) films, characterization, post-deposition treatments 2) Time-resolved photoluminescence at elevated temperatures and thermographic sensitivity

36 MATERIALS SCIENCE↗
DOE OSTI2024

RELAP5-3D validation studies based on the High Temperature Test facility

In the spring and summer of 2019, experiments were conducted at the High Temperature Test Facility (HTTF) that form the basis of an upcoming high-temperature gas-cooled reactor (HTGR) thermal hydraulics (T/H) benchmark. HTTF is an integral effects test facility for HTGR T/H modeling validation. This paper presents RELAP5-3D models of two of those experiments: PG-27, a pressurized conduction cooldown (PCC); and PG-29, a depressurized conduction cooldown (DCC). These models used the RELAP5-3D model of HTTF originally developed by Paul Bayless as a starting point. The sensitivity analysis and uncertainty quantification code, RAVEN was used to perform calibration studies for the steady-state portion of PG-27. Here we developed four PG-27 calibrations based on steady-state conditions. These calibrations all used an effective thermal conductivity equal to 36 % of the measured thermal conductivity, but they differed with respect to the frictional pressure drops and radial conduction models. These models all captured the trends in steady-state temperature distributions and transient temperature behavior well. All four calibrations show room for improvement in predicting the transient temperature rise. The smallest error in temperature rise during the transient was a 21 % underprediction, and the largest was a 48 % underprediction. The errors in transient temperature rise are largely a result of a mismatch in power density between the RELAP5-3D model and the experiment due to the location of active heater rods along the boundary between heat structures in the model. The best of these calibrations was applied to PG-29 to model the DCC. Once again, temperatures during the transient were underpredicted but trends in temperature were captured. The RELAP5-3D model captured trends in the data but could not reproduce measured temperatures exactly. This result is not attributed to deficiencies in the experimental data or to RELAP5–3D itself. Rather, this result likely arises due to the some of the assumptions and decisions made when the RELAP5-3D model was first developed, prior to the execution of HTTF experiments. An agreement in prediction of temperature trends but challenges reproducing HTTF temperatures within measurement uncertainty is consistent with previous analyses of HTTF in the literature. Future RELAP5-3D validation activities centered around HTTF may be able to provide greater insight into the code’s capabilities for HTGR modeling with a more finely nodalized model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

International database of reference gamma spectra for nuclear safeguards applications

Nuclear safeguards missions use gamma spectroscopy as a non-destructive measurement technique for examining nuclear materials. Despite advances in the development of detection equipment as well as software codes, one of the concerns is the lack of well-documented spectra needed to test and validate isotopic analysis codes for their applicability. To address this need, this work introduces IDB, an international database of the reference gamma spectra of uranium, plutonium and mixed oxide nuclear material samples. IDB provides access to well-characterized sets of gamma spectra described by rich metadata, including information on the sample, measurement configuration and detector specifications. These spectra are accessible in different formats, also compatible with analysis code standards, thus promoting their sustainability and maintenance.

None, Dipti [International Atomic Energy Agency (I↗
DOE OSTI2024

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Basic Physical Processes Involving Dust in Fusion Plasmas

This report presents the main outcomes of our research focused on understanding the behavior and effects of dust particles in fusion plasmas, particularly in the edge regions of tokamaks and stellarators. We conducted computer modeling studies of burst injections of carbon and tungsten dust particles in DIII-D like divertor plasmas. The studies investigated effects of transient influx of the low- and high-Z material plasma contaminants in form of dust on the edge plasma dynamics in a modern mid-size tokamak. We also performed theoretical and computational studies of the forces acting on non-spherical dust grains in magnetized and nonmagnetized plasmas. In addition, we cooperated with General Atomics slag management group on development of DIII-D dust injection experiments to measure trajectories of dust grains in the divertor plasma and compare them with DUSTT code predictions for validation of the dust modeling capabilities. We collaborated with JET experimentalists on evaluation of radiative losses induced by injection of mixed neon-deuterium ice pellets for disruption and run-away electron generation mitigation during thermal and current quench phases. We also cooperated with experimentalists at LHD fusion device (Japan) on modeling support of experiments on injection of boron granules in fusion plasma discharges to assess their dynamics and effectiveness for in situ dynamic wall conditioning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Glovebox Fires Phase 1: Simulations of Open Burner Experiments

Understanding the duration of a fire required to breach a glovebox via glove ports and/or windows within a facility is important in developing mitigation strategies for safety. To facilitate this understanding, both simulation and experiments are utilized. The approach is to first validate the Computation Fluid Dynamics (CFD) code, Fire Dynamics Simulator (FDS), with experimental data collected by New Mexico Tech (NMT) and then perform simulations of full-scale rooms containing gloveboxes to assess numerous scenarios which would otherwise be cost-prohibitive experimentally. This report provides comparison to first-phase experiments involving a fire without a glovebox. A subsequent report will provide comparison to experiments involving a glovebox. The results indicate good agreement with FDS tending to over predict the pre-mixed and diffusion-mode tests by 6% and 10%, respectively.

42 ENGINEERING↗
DOE OSTI2026

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Polarized Deep-Inelastic Scattering with Spin Correlations in Herwig 7

This repository is the research software and reproducibility companion for the HerwigPol polarized deep-inelastic scattering implementation developed for Herwig 7. It brings together the modified Herwig and ThePEG source snapshots, the curated POLDIS fixed-order reference code, the custom Rivet analyses, the DIS validation workflow, and the paper source in a single formal repository layout. The repository is intended to preserve the source-level ingredients needed to rebuild and re-run the validated DIS studies. It therefore tracks code, input cards, workflow drivers, and technical notes, while intentionally excluding generated artifacts such as build products, campaign outputs, merged YODA files, plots, and rendered paper outputs.

Papaefstathioou, Andreas [Kennesaw State Universit↗