Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CGMF”

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

Energy-dependent optimization of the prompt fission neutron spectrum with CGMF

Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

$\overline{TKE}$ Parameterization and $\bar{v}$ Uncertainty Analysis for CGMF

Previous work was performed on tuning CGMF parameters for 235 U, 238 U, and Plutonium isotopes. Now work is being done to tune minor uranium isotopes. However, uranium isotopes like 232 U and 236 U have almost no experimental data. We are applying cross-isotope models to extrapolate and tune CGMF on isotopes that lack experimental data. There exist several internal CGMF physics quantities that affect the output of CGMF—multi-chance fission probability, excitation energy sharing, spin-cutoff factor, spin scaling, and fragment total kinetic energy to name a few. The mean fragment total kinetic energy, $\overline{TKE}$, is particularly interesting because of its strong anti-correlation with $\bar{v}$. We are most interested in the mean fragment total kinetic energy before neutron emissions. $\overline{TKE}$ is assumed to be pre-neutron emission unless otherwise stated. Currently in CGMF, the $\overline{TKE}$ model for 233,234,235,238 U are tuned independently to reproduce ν for the associated isotopes. In this report, we will tune a cross-isotope $\overline{TKE}$ model to experimental $\overline{TKE}$ data for 232,233,234,235,236,238 U. Because of the unreliable and sparse nature of $\overline{TKE}$ experimental data, future work will use more reliable experimental $\bar{v}$ data to infer the $\overline{TKE}$ model (and likely other internal CGMF parameters) for uranium isotopes. Such work has been performed previously using a sensitivity analysis and Kalman filter methods.

07 ISOTOPE AND RADIATION SOURCES↗

Open-source release of CGMF 1.1 and Integration into the MCNP6.3 ® Code [Slides]

As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Consistent $\overline{ν}$ evaluation for minor Pu isotopes with CGMF

This report follows up on the evaluation work done in 2023 where we developed a procedure and performed the first consistent evaluation of the average prompt neutron multiplicity for minor Pu isotopes, including 238,240-242 Pu, using CGMF. Details on the necessary updates to the release version of CGMF, the experimental data and experimental uncertainty quantification that went into the evaluation, and this first evaluation effort are documented in and will not be repeated in this report. Instead, we discuss here the further investigations into the CGMF model space to improve the evaluation results from FY23. Additionally, we note that these evaluation results have been transformed into the ENDF format for mean values and covariances, and validation with critical assemblies has been performed; that work is documented elsewhere. In this short report, we discuss the updated evaluation efforts performed during FY24 in Sec. 2, show results from CGMF for other prompt fission observables in Sec. 3, and then briefly conclude in Sec. 4.

07 ISOTOPE AND RADIATION SOURCES↗

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↗

Evaluated 238 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 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here 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 validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 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 for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗

Producing ENDF/B-quality Evaluations of 239 Pu(n,f) and 235 U(n,f) Average Prompt Neutron Multiplicities using the CGMF Model

This report documents evaluations of the neutron-induced 239 Pu and 235 U average prompt neutron multiplicities, $\overline{v}_p$, using the CGMF model developed at LANL. These evaluations are not updates of existing ENDF/B-VIII.0 nuclear data, but were both re-done from "scratch". That is, all experimental data were extracted from EXFOR, re-normalized to the newest nuclear data representing monitor observables, and detailed uncertainties were estimated from information in EXFOR as well as from templates of expected measurement uncertainties. We also included experimental data that were not yet available for ENDF/B-VIII.0, for instance, the data of Marini et al. for 239 Pu(n,f) $\overline{v}_p$ or the data of Khoklov et al. for 235 U(n,f) $\overline{v}_p$. Finally, we include the CGMF model; it computes $\overline{v}_p$ based on model parameters that at the same time calculate fission fragments as a function of mass, Y (A), the average total kinetic energy as a function of incident-neutron energy, $\langle TKE \rangle$(E inc ), etc. Such a detailed fission model that ties together many fission quantities has not been used to date for any $\overline{v}_p$ evaluation in ENDF/B. Here, we show that we can get evaluated 235 U(n,f) and 239 Pu $\overline{v}_p$ that not only correspond well to experimental $\overline{v}_p$, but that the associated evaluated model parameters also yield parameterizations of Y (A), $\langle TKE \rangle$(E inc ), etc., that correspond well to their respective experimental data. The new evaluated 239 Pu(n,f) $\overline{v}_p$ was also combined into a 239 Pu test file with the newest nuclear data for the 239 Pu prompt-fission neutron spectrum and fission cross sections. This new 239 Pu file performed reasonably well in predicting PU-MET-FAST assemblies, reaction rates in Jezebel and Flattop and three LLNL pulsed spheres. Due to that, the new 239 Pu(n,f) $\overline{v}_p$ evaluation presented here is being considered for ENDF/B-VIII.1. Hence, we can show here that CGMF is able to produce ENDF/B quality $\overline{v}_p$ nuclear data.

239Pu, 235U, Average Prompt-fission Neutron Multip↗

Baseline CGMF Calculations for Pu Isotopes

Updates were made to CGMF to run neutron-induced fission calculations on Pu isotopes of A =238, 239, 240, 241, and 242 (the release version of CGMF can calculate 239,241 Pu(n,f) and 238,240,242 Pu(sf) only). Here, we show the baseline calculations from CGMF for these neutron-induced fission reactions, including updates made for total kinetic energy systematics and fission barriers fit to be more consistent with the neutron-induced fission cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Consistent $\overline{ν}$ evaluation for minor U isotopes with $\tt{CGMF}$

Following several successful prompt $\overline{ν}$ evaluations using $\tt{CGMF}$, including consistent evaluations for minor Pu isotopes, we detail in this report our efforts to perform a consistent $\overline{ν}$ evaluation for minor U isotopes during FY25. Although we have not yet produced a finalized evaluation, we present the progress that we have made towards such an evaluation for 232,233,234,236,237,239 U prompt $\overline{ν}$. Our milestone explicitly calls out evaluations for 233 U, 234 U, and 236 U, however, to better constrain the model with reliable experimental $\overline{ν}$ data, we also include 235 U and 238 U in the evaluation procedure. Then, we additionally produce evaluations for 232 U, 237 U and 239 U $\overline{ν}$ as a byproduct. Elsewhere, we will report our efforts on a stand-alone 233 U $\overline{ν}$ evaluation. This report is organized in the following manner. In Sec. 2, we briefly outline the updates to CGMF that were needed to be able to calculate all of these minor U fission reactions. The experimental data overview is given in Sec. 3. The evaluation methodology and results are presented in Secs. 4 and 5, respectively. Finally, we conclude and outline work for FY26 in Sec. 6.

07 ISOTOPE AND RADIATION SOURCES↗

Fitting $\overline{ν}$ for Minor Pu Isotopes

After successful fitting of prompt ν for 235 U, 238 U, and 239 Pu(n,f) using CGMF, we have moved on to the minor plutonium isotopes. Minor isotopes pose a greater challenge both because there is less data available and we do not, by default, have parametrizations in CGMF already. Initially, we had planned to mitigate these challenges—and provide consistency within CGMF—by taking a stepped approach to the optimization. Continuing from the 239 Pu work, we would then fit 241Pu(n,f) $\overline{ν}$, where there is also a number of experimental measurements, keeping consistency between the parameters that are included in the calculation for 239 Pu and 241 Pu. With these parametrizations settled, we could consistently optimize 240 Pu(n,f) $\overline{ν}$. Following that step, we would move on to 242 Pu(n,f) (and increasing neutron number) and 238 Pu(n,f) (and decreasing neutron number). By including the fissioning systems in this manner, we should be able to minimize the unknown parameters in CGMF. We were also able to develop systematics for the CGMF input parameters along the Pu isotopic chain. However, after studies in the beginning of FY23, where we found reasonable agreement between CGMF and experimental prompt neutron multiplicities for 238−242 Pu, using a compound mass dependent parametrization from, and from discussions with I. Stetcu and P. Talou, we instead performed the evaluation of all five isotopes simultaneously. For all parameters in CGMF, they are either the same for each compound nucleus, or have a dependence on the compound mass. The details of these updates are given in. This report builds upon those details. In this short report, we first give an example of how we have updated CGMF to include 240 Pu(n,f) and 242 Pu(n,f), keeping consistency with the current 239 Pu and 241 Pu calculations, but without rigorous optimization (Sec. 2). Then, we will show in Section 3 what experimental data exist for the various isotopes to provide some insight into why 241 Pu and 239 Pu are used as anchor points. The evaluation is detailed in Section 4, and the comparison between CGMF calculations using the evaluated parameters and other prompt observables besides average neutron multiplicity are discussed in Section 2.2. Finally, we present conclusions and future work in Section 5.

07 ISOTOPE AND RADIATION SOURCES↗

Fitting $\overline{\nu}$ for minor Pu isotopes

After successful fitting of prompt $\overline{\nu}$ for 235 U, 238 U, and 239 Pu(n,f) using CGMF, we will move on to the minor plutonium isotopes. Minor isotopes pose a greater challenge both because there is less data available and we do not, by default, have parametrizations in CGMF already. To mitigate these challenges— and provide consistency within CGMF—we will take a stepped approach to the optimization. Continuing from the 239 Pu work, we will then fit 241 Pu(n,f) $\overline{\nu}$, where there is also a number of experimental measurements, keeping consistency between the parameters that are included in the calculation for 239 Pu and 241 Pu. With these parametrizations settled, we can consistently optimize 240 Pu(n,f) $\overline{\nu}$. Following that step, we will move on to 242 Pu(n,f) (and increasing neutron number) and 238 Pu(n,f) (and decreasing neutron number). By including the fissioning systems in this manner, we should be able to minimize the unknown parameters in CGMF. We will possibly also be able to develop systematics for the CGMF input parameters along the Pu isotopic chain. This work can serve as a guide to broadening the reactions available in CGMF. In this short report, we first give an example of how we have updated CGMF to include 240 Pu(n,f) and 242 Pu(n,f), keeping consistency with the current 239 Pu and 241 Pu calculations, but without rigorous optimization (Sec. 2). Then, we will shown in Section 3 what experimental data exist for the various isotopes to provide some insight into why 241 Pu and 239 Pu are used as anchor points.

07 ISOTOPE AND RADIATION SOURCES↗

Correcting the PFNS for more consistent fission modeling

For FY20, we had a deliverable to write a report detailing efforts to simultaneously evaluate both the prompt neutron multiplicity, $\overline{ν}$, and the prompt neutron fission neutron spectrum, PFNS, using CGMF. CGMF is the LANL-developed fission fragment decay code that consistently evaporates prompt neutrons and γ rays using the Hauser-Feshbach statistical theory of compound reactions. The decay begins by constructing the initial conditions of the fission fragments, then decaying each one from the excited state by neutrons and γ rays, conserving energy, momentum, spin, and parity in each step of the emission. The initial conditions of the fragments, along with the multiplicity, energy, and direction of each emitted neutron or γ ray, are recorded, allowing for the full reconstruction of the fission event. These event histories allow us to reconstruct average quantities, as well as correlations between observables, that can be compared with experimental or evaluated data. In that initial report, although there was already a favorable comparison between $\overline{ν}$ from CGMF, experiment, and the current ENDF/B-VIII.0 evaluation, we showed that there was still significant work to be done to improve the PFNS from CGMF. Historically, the PFNS is calculated too soft by Hauser-Feshbach fission models, and CGMF is no exception. The incorrect shape presents a significant challenge in fission modeling, including for our understanding of the fission process and for our ability to consistently calculate and predict a variety of prompt fission observables (such as fission fragment initial conditions, neutron and γ-ray multiplicities and energies, and the correlations between all observables). In our companion report, we detail our success in using CGMF to evaluate $\overline{ν}$. Although not included in the optimization explicitly, we also keep the initial conditions of the fission fragments physical, along with reproducing reasonably well the neutron multiplicity distribution. As we would expect from the sensitivities calculations from, the average neutron energies change very little from the $\overline{ν}$ optimization along with the PFNS (as will be shown in Sec. 2.6). The conclusion was that the global and statistical models would have to be investigated instead of just the fission fragment initial conditions (as is sufficient for $\overline{ν}$). This report details those efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Compare $\overline{ν}$ parameterizations for 235 U(n,f) with fitted fission barriers

In this document, we compare observables computed with CGMF to some available data for 235 U(n,f), to investigate the multi-chance fission probabilities (fission barriers from CoH were refit) and two TKE parametrization. The four parametrizations include i) default CGMF multi-chance fission probabilities with bends in TKE for all fission chances, ii) default CGMF multi-chance fission probabilities with a TKE bend in only first-chance fission iii) refitted multi-chance fission probabilities with bends in TKE for all fission chances, and iv) refitted multi-chance fission probabilities with a bend in TKE for only first-chance fission. In all figures where the four parametrizations are on the same plot, the ones with bends in TKE for all chance fission are in black, with only first-chance bend in red, the default multi-chance fission probabilities are dashed and the new multi-chance fission probabilities are the solid lines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Producing Evaluation-quality 239 Pu Average Prompt Fission Neutron Multiplicities using a Correlated Fission Model

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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and evaluating 239 Pu and 235 U PFNS and average prompt-neutron multiplicity [Slides]

The following are currently in VIII.1 LANL and IAEA test files: 239 Pu nu-bar including CGMF modeling and CEA data, 239 Pu PFNS including Chi-Nu & CEA data, 239 Pu(n,f) cross section including fission TPC data. The following are currently being tested: 235 U nu-bar including CGMF modeling and 235 U PFNS including Chi-Nu data. Upcoming tasks include the correction of 235 U Chi-Nu PFNS at higher E inc , benchmarking 235 U PFNS and nu-bar evaluations, and getting 238 U nu-bar.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shedding light on the 239Pu fission source term with new high-precision experiments and advanced fission modeling

In the last decade, there has been a renaissance of fission research resulting in new high-precision experiments and advanced fission modeling. For instance, the Chi-Nu and CEA teams supplied, for the first time, the 239 Pu prompt fission neutron spectrum (PFNS) for broad ranges of incident and outgoing neutron energies. The CEA team also measured 239 Pu average prompt neutron multiplicities, ν ̄ p , with lower statistical uncertainties and a technique significantly different than the one used in the past. The NIFFTE collaboration provided 239 Pu( n ,f)/ 235 U( n ,f) cross section shape ratios with uncertainties below 1% utilizing a novel detector type. Advanced fission event generators were developed, among them CGMF, FIFRELIN, FREYA, and GEF, which calculate post-scission fission observables in a correlated manner. These new experimental data and more consistent fission models change the evaluated PFNS, ν ̄ p , and ( n ,f) cross sections, some only modestly, compared to ENDF/B-VIII.0. In turn, the individual new nuclear data distinctly change simulated effective neutron multiplication factors of fast critical assemblies, but their combined impact is small, while affecting the prediction of LLNL pulsed sphere neutron leakage spectra and reaction rates only within experimental uncertainties. Also, the parameters obtained from fitting to ν ̄ p reproduce various post-scission fission observables within the uncertainties of experimental data. This indicates that new differential experiments and consistent fission modeling reduce compensating errors present in ENDF/B-VIII.0.

fission cross section↗