Search NASA⌕ Search

Engineering topics

Favorite, Jeffrey A.

Publications and source records attributed to Favorite, Jeffrey A..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

SOURCES4D

SOURCES is a code for computing neutron source rates and spectra from spontaneous fission (including delayed neutrons) and (alpha,n) reactions in homogeneous materials and (alpha,n) reactions in single-interface and two-interface geometries. SOURCES is a Los Alamos National Laboratory (LANL) code that is written in FORTRAN and distributed through the Radiation Safety Information Computation Center (RSICC). LANL’s last release of SOURCES to RSICC was SOURCES4C in 2002. This disclosure covers the latest version of SOURCES, SOURCES4D. This version adds sensitivity capabilities for (alpha,n) sources in homogeneous materials. Specifically, SOURCES4D writes new output that can be used to calculate, in post-processing, first and second derivatives of the (alpha,n) source rate density and spectrum with respect to nuclide densities in a homogeneous material and first derivatives of the (alpha,n) source rate density and spectrum with respect to nuclide stopping powers and (alpha,n) cross sections (nuclear data). These derivatives are useful for uncertainty quantification, predictive modeling, and other applications in neutron transport problems.

Favorite, Jeffrey A.↗
DOE OSTI2024

(U) A Code System for Cross-Section Uncertainty Propagation for PARTISN

A rudimentary code system has been developed for propagating neutron cross-section uncertainties to response uncertainties using the PARTISN multigroup discrete-ordinates neutron transport code. The code system uses the first-order sandwich rule with sensitivities computed using SENSMG and covariances computed using NJOY. Databases of reaction cross sections and covariances using ENDF/B-VIII.0 data in the LANL standard 30-group structure have been precalculated. The final uncertainty in a k eff example problem compares well with the result from TOFFEE, a similar code that propagates neutron cross-section uncertainties to response uncertainties using the MCNP6 Monte Carlo neutron transport code. Presently, neither code system propagates uncertainties in $\overline{v}$, $\overline{μ}$, or χ, and neither code system propagates nuclide-nuclide covariances

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Chlorine Worth Study Nuclear Data

The Chlorine Worth Study (CWS) was a series of experiments that took place at the National Criticality Experiments Research Center (NCERC), operated by Los Alamos National Laboratory (LANL). The focus of the experiments was to develop new integral benchmarks for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook with high sensitivity to chlorine in the thermal neutron energy region and which match sensitivities of aqueous chloride operations at LANL. This work discusses the experiment and how the experimental and simulated results using different nuclear data libraries compare to each other.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Scientific Computing [Slides]

Presentation by the Los Alamos National Laboratory Radiation Transport Applications Group (XTD-RTA) for the CONNECT2 Undergraduate Research Experience Seminar at the University of Texas, San Antonio.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

(U) Release of SOURCES4D Featuring Sensitivity Capabilities

SOURCES is a code for computing neutron source rates and spectra from spontaneous fission (including delayed neutrons) and (α,n) reactions. The first release since 2002, SOURCES4D, is now available. The main updates include new output that allows the user to calculate, in post-processing, first and second derivatives of the (α,n) source rate density and spectrum with respect to nuclide densities in a homogeneous material and first derivatives of the (α,n) source rate density and spectrum with respect to nuclide stopping powers and (α,n) cross sections (nuclear data). The value of π was made more accurate in order to eliminate negative spontaneous-fission sources. Otherwise, SOURCES4D has made no changes to the physics or data of neutron source calculations used in SOURCES4C. First and second derivatives are calculated in an example problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

(U) A Tool to Set Up an MCNP6 Input to Use Correlated Sampling with Batch Statistics

Correlated sampling using batch statistics with MCNP6’s tally fluctuation charts (TFCs) is a powerful means of reducing the statistical uncertainties of various combinations of tallies. Because MCNP6 prints a maximum of 20 entries in the TFCs, at least five inputs must be run to obtain the recommended number of at least 100 batches. A new software tool, MAKE_COSUBS, reads an MCNP6 input file and sets up a sequence of input files to be run. The user runs MAKE_COSUBS for the unperturbed and all perturbed MCNP6 inputs, runs MCNP6 for the inputs, and finally analyzes the outputs using COSUBS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

(U) Correlated Sampling Using Batch Statistics to Reduce the Uncertainty of Combinations of KSEN Outputs with MCNP6

The relative sensitivity of k eff to the densities of nuclides in a material are combined to compute relative sensitivities to other inputs. When computed in a single Monte Carlo run, the nuclide density sensitivities are correlated, and the statistical uncertainties propagated to other inputs will be incorrect unless those correlations are accounted for. Equations are presented to apply correlated sampling using batch statistics for the sum of an arbitrary number of random tallies and the difference of two random tallies when each is multiplied by a different constant. When correlated sampling is used on a recent benchmark evaluation, the correct statistical uncertainties for certain combinations of sensitivities are dramatically smaller than the incorrect uncertainties. A user-controlled, regular output of MCNP6’s KSEN sensitivities that allows batch statistics to be applied to combinations would be extremely valuable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

(U) Average Scattering Angle Cosine Sensitivities

The Evaluated Nuclear Data File (ENDF) provides covariances of the average scattering angle cosine $\overline{μ}$ for 117 nuclides, and the NJOY nuclear data processing code can process these covariances. To propagate them to the uncertainty of a response due to uncertainties in $\overline{μ}$ requires sensitivities of the response with respect to $\overline{μ}$. This report uses NJOY’s definition of $\overline{μ}$ and derives a relationship between the sensitivity of a response to $\overline{μ}$ and the sensitivity of the response to the anisotropic P 1 elastic scattering cross section. This relationship has been implemented in the multigroup neutron sensitivity code SENSMG, which relies on the discrete ordinates code PARTISN for the neutron transport calculations. Kodeli has also calculated sensitivities with respect to $\overline{μ}$ and associated uncertainties. However, the formal derivation of this paper shows that the definition of the sensitivity with respect to $\overline{μ}$ is ambiguous because there are multiple ways to change $\overline{μ}$ by changing scattering cross sections. Here we derive two such ways. Both are implemented in SENSMG. This report is organized as follows. Section II discusses $\overline{μ}$ and its significance. Section III presents the equations used for the sensitivity of a response with respect to $\overline{μ}$, and Section IV discusses how to compute these sensitivities with SENSMG and MCNP6.3. Section V is a fast-spectrum test problem, and Sec. VI is a thermal-spectrum problem. Section VII is a summary. The example problem inputs are listed in the appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Preliminary Chlorine Worth Study Benchmark Evaluation

The Chlorine Worth Studies (CWS) experiments with polyvinyl chloride (PVC with chemical formula (C 2 H 3 Cl) n ), chlorinated polyvinyl chloride (CPVC with chemical formula (C 9 H 11 Cl 7 ) n ), and high density polyethylene (HDPE with chemical formula (CH 2 ) n ) were a series of measurements performed at the National Criticality Experiments Research Center (NCERC). The purpose of the CWS experiments was to perform integral experiments that were highly sensitive to the thermal 35 Cl(n,γ) reaction and matched the sensitivities of aqueous chloride operations at the plutonium facility at Los Alamos National Laboratory (LANL). The CWS experiments were performed on the Planet critical assembly machine at NCERC and utilized weapons grade plutonium (WGPu) plates as fuel. The design process and design of the CWS experiment were discussed previously. This paper discusses the benchmark evaluation of the experiment, intended for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). Criticality calculations were performed with MCNP version 6.3. Results presented here are preliminary, as the benchmark has not yet been submitted to the ICSBEP.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

IER 555: Godiva Benchmark Update CED-2 (Final Design Report)

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) evaluation of the Godiva IV critical assembly, HEU-MET-FAST-086: GODIVA-IV DELAYED-CRITICAL EXPERIMENTS (HMF-086), was completed by Russ Mosteller. Five critical experiment configurations performed at the Los Alamos National Laboratory (LANL) Technical Area (TA)-18 were evaluated as acceptable benchmark cases. The five cases consist of four delayed critical configurations which differ in control rod positions and one prompt critical configuration. All cases calculated a lower $k_{eff}$ than measured by experiment. This data is referred to as the TA-18 Godiva IV benchmark in this report. In 2005, Godiva IV was disassembled for relocation to the Nevada Test Site (NTS), now Nevada National Security Site (NNSS), at the National Criticality Experiments Research Center (NCERC). Following the disassembly and subsequent reassembly and startup of Godiva IV at NCERC, additional information about the Godiva IV components was obtained. An errata note was added to the HMF-086 evaluation in the ICSBEP handbook to provide this new information until a revision to the benchmark evaluation could be performed. In addition to those corrections, there are differences between the Godiva IV assembly at TA-18 and the Godiva IV assembly at NCERC. These differences include both assembly-specific differences (differences in the safety block gap, differences in the control rod positions, a new NCERC Top Hat and contamination shield) as well as environmental differences, such as the size and shape of the experimental building where the assembly is located. An additional model with similar cases, referred to as the NCERC Godiva IV benchmark in this report, will be added to the revised HMF-086 to capture these additional differences. This will provide the best benchmark model of Godiva for use by those performing experiments at NCERC. The IER 555 CED-2 report documents the information that will be updated in the HMF-086 revision, both the corrections to the TA-18 Godiva IV benchmark and the subsequent changes to create a NCERC Godiva IV benchmark. It describes the measurements that will be performed for cases similar to those performed at TA-18. It also describes measurements that will be included in the evaluation as additional data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

(U) SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , α , and Subcritical Multiplication Using PARTISN

SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff , α, and subcritical multiplication using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, a nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every nuclide and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations and outer boundaries. It computes sensitivities to user specified reactions whose cross sections are available in a user-supplied NJOY output file. The tool can be used for one-dimensional spherical and slab (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used for fixed-source and eigenvalue problems. For most responses, the tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. The tool is thus limited to computing sensitivities only for GPT-allowable responses. For subcritical multiplication, the tool implements sensitivities derived by O’Brien and Clark. SENSMG has a similar role as the old SWANLAKE (Ref. 8), FORSS (Ref. 9), and SENSIT (Ref. 10) codes. It has capabilities similar to those of SUSD3D (Refs. 11 and 12), which also uses PARTISN. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI presents some sample problems and discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists most of the input files whose results are presented in Sec. VI. Appendix B provides some useful details on one of the cross-section libraries that SENSMG supports.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Updated Godiva-IV Benchmark Preview

A note of errata prepended to the Godiva-IV delayed-critical benchmark (HEU-MET-FAST- 086) identifies two corrections to be made to the model: The glory hole in the spindle should be made larger, and the height of the safety block should be made smaller (and therefore its density made larger). In addition, the safety block at its full-in position is closer to the inner subassembly plate than was modeled in the benchmark. These changes have been made to HEU-MET-FAST-086 Case 4 in order to estimate the effect on $\kappa$ eff and on the neutron flux spectrum. Using smaller separation, a smaller safety block, and a larger glory hole caused keff to increase by 453 ± 1 pcm from the benchmark. The latest nuclear data, ENDF/B-VIII.0, have also been used, causing $\kappa$ eff to increase another 37 ± 1 pcm. Flux spectra were compared in a modeled fission foil (near the center of Godiva-IV in the glory hole) and at three external (point) detectors. Within the fission foil, using ENDF/B-VIII.0 induced changes in the flux spectrum similar in size to the changes due to using smaller separation, a smaller safety block, and a larger glory hole. At the detectors, using ENDF/B- VIII.0 induced changes in the flux spectrum much larger than those due to changing the model. In other words, the corrections to the benchmark model cause a large increase in $\kappa$ eff , but the changes to the neutron flux spectrum are small compared to those caused by using the latest nuclear data. This study presents a preview of results expected during the reevaluation of the Godiva IV benchmark, but it is not a substitute for the full reevaluation.A note of errata prepended to the Godiva-IV delayed-critical benchmark (HEU-MET-FAST- 086) identifies two corrections to be made to the model: The glory hole in the spindle should be made larger, and the height of the safety block should be made smaller (and therefore its density made larger). In addition, the safety block at its full-in position is closer to the inner subassembly plate than was modeled in the benchmark. These changes have been made to HEU-MET-FAST-086 Case 4 in order to estimate the effect on $\kappa$ eff and on the neutron flux spectrum. Using smaller separation, a smaller safety block, and a larger glory hole caused $\kappa$ eff to increase by 453 ± 1 pcm from the benchmark. The latest nuclear data, ENDF/B-VIII.0, have also been used, causing $\kappa$ eff to increase another 37 ± 1 pcm. Flux spectra were compared in a modeled fission foil (near the center of Godiva-IV in the glory hole) and at three external (point) detectors. Within the fission foil, using ENDF/B-VIII.0 induced changes in the flux spectrum similar in size to the changes due to using smaller separation, a smaller safety block, and a larger glory hole. At the detectors, using ENDF/B- VIII.0 induced changes in the flux spectrum much larger than those due to changing the model. In other words, the corrections to the benchmark model cause a large increase in $\kappa$ eff , but the changes to the neutron flux spectrum are small compared to those caused by using the latest nuclear data. This study presents a preview of results expected during the reevaluation of the Godiva IV benchmark, but it is not a substitute for the full reevaluation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

(U) Calculating First-Order Sensitivities to Material Density and Composition in Fixed-Source Problems Using the MCNP6 Perturbation Capability

The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first-order sensitivity analyses in fixed-source problems.2 Sensitivities of a response with respect to material and nuclide densities can be used for efficient first-order uncertainty quantification. This report discusses how to compute sensitivities of material densities and compositions in order to propagate uncertainties using the standard uncertainty propagation formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Rossi-alpha Analysis of CURIE Experiment

Critical assembly measurement and operations are crucial to the development of benchmark data to support research into criticality safety, radiation-detection development, and the overall application of nuclear technologies. Accurate nuclear data are needed for accurate predictive simulations, and validation using critical experiments is an important part of the nuclear data pipeline. The accuracy of nuclear data are improved using more robust and targeted measurements. In particular, the intermediate energy range of uranium is of great interest. LANL has successfully performed the Zeus series of experiments on the Comet assembly to investigate the intermediate energy range for HEU with various moderators. A successor to these experiments is the benchmark for the Critical Unresolved Region Integral Experiment (CURIE), designed to be sensitive to the unresolved resonance region (URR). This experiment is designed using polytetrafluoroethylene, more commonly known as Teflon, moderators. The CURIE experiment was successfully conducted at the National Criticality Experiments Research Center (NCERC). The Rossi-alpha method was used to evaluate the propensity of the CURIE configurations to sustain fission chains by estimating the prompt neutron decay constant α. This work evaluates the α at delayed critical using Rossi-alpha for different Teflon moderator thicknesses in the CURIE experiment. These results will help improve understanding of the CURIE benchmark experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

(U) Using the MCNP6 Perturbation Capability for Source Nuclide Density Sensitivities

The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first- and second-order sensitivity analyses in fixed-source problems. Unfortunately, the PERT capability does not yet work with point-detector (F5) or pulse-height (F8) tallies. Also, the second-order Taylor series term lacks the cross terms that would allow the calculation of mixed sensitivities. Finally, the PERT capability does not compute any sensitivities involving the source in a fixed-source problem. When neutrons or gamma rays are emitted due to radioactive decay, the source rate density depends intrinsically on the nuclide densities in the source material. When the PERT capability in MCNP6 is used to compute sensitivities of a response to nuclide and material densities, those sensitivities do not include the effect of the source.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

(U) Updated Godiva-IV Benchmark Preview

A note of errata prepended to the Godiva-IV delayed-critical benchmark (HEU-METFAST-086) identifies two corrections that need to be made to the model: The glory hole in the spindle should be made larger, and the height of the safety block should be made smaller (and therefore its density made larger). In addition, the safety block at its full-in position is closer to the inner subassembly plate than was modeled in the benchmark. These changes have been made to HEU-METFAST-086 Case 4 in order to estimate the effect on k eff and on the neutron flux spectrum. Using smaller separation, a smaller safety block, and a larger glory hole caused k eff to increase by 453 ± 1 pcm from the benchmark. The latest nuclear data, ENDF-B/VIII.0, have also been used; this caused k eff to increase another 37 ± 1 pcm. Flux spectra were compared in a modeled fission foil and at three external point detectors. Within the fission foil, using ENDF/B-VIII.0 for Godiva-IV induced changes in the flux spectrum similar in size to the changes due to using smaller separation, a smaller safety block, and a larger glory hole. At the point detectors, using ENDF/B-VIII.0 induced changes in the flux spectrum much larger than those due to changing the model. In other words, the corrections to the benchmark model cause a large increase in k eff , but the changes to the neutron flux spectrum are small compared to those caused by using the latest nuclear data. This study presents a preview of results expected during the reevaluation of the Godiva IV benchmark, but it is not a substitute for the full reevaluation

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗