Generation of high-resolution thermal scattering laws for solid moderators using fast Fourier transforms
Not provided.
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.
Not provided.
The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.
Accurate thermal scattering laws (TSLs) are essential for reliable neutron transport simulations, particularly for systems with low-energy neutrons. Pulsed-neutron die-away (PNDA) experiments provide a highly sensitive platform for TSL validation but face challenges from fast neutron leakage and material constraints in nonhydrogenous targets such as beryllium. This work investigates the use of high-density polyethylene (HDPE) as a premoderator to reduce leakage, improve thermalization, and shorten experimental run times while preserving sensitivity to the Be TSL. Here, three premoderated configurations, encapsulation, slab, and interstitial, were evaluated using MCNP6.3® simulations with ENDF/B-VIII.0 nuclear data. The effects of the Be TSL was quantified by comparing decay constants from simulations with and without 𝑆(𝛼,𝛽) treatments for both beryllium and HDPE. The results show that premoderation enables the use of as little as 5.20% of the beryllium volume required for an unreflected geometry, achieving uncertainties below 0.5% and reducing run time by up to 72.2% in the studied configurations. The interstitial configuration achieved the best balance between low statistical uncertainty and high TSL sensitivity, outperforming the slab and encapsulation designs. These findings demonstrate that carefully optimized premoderation can significantly enhance the efficiency, feasibility, and precision of PNDA experiments for Be TSL validation.
Not provided.
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.
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.
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).
A collaborative effort between Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) is underway to provide a technical basis and methodology for the criticality safety community to use the sum of fractions (SoF) method for generating limits for mixtures of “selected actinide nuclides” included in the ANSI/ANS 8.15 standard. The PNNL scope in this project is to define a range of mixtures of 233 U, 235 U, and 239 Pu moderated with either light water or polyethylene and to examine the critical masses for these mixtures. The ORNL scope is primarily to provide validation support for these studies. More complete discussion of the project and its validation aspects will be presented at the upcoming International Conference on Nuclear Criticality Safety (ICNC) this Fall in Sendai, Japan. Clear differences in benchmark similarity to application systems as assessed by the integral parameter ck are noted in the validation studies performed as part of this project as a function of moderator. The c k value is a correlation coefficient that represents that amount of shared uncertainty in k eff due to cross sections between two systems. Individual nuclide-reaction contributions between the two systems can be simply summed to arrive at the total c k value. Specifically, the c k values for light-water–moderated solution experiments are higher for a water-moderated application than for a polyethylene-moderated application. This result is neither totally unexpected nor surprising, but the magnitude of the difference was difficult to anticipate. The TSUNAMI sequence, in the SCALE 6.2.4 code package developed by ORNL, was used to generate eigenvalues and reactivity effects with perturbation-theory based approach through sensitivity coefficients for all nuclides in the system with all reactions and energy groups. The TSUNAMI-Indices and Parameters (IP) sequence then uses the sensitivity data generated through TSUNAMI to generate relational parameters (i.e., c k ) to determine the degree of similarity between systems. One detail of the SCALE material and data implementation must be discussed at this point. Several thermal scattering laws (TSLs) are available for 1 H. SCALE uses a different nuclide ID number for each TSL; essentially, each version of 1 H is treated as a unique nuclide. For example, 1 H bound in water ( 1 H-H2O) is assigned the nuclide ID 1001, whereas 1 H bound in polyethylene (h-poly) is assigned the nuclide ID 9001001. The same cross section data are used for all reactions in 1 H, regardless of TSL, except for scattering below the TSL cutoff energy. TSUNAMI-IP treats different nuclide IDs as different nuclides; thus, no uncertainty is shared between 1 H-H2O and h-poly, despite much of the same data, including covariance data, being used for both nuclides. This presents a question: how much of the difference in assessed similarity between water- and polyethylene-moderated systems is due to the moderators, and how much is caused by the treatment of 1 H-H2O and h poly with cross section and covariance data. The extended edits generated by TSUANMI-IP allow for an investigation of this issue specifically, as well as a demonstration of the general techniques available within TSUNAMI to understand the results of the similarity assessment. This paper presents and analyzes the similarity assessment of both water- and polyethylene-moderated systems for a single benchmark: PST-002-001.
This report summarizes the experiments performed for IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments were conducted in two separate two-day campaigns: one campaign for high-density polyethylene and another for Lucite. All experiments were performed at LLNL. They will become high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high-priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). The report presents the PNDA design and its equipment, and it documents the experimental die-away curves. It presents the characterization measurements that have been performed to-date. Importantly, mass spectrometry measurements to determine sample impurities were not yet completed. These will be included in the final benchmark. The report gives the fitted decay constants for the die-away curves of each target sample. It also provides the data for the die-away curves in the appendices. The HDPE experiments examined twelve targets of varying size. The Lucite measurements included ten targets of varying size.
This report summarizes the experiments done for IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments were conducted in two separate two-day campaigns: one campaign for high-density polyethylene and another for Lucite. All experiments were performed at Lawrence Livermore Laboratory. They will become high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). The report presents the PNDA design and its equipment, and it documents the experimental die-away curves. It presents the characterization measurements that have been performed to-date. Importantly, mass spectrometry measurements to determine sample impurities were not yet completed. These will be included in the final benchmark. The report gives the fitted decay constants for the die-away curves of each target sample. It also provides the data for the die-away curves in the appendices. The HDPE experiments examined twelve targets of varying size. The Lucite measurements included ten targets of varying size. Experimenters in the campaign were Daniel Siefman, William Zywiec, and Ruby Araj.
This presentation finds that PMT-004 has four new benchmark cases highly sensitive to PE TSL (2 cases) and PMMA TSL (2 Cases). PE cases were well predicted using MCNP6.2 and ENDF/B-VIII.0. PMMA cases overpredicted by approximately 0.6-0.7% in k eff at 20°C. Accepted into 2023 version of the ICSBEP Handbook. Temperature had a large impact on reactivity of the critical configuration. Implications for validation work for thermal cases- need to adjust TSL data to correct temperature as it can have hundreds of pcm effects for a few °C. Future thermal experiments should try and measure reactivity at multiple temperatures to aid in data testing and benchmark adjustment.
This presentation finds that PMT-004 has four new benchmark cases highly sensitive to PE TSL (2 cases) and PMMA TSL (2 Cases). Additionally discovered was that the PE cases were well predicted using MCNP6.2 and ENDF/B-VIII.0. PMMA cases overpredicted by approximately 0.6-0.7% in k eff likely due to the lack of integral benchmarks sensitive to PMMA. The NCSU ENDF/B-VIII.1 PMMA TSL is an improvement over VIII.0. This was Accepted at last ICSBEP meeting pending comment resolution, should be published in 2024 version of the Handbook. It was found that temperature had a large impact on reactivity of the critical configuration. This has Implications for validation work for thermal cases- need to adjust TSL data to correct temperature as it can have hundreds of pcm effects for a few °C. Future thermal experiments should try and measure reactivity at multiple temperatures to add in data testing and benchmark adjustment.
This presentation found validation efforts revealed noticeable differences in c k for light-water and polyethylene moderated systems when compared to PST-002-001. Additionally, although 1 H-H 2 O and 1 H-poly use the same cross section and covariance data, the treatment of TSLs in SCALE lead to no contribution to ck between polyethylene application and benchmark experiment. When H 2 O contribution were removed the difference in c k appears to come from 239 Pu chi. Sensitivity profiles and data-induced uncertainty confirmed that the polyethylene application was more closely resembling the benchmark experiment. Specific differences between similar systems can be accounted for and examined to help understand differences in similarity (c k ) between systems was a conclusion reached by this presentation.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.