Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ENDF”

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 199 records · Page 11

Open Source Release of NJOY2016 and NJOY21

NJOY has been the standard for nuclear data processing for many decades. NJOY is a comprehensive computer code used to produce point-wise and multigroup nuclear data application libraries from cross sections and related quantities obtained from evaluated nuclear data in the Evaluated Nuclear Data Format (ENDF). NJOY was first publicly released in 1977 and distributed by what is now the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory and the National Energy Software Center (NESC) at Argonne National Laboratory. From that time until the release of NJOY2012, NJOY has been distributed by RSICC and made available to researchers outside the United States by agreements between RSICC and the NEA Databank. NJOY2012 was distributed by the Feynman Center at Los Alamos National Laboratory which did not have agreements with the NEA Databank, thus making it more difficult for non US citizens to obtain copies of NJOY2012. Furthermore, NJOY has had export control restrictions preventing some researchers from obtaining NJOY at all. Two new versions of NJOY have been released, NJOY2016 and NJOY21, in January of 2017. Both NJOY2016 and NJOY21 have been released with no export control restrictions and are distributed with an open-source license. This is a big departure from previous releases and impacts not only who has access to NJOY, but also how updates are made available to users and how collaboration is facilitated among interested developers. In addition to distributing NJOY2016 and NJOY21 with an open-source license, these two new versions provide new capabilities for NJOY which will maintain NJOY as the standard for nuclear data processing into the foreseeable future. We report in this summary about the changes that come with these two new versions of NJOY and how this impacts the users and long-time contributors to NJOY.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Speed and Memory Improvements to MCNP6 Delayed-Gamma Line Treatment

The delayed-gamma (DG) line-emission (line) sampling was first introduced into MCNPX in 2005. At the time, storage demands for problems with large numbers of DG lines was a major concern. In order to accommodate large amounts of line data while keeping the memory footprint of such a problem from exceeding memory limits, complicated logic was included to sift through various storage options. In mid-2015 concerns arose over the adequacy of the MCNP6 DG source option (inherited from MCNPX) with line sampling for modeling high-fidelity problems. One indication that there was a problem came when studying the MCNP6 produced DG spectra of U-238 and its progeny using high-fidelity spectroscopic tallies when compared to DG line amplitudes directly from ENDF/B-VII.1 emission data. Despite producing a correct integral number of gammas, the comparison to ENDF data showed a broadening in the energy lines and a reduction in the amplitude of more prominent lines. It was later found that the cause of this broadening was one of the memory-saving techniques used to store DG line data when problems required many DG lines. Specifically, the DG lines were integrated into fine energy bins for sampling. When a DG line was to be sampled, a random number would be used to choose which energy bin to sample based on emission probability. Once chosen, a second random number was used to perform nearest-line sampling in that bin. The integration of lines led to an oversampling of low-probability lines, especially when in close proximity to high-probability lines. This treatment also caused a similar issue with the time-dependent behavior of DG line emission which occurred when low-probability lines were sampled near high probability lines with the former decay constant being used. The effect on time-dependent gamma production was demonstrated in analytic benchmarks performed by Weldon et al.. To address these issues a new gamma line storage and sampling treatment was introduced in 2015. This new treatment, stores gamma emission data in a line-by-line CDF which is sampled directly, thereby resolving both the line-broadening issue, and the time-dependence issue as was demonstrated by Tutt el al.. In addition, the new DG sampling treatment was written to use dynamically dimensioned storage arrays as opposed to large fixed-dimensioned arrays, as in the former treatment, and also makes more efficient use of line data that is stored by reducing or eliminating repeated data. As a result, storage arrays are generated only for DG data relevant to the problem being executed. The result should be a reduction in both memory requirements and runtime. However, these reductions could not be realized until the logic for the old treatment could be removed. This paper provides an update on recent work which has focused on the removal of the old DG sampling logic and attempts to quantify the memory reduction and speed up of selected problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainties in secondary angular distributions and their impact on shielding, criticality and kinetics analysis - 25556

The availability of the covariance data improved over the last decades. Use of covariance matrices of different origin (such as JEFF, JENDL, ENDF, TENDL, SCALE, etc.) and types (files MF31, MF33, MF34, MF35) in sensitivity and uncertainty analysis of vast list of cases ranging from fission to fusion and from criticality, kinetics and shielding to adjustment applications was already studied in the past. The present work focuses on the status and the availability of the covariance of the angular distributions of scattered neutrons (secondary angular distributions - SAD covariances). SAD uncertainties are usually neglected in classical uncertainty propagation codes, even though the methods and codes for calculating the sensitivities are available for decades and excellent agreement was demonstrated between different deterministic and Monte Carlo codes. Of particular concern is still the availability, quality and completeness of covariance matrices. Indeed, in most evaluations the SAD covariances are limited to P1 term, if available at all. Furthermore, some differences were observed between different evaluations, e.g. in the case of iron. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗

Testing of ENDF/B-VIII.0 in the GNDS format with LLNL transport codes - 25606

ENDF/B-VIII.0 was released in the Evaluated Nuclear Data Format and the Generalized Nuclear Data Structure format (GNDS), a new international standard for storing nuclear data developed under the auspices of the Nuclear Energy Agency Working Party on Evaluation Cooperation Subgroup 38. GNDS relies on modern coding practices, and offers a structured, flexible and physics-based standard for storing and exchanging nuclear data to address limitations in legacy nuclear data formats such as ENDF-6. GNDS supports the storage of evaluated nuclear data as well as processed data suitable for both deterministic (i.e., multi-group) and Monte Carlo transport codes. Lawrence Livermore National Laboratory has implemented a complete suite of codes to test GNDS data. Simulations for the ENDF/B-VIII.0 cross-section library in GNDS format are presented. Deterministic and Monte Carlo radiation transport simulations with the Ardra and Mercury codes respectively, are validated against reactivities in fast criticality safety benchmark assemblies. Finally, Mercury simulations of {sup 238}U(n,f), {sup 237}Np(n,f) and {sup 239}Pu(n,f) reaction rates in Godiva, Jezebel and Flattop are compared to MCNP results and measured data. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Importance of 3-D S{sub N} depletion in non-proliferation using BSOLVE

We present a single Pressurized Water Reactor (PWR) 3-D fuel rod design for depletion analysis using BSOLVE, our newly developed Runge-Kutta-Fehlberg based depletion code. BSOLVE is coupled with the deterministic 3-D S{sub N} particle transport code, PENTRAN, applied here with a 4-neutron energy group comparison to Continuous Energy (C/E) SERPENT2 Monte Carlo results. Differences are expected, as PENTRAN+BSOLVE retains full (multi-group) energy information for reactions, nuclide specific fission contributions, and energy dependent fission yields, using the latest available ENDF-BVIII data, important to retain accurate burned fuel inventories; SERPENT2 collapses burnup reactions to a single energy value. For depletion times up to ∼ 700 days and typical PWR power densities, relative differences between multigroup 3-D S{sub N} with full energy data and Monte Carlo one group burnup for trans-uranium nuclide concentrations and fission products are up to ∼20%. System eigenvalues are consistent, but with differences early and late in the cycle attributed to multigroup vs. C/E Monte Carlo cross sections. This work highlights the importance of low variance transport driven burnup for non-proliferation concerns, since plutonium quality varies significantly along axial lengths, and is more challenging to converge using Monte Carlo; details of depletion steps with spatial/zone dependent plutonium quality are provided. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗

Analysis of NCERC Critical Experiments with ENDF/B Nuclear Data Libraries

Nuclear data (ND) libraries are the backbone of the nuclear industry, as they are the collections of tabulated probabilities that define sub-atomic particle interactions with matter. In the areas of criticality safety and experiments, no evaluated nuclear data files (ENDF) are more important than those containing neutron cross section data. It is upon these files, and accompanying radiation transport codes, that practitioners are enabled to safely design both subcritical and critical systems. Likewise, in a symbiotic fashion, it is the same critical assemblies which are used primarily to validate that the cross sections are correct. The ENDF/B library, the United States’ national library maintained by the National Nuclear Data Center (NNDC) at Brookhaven National Laboratory (BNL), is soon releasing a new version, ENDF/B-VIII.1. Prior to the official release, several beta versions of the library were prepared and tested in simulation suites. The work presented here are results comparing the newest ENDF/B beta release (ENDF/B-VIII.1b3) and ENDF/B-VIII.0 with recent experiments done at the National Criticality Experiments Research Center (NCERC) as well as correlated experiments from the Los Alamos Critical Experiments Facility (LACEF). These NCERC and LACEF experiments were performed in part to provide validation for various cross sections that were identified as insufficient in the ENDF/B-VIII.0. In particular, lead, copper, fluorine, and chlorine, as well as the major actinides, were targeted from the last decade of critical experiments.

97 MATHEMATICS AND COMPUTING↗

Thermal scattering law of C 2 H 4 n : Integrating experimental data with DFT calculations

Improvements in determination of the thermal scattering law of moderator materials (measuring, calculating and validating) are important for accurate prediction of neutron thermalization in nuclear systems. In this paper a methodology for producing thermal scattering libraries from the experimental data for polyethylene C 2 H 4 n is discussed. Double differential scattering cross section (DDSCS) experiments were performed at the Spallation Neutron Source of Oak Ridge National Laboratory (SNS ORNL). New scattering kernel evaluations, based on phonon spectrum for C 2 H 4 n , are created using the NJOY2016 code. Two different methods were used: direct and indirect geometry neutron scattering at ARCS and SEQUOIA, and VISION instruments, respectively, where the phonon spectrum was derived from the dynamical structure factor S(Q, ω ) obtained from the measured DDSCS. In order to compare and validate the newly created library, the experimental setup was simulated using MCNP6.1. Compared with the current ENDF/B-VII.1, the resulting RPI C 2 H 4 n libraries improved both double differential scattering and total scattering cross sections. A set of criticality benchmarks containing C 2 H 4 n from HEU-MET-THERM resulted in an overall improved calculation of K eff , although the libraries should be tested against benchmarks more sensitive to C 2 H 4 n . The DFT + oClimax method is used and is shown to be most comprehensive method for analysis of moderator materials. Finally, the importance of DFT + oClimax method lies in the fact that it can be validated against all data measured at VISION, ARCS and SEQUOIA, and experimental total scattering cross section measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗

Toward a better thermal scattering law of (C 5 O 2 H 8 ) n : Inelastic neutron scattering and oClimax + NJOY2016

With the advancements in technology (both experimental and computational) the determination of the“true” experimental phonon spectrum became more accessible. Here a methodology for producing thermal scattering libraries from the experimental data (namely the DFT+oClimax method) for lucite (C 5 O 2 H 8 ) n is discussed. Double differential scattering cross section (DDSCS) experiments were performed at the Spallation Neutron Source of Oak Ridge National Laboratory (SNS ORNL). Novel scattering kernel evaluations, based on the phonon spectrum for (C 5 O 2 H 8 ) n ,were created using oClimax and NJOY2016 codes. In order to compare and asses the performance of the newly created library, the experimental setup was simulated using MCNP6.1. Compared to the current ENDF/B-VIII.0, the resulting RPI (C 5 O 2 H 8 ) n library improved the calculation of both double differential scattering and total scattering cross sections. A set of criticality benchmarks containing (C 5 O 2 H 8 ) n from HEU-MET-THERM resulted in an overall improved calculation of K eƒƒ . The DFT + oClimax method is shown to be the most comprehensive method for analysis of moderator materials, due to the fact that it can be verified against all data measured at VISION, ARCS and SEQUOIA neutron spectrometers at SNS ORNL, and experimental total scattering cross section measurements. This technique also provides a new technique for calculating any phonon spectrum-related quantities such as scattering law kernel, specific heat capacity, thermal conductivity, etc. for any solid state material.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On a measurement approach to support evaluation of thermal scattering law data

Inelastic thermal neutron scattering in materials that act as neutron moderators, reflectors, and filters results in shaping the neutron spectrum at low energies. This phenomenon is described using differential scattering cross sections calculated from three components including the bound atom (i.e., nuclear) scattering cross section of the neutron, the ratio of the outgoing and incoming neutron energy, and the thermal scattering law (TSL), i.e., S(α,β), where α and β represent dimensionless momentum and energy exchange variables, respectively. To date, no TSL libraries are generated using measured data. However, valuable information may be derived from measurements and “targeted” experiments that can validate TSL data and the related inelastic scattering cross sections. As a demonstration, a suite of coordinated measurements and experiments is described that was designed and used to support the evaluation of the TSL for “nuclear” graphite. This experimental suite includes neutron powder diffraction (for structure analysis), positron annihilation (for nano porosity assessment), inelastic neutron scattering measurements using a chopper spectrometer, transmission experiments using neutrons with energy below the Bragg cutoff thereby accessing the total (inelastic) cross section, and a slowing-downtime experiment to observe the developing neutron spectrum in the material. We find that this experimental suite was key to understanding the difference in TSL between “nuclear” and “ideal” graphite and for the inclusion of “nuclear” graphite in the ENDF/B-VIII.0 nuclear data library release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of ZPPR-15 experimental data for VTR software validation

As part of the Versatile Test Reactor (VTR) program, extensive validation work is being performed to demonstrate the fidelity of the Argonne Reactor Code (ARC) software package presently used for the core design studies. It is anticipated that this software will be also used for initial operations although additional software may be introduced at a later time to supplement or replace some components. The validation work was focused on accumulating reactor models that are consistent with the VTR and thus can be used to validate the accuracy of predictions made by the software used to design and operate the VTR. The U.S. does not have any active experimental facilities that can neutronically mockup the VTR core design. However, in the past, the U.S. had extensive experimental facilities for fast reactor neutronics. In that regard, the ZPPR-15 experiments carried out at ANL during the 1980s work best because the fuel form is closest to the VTR design, there is a considerable amount of processed data available for use, and existing staff are already familiar with the experimental machine and measurements. Because of previous collaborative work, many of the ZPPR-15 loadings of interest have already been processed into MCNP and preliminary ARC models, and a full uncertainty quantification has already been completed for many of the ZPPR-15 measurements. This manuscript covers the recent analysis work on the ZPPR-15 experiments and compares it with the experimental results. With respect to the ARC software components, this manuscript focuses on validation of the MC 2 , DIF3D and GAMSOR using the ZPPR-15 experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of the Neutron Data Standards

With the need for improving existing nuclear data evaluations, (e.g., ENDF/B-VIII.0 and JEFF-3.3 releases) the first step was to evaluate the standards for use in such a library. This new standards evaluation made use of improved experimental data and some developments in the methodology of analysis and evaluation. In addition to the work on the traditional standards, this work produced the extension of some energy ranges and includes new reactions that are called reference cross sections. Since the effort extends beyond the traditional standards, it is called the neutron data standards evaluation. This international effort has produced new evaluations of the following cross section standards: the H(n,n), 6 Li(n,t), 10 B(n,α), 10 B(n,α 1 γ), nat C(n,n), Au(n,γ), 235 U(n,f) and 238 U(n,f). Also in the evaluation process the 238 U(n,γ) and 239 Pu(n,f) cross sections that are not standards were evaluated. Evaluations were also obtained for data that are not traditional standards: the Maxwellian spectrum averaged cross section for the Au(n,γ) cross section at 30 keV; reference cross sections for prompt γ-ray production in fast neutron-induced reactions; reference cross sections for very high energy fission cross sections; the 252 Cf spontaneous fission neutron spectrum and the 235 U prompt fission neutron spectrum induced by thermal incident neutrons; and the thermal neutron constants. The data and covariance matrices of the uncertainties were obtained directly from the evaluation procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dr. Cecil Lubitz (1925–2021)

Dr. Cecil Lubitz passed away on March 5, 2021 at his home in upstate New York. A graduate of the Class of 1945 from the Naval Academy, Cecil served in the US Navy until 1947 after which he earned a MS in electrical engineering (1949) and a PhD in Physics (1960) from the University of Michigan. Following graduate school, Cecil came to work at Knolls Atomic Power Laboratory (KAPL) on applied nuclear physics in the Naval Nuclear Propulsion Program (NNPP). Cecil worked at KAPL for 54 years and had an impactful career with many substantial contributions to the NNPP, national and global nuclear data communities.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗