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 289 records · Page 16

Modernization and optimization of a legacy open-source CFD code for high-performance computing architectures

Legacy codes remain a crucial element of today's simulation-based engineering ecosystem due to the extensive validation process and investment in such software. The rapid evolution of high-performance computing architectures necessitates the modernization of these codes. One approach to modernization is a complete overhaul of the code. However, this could require extensive investments, such as rewriting in modern languages, new data constructs, etc., which will necessitate systematic verification and validation to re-establish the credibility of the computational models. The current study advocates using a more incremental approach and is a culmination of several modernization efforts of the legacy code MFIX, which is an open-source computational fluid dynamics code that has evolved over several decades, widely used in multiphase flows and still being developed by the National Energy Technology Laboratory. Two different modernization approaches,‘bottom-up’ and ‘top-down’, are illustrated. Here, preliminary results show up to 8.5x improvement at the selected kernel level with the first approach, and up to 50% improvement in total simulated time with the latter were achieved for the demonstration cases and target HPC systems employed.

97 MATHEMATICS AND COMPUTING↗

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of Transport Processes in Liquid-Metal Fusion Blankets: Past, Present, and Future

The successful development of robust breeding blanket systems will strongly rely on computational tools for predicting the complex behavior of the electrically conducting liquid-metal (LM) breeder flowing in the complex-shaped blanket ducts in the presence of a strong plasma-confining magnetic field, volumetric heating, and tritium generation. Associated transport processes involve magnetohydrodynamic (MHD) flows, heat transfer, corrosion, and tritium transport. This paper is an overview of past and present efforts in the development, application, and verification and validation (V&V) of such computational tools. As a result of the ongoing campaign on V&V of computer codes for LM blankets, the international fusion community has identified several candidates that promise to become real blanket design and analysis tools in the near future. Among them are HIMAG, MHD-UCAS, COMSOL Multiphysics, ANSYS FLUENT, ANSYS CFX, and OpenFOAM. The progress, over the last decade, in the application of such codes in blanket studies is tremendous. This is illustrated with two examples for a dual-coolant lead-lithium (DCLL) blanket: (1) integrated computer modeling for the recently designed DCLL blanket in the United States and (2) application of the code MHD-UCAS to the analysis of PbLi flows and heat transfer in a generic DCLL blanket prototype at high Hartmann (Ha ~ 10 4 ) and Grashof numbers (Gr ~ 10 12 ). Here, this paper also presents an approach to the development of a new integrated computational tool called the virtual dual-coolant lead-lithium (VDCLL) blanket, which elaborates the existing U.S. MHD code HIMAG.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Re-evaluation of Spent Nuclear Fuel Assay Data for the Three Mile Island Unit 1 Reactor and Application to Code Validation

Destructive radiochemical assay measurements of spent nuclear fuel rod segments from an assembly irradiated in the Three Mile Island unit 1 (TMI-1) pressurized water reactor have been performed at Oak Ridge National Laboratory (ORNL). Assay data are reported for five samples from two fuel rods of the same assembly. The TMI-1 assembly was a 15 X 15 design with an initial enrichment of 4.013 wt% 235 U, and the measured samples achieved burnups between 45.5 and 54.5 gigawatt days per metric ton of initial uranium (GWd/t). Measurements were performed mainly using inductively coupled plasma mass spectrometry after elemental separation via high performance liquid chromatography. High precision measurements were achieved using isotope dilution techniques for many of the lanthanides, uranium, and plutonium isotopes. Measurements are reported for more than 50 different isotopes and 16 elements. One of the two TMI-1 fuel rods measured in this work had been measured previously by Argonne National Laboratory (ANL), and these data have been widely used to support code and nuclear data validation. Recently, ORNL provided an important opportunity to independently cross check results against previous measurements performed at ANL. The measured nuclide concentrations are used to validate burnup calculations using the SCALE nuclear systems modeling and simulation code suite. These results show that the new measurements provide reliable benchmark data for computer code validation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Compiler-Assisted Detection of Transient Memory Errors

The probability of bit flips in hardware memory systems is projected to increase significantly as memory systems continue to scale in size and complexity. Effective hardware-based error detection and correction requires that the complete data path, involving all parts of the memory system, be protected with sufficient redundancy. First, this may be costly to employ on commodity computing platforms and second, even on high-end systems, protection against multi-bit errors may be lacking. Therefore, augmenting hardware error detection schemes with software techniques is of consider- able interest. In this paper, we consider software-level mechanisms to comprehensively detect transient memory faults. We develop novel compile-time algorithms to instrument application programs with checksum computation codes so as to detect memory errors. Unlike prior approaches that employ checksums on computational and architectural state, our scheme verifies every data access and works by tracking variables as they are produced and consumed. Experimental evaluation demonstrates that the proposed comprehensive error detection solution is viable as a completely software-only scheme. We also demonstrate that with limited hardware support, overheads of error detection can be further reduced.

def-use analysis↗

Extension of Argonne’s Plant Dynamics Code (PDC) Capabilities to TerraPower Direct-Cycle Supercritical CO2 (“Pascal”) Nuclear Reactor Concept. Final CRADA report

Argonne National Laboratory (the Contractor) and TerraPower, LLC, (the Participant), entered into a Cooperative Research and Development Agreement (CRADA) to extend the Plant Dynamics Code (PDC) to model the Pascal reactor plant. PDC is the system-level computer code developed at Argonne, primarily for design and transient analysis of supercritical carbon dioxide (sCO2) energy conversion systems for various power plants. Pascal is a supercritical carbon dioxide (sCO2)-cooled direct-cycle nuclear power reactor under development by TerraPower. New features were added to the code to model components and systems specific to the Pascal plant and to make the code readily available as a tool for design, control, and safety analysis for this and other direct-cycle power plants. TerraPower provided design information of the Pascal reactor and the requirements for the development of the new PDC features. Argonne modified the PDC capabilities to introduce the new features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sockeye Heat Pipe Code Theory Development: Based on the 7-Equation, Two-Phase Flow Model of RELAP-7

Sockeye is a MOOSE-based computer code designed to simulate the thermal-hydraulic two-phase flows of heat pipes in a nuclear micro-reactor. For this application, a significantly modified version of the fully non-equilibrium, compressible, 7-equation, two-phase flow model of RELAP-7 was developed. This summary report documents the significant points of this model development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

MCNP ® Code Version 6.3.0 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.0, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

GALE 3.2 Verification Report

This report is documentation of the verification activities that have been performed on GALE-3.2 to ensure that the codes meet the requirements and specifications set forth by NRC and that the code meets its intended purpose. These activities are done in accordance with the GALE Software Quality Assurance Plan (PNNL-24249) and the GALE Configuration Management & Maintenance Plan (PNNL-24250). These two plans meet the requirements set for by both PNNL software quality assurance requirements and the NRC’s Software Quality Assurance Program and Guidelines (NUREG/BR-0167). The verification activities are documented to conform to Level 2 requirements of NUREG/BR-0167. Verifications were performed for both GALE-3.2 BWR and GALE-3.2 PWR Codes. The GALE 3.2 BWR and PWR codes are computerized mathematical models for calculating the releases of radioactive material in gaseous and liquid effluents. These computer codes are tools which can be used to determine conformance with the requirements of Appendix I to Title 10, Part 50 of the Code of Federal Regulations (10 CFR 50).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hand Calculation Methods for Nuclear Criticality Safety

This primer provides an overview of the most common hand calculation methods used for criticality safety calculations. The most widely used tools available to a nuclear criticality safety (NCS) practitioner are probably the common Monte Carlo or deterministic criticality safety codes, which can be used to model very complex systems. However, use of these codes can obscure the parameters to which a particular fissile system may be sensitive, whereas the hand calculation methods can be used to delve into the ways each parameter may affect the reactivity of a fissile material system. Furthermore, practitioners must avoid using computer codes as devices that take inputs and simply provide outputs (i.e., a “black box”). Many years ago, pioneers such as Joe Thomas, David Smith, and Hugh Paxton, among others in the field of nuclear criticality safety, took the time before the advent of high-speed desktop computers to create simple hand methods for criticality safety analyses. Some of the methods can be used for single fissile units; others are applicable to fissile units arranged into simple array configurations. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the various results. The NCS practitioner will need to spend time to master the methods that could be most useful; however, they can provide the practitioner with fast and accurate answers to criticality safety problems if they are used correctly and if critical data exist for the problem at hand. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information more quickly than using a criticality code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Parallel computation safety analysis irradiation targets fission product molybdenum in neutronic aspect using the successive over-relaxation algorithm

One of the research activities in support of commercial radioisotope production program is a safety research on target FPM (Fission Product Molybdenum) irradiation. FPM targets form a tube made of stainless steel which contains nuclear-grade high-enrichment uranium. The FPM irradiation tube is intended to obtain fission products. Fission materials such as Mo{sup 99} used widely the form of kits in the medical world. The neutronics problem is solved using first-order perturbation theory derived from the diffusion equation for four groups. In contrast, Mo isotopes have longer half-lives, about 3 days (66 hours), so the delivery of radioisotopes to consumer centers and storage is possible though still limited. The production of this isotope potentially gives significant economic value. The criticality and flux in multigroup diffusion model was calculated for various irradiation positions and uranium contents. This model involves complex computation, with large and sparse matrix system. Several parallel algorithms have been developed for the sparse and large matrix solution. In this paper, a successive over-relaxation (SOR) algorithm was implemented for the calculation of reactivity coefficients which can be done in parallel. Previous works performed reactivity calculations serially with Gauss-Seidel iteratives. The parallel method can be used to solve multigroup diffusion equation system and calculate the criticality and reactivity coefficients. In this research a computer code was developed to exploit parallel processing to perform reactivity calculations which were to be used in safety analysis. The parallel processing in the multicore computer system allows the calculation to be performed more quickly. This code was applied for the safety limits calculation of irradiated FPM targets containing highly enriched uranium. The results of calculations neutron show that for uranium contents of 1.7676 g and 6.1866 g (× 10{sup 6} cm{sup −1}) in a tube, their delta reactivities are the still within safety limits; however, for 7.9542 g and 8.838 g (× 10{sup 6} cm{sup −1}) the limits were exceeded.

07 ISOTOPE AND RADIATION SOURCES↗

GALE 3.1 Verification Report

This report is documentation of the verification activities that have been performed on GALE-3.1 to ensure that the codes meet the requirements and specifications set forth by NRC and that the code meets its intended purpose. Additional benchmark tests for GALE-3.1 are documented in Appendix A of NUREG-0016, Revision 2 for GALE-BWR 3.1 and Appendix A of NUREG-0017, Revision 2 for GALE-PWR 3.1. These activities are done in accordance with the GALE Software Quality Assurance Plan (PNNL-24249) and the GALE Configuration Management & Maintenance Plan (PNNL-24250). These two plans meet the requirements set for by both PNNL software quality assurance requirements and the NRC’s Software Quality Assurance Program and Guidelines (NUREG/BR-0167). The verification activities are documented to conform to Level 2 requirements of NUREG/BR-0167. Verifications were performed for both GALE-3.1 BWR and GALE-3.1 PWR Codes. The GALE 3.1 BWR and PWR codes are computerized mathematical models for calculating the releases of radioactive material in gaseous and liquid effluents. These computer codes are tools which can be used to determine conformance with the requirements of Appendix I to Title 10, Part 50 of the Code of Federal Regulations (10 CFR 50).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DEVELOPMENT OF A VERIFICATION BENCHMARK SUITE FOR THE ATR REACTOR PHYSICS UPGRADE

The result of a numerical algorithm is only as good as the effort made toward its verification. Comprehensive computer codes, especially for application to nuclear reactors, require in-depth assessment of the numerical methods used and their computational implementation. This presentation describes a suite of semi-analytical benchmarks for the ATR Reactor Physics Upgrade as part of the ATR Life Extension Project (LEP). Several neutron transport codes including deterministic codes HELIOS, NEWT and probabilistic codes MCNP6, SERPANT and MC21 are upgrade candidates. To properly evaluate code performance capabilities, code verification is an essential element. A credible verification program requires education to guide proper application, training in proper input preparation and output expectation and verification to instill confidence in the coded numerical procedures. The subset of the benchmark suite to be described here is primarily of the 1D transport variety (however, a 2D case does indeed exist) with the intent to provide unquestionably high order accurate semi-analytical numerical solutions to the transport equation. We base the analytical solutions on analytical Fourier series and discrete ordinates representations. Here we present two 1D benchmarks based on Fourier series.

99 GENERAL AND MISCELLANEOUS↗

SNL-NJOY-2016

This is a wrapper around the LANL NJOY-2016 code that interfaces with extended code capabilities and minor modification of the LANL NJOY-2016 code that supports modeling radiation damage to materials. The NJOY Nuclear Data Processing System is a modular computer code designed to read evaluated data in ENDF format, transform the data in various ways, and output the results as libraries designed to be used in various applications. The wrapper provided here permits Sandia-specific control parameters to be used in the input data file. The modifications incorporated here enhance the ability of NJOY to address material damage response functions for many materials, e.g. to include the NRT and arc-dpa forms of the damage energy in addition to the default NJOY-2016 implementation of the sharp-threshold Kinchin-Pease threshold energy treatment. All of the modifications provided here are being made available to the GitHub-based NJOY-2016 code. As useful enhancements found here are incorporated into the baseline NJOY-2016 code, they will be eliminated from this version so as to maintain our compatibility with the baseline NJOY-2016 code. SAND2020-13060 M 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-NA00035

Griffin, Patrick↗

MCNP6 simulation of light and medium nuclei fragmentation at intermediate energies

Fragmentation reactions induced on light and medium nuclei by protons and light nuclei of energies around 1 GeV/nucleon and below are studied with the Los Alamos transport code MCNP6 and with its CEM03.03 and LAQGSM03.03 event generators. CEM and LAQGSM assume that intermediate-energy fragmentation reactions on light nuclei occur generally in two stages. The first stage is the intranuclear cascade (INC), followed by the second, Fermi breakup disintegration of light excited residual nuclei produced after the INC. CEM and LAQGSM account also for coalescence of light fragments (complex particles) up to sup>4 He from energetic nucleons emitted during INC. We investigate the validity and performance of MCNP6, CEM, and LAQGSM in simulating fragmentation reactions at intermediate energies and discuss possible ways of further improving these codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MCNP® Code Version 6.3.2 Theory & User Manual (Rev. 1)

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.2, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

42 ENGINEERING↗

Developing and Implementing the Data Mining Algorithms in RAVEN

The RAVEN code is becoming a comprehensive tool to perform probabilistic risk assessment, uncertainty quantification, and verification and validation. The RAVEN code is being developed to support many programs and to provide a set of methodologies and algorithms for advanced analysis. Scientific computer codes can generate enormous amounts of data. To post-process and analyze such data might, in some cases, take longer than the initial software runtime. Data mining algorithms/methods help in recognizing and understanding patterns in the data, and thus discover knowledge in databases. The methodologies used in the dynamic probabilistic risk assessment or in uncertainty and error quantification analysis couple system/physics codes with simulation controller codes, such as RAVEN. RAVEN introduces both deterministic and stochastic elements into the simulation while the system/physics code model the dynamics deterministically. A typical analysis is performed by sampling values of a set of parameter values. A major challenge in using dynamic probabilistic risk assessment or uncertainty and error quantification analysis for a complex system is to analyze the large number of scenarios generated. Data mining techniques are typically used to better organize and understand data, i.e. recognizing patterns in the data. This report focuses on development and implementation of Application Programming Interfaces (APIs) for different data mining algorithms, and the application of these algorithms to different databases.

97 MATHEMATICS AND COMPUTING↗