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.
Banner Highlight for Structure-Dependent Doppler Broadening Using a Generalized Thermal Scattering Law [DOI 10.3390/jne2020013].
Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).
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.
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)
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.
Scattering of thermal neutrons and Doppler broadening of epithermal neutron resonances in uranium and its compounds may be sensitive to crystal binding. The thermal scattering law (TSL) for uranium dioxide, which captures crystal binding effects, has been reevaluated for ENDF/B-VIII.0. Phonon spectra were generated using ab initio lattice dynamics for the paramagnetic phase and validated against experiment. Improved agreement with the Debye-Waller coefficient as a function of temperature is found relative to the spectrum used for the ENDF/B-VII.1 evaluation. The TSL was generated using the phonon expansion method within the NJOY nuclear data processing package and was found to be in reasonable agreement with inelastic neutron scattering measurements. The present evaluation predicts a reduction in the inelastic scattering cross section relative to ENDF/B-VII.1 and a total scattering cross section consistent with neutron transmission experiments
Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.
The increasing demand for green, low-carbon energy solutions has amplified the focus on advanced nuclear reactor technologies. Among these, Molten Salt Reactors (MSRs) have been spotlighted because of their special features. Historically, Oak Ridge National Laboratory (ORNL) started the Molten Salt Reactor Experiment (MSRE) in 1956. Later, ORNL and the University of California, Berkeley (UCB) jointly developed the MSRE benchmark. This was then reviewed by the International Reactor Physics Experiment Evaluation Project (IRPhEP) committee and added to their 2019 handbook. In the current study, the benchmark model was recreated via the Serpent code. The calculated effective multiplication factor was 1.02087 ± 0.00019, which gives a deviation of about 2000 pcm compared to the benchmark/experimental result. At North Carolina State University (NCSU), recent evaluations were conducted on thermal scattering cross sections for molten salt FLiBe and 20% porous nuclear graphite, anticipated for incorporation in the ENDF/B-VIII.1 library. The influences of these thermal scattering law (TSL) data and geometry changes on the multiplication factor were examined using the ENDF/B-VII.1 and ENDF/B-VIII.0 libraries. Introducing the TSL data for FLiBe and 20% porous graphite resulted in an effective multiplication factor of +270 pcm compared to the reference value obtained with ENDF/B-VII.1 in this study. The adjustment of the reactor vessel dimensions resulted in decrease of -448 pcm in the calculated effective multiplication factor.
The Full Law Analysis Scattering System Hub (FLASSH) is a modern, advanced code which evaluates the thermal scattering law (TSL) along with accompanying cross sections. FLASSH features generalized methods which accommodate any material structure. Historical approximations including the incoherent and cubic approximations have been removed. Instead, the latest release of FLASSH features advanced physics options including distinct corrections (1-phonon contributions) and non-cubic formulations. The non-cubic elastic and inelastic contributions are necessary to accurately evaluate 1-phonon contributions. Both non-cubic and 1-phonon calculations require high-density sampling of the various scattering directions. Optimization and parallelization of these routines were therefore necessary to produce results in a reasonable timeframe. With these notable improvements to the generalized TSL, FLASSH 1.0 meets benchmark requirements, demonstrating noticeable agreement with experiment for both TSLs and the resulting integrated cross sections. Additional features including a graphical user interface (GUI), plotting diagnostics, and formatted output options including ACE files allow users to complete a TSL evaluation with minimal input and maximum flexibility. The user GUI creates input files for FLASSH, reducing user error and also providing built-in error checks. Autofill options and suggested input values help make TSL evaluation accessible to novice users. The FLASSH code is compiled to run on both Windows and Linux platforms with automatic parallelization. (authors)
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.
Precise estimation of neutron thermalization in moderators relies on high-fidelity thermal scattering cross-section data governed by the thermal scattering law (TSL). The Institut de Radioprotection et de Sûreté Nucléaire (IRSN) has been working on the development of improved TSL for light water ice. Many polymorphic phases of light water ice exist depending on the thermodynamic conditions. The most common type of ice at standard pressure and temperature below water freezing point (273.15 K) is ice-Ih. It is essential to have high-resolution experimental double differential data for developing and/or validating TSL for moderator materials.Existing experimental double-differential scattering data for ice-Ih are extremely sparse and of limited quality. New high-quality double-differential measurements for ice-Ih over multiple temperatures and incident neutron energies would directly support the validation and improvement of ice-Ih TSL models for criticality safety applications.Series of time-of-flight (TOF) inelastic neutron scattering experiments on ice-Ih at temperatures starting at 271 K and down to 6 K, have been carried out at the SEQUOIA spectrometer at the Spallation Neutron Source (SNS) at the Oak Ridge National Laboratory (ORNL), United States. The experiments have been performed for incident neutron energies, Ei= 11, 55, 160, 250, and 600 meV, to explore different excitation energies in the vibrational phonon spectrum. This paper presents the thermodynamic conditions and the details of the TOF measurements on ice-Ih and the derived phonon spectrum from the experimentally measured double differential data. A study of the variation of the phonon spectrum of ice-Ih as a function of temperature is highlighted, and a preliminary TSL evaluation is developed based on the experimental phonon spectrum.
Nuclear technologies, such as thermal reactors and critical systems, are driven by reaction rates which are directly dependent on thermal neutrons. The ability to calculate reaction rates is directly linked to the capability to predict and understand the cross sections. For atoms bound within a material structure, the available modes of energy exchange are captured in the Thermal Scattering Law (TSL, i.e. S(α,β)). The TSL is fundamentally a material property and can be defined explicitly to account for the excitation and deexcitation of the lattice energy states. This work describes the development of the Full Law Analysis Scattering System Hub (FLASSH) which calculates the TSL and cross sections. FLASSH is a modern code which has been designed to produce accurate TSL evaluations, remove historical approximations, and improve formulations while making advance physics accessible and simple to use.
Zirconium Carbide (ZrC) is a refractory material currently being investigated for high temperature reactor applications. ZrC was studied as a base of (U,Zr)C solid solution and UC 2 composite fuels for NERVA-type cores during the Rover program. Testing of the material was performed at Nuclear Furnace 1 (NF-1) located at Los Alamos. Thermal Scattering Law (TSL) evaluations generated with FLASSH are a new contribution to ENDF/B-VIII.1.
The presentation discusses the eight new Thermal Scattering Law (TSL) evaluations that have been submitted to the National Nuclear Data Center (NNDC). It explains that evaluations are performed using the NCSU FLASSH code and that NJOY/LEAPR input files are provided to allow for independent testing. The presentation concludes by stating that more evaluations are expected to be submitted before the end of 2020.
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.