Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “I 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗
DOE OSTI2024

Validation of the SCALE/Polaris–PARCS Code Procedure With the ENDF/B-VII.1 AMPX 56-Group Library: Boiling Water Reactor

The SCALE/Polaris–PARCS code procedure has been used in the confirmatory analysis for boiling water reactors by the US Nuclear Regulatory Commission. In this study, the SCALE/Polaris v6.3.0–PARCS v3.4.2 code procedure with the Evaluated Nuclear Data File (ENDF)/B-VII.1 AMPX 56-group library was validated by comparing the simulated results with the measured data for operating boiling water reactors, including Peach Bottom Unit 2 cycles 1–3, Hatch Unit 1 cycles 1–3, and Quad Cities Unit 1 cycles 1–3. The uncertainties and biases of the SCALE/Polaris–PARCS code package for boiling water reactor physics analysis were evaluated in the validation for key nuclear parameters such as reactivity and traversing in-core probe data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

EMPIRE-PIC: A Performance Portable Unstructured Particle-in-Cell Code

In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

code for "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2"

This is all code that used in the journal article "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2", including: Running and post processing E3SM: Run E3SM: run_E3SM.2023.scidac.MMF.amip.py Regrid: regrid.E3SM.py post processing, diagnostics, archive: run_post.2023.scidac.MMF.py TEM calculation, including filter TEM, E3SM: calculate.TEM.v2.py calculate.TEM.monthly_mean.v1.py TEM, ERA5: ERA5_calculate_TEM_90.py merge TEM data: merge_TEM_E3SM_ERA5.py TEM, filter, E3SM: filter_TEM_E3SM.py TEM, filter, ERA5: filter_TEM_ERA5.py WK wave analysis code from (https://github.com/brianpm/wavenumber_frequency): wavenumber_frequency_functions.py calculate the wave spectra so easy to plot: wk_calculation.py Make figures some functions for plot: plot.py code for drawing paper plots: figures_for_paper.ipynb

Hannah, WalterM↗
DOE OSTI2022

Inputs and Outputs for paper "The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations"

Contains the DESC and VMEC input and output files used in the paper, as well as the plotting scripts used to create the figures in the paper (Current with the arXiv Mar 31 2022 version, https://arxiv.org/abs/2203.17173v1): 3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two 3D equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the MHD force error in real space directly. Accuracy as measured by final plasma energy and satisfaction of MHD force balance, as well as other metrics, will be presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. DESC's global Fourier-Zernike basis also yields the solution everywhere in the plasma volume, not just on discrete flux surfaces. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

V2XConnect: Harmonizing the Landscape of Bidirectional Charging Codes, Standards, and Communication Protocols

The utility grid is constantly evolving with new generation sources and diverse loads from backup generators, backup batteries, rooftop photovoltaics (PV), and flexible charging loads. Each of these technologies presents consumers with an opportunity to work with their local utility to promote grid stability, reliability, and affordability. Through a wide range of programs, utilities are incentivizing the use of these resources to benefit the grid. However, especially for devices with inverters that supply power back to the grid, it is important that these assets comply with local grid codes to ensure human safety, power quality, and grid reliability. One of the most promising of these new technologies is the bidirectional electric vehicle (EV) coupled with bidirectional electric vehicle supply equipment (EVSE). This report introduces the technologies that enable bidirectional charging and how these systems operate. It then outlines the necessary interconnection codes and industry standards that govern the operations of grid-tied inverters, including those in bidirectional charging. A key challenge is then outlined from the complex landscape of communication protocols and proprietary systems designed to coordinate these grid assets. Finally, a single gateway solution is proposed to harmonize communications between all assets across an entire site for a cohesive response to grid codes and utility programs that is flexible enough to support a range of market options.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Path-synchronous performance monitoring of interconnection networks based on source code attribution

Examples disclosed herein relate to path-synchronous performance monitoring of an interconnection network based on source code attribution. A processing node in the interconnection network has a profiler module to select a network transaction to be monitored, determine a source code attribution associated with the network transaction to be monitored, and issue a network command to execute the network transaction to be monitored. A logger module creates, in a buffer, a node temporal log associated with the network transaction and the network command. A drainer module periodically captures the node temporal log. The processing node has a network interface controller to receive the network command and mark a packet generated for the network command to be temporally tracked and attributed back to the source code attribution at each hop of the interconnection network traversed by the marked packet.

Chabbi, Milind M.↗
DOE OSTI2022

System and method for identifying and comparing code by semantic abstractions

Certain embodiments of the present invention are configured to facilitate analyzing computer code more efficiently. For example, by conducting a first level abstraction (e.g., symbolic interpretation and algebraic simplification) and a second level abstraction (e.g., generalization) of the computer code, the analysis may more accurately account for variations in the code that may occur as a result of register renaming, instruction reordering, choice of instructions, etc. while minimizing the cost of computations required to perform the analysis.

Lakhotia, Arun↗
DOE OSTI2024

The Impact of the Halden Reactor Project on the Development of the Bison Fuel Performance Code

The Bison fuel performance code has been under development at Idaho National Laboratory by various programs within the U.S. Department of Energy since 2009. As part of any fuel performance code development, validation is necessary. To this end, several Halden experiments have been used. Early on, the focus was on validating the thermal mechanical behaviour, including pellet cladding mechanical interaction and fission gas release models implemented in Bison. Support for design-basis accident capabilities began in approximately 2014. Also in 2014, the development of models for accident tolerant fuel (ATF) concepts also began. This includes dopedUO2 fuel for which Halden has some experiments with temperature and fission gas release measurements. The Bison team has been fortunate to have two senior developers be secondees to Halden for a year: the late Dr. Giovanni Pastore from May 2015 to April 2016 and the now retired Dr. Richard Williamson in 2018. The paper will provide examples of Bison validation using Halden data for standard UO2 fuel under normal operation and loss-of-coolant accident condition, and ATF concepts. A tribute to the legacy that both the late Giovanni and retired Richard left on the Bison code is provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Status of the MARS code

This report describes major features of the most recent version of the MARS code as well as ongoing developments. The list of features includes various options for geometry models, a beam line builder based on MADX code, import of geometry models in GDML format, use of structured and unstructured meshes for scoring purposes, an update to the recent TENDL library for a number of projectiles at low energies (up to 250 MeV), and a recently implemented method to calculate spatial distribution of residual dose in a single computer run without an intermediate source. Examples of the code application to various projects are presented as well.

Rakhno, Igor [Fermilab] (ORCID:0000000265828058)↗
DOE OSTI2025

Using Calibrated Sodium Data for Preliminary Validation of the SRT Code for Advanced Reactors

Various types of non-light water reactors are currently engaged in the U.S. licensing process. Because of inherent differences compared with well-established large light water reactors, appropriate assessment tools are needed. Specifically, source term analysis, which determines environmental dose impacts from potential accident scenarios, is a crucial part of design and licensing. The U.S. Nuclear Regulatory Commission has emphasized the importance of mechanistic source term analysis for advanced reactor deployments. To align with these needs, Argonne National Laboratory has developed the Simplified Radionuclide Transport (SRT) source term analysis code for metal fuel Sodium-cooled Fast Reactors (SFRs) and microreactors. SRT conducts time-dependent radionuclide transport and retention in SFRs for core and ex-core radionuclide source accident sequences. The main objective of SRT is to provide rapid sensitivity and uncertainty analyses, incorporating parametric uncertainties and summarizing probabilistic results. As part of the code validation process, a study focused on the bubble scrubbing module was performed using an experiment recently carried out by the University of Wisconsin-Madison. Based on the analysis, the modeling approach in SRT provides accurate results for small and large aerosols, while slight underprediction of radionuclide aerosol removal are observed for medium sized aerosols. However, the deviation is minor, considering the highly uncertain phenomenon and range of results, and is in the conservative direction. In addition, uncertainty information derived from the experiments is further implemented, reflecting the actual span of parameters, which leads to enhanced agreement with code predictions. The results demonstrate that SRT provides reasonable predictions for the bubble scrubbing process in sodium pool.

Kam, Dong Hoon↗
DOE OSTI2023

Code for multi-shape Gaussian process (GP) fitting with uncertainty quantification (UQ)

This is code associated with the publication “Nonparametric Multi-shape Modeling with Uncertainty Quantification,” authored by Hengrui Luo (Lawrence Berkeley National Laboratory) and Justin Strait (Los Alamos National Laboratory). The code is used to fit multiple-output Gaussian process (GP) models of planar closed curves to collections of ordered point sets, allowing for flexible nonlinear prediction of the underlying curve under dense or sparse point set samplings and with or without noise, as well as tractable uncertainty quantification. To do this, we employ use of a periodic kernel to account for the nonlinear input space of closed curves, and combine with coregionalization models to account for dependence both (i) between curve coordinates, and (ii) between pairs of curves. Functions in the code are capable of fitting these models, as well as performing additional tasks with the fitted curves such as (a) shape registration and alignment, (b) shape averaging, and (c) fitting for curve sub-populations / clusters.

Strait, Justin↗
DOE OSTI2023

Python codes for the paper - "Trade-offs in the latent representation of microstructure evolution"

SAND2024-00946O Python codes used in "Trade-offs in the latent representation of microstructure evolution," a manuscript accepted for publication in Acta Materialia, are part of a repository. The code was developed to perform analysis of microstructure evolution. The repository consists of two main directories: models, which train and test models such as autoencoders and diffusion maps, and analysis, which analyzes microstructures. Source code is used to perform dimensionality reduction of microstructure data for analysis of its evolution in time. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dingreville, Remi↗
DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗