SEARCH · Engineering Papers
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.
This presentation outlines improvements in the Monte Carlo N-Particle Code (MCNP) version 6.3. MCNP is an internationally renowned Monte Carlo particle transport code developed by Los Alamos National Laboratory.
The EOSPAC utility package is a collection of interface routines, which can be used to access the SESAME data library and perform various data adjustments and interpolations on the SESAME data. The SESAME data library contains both thermodynamic (e.g., equation of state) and transport coefficients (e.g., opacity and conductivity). Note, for simplicity, the term EOS (equation of state) used herein includes both thermodynamic variables and transport coefficients. The EOSPAC utility package is designed to be used by physics codes (henceforth ”host codes”) written in multiple languages and on multiple platforms. The remainder of this manual is organized into several sections. Chapter 2 discusses conventions such as data organization and routine names. Chapter 3 provides a general overview of basic theory and models implemented within EOSPAC. Chapter 4 provides a general overview of how to use the EOSPAC interface library. Chapters 5 to 7 describe the public interfaces of EOSPAC in detail. Chapter 8 provides a brief introduction to some related tools, which may be of use to the user. Chapter 9 provides details related to some selected numerical features of EOSPAC. Chapter 10 gives examples for using the interface routines described in chapters 5 to 7. Chapter 11 provides technical support contact information. Chapter 12 contains a brief set of acknowledgments. Chapter 13 contains a list of referenced documents. Finally, chapter 14 lists the “table types: mnemonic conventions”, “table types: grouped by category, sorted by name”, “table types: eospac version 5 cross reference”, “options: setup phase”, “data information parameters”, “meta-data information parameters”, “options: interpolation phase”, and the “error codes”.
Random number generators are integral components to Monte Carlo codes. They provide the pseudorandom number sequence used to actually sample the distributions of interest. As a result, they are one of the most important components to the software. The current recommended MCNP random number generator is a 63-bit linear congruential generator (LCG). This generator is quite fast, but it has some drawbacks. First, it only has a period of 2 63 . Due to the necessarily non-optimal usage of random numbers to ensure parallel reproducibility, this amount is too few to guarantee random number sequences are not reused in all configurations the code runs under. As simulation size increases, users will need to be aware of the limitations of the generator and tune configuration variables to best suit their simulations, or they will need to assume that reuse is not negatively affecting their answers. Neither of these are optimal. Second, small LCGs are fairly weak in bit generation quality, and this can have an unknown impact on the quality of the simulation. This paper is an investigation into whether or not more modern random number generators can supersede the current ones. The goal is to find a generator that is similar or superior in speed to the LCGs, has a state space large enough to make strong guarantees about random number reuse, and passes all modern random number test suites. If such a generator is found, it would eliminate the need for the user to even be aware of the limitations of the random number generator and would simplify the use of the code. This paper will be broken into several parts. Sec. 2 will discuss the evolution of the random number generator within the MCNP code. Sec. 3 will go over what a Monte Carlo code needs from a generator to be reproducible and portable and how the current generator behaves in that light. Sec. 4 goes through how each generator was tested. Finally, Sec. 5 will discuss improvements that could be made to the code.
We develop the “integrated code” FLAG, incorporating a variety of physics/meshing/software packages into one multi-physics code.
REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.
H3D Inc. has introduced a new type of CZT detectors that has resolution several times better than that of a typical, single crystal CZT detector. The better resolution has made it easier for software codes to analyze its spectra for isotopic composition of nuclear materials. The Fixed-energy Response-function Analysis with multiple efficiency (FRAM) is one of such software codes. The good peak resolution of the H3D CZT detectors allows the current FRAM v.6.1 and earlier to successfully analyze uranium spectra with minimal modification of the parameter sets and get reasonable results. (A parameter set governs how an analysis is performed.) The exceptionally large high-energy tail of the peaks makes it difficult to achieve good uranium results and impossible to successfully analyze plutonium spectra. The FRAM code needed to be modified to better fit the high-energy tail of the peaks in order to analyze plutonium data and to get better uranium analysis results. Results of this modification, encoded in beta version FRAM v.7.0i, are reported herein.
We present an overview and description of the atmospheric transport model (ATM) component of the Defense Land Fallout Interpretative Code (DELFIC) and its revision history. In particular, we discuss the hybrid Eulerian-Lagrangian model of the DELFIC diffusive transport module (DTM), its assumptions and constraints, and assess its capability for modeling the transport, dispersion, and deposition of nuclear fallout debris. This hybrid model is an efficient transport model for the wide range of particle sizes found in nuclear fallout debris. Furthermore, we validate the 1979 version of DELFIC using the 1979 test case input and output presented in its documentation. Except for some rounding differences and exposure rate differences from the OPM, the DTM output is in very good agreement with the 1979 results. Overall, the DELFIC DTM provides a reasonable and rapid approximation to the transport and deposition of local nuclear fallout debris.
The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.
For decades, Los Alamos National Laboratory has been at the forefront of neutron transport methods research and code development. One such code is PARTISN, the LANL parallel time-dependent discrete ordinate neutron transport code. In this presentation, we describe the various research efforts currently underway by the PARTISN and other code teams. Some examples of current research are a block automated mesh refinement scheme, the application of tensor trains to the discretized neutron transport equation, and GPU code porting. The block automated mesh refinement scheme uses cross section information to refine and coarsen the solution mesh to improve time to solution and reduce memory. The tensor train approach expresses discretized transport operators as tensor products of vectors and matrices to compress the size of linear systems being solved by transport codes. Rather than relying on matrix-free methods such as the transport sweep, we have access to an operator that can be inverted, reshaped, or manipulated algebraically. Finally, we describe how PARTISN is used, what problems we are looking to solve, and what the future holds for neutron transport at LANL. In addition to this, we briefly describe the various research efforts in other particle transport teams using both deterministic and Monte Carlo methods. In the presentation, we list possible opportunities for collaboration between the laboratory and faculty and students.
Abstract not provided.
The GP-SANS, Bio-SANS and EQ-SANS instruments at ORNL utilize drtsans for data reduction. drtsans is built on Python, and it can be run using python scripts and Jupyter notebooks. The flexibility afforded by Python makes it possible to incorporate additional actions into the scripts used for data reduction, such as analysis and visualization. Here, a new set of tools for visualizing and manipulating SANS data that can be incorporated into the data reduction scripts for the ORNL SANS instruments, or employed during post–processing, is presented that expands the capabilities of the two previously-released tool sets.
An improvement to the handling of contributions from neutron inelastic scattering to next-event estimators has been implemented in the MCNP6 ® software for release with version 6.3.1. The kinematic equations that govern the outgoing energy of inelastic neutron scattering contributions to next-event estimators have two roots. Historically, the implementation in the MCNP coding has only used the upper root to the quadratic equation. This includes all versions predating version 6.3.1 including all versions of MCNP5 and MCNPX software. However, a review of the neutron next-event estimator physics has shown that this does not reproduce the track-length estimator results at low energies. Several examples are presented that test single neutron inelastic scattering reaction types, Level Scattering (Law 3), Tabulated Energy Angle (Law 61), and Kalbach-Mann distribution (Law 44). The test problems compare the track-length estimator (f4 tally) in the MCNP software, with the existing implementation of the neutron next-event estimator (f5 tally), and the modified changes to the neutron next-event estimator implementation. As the MCNP implementation ignores the lower root, the unmodified neutron next-event estimators will generally underestimate the lower energy contribution. However, a second issue with the Kalbach-Mann distribution (Law 44) implementation allows contributions to backward scattering in the center-of-mass frame that is not kinematically possible, thus overestimating backward scattering contributions. A third issue with the way the MCNP implementation handles floating point comparison for scattering directly ahead or directly backward in the center-of-mass frame generally leads to underestimation (except for backwards scattering for Law 44).
An evaluation of n+ 88 Sr reactions up to 950 keV was compiled from a publicly available set of resonance parameters. Starting from the repository of the ENDF/B-VIII.0 library, an evaluated nuclear data file was assembled and submitted for inclusion to the Evaluated Nuclear Data File (ENDF)/B-VIII.1 release. This evaluation features an external function representation defined by pole strengths and R ∞ parameters over the entire resolved resonance range. This work improves the modeling of the R-matrix parameterization from multi-Level Breit-Wigner resonance parameters to Reich-Moore formalism, and it also extends thrice the upper energy range of the current ENDF/B-VIII.0 resolved resonance evaluation.
This presentation touches on thermal scattering cross section and scattering lengths including criticality benchmarks (fission). Secondly, development and completeness of charged particle libraries: criticality benchmarks (fission), neutron source (fission/fusion), simulation benchmarks and light elements generation (fusion). Additionally, this presentation details on the covariances (methodology and completeness) with adjustment and optimization of nuclear data libraries (fission/fusion). Finally nuclear data reproducibility and conclusions are discussed.
This presentation finds that INS measurements are not as variable based on sample composition, because phonon spectra is a properly of the bulk material. Additionally, INS and transmission measurements are and should be the most significant part of TSL validation, as they are the fundamental quantities used by Monte Carlo neutronics codes. Relaying solely on theoretical atomistic calculations for TSL evaluation is not enough, they are highly dependent on inter-atomic potentials (which are not the most accurate depending on the quantity trying to be reproduced) and evaluators knowledge and understanding of the material being studied. When dealing with Crit. benchmarks - While extremely useful for validation of nuclear data at all energies, critical benchmarks are not the best tool to provide a definitive answer on conflicting TSLs. Finally, INS and transmission measurements need to be the basis of validation of TSLs.
Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem
Correlated sampling using batch statistics with MCNP6’s tally fluctuation charts (TFCs) is a powerful means of reducing the statistical uncertainties of various combinations of tallies. Because MCNP6 prints a maximum of 20 entries in the TFCs, at least five inputs must be run to obtain the recommended number of at least 100 batches. A new software tool, MAKE_COSUBS, reads an MCNP6 input file and sets up a sequence of input files to be run. The user runs MAKE_COSUBS for the unperturbed and all perturbed MCNP6 inputs, runs MCNP6 for the inputs, and finally analyzes the outputs using COSUBS.
Many thermal radiative transfer (TRT) test problems have been introduced over the years. Here, we will combine the good features of several tests into one new problem. Our goals are to have a multi-dimensional, multi-group TRT problem that stresses the physics and numerics of codes in a realistic way, but the geometry and materials are idealized just enough to be simple to set up in any code. The proposed problem has no connection to any physical experiment; it is simply meant to have a combination of properties that stress the solvers with both the optically thick and thin limits. To achieve this, we are willing to give up any form of analytic solution.