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.
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.
This report documents an evaluation of 238 U prompt fission neutron spectra (PFNS) which is a deliverable for a FY2024 Q4 NCSP (Nuclear Criticality Safety Program) milestone. This evaluation is new; its prior input is based on extended Los Alamos and exciton models implemented in the code CoH. Experimental covariances were estimated for five experimental data sets. One of these data sets that was measured by the Chi-Nu team of LANL and LLNL. It covers the 238 U PFNS for continuous incident-neutron energies of 1–20 MeV and outgoing-neutron energies from 10 keV– 10 MeV with high precision. Contrary to Chi-Nu data, previous data sets were measured in a limited energy range. The resulting evaluated data correspond well to the experimental PFNS taken into account for the evaluation. The evaluated PFNS also produce average mean energies in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of the Flattop, Flattop-Pu and BigTen ICSBEP critical assemblies (which all have thick reflectors with high percentages of 238 U), the differences of simulated values compared to those using ENDF/B-VIII.1β3 is modest (less than 25 pcm). In addition to that, the new PFNS predict on average 238 U LLNL pulsed-sphere neutron-leakage spectra slightly better than ENDF/BVIIII.0 and ENDF/B-VIII.1β3 PFNS. The differences are, however, well within the experimental uncertainties.
This report documents an evaluation of 240 Pu prompt fission neutron spectra (PFNS) including Chi-Nu experimental data. This work is a FY2025 Q4 NCSP (Nuclear Criticality Safety Program) milestone and named a “strategic priority”. This evaluation is new; it differs substantially from the most recent evaluation included in the ENDF/B library, which was for ENDF/B-VII.1 and then carried over unchanged to ENDF/B-VIII.1. The key difference between ENDF/B-VIII.1 and the new evaluation is that we have, for the first time, realistic experimental 240 Pu PFNS covering a broad incident and outgoing neutron energy range.
63,65 Cu angular distributions resolved for ENDF/B-VIII.1, consistent with both differential data and integral benchmark performance. 139 La RRR and URR evaluations will be merged with the LANL high energy evaluation, and to be submitted to the future ENDF-B release (post-VIII.1). 14 N RRR evaluation is planned to produce n+ 14 N, p+ 14 C, and a+ 11 B contributions to be submitted to the future ENDF-B release (post-VIII.1). 181 Ta covariances repaired and reported as intended in an errata to ENDF/B-VIII.1.
The ENDF/B-VIII.1 evaluation library has seen a great growth in the thermal neutron scattering sub-library. The SCALE code system has traditionally approached CE transport by assuming that the CE library on disk represented the fully expanded cumulative probability distributions, conditional on exiting angle and marginal on exiting energy. While this is a complete description of the data, it comes at the potential cost of large amounts of on-disk storage. This approach was strained by several TSL files in ENDF/B-VIII.1, such as graphite, which contained data for a large number of Bragg edges. In the fully expanded probability distributions, this was found to be a disproportionately large fraction of the SCALE CE library.
My Uncle Willie circa 1600 wrote “What’s in a name; a rose by any other name would smell as sweet.” I fear in this case we have a somewhat similar problem in that we may be using the same word but are not using the same definition; specifically, the word Unresolved. The simplest physics definition as it applies to neutron resonances, is the energy point where we can no longer see/measure ALL – let me repeat that – ALL - of the individual resonances. That seems simple and clear, but the question is: how to represent resonances beyond this point in order to accurately reproduce the effects we have seen in measurements and expect/need to reproduce in our applications. We know there are more, unseen resonances, otherwise we wouldn’t say Unresolved. The ENDF approach is well defined in ENDF-102 and simple: for ENDF data the only way to represent Unresolved data is by using a theoretical model to define the distribution of resonances, including those that are too narrow to measure (i.e., are unresolved). It is important to note that in ENDF this is the one and only Unresolved model, e.g., there is no provision in ENDF to accurately define individually ALL resonances above the Resolved energy range – by ALL here I mean both those that we can measure and those that we cannot individually measure, but that theory and integral measurements tells us are present. An alternative approach, which would appear to be equally valid, would be to include the latest measured data as tabulated energy expendent data extending upwards in energy above the Resolved energy range. In this approach the evaluation would not include an ENDF style Unresolved energy range; it would only include a Resolved resonance region, followed by tabulated higher energy points, representing the resonances that could be measured beyond the Resolved range. But an important point to note: By listing these resonances above the resolved energy one admits that at least some resonances in this energy range are missing as Unresolved; i.e., they are too narrow or overlapping to measure. The purpose of this paper is to illustrate that the later approach, while done with good intentions, and appearing to be valid/adequate in plots, does not meet the need of our engineering applications. Why? As we will see below, of these two possible approaches, only the ENDF use of a model to statistically include the missing, i.e., unresolved, resonances, can meet our engineering needs to reproduce the integral effects we have measured and understand. Only with this statistical model can we predict and include in our calculated results the important effects of temperature (Doppler broadening), and energy integrals (self-shielding). Below I will first present results using two ENDF/B-VIII.1 evaluations, U235 and U238, that use the correct ENDF-102 definition of an Unresolved resonance region, using a statistical model to include the effects of resonances that theory predicts are present, but are too narrow to measure. These two evaluations reproduce the expected temperature (Doppler) and energy integral (self-shielding) effects that we expect. Next I will present results using one ENDF/B-VIII.1 evaluation, 26-Fe-56, that does not use an ENDF-102 Unresolved resonance region; instead above its Resolved energy range it lists many tabulated energy points, that look like measured data, but by definition, since they are included above the ENDF Resolved energy range there are missing Unresolved resonances, i.e., there are missing the resonances that are too narrow to resolve, i.e., are unresolved. My conclusion, and I hope yours, is that the below figures illustrate that this approach does not reproduce the temperature and energy integrals that we expect and need to accurately calculate results for our fission reactor calculations. As such this approach should not be used in ENDF formatted evaluations.
12 new/revised TSL evaluations were contributed to ENDF/B-VIII.1. The thermal diffusion length $L$ is an integral property of a single material’s absorption and scattering cross sections (both differential and integral). No other neutron reactions or materials are involved. The MC21-calculated $L$ for water is consistent with the spread of experimental data and is sufficiently sensitive to different H-H 2 O TSL physics models to use the method as a TSL integral performance benchmark. Modern high-quality diffusion experiments at elevated $T$ would allow direct low-cost physics benchmarking of water TSLs when public elevated-$T$ critical benchmarks are limited. NNL is working with LLNL and RPI to develop experimental PNDA capability. An ICSBEP Fundamental Physics benchmark of the Nassar and Murphy PNDA experiment is being developed to provide an example of how to evaluate these type of benchmarks.
In this presentation 6 Li Distinct changes are observed. 235 U Discussion shows Improvements (seen in VIII.1beta0) are coming from new PFNS informed by Chi-Nu experimental data. Additionally, 239 Pu large changes coming from inelastic scattering. 9 Be ENDF/B-VIII.0 shows to be better right after elastic peak while 16 O had slight changes. 19 F had small improvements in ENDF/B-VIII.1 and 56 Fe had noticeable changes in thicker spheres. Examples of negligible changes for the sake of completeness are given.
This report is in answer to the Nuclear Criticality Safety Program FY24 quarter 4 milestone that requires: “Format and V&V nu-bar means and covariances” for 238,240-242 Pu average prompt fission neutron multiplicities, $\overline{ν}$ p , that were obtained by a consistent evaluation leveraging the fission-event generator CGMF and a detailed uncertainty quantification of experimental data. It is described how nuclear data mean values and covariances were formatted using ENDFtk. Implementing the new 238,240-242 Pu $\overline{ν}$ p into the ENDF/B-VIII.1β 4 library leads to only small overall changes in criticality values of the Jezebel, Dirty Jezebel, Jupiter-001, Jupiter-002, EUCLID 3x2 and EUCLID 8x1 critical assemblies. Simulated k eff uncertainties due to $\overline{ν}$ p covariances change only little if cross-isotope covariances are considered or not for those assemblies with low percentage content of minor Pu isotopes. However, for the Dirty Jezebel critical assembly, that has a sizeable 240 Pu and non-negligible 241 Pu content, the simulated k eff uncertainties due to considering or neglecting cross-isotope $\overline{ν}$ p covariances is 443 versus 374 pcm.
The R -matrix analysis of neutron-induced reactions for two stable chlorine isotopes ( 35 Cl and 37 Cl) was performed in the energy range of thermal up to 1.2 MeV. Starting from the repository of the ENDF/B-VIII.0 library and following recent measurement series, this work represents a significant improvement, particularly in the evaluation of the ( n , p ) reaction channel. The evaluation methodology used the R -matrix code SAMMY to generate a set of Reich–Moore resonance parameters for both stable chlorine isotopes. Consistent with recently measured 35 Cl( n , p ) data, the presented evaluation features a dramatic increase in the magnitude of the ( n , p ) reaction channel over the ENDF/B-VIII.1 nuclear data library and previous ENDF/B libraries.
This presentation focus on AMPX development for 2024 it also touches on ENDF/B-VIII.1. The slides touches on Thermal Neutron Scattering Library, AMPX Development Highlights, and Transport libraries for SCALE 7.0.0.
This presentation finds that as it was designed, Cerberus is extremely sensitive to copper nuclear data. Just changing copper nuclear data to another library can swing simulated keff by >1000 pcm. ENDF/B-VIII.1 and JENDL-5.0 both show large downward trends in keff as copper interstitial thickness is increased
Westcott 𝑔 factors are used in Neutron Activation Analysis (NAA) and Prompt Gamma-ray Activation Analysis (PGAA) to evaluate the impact of non-1∕𝑣 behavior in the neutron-capture cross sections of certain nuclei on activation product yields. This non-1∕𝑣 behavior arises from the presence of neutron resonances in the neutron- capture cross sections that overlap with the source neutron spectrum at low (< 5 eV) energies. Historically, Westcott 𝑔 factors that have been cataloged for NAA and PGAA applications are the result of calculations that assume a Maxwellian neutron flux distribution with a given temperature. In this work, we use this approach with updated neutron-capture cross sections from the Evaluated Nuclear Data File, version VIII.1 (ENDF/B-VIII.1) to tabulate Westcott 𝑔 factor values for a broad range of Maxwellian distribution temperatures, comparing the results against currently-available 𝑔 factors from International Atomic Energy Agency tables and other sources. Here, it was discovered during this analysis that the use of guided thermal and cold-neutron beams at certain facilities necessitates an approach for evaluating Westcott 𝑔 factors based on arbitrary non-Maxwellian spectra. In this paper, we present an approach for calculating 𝑔 factors with user-specified neutron spectra, and we demonstrate these methods to obtain Westcott 𝑔-factors for guided- and cold-neutron beams at the Budapest Research Reactor and the Forschungsreaktor München II reactor. As part of this work, open-source software has been developed that can be used to perform these calculations for applications in PGAA and NAA experiments.
Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.
An evaluation of the average prompt fission neutron multiplicity, $\bar{ν}_p$, of 239 Pu(n,f) is shown. This evaluation includes (a) the correlated fission model CGMF, and (b) a detailed analysis of past and recently published experimental data. Using CGMF-calculated $\bar{ν}_p$ as prior enables to link, through the use of evaluated model input parameters, $\bar{ν}_p$ to other fission observables such as the prompt fission neutron spectrum (PFNS), preneutron emission fission yields as a function of mass, and the average total kinetic energy of the fragments. These evaluated parameters produce realistic predictions of many fission observables, while the evaluated $\bar{ν}_p$ agrees well (χ 2 ≈ 1) with data. Moreover, with the new evaluated $\bar{ν}_p$, the effective neutron multiplication factor of fast Pu ICSBEP critical assemblies are predicted with a mean bias of 58 pcm compared to 18 pcm with ENDF/B-VIII.0, when paired with a new 239 Pu PFNS and fission cross section. Due to these encouraging validation results, the evaluated $\bar{ν}_p$ is currently part of a release candidate for the 239 Pu ENDF/B-VIII.1 file. Hence, a correlated fission model was used for the first time for evaluating $\bar{ν}_p$ that is of evaluation quality. This is an important step towards consistent evaluations of prompt fission observables.
The inverse-kinematics p( 7 Li,n) 7 Be reaction produces forward-focused neutron emission, offering enhanced usable flux and reduced shielding requirements. Reliable simulation of such neutron fields is essential for the development of compact accelerator-based neutron sources. In this study, a PHITS-based simulation framework for the reaction was experimentally assessed using fast-neutron measurements. Forward-directed neutrons were measured with a diamond neutron detector and quantitatively compared with simulations with newly prepared IK-Frag cross-section file based on the proton-induced reaction data in ENDF/B-VIII.1, TENDL-2023, and JENDL-4.0/HE. Measurements and simulations were performed for incident 7 Li 3+ energies ranging from 15.0 to 25.0 MeV using a 50 μm-thick polypropylene target. For all conditions, the PHITS-based simulation framework reproduced the deposited energy spectra at the correct order of magnitude. The comparison of deposited energy spectra in the diamond detector showed high correlation coefficients across all investigated energies, indicating reasonable agreement in spectral shape between simulations and measurements. This work represents an initial step toward establishing a benchmark for PHITS simulations of the inverse kinematic reaction between an incident lithium-ion and a proton target.
Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.
Zirconium hydride (ZrH x ) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high- enriched uranium (HEU) and ZrH x in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of 235 U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH x .