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 deals with AMPX development in FY2023 focused on preparing for the release and testing of ENDF/B-VIII.1. It includes updates for Thermal Scattering Law in ENDF/B-VIII.1. The presentation finishes with touching on resolution of a SCALE Issue with thermal kernels and photonuclear physics support.
Recently, a cold moderator was designed and developed for use at the Rensselaer Polytechnic Institute Linear Accelerator (RPI LINAC). This cold moderator proved to easily and safely couple to an existing neutron producing target, while enhancing neutron flux below 0.02 eV by up to a factor of 8 by cooling polyethylene down to 29 K. This cold moderator capability allowed for significantly improved counting statistics below 0.02 eV not previously possible due to a poor signal to background ratio. Additionally, testing was performed to characterize the energy resolution of the new cold moderator system and found the system easily capable of resolving resonances in Ta-181 at 4 and 10 eV, while also clearly resolving the Bragg edges found in Be metal below 0.01 eV. Furthermore, following the design and development of a cold polyethylene moderator, a series of thermal total cross section measurements were performed for polyethylene, polystyrene, Plexiglas and yttrium hydride in the thermal region. These measurements serve to help validate thermal scattering law (TSL) evaluations in the 0.0005 – 1 eV energy range. For polyethylene and polystyrene, two sets of experiments were performed – one with the Enhanced Thermal Target (ETT) and another with the ETT plus the new cold moderator capability (ETTC). The yttrium hydride and Plexiglas measurements were only performed with the ETTC. The measurements for polyethylene help to validate the data processing methodology when using the ETTC, while extending the measured range of polyethylene down to 0.0005 eV. Two different Plexiglas, Plexiglas G and Plexiglas G-UVT, and two different concentrations of yttrium hydride, H/Y = 1.85 and 1.68, were measured. Overall, all materials had generally good agreement with their ENDF/B-VIII.0 TSL evaluations, though some discrepancies were noticed. In the case of the yttrium hydride, the high energy oscillations in the hydrogen cross section and the low energy Bragg edges in the yttrium cross section were clearly seen. These measurements represent the first total cross section measurements that encompass the entire thermal region from 0.0005 – 1 eV for polystyrene and yttrium hydride.
Recently, a cold moderator was designed and developed for use at the RPI LINAC. This cold moderator proved to easily and safely couple to an existing neutron producing target, while enhancing neutron flux below 0.02 eV by cooling polyethylene down to 29 K. This cold moderator capability allowed for significantly improved counting statistics below 0.02 eV not previously possible due to a poor signal to background ratio. Additionally, testing was performed to characterize the energy resolution of the new cold moderator system and found the system easily capable of resolving resonances in Ta-181 at 4 and 10 eV, while also clearly resolving the Bragg edges found in Be metal below 0.01 eV. Following the design and development of a cold polyethylene moderator, a series of thermal total cross section measurements were performed for polyethylene, polystyrene, Plexiglas and yttrium hydride in the thermal region. These measurements serve to help validate thermal scattering law (TSL) evaluations in the 0.0006 – 20 eV energy range. For polyethylene and polystyrene, two sets of experiments were performed – one with the Enhanced Thermal Target (ETT) and another with the ETT plus the new cold moderator capability (ETTC). The yttrium hydride and Plexiglas measurements were only performed with the ETTC. The measurements for polyethylene help to validate the data processing methodology when using the ETTC, while extending the measured range of polyethylene down to 0.0007 eV. Two different Plexiglas, Plexiglas G and Plexiglas G-UVT, and two different concentrations of yttrium hydride, H/Y = 1.85 and 1.68, were measured. Overall, all materials had generally good agreement with their ENDF/B-VIII.0 TSL evaluations, though some discrepancies were noticed. In the case of the yttrium hydride, the high energy oscillations in the hydrogen cross section and the low energy Bragg edges in the yttrium cross section were clearly seen. These measurements represent the first total cross section measurements that encompass the entire thermal region from 0.0006 – 3 eV for polystyrene and yttrium hydride.
In order to test the performance of new neutron thermal scattering law (TSL) evaluations it is desirable to have experimental data that is highly sensitive to the TSL and provides high fidelity information on the energy dependent performance of TSL evaluations. Three relevant experiments are discussed including: accurate thermal total cross section measurements, thermal neutron die-away experiments, and neutron leakage experiments. The experimental setups and results are reviewed, and examples provided for some moderators including polyethylene, Plexiglas, and YHx. For the experiments preformed thus far, there is generally good agreement between the measured total cross section and simulations using current TSL evaluations, however in certain energy ranges differences were observed. Similarly, neutron die-away and leakage measurements for samples at room temperature are in good agreement with data computed from TSLs, however leakage measurements for polyethylene at 29K show discrepancies with TSL evaluations.
There is a strong need for new benchmarks to validate neutron thermal scattering laws (TSLs). Lawrence Livermore National Laboratory (LLNL) has designed a Pulsed-Neutron Die Away (PNDA) testbed for this purpose. The experiment has a deuterium-tritium (D-T) neutron generator that impinges a 10 -4 s, mono-energetic pulse of 14.1 MeV neutrons on a target sample. After the pulse, the neutron population moderates and establishes a thermal equilibrium within the sample, with a fundamental spatial mode and characteristic decay-time eigenvalue, ⍺. The ⍺ eigenvalue can be extracted from the experimental measurements of the time-dependent neutron flux coming off the surface of the sample and can then be used as an integral parameter (similar to k eff ) to validate nuclear data involved with neutron migration, thermalization, and absorption. For moderating materials and geometric configurations, the ⍺ eigenvalue is heavily dependent on thermal neutron scattering of the target material. For that reason, a PNDA experiment can have a higher sensitivity to TSLs than is commonly available with the k eff parameter in critical experiments.
Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.
This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.
The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)
To address the mounting need for below room temperature nuclear data validation, the Low-Temperature Thermal Epithermal eXperiments (LT-TEX) have been designed. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University.
Anhydrous Hydrogen Fluoride (HF) at high temperatures and pressures is used to process and manufacture nuclear fuel. As HF is often used directly with uranium, correct neutron thermal scattering cross sections are crucial to criticality safety applications. Classical molecular dynamics (CMD) simulation of the flexible HF system was used to create the thermal scattering law (TSL) and cross sections. The initial 2-site model is used in LAMMPS, and it can not capture the H-bond. To correctly represent the H-bond effects, a second, 3-site model was constructed in GROMACS. The 3-site model handled H-bonds by connecting a massless charge to the molecule. Key model parameters were compared to experimental data to verify the approach and models. To get the normalized VACF, the model was compared using hydrogen and fluorine bond length, density, potential energy, and diffusion coefficient. The phonon DOSs for both models were derived from the normalized VACF. DOSs were used to estimate the TSL ( S ( α, β )) and neutron thermal scattering cross sections for hydrogen in HF. The TSLs were evaluated using the FLASSH code with the Schofield diffusion model. It was observed that the representation of the hydrogen bonding changes the TSL's diffusional contributions. This is represented in the low energy scattering cross section, where intermolecular binding effects shift the cross section.
This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.
The Sandia Nuclear Criticality Safety (NCS) program’s benchmark suite was recently updated. This suite is used to ensure that NCS calculations using computer-based neutron transportation codes have an established baseline comparison of calculated versus known experimental results. The Evaluated Nuclear Data File (ENDF) version used for the MCNP models in the suite was changed from ENDF/B-VII.1 to ENDF/B-VIII.0. Additionally, relevant thermal scattering law data libraries (TSLs) were updated. The sensitivity of the calculational bias of each benchmark model to these changes is discussed. Implementation of the ENDF/B-VIII.0 library and updated TSLs results in improvements to bias distribution in the intermediate enriched uranium, plutonium, and mixed uranium–plutonium (IEU, PU, and MIX) fissionable material benchmark categories, but a small bias increase in low- and high-enriched uranium categories (LEU and HEU, respectively). The results also highlight the sensitivity of the benchmarks, with average lethargy of neutrons causing fission energies (EALF) in the intermediate energy range to ENDF/B library changes. The most numerous bias changes were observed in the thermal energy region when transitioning from ENDF/B-VII.1 to ENDF/B-VIII.0. In conclusion, most of the unique bias changes observed in MCNP 6.3.0 between the two nuclear data libraries were in the LEU-COMP-THERM evaluation subset.
The goal of IER-479 is to design uranium critical experiments that can be used to validate low temperature cross sections and criticality safety analyses over multiple neutron energy regimes. Currently, there are no benchmarks in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook at temperatures lower than room temperature (International Criticality Safety Benchmark Evaluation Project Handbook, 2019). However, there are many needs for validation of criticality safety analysis at lower temperatures, including meeting transportation requirements and operations conducted outside or in unheated facilities. Additionally, NCSP has funded North Carolina State (NCSU) to generate new thermal scattering laws, including at lower temperatures, and the lack of integral benchmarks impedes data testing of these new cross sections. To address these needs, this report will present a critical experiment design covering various fission energy regimes with a goal temperature of -40°C (-40°F), which is based on the lower bound of expected non-cryogenic operational temperatures. The goal of the U.S. Nuclear Criticality Safety Program’s (NCSP) Thermal/Epithermal eXperiments (TEX) is to design and conduct new critical experiments to address high priority nuclear data needs from the nuclear criticality safety and nuclear data communities. The TEX program includes two series of baseline experimental configurations, one based on plutonium fuel (plutonium-aluminum Zero Power Physics Reactor (ZPPR) plates) and the other based on uranium fuel (highly enriched uranium (HEU) plates), that are moderated with varying thickness of polyethylene to create assemblies which span the thermal, intermediate, and fast fission energy regimes. The configurations are designed to be easily modified (for example, to add diluent materials of interest) to allow for efficient generation of additional benchmark configurations and allow for added nuclear data testing utility when comparing modified configurations to baseline configurations. The goal of IER-479 is to use the TEX-HEU concept (stack of HEU plates and polyethylene moderators) to design a critical experiment that can be used to validate low temperature cross sections and criticality safety analyses.
For advanced reactor applications, Neural Thermal Scattering (NeTS) modules were developed to predict the thermal scattering law (TSL or $S(α, β, T)$) of a nuclear graphite neutron moderator. NeTS are multi-layer, feedforward artificial neural networks, which act as universal function approximators designed for TSL datasets. In this case, a 4-layer neural network with 164 neurons per layer is trained using FLASSH evaluated data in PyTorch and serialized as a torchscript dictionary to predict $S(α, β, T)$ on-the-fly. Relative, absolute and maximum percent deviations of NeTS from File 7 data generated using the FLASSH code are on the order of 0.01%, 0.1% and 1%, respectively, with low inference latencies of 0.000172 s per $S(α, β, T)$ at a given temperature. Capturing the full dimensionality of possible inelastic neutron-lattice interactions, NeTS functionality is embedded in the Serpent Monte Carlo code, where $S(α, β, T)_{NeTS}$ sampling is conducted on-the-fly and compared to ACE look-up-tables for predicting TREAT criticality. k-eff differences between sampling algorithms of 6 pcm are observed and are within the order of Monte Carlo uncertainty. Compared to discrete and continuous-energy ACE files (30 MB and 131 MB per temperature), the NeTS format is on the order of 200–300 kB for a continuous-temperature, interpolation-free representation of $S(α, β, T)$ and cross sections. NeTS-in-Serpent runtimes comparable with ACE look-up tables are achieved by scaling NeTS for high performance computing architectures with hybrid OpenMP + MPI parallelization. This work validates a novel, self-contained reactor physics framework for predictive cross sections, and demonstrates a general methodology for embedding modern machine learning libraries within existing neutronic analysis frameworks.
This presentation provides overviews of ENDF and GNDS reading in AMPX, AMPX covariance processing updates, thermal scattering law updates, and data libraries for SCALE 6.3.0.
On March 1 st , 2023, the National Nuclear Data Center (NNDC) released the ENDF/B-VIII.1β1 nuclear data library, which contained neutron, spontaneous and neutron-induced fission yield, and alpha sublibraries. The library was released in the standard Evaluated Nuclear Data File (ENDF) format. The files can be accessed on Brookhaven National Laboratory’s GitLab. The thermal scattering law (TSL) sublibrary will be released at a later date. The files provided in the neutron sublibrary were processed into A Compact ENDF (ACE)-formatted files, verified, and validated by the XCP-5 Nuclear Data Team. The resulting collection of ACE files will be referred to as the “e81b1” library. This report details the processing of these files and the quality assurance approach taken. This is not intended to be a full validation effort; rather, this library is intended to simply reproduce the released files for further validation testing by the community. The e81b1 ACE library will be released to the NNDC for access through GitLab and hosted locally on LANL computing resources. This library is only intended for analysis to support the future ENDF/B-VIII.1 release and should not be used for production results.
A new cold neutron moderator system was developed and constructed at the Rensselaer Polytechnic Institute (RPI) Gaerttner Linear Accelerator. This system was then used to measure the neutron total cross section in the thermal energy range for polyethylene, polystyrene, Lucite, and yttrium hydride. In tandem with these measurements, a fitting procedure was developed at Oak Ridge National Laboratory to vary computationally simulated phonon density of states until a fit with differential and integral data was achieved. This model, combined with the new RPI total cross section data, was used to generate thermal scattering law files for the materials measured at RPI.
The low-temperature variant of the TEX HEU (called Low-Temperature TEX or sometimes LT TEX) campaign is a highly anticipated and necessary experimental series by the greater nuclear science community. Fundamentally, the need for low-temperature integral experiments is required to perform validation of cross sections below room temperature. There has been substantial international interest in low-temperature benchmarks to validate below room temperature cross sections, namely talks given at the 2019 International Conference on Nuclear Criticality (ICNC): UK (Watson, 2019), France (Milin, 2019), and UK (Gan & Wilson, 2019). Additionally, NCSP funded thermal scattering laws (TSLs) were produced by North Carolina State University and require low-temperature benchmarks to validate them. Validation of low-temperature cross sections is also necessary for criticality safety applications. One particularly important application is to ensure that during transportation, fissile materials must remain subcritical under normal ambient conditions which is defined as temperatures down to -40°C/°F by the United States 10 CRF 71 as well as a regulation put forward by the International Atomic Energy Agency (IAEA).