Search NASA⌕ Search

SEARCH · Search NASA

Results for “keff”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Preliminary Investigation of Target keff Search for Fuel Feed in MSRs

Fluid fuel molten salt reactor (MSR) concepts may achieve continuous-at-power operation with continuous or batch fresh fuel feed to counter reactivity loss from fuel burnup and fission product buildup. Accurate modeling and simulation of MSRs requires consideration of online fuel feed during depletion as well as a capability to determine the required feed amount to maintain a target multiplication factor (keff). As a first step supporting a full implementation of a new search capability into the SCALE code system’s TRITON reactor physics and depletion sequence, a preliminary investigation was performed to identify the best approach for a critical feed search capability considering computation time and accuracy. Using a simple iterative approach in which initial feed amounts were guessed based on previous depletion steps, the main drivers of computation time were found as the Monte Carlo neutron transport settings and the applied tolerance for the target keff. In the applied model, with tight convergence criteria (keff convergence to 8 pcm and a tolerance of ±15 pcm), the critical search calculation required twice the number of neutron transport calculations compared to a depletion calculation without critical search. In contrast, with a relaxed tolerance of ±50 pcm, the overall computation time for the critical search for merely increased by 30-40%. Further studies will be performed before implementing the critical search capability into TRITON.

Faulkner, Jonathon [ORNL]↗
DOE OSTI2023

How can a diverse set of integral and semi-integral measurements inform identification of discrepant nuclear data?

Nuclear data are used for a variety of applications, including criticality safety, reactor performance, and material safeguards. Despite the breadth of use-cases, the effective neutron multiplication factor, keff, of ICSBEP critical assemblies are primarily used for nuclear data validation; these are sensitive to specific energy regions and nuclides and are unable to uniquely constrain nuclear data. As a consequence, general-purpose nuclear data libraries, such as ENDF/B-VIII.0, may have deficiencies that, while not apparent in criticality applications, negatively impact other applications, such as non-destructive analysis of special nuclear material and neutron diagnosed subcritical experiments. Recent work by the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project developed a machine learning tool, RAFIEKI, which uses random forests and the SHAP metric to determine which nuclear data contribute most to predicted bias between measured and simulated responses (e.g. keff). This paper contrasts RAFIEKI analysis applied to keff only against RAFIEKI analysis with keff paired with either LLNL pulsed sphere measurements or subcritical benchmarks. Two examples show that a) including pulsed sphere measurements substantially increases 9Be nuclear data importance to bias between 2 and 15 MeV, and b) including subcritical benchmarks has the potential for disentangling compensating errors between 240Pu (n,el) and (n,il) cross-sections between 0.1 and 10 MeV. These results show that RAFIEKI analysis applied to response sets that include, but go beyond, keff can aid nuclear data evaluators in identifying issues in nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗
DOE OSTI2026

Nuclear Data Libraries Sensitivity Studies for ORSA Using SCALE

Subcritical assemblies offer valuable training capabilities in nuclear criticality safety (NCS) for individuals handling fissile material. At Oak Ridge National Laboratory(ORNL), the Oak Ridge Subcritical Assembly (ORSA), a new experimental facility, is being established to provide hands-on training for the Nuclear Criticality Safety Program (NCSP).It is essential to accurately determine the neutron multiplication factor (keff) to ensure that ORSA remains subcritical and safe during operations. This study investigated the sensitivity of keff to variations across Evaluated Nuclear Data File (ENDF/B) libraries, consisting of ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1. The analysis was conducted using the CSAS6 sequence in the SCALE-6.3 code system. Individual isotopes in the ORSA model were replaced one at a time with ENDF/B-VII.1 as the base library and changing to ENDF/B-VIII.0or ENDF/B-VIII.1. The results demonstrated that the changes in keffof the nuclides associated with the ORSA model were mostly within the uncertainty of the base model (~24 pcm), except for primary nuclides like Uranium-235and H-poly with few other nuclides. The relative delta keff values of Uranium-235 and H-poly, expressed in pcm, were +280 and -212 in ENDF/B-VIII.0 and +332 and -275 in ENDF/B-VIII.1, respectively, which were notable changes in reactivity. These results demonstrated that ORSA was largely insensitive to variations across these nuclear data libraries.

Hong, Evan [North Carolina State University]↗
DOE OSTI2025

The HTR-Proteus Benchmark: Analysis and Use as a Verification and Validation Case

This presentation outlines the evaluation and application of the HTR-Proteus benchmark as a verification and validation (V&V) case for advanced reactor modeling tools. The work supports the U.S. Department of Energy’s HALEU Availability Program (HAP) and the joint DOE/NRC DNCSH project, which aims to reduce criticality safety uncertainties in commercial-scale HALEU fuel cycle and transportation systems. The HTR-Proteus experiments, conducted at the Paul Scherrer Institute, provide high-fidelity data for TRISO-fueled, graphite-moderated pebble bed reactors with high neutron leakage—conditions relevant to HALEU transport scenarios. This study focuses on Cores 4.2 and 4.3 of the HTR-Proteus benchmark, analyzing key sources of uncertainty including pebble packing, TRISO particle positioning, and core height. Using Project Chrono for realistic pebble geometries and Serpent, SHIFT, and MCNP for neutronics simulations, the study quantifies the impact of these uncertainties on the effective multiplication factor (keff). Results show that a sample size of 110 pebble configurations is sufficient to converge keff, with ±30 pcm uncertainty due to packing randomness. TRISO positioning and core height variations also significantly influence keff, highlighting the importance of detailed modeling in V&V efforts. The benchmark serves as a valuable test case for validating the Griffin reactor physics code and improving confidence in HALEU system simulations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

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.

Critical benchmarks↗
DOE OSTI2026

Effects of Uranium Impurities in Downblended HEU on HTGR Performance

Many advanced reactor designs require fuel enriched between 5% and 20% 235 U. To assist in producing fuel at these enrichment levels, government-owned inventories of highly enriched uranium can be downblended. However, fuel produced from these inventories contain uranium impurities that are not often found when enriching natural uranium or accounted for when modeling reactor cores. To address this concern, this work models reactor designs like the X-energy Xe-100 and the Ultra Safe Nuclear Company’s Micro Modular Reactor, and compares their performance with fuel from enriching natural uranium to fuel from downblended highly enriched uranium. This paper evaluates the models based on the effective neutron multiplication factor, k eff , effective delayed neutron fraction, β eff , and energy- and spatially dependent neutron flux, ϕ, as well as the fuel, coolant, moderator, and total reactivity temperature feedback coefficients, α F , α C , α M , and α T . The results show that the fuel from downblended highly enriched uranium inventories leads to differences in each of the metrics, especially in the keff values. In the Xe-100–like and Micro Modular Reactor–like models, k eff changes by about 1400 pcm and up to 1200 pcm, respectively. Total reactivity feedback coefficients α T are negative with the impure fuels and the keff values remain above 1 for each core configuration and fuel composition. These results show that the impure fuel compositions do not necessarily prevent achieving key design parameters, such as cycle length, or from operating in a safe condition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Cerberus and the Zeus Suite of Critical Experiment Benchmarks [Slides]

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

ENDF/B-VIII.1↗
DOE OSTI2023

MCNP Neutronic Design Optimization of an Accelerator-Driven Subcritical Assembly for Mo-99 Production

Molybdenum-99 (99Mo) has been in short supply in the recent decades. The most common productionmechanism, irradiation in a reactor fueled by highly enriched uranium (HEU), raisesproliferation concerns. Novel efforts are being pursued to meet the demand for 99Mo withoutHEU irradiation, relying on the processing of low-enriched uranium or natural uranium. Onesuch concept is to use an electron linear accelerator to generate bremsstrahlung that, through aphotonuclear reaction, can help multiply neutrons in a subcritical assembly. The goal of this workwas to neutronically optimize the subcritical assembly for keff and 99Mo production. This optimizationis unique for The design was optimized in terms of reflector type, pitch, and LEU andNU mass. Results of the paper show the effects of fuel mass, pitch, and reflector material on keffand the NU fission fraction. keff was found to be mostly driven by LEU mass, and the NU fissionfraction was dependent on several parameters, primarily pitch and mass.

Conant, Andrew↗
DOE OSTI2025

Automated Direct Perturbation Calculations with SCALE TSUNAMI [Abstract]

In nuclear criticality safety analysis, the sensitivity of the eigenvalue keff to uncertainties in nuclear data and its evaluation are crucial. The TSUNAMI sequences within the SCALE code system offer users various options with both multigroup (MG) and continuous-energy (CE) 3D Monte Carlo (MC) transport capabilities for calculating keff sensitivity coefficients and storing them in a sensitivity data file (SDF). Each methodology available in TSUNAMI offers distinct advantages and limitations, and its effectiveness can vary based on the specific problem being solved. As a best practice, practitioners typically use the direct perturbation (DP) method as a confirmatory step alongside their sensitivity calculations to verify the accuracy of the sensitivity data generated. In this process, DP calculations are usually performed on select nuclides, those considered most important for validating their total sensitivities. However, because of code limitations, analysts use a workaround method when conducting DP calculations for a single nuclide: rather than perturbing the nuclide's microscopic cross section, an equivalent number density for this nuclide is calculated to reflect the effect of a change in the macroscopic cross section due to a perturbation in the microscopic cross section. The current approach requires rerunning the CSAS criticality calculation several times with model changes. Although this method can yield results with acceptable accuracy, it is labor-intensive and prone to errors.

AZURE: SAMMY↗
DOE OSTI2022

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation [Slides]

During the evaluation of Oak Ridge National Laboratory Health Physics Research Reactor operation data for critical benchmark creation, a detailed, functional SCALE model was built, and experimental data was evaluated for the creation of a critical benchmark. The estimated experimental uncertainty is about 3800 pcm, which is considered very high. The relative difference between expected and calculated keff values is about 1.5 %, also considered very high. It was concluded that a good quality critical benchmark worthy of the ICSBEP standards cannot be created from HPRR data in the present conditions. Locating the HPRR fuel to obtain an uncertainty on the density would aid in solving this issue.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

SCALE inventory and reactivity analysis as part of the Hermes 2021 PSAR review

The readiness of SCALE for comprehensive studies of pebble-bed reactors has been demonstrated through detailed analysis of a fluoride salt–cooled, high-temperature pebble-bed reactor (PB-FHR). The methods developed for pebble-bed reactor modeling in SCALE, particularly for inventory generation, have proven effective in gaining insights into the reactor physics of this advanced reactor. Excellent agreement with another code package has been observed, further highlighting SCALE’s strong performance. The SCALE results supported the US Nuclear Regulatory Commission’s construction permit application review of the Hermes low-power PB-FHR demonstration reactor. A SCALE model of the Hermes reactor was developed at Oak Ridge National Laboratory using information from the Preliminary Safety Analysis Report (PSAR) and supplemented with publicly available data. SCALE reactivity coefficient simulations reproduced PSAR results within 1σ statistical uncertainties. Sensitivity studies emphasized the importance of graphite specifications for accurate keff predictions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Resolved Resonance Region Analysis of 206 Pb, 207 Pb, and 208 Pb for Next Generation Lead-Cooled Fast Systems

Recently, great interest has been generated in using lead as a coolant for fast neutron systems and as a result it is important to investigate the ENDF/B-VIII.0 isotopic evaluations that comprise stable le ad. To this end, resonance parameters for 206 Pb, 207 Pb, and 208 Pb were re-evaluated because their resolved reso­nance regions extend beyond 0.5 MeV meaning resonance parameters used in reconstructing cross sections and elastic scattering angular distributions impact fast systems. The impact of resonance parameters is demon­strated by the differences between the evaluations in predicting experimental results from the fully modeled RPI Quasi-Differential Scattering Experiment via MCNPv6.2. In addition, MCNP KCODE calculations of lead-sensitive fast spectra critical benchmarks showed variations of keff on average of 400 pcm, caused solely from differences in elastic scattering angular distributions in 208 Pb. Re-evaluation entailed fitting data with the program SAMMY.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Continuous-Energy ENDF/B-VIII.0 Cross Section and SCALE 6.2.4 Performance for Nuclear Criticality Safety Applications: 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238 U, 239 Pu

As part of the nuclear data evaluation and validation cycle, the ENDF/B-VIII.0 cross-section library released in 2018 requires testing to determine areas of improvement and deterioration. Previous work by the authors investigated the performance of 16 O, 56 Fe, and 63,65 Cu cross sections, with this study acting as an extension of the prior work. In addition to the isotopes and nuclear criticality safety benchmarks of interest to the prior work, benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook were selected for their k eff sensitivity to 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238U, or 239 Pu cross sections and were modeled in the SCALE code system maintained by Oak Ridge National Laboratory. In total, 253 benchmark configurations were selected for their sensitivities and modeled using SCALE 6.2.4 Criticality Safety Analysis Sequences (CSAS) continuous-energy Monte Carlo k eff calculations. This collection includes and expands upon the 99 benchmarks in the prior work. The AMPX-processed ENDF/B-VIII.0 library was decomposed into individual ENDF/B-VIII.0 datum libraries for each isotope of interest. Doing so allowed for the individual substitution of an ENDF/B-VIII.0 cross section in the place of ENDF/B-VII.1, determining isotope-specific effects of ENDF/B-VIII.0 relative to ENDF/B-VII.1. Full library calculations with entirely ENDF/B-VII.1 data or entirely ENDF/B-VIII.0 data were also executed. As a measure of performance, the average relative deviation was determined as the ratio of the deviation between calculated and experimental keff to the propagated calculational and experimental uncertainty. With calculated full library and isotope-specific ENDF/B-VIII.0 k eff ’s, an optimized combination of data libraries was estimated and confirmed with SCALE calculations. This showed that reverting 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 resulted in improved performance relative to the full ENDF/B-VIII.0 library. Across all 253 benchmarks, the average relative deviation was 1.29σ for the full ENDF/B-VII.1 library, 1.17σ for the full ENDF/B-VIII.0 library, and 0.97σ for the optimized combination. The reversion of 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 in the 99 benchmarks of the prior work resulted in further improved experimental agreement compared to the previously reported improvement from 16 O and 65 Cu alone. Therefore, it is suggested that applications with significant sensitivities to 239 Pu, 58 Ni, 16 O, and 65 Cu consider their choice of nuclear data library.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Release of Evaluated 235 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 235 U from 200 keV to 15 MeV that is a potential release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation had to be re-done from "scratch", as the input to the $\overline{v}_p$ evaluation of the previous library, ENDF/B-VIII.0, was lost. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. Another major difference to ENDF/B-VIII.0 is that this evaluation includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to perform two types of validation for the new 235 U $\overline{v}_p$: On the one hand, one can employ evaluated CGMF parameters obtained from fitting to experimental 235 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data indicating the evaluated 235 U(n,f) $\overline{v}_p$ are physical. On the other hand, one can validate 235 U $\overline{v}_p$ with respect to integral responses such as fast ICSBEP critical assemblies or LLNL pulsed spheres. LLNL pulsed-sphere neutron-leakage spectra are minimally impacted by the new 235 U $\overline{v}_p$ as these experimental data are shape data and the $\overline{v}_p$ would mostly lead to a change in normalization of the data as the spheres are relatively thin (0.7 and 1.5 mean-free path) and, thus, mostly depend on 235 U $\overline{v}_p$ from 12-15 MeV. The change in the predicted effective neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, however, is large compared to values using ENDF/B-VIII.0 and experimental k eff : The average bias is 108 pcm across all studied k eff values versus 12 pcm for ENDF/B-VIII.0. A reasonable performance in simulating keff (mean bias of 14 pcm) can be retained by tweaking 235 U $\overline{v}_p$ from 3-5 MeV, and combining it with a recent 235 U PFNS evaluation that is also a ENDF/B-VIII.1 release candidate.

235U↗
DOE OSTI2023

Parametric Study of 239 Pu and 235 U Mass Equivalency in Spherical Geometries

In this technical document, parametric studies were performed to determine if a bounding mass equivalency can be used to simplify multi-material models in MCNP6 for criticality evaluations. The studies explored the simplification of MCNP6 models by substituting materials, particularly 239 Pu for 235 U. The initial hypothesis for the material simplification is based on the subcritical mass limits from the ANSI/ANS- 8.1 standard which states that for water reflected metal spheres the limits are 20.1 kg 235 U and 5.0 kg 239 Pu. Using those masses, it was hypothesized that a 20.1 kg sphere of 235 U and a 5.0 kg sphere of 239 Pu could be considered to have the same multiplication. Given that equivalency, 4kg 235 U could be modeled by 1kg 239 Pu (i.e., a mass ratio of 4:1). The parametric studies compare the materials 239 Pu and 235 U independently and in multi-material spherical models. Critical mass, critical thickness, critical mass ratio and critical thickness ratio values are presented as functions of inner radius. keff values are presented as a function of inner radius at different mass ratios and total masses. Results indicate that the conservative bounding mass equivalency is one. That is, if substituting 239 Pu for 235 U it should be done on a gram for gram basis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗