Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer 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 37 records · Page 2

MSTor 2023 : A new version of the computer code for multistructural torsional anharmonicity, now with automatic torsional identification using redundant internal coordinates

The MSTor program is a computer program for calculating partition functions and thermodynamic functions of complex gas-phase molecules with multiple torsions; the multi-structural approximation with torsional anharmonicity (MS-T) can be used based on either a coupled torsional potential or an uncoupled torsional potential. The program can also carry out calculations in the multiple-structure local harmonic (MS-LH) approximation or multi-structural local quasiharmonic (MS-LQ) approximation or by the dual-level MS-T method. Furthermore, the program package includes eight utility codes that can be used as stand-alone programs to help the user to generate the input files for the code or to generate comparison results. The 2023 version of MSTor includes a new capability, namely automatic identification of torsional modes with redundant internal coordinates.

97 MATHEMATICS AND COMPUTING↗

Interface design of VSOP'94 computer code for safety analysis

Today, most software applications, also in the nuclear field, come with a graphical user interface. VSOP'94 (Very Superior Old Program), was designed to simplify the process of performing reactor simulation. VSOP is a integrated code system to simulate the life history of a nuclear reactor that is devoted in education and research. One advantage of VSOP program is its ability to calculate the neutron spectrum estimation, fuel cycle, 2-D diffusion, resonance integral, estimation of reactors fuel costs, and integrated thermal hydraulics. VSOP also can be used to comparative studies and simulation of reactor safety. However, existing VSOP is a conventional program, which was developed using Fortran 65 and have several problems in using it, for example, it is only operated on Dec Alpha mainframe platforms and provide text-based output, difficult to use, especially in data preparation and interpretation of results. We develop a GUI-VSOP, which is an interface program to facilitate the preparation of data, run the VSOP code and read the results in a more user friendly way and useable on the Personal 'Computer (PC). Modifications include the development of interfaces on preprocessing, processing and postprocessing. GUI-based interface for preprocessing aims to provide a convenience way in preparing data. Processing interface is intended to provide convenience in configuring input files and libraries and do compiling VSOP code. Postprocessing interface designed to visualized the VSOP output in table and graphic forms. GUI-VSOP expected to be useful to simplify and speed up the process and analysis of safety aspects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The ocean model for E3SM global applications: Omega version 0.1.0 – a new high-performance computing code for exascale architectures

This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.

54 ENVIRONMENTAL SCIENCES↗

Debugging Computer Code [Slides]

Debugging? Klingons do not debug. Our software does not coddle the weak. For the rest of us, finding and removing defects in code is a fundamental programming skill. Besides preventing code from performing its intended function, bugs destroy users’ confidence in the reliability of code even when no errors are obvious. A through discussion of debugging techniques could easily fill a semester course in computer science. We have an hour here today. I will take that time to discuss broad principles rather than fine points of debugging.

97 MATHEMATICS AND COMPUTING↗

PROBLEM DEPENDENT DOPPLER BROADENING OF CONTINUOUS ENERGY CROSS SECTIONS IN THE KENO MONTE CARLO COMPUTER CODE

For many Monte Carlo codes cross sections are generally only created at a set of predetermined temperatures. This causes an increase in error as one moves further and further away from these temperatures in the Monte Carlo model. This paper discusses recent progress in the Scale Monte Carlo module KENO to create problem dependent, Doppler broadened, cross sections. Currently only broadening the 1D cross sections and probability tables is addressed. The approach uses a finite difference method to calculate the temperature dependent cross-sections for the 1D data, and a simple linear-logarithmic interpolation in the square root of temperature for the probability tables. Work is also ongoing to address broadening theS (alpha , beta) tables. With the current approach the temperature dependent cross sections are Doppler broadened before transport starts, and, for all but a few isotopes, the impact on cross section loading is negligible. Results can be compared with those obtained by using multigroup libraries, as KENO currently does interpolation on the multigroup cross sections to determine temperature dependent cross-sections. Current results compare favorably with these expected results.

Monte Carlo↗

Modernization of ORIGEN Library Creation in the SCALE Computer Code

The modernization of SCALE features has been a target for recent SCALE releases. ORIGEN, the activation and depletion code in SCALE, is one such legacy code that has been targeted for modernization. A part of these modernization efforts has been to develop capabilities within ORIGEN to directly read ENDF-formatted data files. In addition, the input format for the COUPLE library creation utility was rewritten to a format more accessible to newer SCALE users. Finally, the ability of COUPLE to create ORIGEN libraries without sublibrary data has been added.KEYWORDS: ORIGEN, depletion, SCALE, API

Hart, Shane↗

React-compact: a computer code for modeling graphite corrosion and fuel hydrolysis - 18526

REACT-COMPACT is used to estimate localized graphite corrosion, fuel compact corrosion, and fuel hydrolysis by water vapor for a unit cell geometry within a High Temperature Gas Cooled Reactor (HTGR) fuel block. REACT-COMPACT may also be used to model corrosion in regions with graphite only, including core support posts. REACT-COMPACT solves the species-conservation equation for water vapor in a porous carbon medium for cylindrical geometry. Boundary-layer mass transfer from a helium coolant is modeled at either the inner or outer radius. Within the porous carbon medium, both diffusive and convective transport mechanisms are modeled. The reaction of water vapor with carbonaceous materials is modeled according to a Langmuir-Hinshelwood mechanism. The fuel hydrolysis reaction rate is based on a semi-empirical model derived from experimental data. Model parameters, including densities, void fractions, burn off factors for the intrinsic reaction rates, and diffusivity of water vapor in porous carbon, are evaluated as functions of local burn off. The effects of void-fraction gradients on transport processes are also modeled. REACT-COMPACT also solves the mass conservation equations for reacting fuel and carbonaceous materials, which are coupled in a nonlinear manner to the water-vapor conservation equation because of burn off effects. For this reason, a numerical finite-difference solution is used. This paper describes the theoretical basis for REACT-COMPACT and results for sample problems under representative HTGR conditions. (author)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Q-value (MeV/fission) Determination for the Advanced Test Reactor

This report documents a step-by-step procedure that can be used to calculate a Q-value for the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). The Q-value calculation here is an estimate of the total energy released per fission event in the ATR driver fuel in units of MeV/fission and is recognized as a reactor-specific Q-value for the ATR. Such a Q-value is needed INL ATR reactor physicists to normalize calculated fluxes and reaction rates to measured ATR total core and lobe powers. The calculated ATR Q-value and the procedure used to calculate this value are presented in a detailed manner so as to establish a credible basis, such that this single Q-value can be used consistently among all ATR reactor physicists. A Q-value is essentially the summation of all the kinetic energy from all radiation components released from a single fission event that ends up as heat in the reactor coolant. Estimation of a Q-value for a given reactor can theoretically be very complex because it is dependent on several variables, which include (1) the materials in the reactor core, (2) the neutron energy spectrum, and (3) the mix of actinide isotopes that are being simultaneously fissioned. Since these three variables can potentially change over time or with reactor burnup, this potentially makes the Q-value a time-dependent quantity. However, for ATR, we will invoke some simplifying assumptions, based on the unique features of the ATR core and supported by neutronic calculations, which will reduce the complexity of the problem and lead to a single time-independent Q-value. Of the three variables, the materials in the core have the greatest impact on the ATR Q-value. Since the materials in the core are most important, an array of components (materials) have been carefully selected to compose a typical ATR core configuration. It should be noted that the step-by-step procedure here for calculating a reactor-specific Q-value can be applied to any nuclear reactor. In a more general procedure, if one needed to include time-dependent spectral and actinide composition changes as a function of burnup for a specific reactor, this procedure can still be used, but at discrete time steps throughout the burnup. The result will be a set of discrete Q-values as a function of burnup. However, here we will apply the procedure to the ATR, as an example, and obtain a single Q-value that is valid for the ATR reactor as a function of burnup. In the evaluated nuclear data file (ENDF), Q-values for specific actinides are known for a reasonably high degree of accuracy, especially the major uranium and plutonium isotopes. The isotopic Q-values are composed of the kinetic energy of the fission fragments and the energies of the various radiation components emitted from the fission event. One complication is the fact that some of these energy components are dependent on the incident neutron energy (EINC) that induces the fission event (neutron energy spectrum dependence). For ATR, we will estimate an average neutron energy based on the neutron energy flux weighted by the fission cross sections and use this average neutron energy value in the ENDF formulas. This neutron energy dependence is one of two dependencies that make the final ATR Q-value a reactor-specific value. The second dependency, and the most important variable for ATR, involves the materials in the core and the de-excitation radiations produced from neutron radiative capture radiation; radiation energy that becomes heat in the reactor coolant. This important dependency is discussed in detail below. The third dependency, or the mix of fissioning actinides as a function of ATR burnup, is ignored because the bulk of the fissions over a cycle are always due to U-235 in the high-enriched ATR core. In a fission event, the fission energy is partitioned between the fission fragments and the radiation components that include neutrons, gamma rays, and beta particles. The main variable that makes the calculated Q-value reactor-specific is the amount of radiative capture energy produced in the core due to prompt and delayed neutron capture in the reactor core materials and the subsequent release of nuclear binding energy in the form of de-excitation gamma radiation energy in the reactor core. This energy component, neutron radiative capture energy, is not in the ENDF, because it depends on the reactor and its core materials. The radiative capture gamma energy release can be significant, or approximately 3–12 MeV/fission depending on the materials in the reactor core (fuel materials, cladding, structural components, coolant, targets, reflector, etc.). For the ATR, the energy released from radiative capture gamma radiation is calculated to be near the middle of the 3–12 MeV range. This makes the radiative capture radiation component an important and significant contribution to the total ATR Q-value. The ATR has unique reactor features that will allow us to simplify the complex procedure for determining an ATR reactor-specific Q-value. For example, the ATR is a thermal reactor and contains high-enriched fuel (93.15 wt% U-235). Consequently, the bulk of the fissions in the ATR core from beginning-of-cycle (BOC) to end-of-cycle are from U-235. Even though U-238 makes up approximately 7% of the uranium atoms in the fuel, the number of fissions relative to U-235 is insignificant, since U-238 is fertile and ATR is a thermal reactor. The same is true for U-234 and U-236 (both fertile). The ATR’s high-enriched uranium fuel also inhibits the buildup of large quantities of plutonium isotopes. So, even though there is some buildup of U-236, Np-237, and Pu-239 in the ATR driver fuel, the impact or contribution of these isotopes fissions on the ATR Q-value is minor compared to U-235. Nonetheless, we will include fission contributions from all the major uranium isotopes (U-234, U-235, U-236, and U-238) in the determination of an average isotopic Q-value that is used in the ATR Q-value determination. Although the ATR Q-value calculated here is based on BOC conditions (fresh fuel with no higher order actinide buildup), it is believed that the ATR Q-value is relatively constant with burnup or during any given ATR power cycle (14–60 days). The reason is that the vast majority of the core driver fuel fissions are always from U-235 during the cycle. Also, since we effectively assume instantaneous energy release for the fission product gammas and beta particles and delayed neutrons from a fission event, this corresponds to equilibrium conditions for all fission product concentrations in the ATR driver fuel too. The goal of this effort was to determine a specific numerical and core-average Q-value for the recoverable energy released by a single fission event (MeV/fission) for the ATR. This calculated value has the potential to help standardize and add credibility to not only the Monte Carlo tally normalization process used by many ATR core physics analysts, but also to the normalization process for all reactor physics calculation methods involving ATR analyses. Currently, the ATR core physics analysts use slightly different Q-values. Some values are known to be conservative for specific applications. However, sometimes a more accurate Q-value is needed for the calculation in which ATR test specimens require non-conservative and as-accurate-as-possible Q-values to estimate burnups, heat rates, neutron fluxes, nuclear reaction rates, etc., where corresponding measured post-irradiation examination data provides the ultimate validation for the calculated results. A detailed Q-value procedure was developed and is outlined in this document along with numerical values calculated specifically for the determination of the ATR Q-value. This report attempts to document in some detail the calculation methodology, nuclear data, computer code, computer model, code calculated results, hand calculations, assumptions, and logic that went into the determination of the ATR Q-value.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nyx: A MASSIVELY PARALLEL AMR CODE FOR COMPUTATIONAL COSMOLOGY

We present a new N-body and gas dynamics code, called Nyx, for large-scale cosmological simulations. Nyx follows the temporal evolution of a system of discrete dark matter particles gravitationally coupled to an inviscid ideal fluid in an expanding universe. The gas is advanced in an Eulerian framework with block-structured adaptive mesh refinement; a particle-mesh scheme using the same grid hierarchy is used to solve for self-gravity and advance the particles. Computational results demonstrating the validation of Nyx on standard cosmological test problems, and the scaling behavior of Nyx to 50,000 cores, are presented.

79 ASTRONOMY AND ASTROPHYSICS↗

Scientific Computational Imaging Code (SCICO)

Scientific Computational Imaging Code (SCICO) is a Python package for solving the inverse problems that arise in scientific imaging applications. Its primary focus is providing methods for solving ill-posed inverse problems by using an appropriate prior model of the reconstruction space. SCICO includes a growing suite of operators, cost functionals, regularizers, and optimization routines that may be combined to solve a wide range of problems, and is designed so that it is easy to add new building blocks. SCICO is built on top of JAX rather than NumPy, enabling GPU/TPU acceleration, just-in-time compilation, and automatic gradient functionality, which is used to automatically compute the adjoints of linear operators. An example of how to solve a multi-channel tomography problem with SCICO is shown in Figure 1. The SCICO source code is available from GitHub, and pre-built packages are available from PyPI. It has extensive online documentation, including API documentation and usage examples, which can be run online at Google Colab and binder.

97 MATHEMATICS AND COMPUTING↗

On the Green's function of the partially diffusion-controlled reversible ABCD reaction for radiation chemistry codes

Several computer codes simulating chemical reactions in particles systems are based on the Green's functions of the diffusion equation (GFDE). Indeed, many types of chemical systems have been simulated using the exact GFDE, which has also become the gold standard for validating other theoretical models. In this work, a simulation algorithm is presented to sample the interparticle distance for partially diffusion-controlled reversible ABCD reaction. This algorithm is considered exact for 2-particles systems, is faster than conventional look-up tables and uses only a few kilobytes of memory. The simulation results obtained with this method are compared with those obtained with the independent reaction times (IRT) method. This work is part of our effort in developing models to understand the role of chemical reactions in the radiation effects on cells and tissues and may eventually be included in event-based models of space radiation risks. However, as many reactions are of this type in biological systems, this algorithm might play a pivotal role in future simulation programs not only in radiation chemistry, but also in the simulation of biochemical networks in time and space as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗