SEARCH · Search NASA
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.
A framework for developing scientific computing codes for projects was developed that accommodates both explorations as well as, simultaneously, software engineering rigor. This is accomplished by organizing modular yet composable and related computational environments. Also known as Management of Multiple Computational Environments
The Second-Level Adjoint Sensitivity System (2nd-LASS) that yields the second-order sensitivities of a response of uncollided particles with respect to isotope densities, cross sections, and source emission rates is derived in Refs. 1 and 2. In Ref. 2, we solved problems for the uncollided leakage from a homogeneous sphere and a multiregion cylinder using the PARTISN multigroup discrete-ordinates code. In this memo, we derive solutions of the 2nd-LASS for the particular case when the response is a flux or partial current density computed at a single point on the boundary, and the inner products are computed using ray-tracing. Both the PARTISN approach and the ray-tracing approach are implemented in a computer code, SENSPG. The next section of this report presents the equations of the 1st- and 2nd-LASS for uncollided particles and the first- and second-order sensitivities that use the solutions of the 1st- and 2nd-LASS. Section III presents solutions of the 1st- and 2nd-LASS equations for the case of ray-tracing from a detector point. Section IV presents specific solutions of the 2nd-LASS and derives the ray-trace form of the inner products needed for second-order sensitivities. Numerical results for the total leakage from a homogeneous sphere are presented in Sec. V and for the leakage from one side of a two-region slab in Sec. VI. Section VII is a summary and conclusions.
The objective of this project is to develop, and integrate, high-performance simulation tools capable of predicting plasma-facing component (PFC) operating lifetime and the impact of the evolving surface morphology of tungsten-based PFCs on plasma contamination, including the dynamic recycling of fuel species and tritium retention, in future magnetic fusion devices. Establishing a fundamental physical understanding and developing predictive capabilities of plasma-surface interactions (PSI) requires simultaneously addressing complex and diverse physics occurring over a wide range of length (Angstroms to meters) and time (femtoseconds to years) scales, as well as integrating extensive physical processes across the plasma–surface interface. This requires development of not only detailed physics models and computational strategies at each scale, but also algorithms and methods to couple them effectively in a way that can be robustly validated. Deploying these tools requires the continued development and coupling of leadership-scale computational codes to describe the boundary plasma and the evolving PFC surface, as well as a host of simulations that bridge disparate scales to address complex physical and computational issues at the plasma–surface interface in multi-component materials systems for magnetic fusion energy development beyond ITER.
Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE-HDBK-3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms from postulated accident scenarios. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this fourth year, we improved existing computational capabilities to better model fragmentation situations to capture small fragments during an impact accident. In addition, we have provided additional new information for various sections of Chapters 4 and 5 of the Handbook on free fall powders and impacts of solids, and have provided the damage ratio simulations for containers (7A drum and standard waste box) for various drops and impact scenarios. Thus, this work provides a low-cost method to establish physics-justified safety bounds by considering specific geometries and conditions that may not have been previously measured and/or are too costly to perform during an experiment.
This report documents comparisons between Fuel Model (PARFUME) model predictions versus Bison for selected compacts from the second irradiation test of the Advanced Gas Reactor (AGR) program that occurred from June 2010 to October 2013 in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). PARFUME is a fuel performance analysis and modeling code, used for evaluating gas-reactor tristructural isotropic (TRISO) coated particle fuel for prismatic, pebble bed, plate, and cylindrical type fuel geometries. PARFUME is an integrated mechanistic computer code that evaluates the thermal, mechanical, and physico-chemical behavior of TRISO coated-fuel particles and the probability for fuel failure given the particle-to-particle statistical variations in physical dimensions and material properties that arise during the fuel fabrication process. Bison is a nuclear fuel performance application built using the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library (Permann 2020). Bison is capable of modeling multiple fuel forms in a wide variety of dimensions and geometries. It solves coupled nonlinear partial differential equations, including heat conduction, mechanics, fission product species transport etc., in a fully implicit manner. Comparisons between PARFUME and Bison were performed using four compacts from the AGR-2 experiment. Selected outputs were chosen based on their impact on the probability of the SiC layer failing. In general there was good agreement between PARFUME and Bison with the exception of predicting the gap formed between the buffer and IPyC. Further comparisons between PARFUME and Bison are planned to further develop Bison’s capabilities.
To support the development and utilization of the SAM code for molten-salt-fueled reactor (MSR) safety analysis and licensing, an effort was devoted to enhancing code capabilities and developing reference models for the MSR primary loop. A reference standard problem of a prototypical reactor design is foundational to NRC to verify the adequacy of computer codes and evaluation models for a specific reactor type. A thermo-fluid model of the Molten Salt Reactor Experiment (MSRE) has been developed. The MSRE primary loop model consists of a 2-D core region and external core components in 1-D or 0-D. Both the steady-state and a transient scenario are simulated. The porous medium model is utilized for both the molten salt fluid channels and the moderator matrix in the MSRE core. The use of the porous medium model for the MSRE core is verified with the higher-fidelity simulation results using 1-D representations of the fluid channels and 3-D modeling of core structures. To further enhance SAM multi-scale simulation capabilities for MSRs, a mass transport model is developed and implemented in SAM multi-dimension flow module for modeling of species transport such as delayed neutron precursors in MSRs.
Research during this project at the University of Wisconsin Fusion Technology Institute (UW FTI) on ion and neutral flow through an arbitrary, monotonic potential difference created by nearly transparent electrodes accomplished the following: (1) developed and implemented an integral equation approach for atomic physics effects in helium plasmas; (2) extended the analysis to coupled integral equations that treat atomic and molecular deuterium ions and neutrals; (3) implemented the key deuterium and helium atomic and molecular cross sections; (4) added negative ion production and related cross sections; and (5) benchmarked the code against experimental results. The analysis and codes treat the species D 0 , D 2 0 , D + , D 2 + , D 3 + , D - and, separately at present, He 0 and He + . Extensions enhanced the analysis and related computer codes to include He ++ ions plus planar and cylindrical geometries.
This repository contains example implementations of phase-field codes primarily for educational purposes. These codes demonstrate how a simple phase-field models can be implemented as an introduction for researchers new to the area (students, new staff, collaborators from different domains, etc.). The computer code includes highly documented Jupyter notebooks that implement phase-field models and analyze the results.
We obtain here a solution to a zone-centered discretization of the one dimensional time-dependent diffusion equation with arbitrary initial conditions and source, constant absorption and scattering opacity, and constant zone size and time step. The solution is obtained using the discrete Green’s functions of the discretized equation. The solution of the discretized equation is useful in the testing of computer codes, because the code can be expected to agree with the solution to the discrete equation, to within small errors caused by roundoff. This is in contrast to solutions of the differential equation, with which the code results only approximately agree. Finally, the usefulness of the solution for tests of an inertial confinement fusion code is demonstrated.
The AXIOM-Unfold application is a computational code for performing spectral unfolds along with uncertainty quantification of the photon spectrum. While this code was principally designed for spectral unfolds on the Saturn source, it is also relevant to other radiation sources such as Pithon. This code is a component of the AXIOM project which was undertaken in order to measure the time-resolved spectrum of the Saturn source; to support this, the AXIOM-Unfold code is able to process time-dependent dose measurements in order to obtain a time-resolved spectrum. This manual contains a full description of the algorithms used by the method. The code features are fully documented along with several worked examples.
The SCALE code system developed at Oak Ridge National Laboratory is widely used and accepted around the world for criticality safety analyses. The well-known KENO V.a three-dimensional Monte Carlo criticality computer code is one of the primary criticality safety analysis tools in SCALE. The KENO V.a primer is designed to help a new user understand and use the SCALE/KENO V.a Monte Carlo code for nuclear criticality safety analyses. It assumes that the user has a college education in a technical field. There is no assumption of familiarity with Monte Carlo codes in general or with SCALE/KENO V.a in particular. The primer is designed to teach by example, with each example illustrating two or three features of SCALE/KENO V.a that are useful in criticality analyses. The primer is based on SCALE 6.2 and 6.3, which includes the Fulcrum graphical user interface (GUI). Each example uses Fulcrum to provide the framework for preparing input data and viewing output results. Starting with a Quickstart section, the primer gives an overview of the basic requirements for SCALE/KENO V.a input and allows the user to quickly run a simple criticality problem with SCALE/KENO V.a. The sections that follow Quickstart include a list of basic objectives at the beginning that identifies the goal of the section and the individual SCALE/KENO V.a features that are covered in detail in the sample problems in that section. Upon completion of the primer, a new user should be comfortable using Fulcrum to set up criticality problems in SCALE/KENO V.a. The primer provides a starting point for the criticality safety analyst who uses SCALE/KENO V.a. Complete descriptions are provided in the SCALE/KENO V.a manual. Although the primer is self-contained, it is intended as a companion volume to the SCALE/KENO V.a training and documentation. The SCALE manual and training schedule are available at https://scale.ornl.gov. The primer provides specific examples of using SCALE/KENO V.a for criticality analyses; the SCALE/KENO V.a manual provides information on the use of SCALE/KENO V.a and all its modules. The primer also contains an appendix with sample input files. In addition, this primer, its errata and sample inputs are also available at https://code.ornl.gov/scale/primers/kenova.
High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems for design, optimization, and licensing calculations requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the OECD/NEA Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). This benchmark contains problems for systems code, computational fluid dynamics (CFD), and coupled systems code/CFD modeling representing lower plenum mixing and both the depressurized and pressurized conduction cooldown (DCC and PCC respectively) transients. Benchmark problems include exercises for code-to-code and code-to-data comparisons as well as an exercise for error scaling between HTTF and the Modular High Temperature Gas-Cooled Reactor, which serves as the basis for the HTTF design. Previous analysis as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. The new model provides a finer nodalization to assess the impact of radial nodalization and allows for asymmetric heating within the core, which was a feature of multiple HTTF experiments. In this paper, we present the new RELAP5-3D model of HTTF. In addition to describing the new model, this paper compares the new and old models and provides results for a full-power steady state, a DCC, and a PCC in HTTF. These analyses are based on the code-to-code comparison exercises for the DCC and PCC problems of the HTGR T/H benchmark. We present the results of these exercises from the new model and compare them to the results of the old model.
The NJOY Nuclear Data Processing System, version 2016, is a comprehensive computer code package for producing pointwise and multigroup cross sections and related quantities from evaluated nuclear data in the ENDF-4 through ENDF-6 legacy card-image formats. NJOY works with evaluated files for incident neutrons, photons, and charged particles, producing libraries for a wide variety of particle transport and reactor analysis codes.
The SCALE code system developed at Oak Ridge National Laboratory is widely used and accepted around the world for criticality safety analysis. The well-known KENO-VI three-dimensional Monte Carlo criticality computer code is one of the primary criticality safety analysis tools in SCALE. The KENO-VI primer is designed to help a new user understand and use the SCALE/KENO-VI Monte Carlo code for nuclear criticality safety analysis. It assumes that the user has a college education in a technical field. There is no assumption of familiarity with Monte Carlo codes in general or with SCALE/KENO-VI in particular. The primer is designed to teach by example, with each example illustrating two or three features of SCALE/KENO-VI that are useful in criticality analysis. The primer is based on SCALE 6.2 and 6.3, which includes the Fulcrum graphical user interface. Each example uses Fulcrum to provide the framework for preparing input data and viewing output results. Starting with a Quickstart section, the primer gives an overview of the basic requirements for SCALE/KENO-VI input and allows the user to quickly run a simple criticality problem with SCALE/KENO-VI. Each following section begins with a list of basic objectives identifying the goal of the section and the individual SCALE/KENO-VI features covered in detail in the section’s sample problems. Upon completion of the primer, a new user should be comfortable using Fulcrum to set up criticality problems in SCALE/KENO-VI. The primer provides a starting point for the criticality safety analyst who uses SCALE/KENO-VI. Complete descriptions are provided in the SCALE/KENO-VI manual. Although the primer is self-contained, it is intended as a companion volume to the SCALE/KENO-VI training and documentation. The SCALE manual and training schedule are available at https://scale.ornl.gov. The primer provides specific examples of using SCALE/KENO-VI for criticality analysis; the SCALE/KENO-VI manual provides information on the use of SCALE/KENO-VI and all its modules. The primer also contains an appendix with sample input files. In addition, this primer, its errata, and sample inputs are also available at https://code.ornl.gov/scale/primers/kenovi/.
This paper describes recent work to develop and incorporate a utility program into the MCNPTools package that enables LNK3DNT file visualization & post-processing. A LNK3DNT file is a binary Cartesian, cylindrical, or spherical structured mesh file that describes material and material relative density distributions by mesh voxel. This mesh file format is principally used by the PARTISN discrete-ordinates radiation-transport code. However, the MCNP ® Monte Carlo radiation-transport code can use LNK3DNT files to specify geometry, and the MCNP code can produce LNK3DNT files representing a user-defined region of the MCNP constructive solid geometry (CSG). Because of the use of LNK3DNT files as MCNP input and output and with other computer codes, improving the ability to visualize and/or quantitatively post-process LNK3DNT files benefits end users. In particular, it is helpful to visually verify that the material and density distributions are as expected. Also, an analyst can orient viewers with respect to other in formation, e.g., by co-plotting a voxelized representation of MCNP CSG with mesh tally and/or particle track results. However, until now, the only capability to interactively interrogate LNK3DNT files is through a 2-D slice plotter provided by the Monte Carlo Application ToolKit (MCATK), which is not currently distributed outside of Los Alamos National Laboratory. The utility created as part of this work converts LNK3DNT binary files to an ASCII XML-based StructuredMesh VTK (.vts) file. The geometries supported are 1-D spherical and cylindrical ( r ), 1-D Cartesian ( x ), 2-D cylindrical ( r-z and r-θ ) and Cartesian ( x-y ), and 3-D cylindrical ( r-z-θ ) and Cartesian ( x-y-z ). This file format is chosen to minimize external dependencies (e.g., on an HDF5 library) and to provide a file that can be interrogated by the user in a text editor or standard visualization applications (e.g., ParaView or VisIt). However, other file formats may provide more utility such as XDMF, which can use HDF5 and is discussed further as something to explore as future work. This paper is organized as follows: a brief introduc tion to the MCNPTools library and utility programs is given, the LNK3DNT-to-VTK (l3d2vtk) utility functionality is de scribed, and then an example post-processing workflow using the ParaView application is given.
We report over the past decades, several computer codes have been developed for simulation and analysis of thermal-hydraulics and system response in nuclear reactors under operating, abnormal transient, and accident conditions. However, simulation errors and uncertainties still inevitably exist even while these codes have been extensively assessed and used. In this work, a data-driven framework (Optimal Mesh/Model Information System, OMIS) is formulated and demonstrated to estimate simulation error and suggest optimal selection of computational mesh size (i.e., nodalization) and constitutive correlations (e.g., wall functions and turbulence models) for low-fidelity, coarse-mesh thermal-hydraulic simulation, in order to achieve accuracy comparable to that of high-fidelity simulation. Using results from high-fidelity simulations and experimental data with many fast-running low-fidelity simulations, an error database is built and used to train a machine learning model that can determine the relationship between local simulation error and local physical features. This machine learning model is then used to generate insight and help correct low-fidelity simulations for similar physical conditions. The OMIS framework is designed as a modularized six-step procedure and accomplished with state-of-the-art methods and algorithms. A mixed-convection case study was performed to illustrate the entire framework.
Sandia and Cray Inc. co-developed Red Storm, a distributed memory, massively paralleled highperformance supercomputer modeled on ASCI Redl, to run computer codes used for conducting materials science simulations for national security. Supercomputers have some of the fastest highperformance systems available and are used primarily for scientific and engineering work requiring exceedingly high-speed computations. Unlike conventional computers, supercomputers have large storage capacity; more than one central processing unit to rapidly retrieve stored data and program instructions; and input/output capability.
In the frame of their research programs on fuel safety, the French “Institut de Radioprotection et de Sûreté Nucléaire” (IRSN) and the Idaho National Laboratory (INL) have developed respectively SCANAIR and BISON codes to describe the thermo-mechanical behavior of irradiated fuel rods during Reactivity Initiated Accidents (RIA) in LWR. A RIA, characterized by a very rapid increase of reactivity and power in some rods of the reactor, can be schematically represented in two main phases. First, the energy deposition leads to a rapid rise of the fuel temperature which induces thermal swelling of the fuel pellets. During this phase, the Pellet Clad Mechanical Interaction (PCMI) leads to clad deformation and potentially to failure depending on fuel enthalpy increase and on the level of clad embrittlement. After the PCMI phase, the increase of clad temperature can lead to the boiling crisis of the water surrounding the rod. During the film boiling phase, the clad to coolant heat transfer becomes very low and the clad can reach high temperature (>700°C). Depending on the internal gas pressure, the ductile clad can undergo large deformation and possible failure. The RIA codes have to predict properly these two phases of the transient in order to be used, in particular, to assess current fuel safety criteria but also to be able to predict the behavior of upcoming fuel (such as accident tolerant fuel). In this paper, codes computation comparison against experimental data obtained on CIP0-1 test is done for the first phase of the transient, the PCMI phase, in order to point out the strengths and weaknesses of the codes. The CIP0-1 experiment was performed with UO2 high burnup fuel in the CABRI sodium loop facility, at 280°C and low pressure (~3 bars) without boiling crisis and no failure of the fuel rod. The initial state of the rod calculated by irradiation codes are compared with measurements before the transient. Then global rod thermal behavior during the CIP0-1 test is assessed with the sodium coolant temperature measurements. The rod mechanical behavior is analyzed with the fuel and clad elongation and the clad residual hoop strain. Finally comparison of the fission gas behavior between codes is performed and compared to the measurements of the fission gas released in the experiment.