Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 901 records · Page 50

Updated Draft ASME Boiler and Pressure Vessel Code Case for Use of Alloy 617 for Construction of Nuclear Components for Section III Division 5

Early in the United States’ very high temperature reactor program, several candidate nickel alloys were considered for use in construction of the intermediate heat exchanger. Based primarily on technical maturity, a downselection was made to focus on Alloy 617. After this downselection, the primary goal of the research and development program was to develop sufficient information on the high temperature properties of the material to qualify it for construction of high temperature nuclear components in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel (BPV) Code. Alloy 617 is a solid solution strengthened material with nominal Ni Cr Co Mo composition, originally developed for aerospace applications such as burner can liners for turbine engines by Huntington Alloys. (Huntington Alloys is now Special Metals Division of Precision Castparts, Inc.) The ASME BPV Code allows use of Alloy 617 for construction of non nuclear pressure vessels, and Alloy 617 is used in fossil fired power plants. Section III, Division 1, Subsection NB of the ASME BPV Code was developed for construction of nuclear components in light water reactors and allows use of ferritic materials up to 700°F and austenitic alloys up to 800°F. Subsection NH of Section III Division 1 was written to allow higher temperature construction with a primary focus on sodium cooled reactors. Recently, a new Division 5 of Section III was published specifically for high temperature reactors (regardless of the primary working fluid) and incorporates both Subsections NB and NH. There are only five alloys currently allowed for use in high temperature nuclear components: 2.25Cr 1 Mo and V modified 9Cr 1Mo steels, Types 304 and 316 stainless steels, and the high nickel austenitic Alloy 800H. This very sparse set of allowed materials is in contrast with the collection of more than 150 materials that are allowed for use in non nuclear pressure vessel construction. A draft Code Case to add Alloy 617 to the list of qualified alloys for use in high temperature nuclear design was submitted to ASME in the early 1990s, but it was withdrawn prior to formal Code action.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Develop SRNL's ALGE3D code for use as national response asset

ALGE3D is a three-dimensional hydrodynamic code which solves the momentum, mass and energy conservation equations. Its original purpose was to predict movement and dissipation of thermal plumes. To fill a gap in existing emergency response technology this code is being modified to predict the movement and dissipation of tracers as they interact with sediments and move between the dissolved and particulate phase. To increase the usability of the code a pre- and post- processing GUI have been developed to aid the user in developing the inputs for the code and processing the code output into understandable figures for decision makers.

97 MATHEMATICS AND COMPUTING↗

Material Model Testing in the DYNA3D Code: Material Models 1 through 5

This report is the first volume in an effort to study the material model implementations in DYNA3D, using the latest developmental version of the code. Input files were derived from a template supplied by Jerry Lin to perform an equation-of-state study. Most of the simulations use 10 cubes, each 1 cm on a side with a density of 8.03 (steel). Each of the cubes tests a different aspect of the material model. The first three cubes are testing normal stresses and normal strains (uniaxial stresses and strains) placed on the cube parallel to the x-axis, y-axis and z-axis, respectively. The x (or y or z) = 1 face of the cube is moved along the x-axis (or y-axis or z-axis) to test the model, while the opposite side remains in the x (or y or z) = 0 plane. For some of the tests, the nodes are free to move perpendicular to the applied strain to allow the cross-sectional area of the cube to change in response to being placed under pressure or tension (to test Poisson’s Ratio), while in others the nodes are frozen to eliminate dynamic effects resulting from the code. The cube is placed in tension along the axis being tested and then released, followed by a cycle of compression and release. This tension/compression cycle is repeated 3 more times. The additional cycles allow for the model to be tested under different loads. For these calculations, the amount by which the cube is stretched and compressed is varied from cycle to cycle to allow testing of model features. Cubes number 4 and 5 are intended to test xy shear stresses/strains. This is done by sliding one face of the cube in its plane parallel to the axis of interest. For example, for cube 4, the upper face of the cube is slid in the y=1 plane in a direction parallel to the x-axis (so there are no changes in the y and z coordinates of the points in the plane) and then returned to its starting position. This preserves the volume of the cube, so there is no change in density or pressure. The cycle is repeated 8 times. For cube 5, the plane of the cube satisfying x=1 is slid in this plane in a direction parallel to the y-axis and then returned to its starting position. This preserves the volume of the cube, so there is no change in density or pressure. The cycle is repeated 8 times. Cubes number 6 through 9 repeat these tests for yz shear stresses/strains (Cubes 6 and 7) and zx shear stresses and strains (Cubes 8 and 9). For reasons discussed in the next paragraph, it is unclear if these tests represent a state of pure shear stress and shear strain in the code. Application of shear stresses and strains should not, in general, result in the appearance of normal stresses and strains in Cubes 4 - 9. In the infinitesimal limit, this result might be obtained. However, DYNA3D is a finite deformation code and in the real deformations used in the code, this result is not achieved and the individual components of the normal stresses and strains are non-zero during the cycling. However, the sums of the normal strains are zero in these tests and, except in the case of thermally dependent models, so are the sums of the normal stresses. Cube 10 tests the response of the material model to a hydrostatic loading of the cube. The load curve used for the first three tests is used again and applied to each of the 8 nodes in a manner causing the cube to be either compressed or expanded. This compression/tension cycle is repeated 4 times. The additional cycles allow for the model to be tested under different loads. The amount by which the cube is stretched and compressed is varied from cycle to cycle to allow testing of model features. This document will discuss the testing of Material Models 1 through 5. Relevant results from the test cubes will be discussed. A discussion of the results from all 10 test cubes will be presented only if necessary.

36 MATERIALS SCIENCE↗

Residential Energy Code Compliance in The Southeast: Georgia and Arkansas

The U.S. Department of Energy (DOE) Building Energy Codes Program launched a series of research studies in 2014 across several U.S. states to investigate typical energy efficiency practices in residential buildings and identify opportunities to reduce household utility bills. A primary goal of the study was to help document baseline practices in new single-family construction, target areas for improvement relative to the state energy code, and quantify related savings potential. This information was intended to assist states in measuring energy compliance with their codes and to identify areas of focus for future education and training initiatives. Following the original baseline studies, project teams in each state conducted a series of education, training and outreach activities aimed at improving compliance with the state residential energy code. A second baseline study was then implemented in each state with the objective of determining whether the investment in education, training, and outreach programs had produced a significant and measurable change in statewide residential energy use. Through the program, DOE planned to establish a sufficient data set to represent baseline construction trends across states and detect significant changes in energy use resulting from education and training activities. The three basic phases employed in the study were as follows: I. Pre-Study : A baseline study to identify the energy use in typical single-family residential buildings in a given state and opportunities for improving energy efficiency II. Education & Training : Education, training, and outreach activities targeting issues identified through the baseline study III. Post-Study : A second study to identify the change in energy use following the education and training activities This report outlines the application of the study in the states of Arkansas, Georgia, Tennessee, and Virginia with a focus on the experiences of the project team across the “Pre-Study” and “Education and Training” phases. Data gathered through Phase I has been submitted to Pacific Northwest National Laboratory (PNNL) for analysis. Findings available to date are published as PNNL technical reports on the DOE Building Energy Codes Program website with Post-Study analysis pending publication in 2019.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from high enrichment uranium fuel to low enrichment uranium fuel. The key capabilities of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilities are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL_4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow calculations are involved. This problem can be corrected in a later version of the code, however, does not affect the application of STAT7 code as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

MELCOR Code Change History: Revision 11932 to 14959 Patch Release Addendum

This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 11932 and 14959. Revision 11932 represents the last official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 14959. Along with the newly updated MELCOR Users' Guide and Reference Manual, users will be aware and able to assess the new capabilities for their modeling and analysis applications. Following the official release an addendum section has been added to this report detailing modifications made to the official release which support the accompanying patch release. The addendums address user reported issues and previously known issues within the official code release which extends the original Quick look document to also support the patch release. Furthermore, the addendums section documents the recent changes to input records in the Users' Guide applicable to the patch release and corrects a few issues in the revision 14959 release as well.

97 MATHEMATICS AND COMPUTING↗

RAZORBACK - A Research Reactor Transient Analysis Code Version 1.0 - Volume 1: User's Manual

Razorback is a research reactor transient analysis computer code designed to simulate the operation of a research reactor (such as Sandia National Laboratories Annular Core Research Reactor (ACRR)). The code provides a coupled numerical solution of the point reactor kinetics equations, the energy conservation equation for fuel element heat transfer, the equation of motion for fuel element thermal expansion, and the mass, momentum, and energy conservation equations for the water cooling of the fuel elements. This input manual describes how an input file is composed, and facilitates an understanding of the various code input parameters. The makeup of the various code output files is also described. This manual also provides instructions for the installation and setup of the code, and how to report bugs and/or errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.0 code which has been developed to perform the steady-state thermal hydraulic safety calculation of the MIT Research Reactor required upon conversion of the reactor from high enrichment uranium fuel to low enrichment uranium fuel. The key capabilites of the code including the statistical processing, water property generation, and thermal hydraulic solution are verified independently. While the first two capabilites are checked using statistical tests and NIST data tables, each component of the thermal hydraulic solution is verified using simple geometry test problems. The comparison between results obtained by the code and hand calculations on these problems shows that each of these components work correctly. Additionally, a comprehensive comparison between STAT7 Version 1.0 and PLTEMP/ANL Version 4.2 calculations for a five-channel test problem is performed for the validation. The good agreement between the obtained results for all thermal hydraulic parameters belonging to the fuel core region is achieved. An issue in bypass flowrate calculation was detected when both fin-friction and bypass flow claculations are involved. This problem can be corrected in a later version of the code, however, does not affect the application of STAT7 code as it has been used to date where the flow distribution in the core region and bypass channel is specified as an input parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Parallelization and Performance Portability in Hydrodynamics Codes

With the eve of Exascale computing, performance and portability are at the forefront of all scientific codes. Adding more cores and more energy to a system is no longer a sustainable way to achieve performance, and extra effort must now be made to improve performance in all areas of code and code development. Using hydrodynamic codes as a basis, this work explores numerous techniques to achieve performance in different ways. Adaptive mesh refinement (AMR) is a necessary technique to improve memory optimization in mesh-based simulations. However it is invasive and conventionally difficult to integrate into existing applications, so we present a new branch of AMR to create a smooth transition to these optimizations, which not only improves performance, but also greatly reduces developer effort. We introduce the concept of this improvement as Phantom-Cell AMR, and assess theoretically the improvements, as well as present an application of its use. Other work included involves and investigation into an efficient data structure that ensures optimal memory layout for cache performance, with a target of making codes performant and portable across all architectures. All of the work targets both performance and portability, not just on CPU hardware, but specifically across GPU architectures. Parallel performance is key to all of the methods presented, but the research makes a great effort to improve the portability of all applications to prepare for current high performance computing systems and those on the horizon.

97 MATHEMATICS AND COMPUTING↗

ESS-Related Proposed Changes to the Group B 2018 ICC International Codes

The purpose of this special briefing paper is to support the objective by providing information about current and upcoming efforts being conducted by U.S. standards and model code developing organizations, specifically code changes associated with the 2018 Group B International Codes (I-Codes) of the International Code Council (100)2 that are relevant to Energy Storage System Safety (ESS)s .

42 ENGINEERING↗

Using Verified Lifting to Optimize Legacy Stencil Codes (Final Project Report)

This project investigated new techniques for compiling stencil and stencil-like computations. Stencils computations are commonly found in applications such as image processing, physical simulations, image processing, and machine learning. In recent years, many high-performance domain-specific languages (DSLs) have been proposed to optimize stencil computations. To leverage such DSLs, however, existing codes often need to be rewritten. Such rewriting is manual, labor intensive, and error-prone. To alleviate such issues, this project investigated the application of program synthesis and artificial learning techniques to enable stencil computations to automatically leverage new high-performance DSLs. Rather than constructing syntax driven rules, verified lifting uses program synthesis to search for a target code fragment to compile the given input code into. In addition, it also searches for a proof that validates how the found target code fragment preserves the semantics of the original input. Thus, the target code fragment is guaranteed to be semantically equivalent to the input.

97 MATHEMATICS AND COMPUTING↗

CodeFlow: A Code Generation System for Flash-X Orchestration Runtime

We propose the CodeFlow toolchain for Flash-X that realizes the “recipe-to-source” code transformation for Flash-X simulations and that is necessary to achive performance portability. We design a high-level language to express operations of simulations in so-called recipes, which are given as input to the toolchain. The tools of the CodeFlow pipeline include code transformation with tree-based source code representation techniques and code orchestration and generation based on control flow graphs. The generated source code utilizes a new runtime, developed for Flash-X, that orchestrates dynamic and asynchronous data movement and task execution. The functionality of CodeFlow is demonstrated using a hydrodynamic problem with a strong shock.

97 MATHEMATICS AND COMPUTING↗

Numerical Codes for the DESC-LSST Analysis Pipeline: Core Cosmology Library Standard Modules and Beyond wCDM Modules (Final Technical Report)

The overall objective of the project is to investigate and develop specific software modules and analysis components for the software pipeline of LSST Dark Energy Science Collaboration (DESC). Following the key projects of DESC Science Roadmap (SRM), we will write and test computer codes for the Core Cosmology Library (CCL) in order to complete its modules, functionalities, and interface to work with the analysis pipelines from the five science probes of DESC (parts of SRM deliverables CX4.2TJP, CX6.2CS). We will also code modules for CCL to test models beyond w-Cold-Dark-Matter (wCDM) and modification to gravity (MG). Interfaces for MG models will also be developed for the TJPCOSMO software which is the main pipeline of the Theory and Joint Probe (TJP) working group (deliverable TJP2.3). In order to use the full power of LSST data to constrain MG models, we will also work on constraints from nonlinear regime by running and analyzing MG N-Body simulations using a Parameterized-Post-Friedmann framework into the Gadget-2 simulation package (deliverable TJP2.2). Preliminary results for the simulations were obtained in the past. In collaboration with other DESC groups, we plan to make these simulations feedable to cosmic emulators that are practical for likelihood analyzes (deliverable TJP2.2, parts of CX6.2CS). We will also modify and integrate our current codes for consistency tests between data sets and probes into the pipeline (parts of deliverables CX8.2TJP, TJP2.3). We understand that other groups will contribute to some of these objectives but our team will focus and collaborate with others on the particular part of testing MG and models beyond wCDM and refine the DESC pipeline for this purpose. PI has been coordinating his work with the TJP and CS working groups and the DESC management team. PI is a full member of DESC since June 2013. He and his students have been contributing to LSST-DESC activities and work including TJP telecons, collaboration meetings, hack-weeks, and workshops. PI chaired or co-chaired sessions at collaboration meetings and hack-weeks about testing gravity and models beyond wCDM using LSST. He is coordinating the TJP2 projects for testing models beyond wCDM including the writing of DESC-research-note, development of code for pipeline, and N-Body simulations for MG and beyond wCDM models testable with LSST analyses. As stressed in the DESC white paper, SRM, and P5 report, one of the important questions in understanding cosmic acceleration and dark energy is to be able to distinguish whether the acceleration is due to a dark energy component in the universe or a modification to gravity. Answering these questions will have a significant impact on the question of cosmic acceleration and dark energy. The methods that we will use include analytical work, numerical code, and N-Body simulations. A first approach that that we will use consists of using growth rate parameters that enter the perturbed dynamics equations. These parameters take distinctive values for distinct gravity theories and have potential to distinguish between Dark Energy and Modified Gravity. The second method is to look for inconsistencies in Dark Energy parameter spaces using specific combinations of cosmological data sets. Our investigation addresses the Dark Energy problem that is relevant to the mission of the HEP program to understand how our universe works at its most fundamental level. It will allow us to make progress on the HEP mission to explore the nature of Dark Energy and the basic nature of space and time using future surveys such as LSST. The investigation supports the DOE HEP program Cosmic Frontier as it will contribute to the study and understanding of dark energy and fundamental properties of the universe. The investigation contributes directly to LSST-DESC key projects and their deliverables as described in the Science Road-map document to build analysis pipeline and to test dark energy and beyond wCDM models using LSST.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Recent MCNP6 ® Code Developments and Improvements for Nuclear Engineering Applications [Slides]

The Los Alamos MCNP Monte Carlo radiation transport code has been the international gold standard for particle transport applications for over three decades. Many developments to the code have taken place with several significant new feature additions, major improvements, and enhancements to existing features. With significant institutional and programmatic investment in the code since the time of the last public release in 2018, important code development and infrastructure modernization has taken place and remains a high priority for all ongoing efforts across the code development team.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational Modeling and Simulation for Nonproliferation: The History (and Present) of Monte Carlo and the MCNP(R) Code at Los Alamos [Slides]

The emergence of the Monte Carlo method as a research tool springs from work done at Los Alamos in the 1940s.Monte Carlo and the MCNP Code have been and continue to be developed at Los Alamos for many decades. From basic science in support of understanding nuclear interaction physics to global security applications, application uses of the code are extensive. Recent R&D projects and code modernization efforts make the MCNP code a great tool for nuclear nonproliferation applications. In collaboration with nuclear safeguards experts, new training has just recently been developed to help new practitioners learn how to use the code for nuclear safeguards applications.

97 MATHEMATICS AND COMPUTING↗

Interim Mechanical Properties Data from FY2022 ORNL Testing of A709 with Precipitation Treatment for ASME Code Case Data Package

A collaborative research and development effort in support of the Alloy 709 Code Case qualification in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code, Section III, Division 5, High Temperature Reactors is being conducted at the US Department of Energy’s Oak Ridge National Laboratory, Idaho National Laboratory, and Argonne National Laboratory. A recent assessment on the Alloy 709 development and testing effort concluded that the mechanical properties of Alloy 709 with the precipitation treatment continued to outperform those of Type 316 stainless steel. The assessment also affirmed the recommendation for its Code qualification. This report provides the detailed integrated data generated to date on the two commercial heats of Alloy 709 in plate product form with precipitation treatment condition that were used as the basis for this assessment. Oak Ridge National Laboratory has continued to perform a subset of the Code Case testing for tensile properties, creep rupture, fatigue, and creep-fatigue. This report also updates the key Alloy 709 Code Case testing status and results in FY 2022 at Oak Ridge National Laboratory.

36 MATERIALS SCIENCE↗

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗

Extension of Plant Dynamics Code Capabilities for Simulation of TerraPower Pascal Reactor

Argonne National Laboratory has been developing the Plant Dynamics Code (PDC) for design and transient analysis of supercritical carbon dioxide (sCO 2 ) Brayton cycles. In previous analyses with PDC, only indirect sCO 2 cycles, where the heat is being added through a heat exchanger, such as sodium-to-CO 2 HX, were analyzed. Under the U.S. Department of Energy Technology Commercialization Fund (TCF), Argonne cooperated with TerraPower to bring the Plant Dynamics Code to commercial market. The main focus of the TCF project is to extend the application base and the code usability by developing the capabilities to be able to simulate reactor systems with direct sCO 2 cycles. These new PDC capabilities have been developed in application to the TerraPower’s Pascal reactor concept. This report documents the Pascal reactor modeling with the PDC, including simulation of the Pascal split-expansion cycle, the development of the reactor module in the PDC, modeling of the Pascal’s shutdown heat removal system in PDC, and other updates to the code. The report also describes the results of the steady state and transient demonstration of the newly developed code features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗