Dorothea Wiarda

Isotopic Depletion and Decay Methods and Analysis Capabilities in SCALE

Abstract The calculation of fuel isotopic compositions is essential to support design, safety analysis, and licensing of many components of the nuclear fuel cycle—from reactor physics and severe accident analysis to back-end fuel cycle issues, including spent-fuel storage and transportation, reprocessing, and radioactive waste management. Versions of the ORIGEN code, developed by Oak Ridge National Laboratory, have been used worldwide for isotopic depletion and decay analysis for more than three decades. The supported version of ORIGEN, maintained as the depletion analysis module for SCALE 6, performs detailed time-dependent isotopic generation and depletion for 1946 nuclides for reactor fuel and activation analysis. Stand-alone ORIGEN calculations can be performed using cross-section libraries developed for a wide range of reactor types and fuel designs used worldwide, including light water reactors UO2 and MOX, CANDU, VVER 440 and 1000, RBMK, and graphite reactors. Alternatively, within SCALE 6, ORIGEN can be automatically coupled to two-dimensional discrete ordinates or three-dimensional Monte Carlo transport solvers that provide problem-dependent cross sections for use in the ORIGEN depletion calculation. The hybrid ability to function as either a stand-alone or coupled depletion code provides ORIGEN advanced capabilities to simulate a broad range of applications for various reactor systems. The nuclear data libraries in ORIGEN have been significantly improved recently, using modern ENDF/B nuclear data evaluations. The most recent developments in SCALE 6.1 include the addition of ENDF/B-VII decay data, energy-dependent fission yields, and fine-group ORIGEN neutron cross sections based on the JEFF-3.0/A special purpose activation files. Advanced methods and data for neutron and gamma source energy spectral analysis are also available in the current version of the code. The ORIGEN code and associated nuclear data libraries have been extensively validated against experimental data that include spent nuclear fuel isotopic assay data for actinides and fission products, radiation source spectra, and decay heat measurements.

Sensitivity and Uncertainty Analysis Capabilities and Data in SCALE

Abstract In SCALE 6, the Tools for Sensitivity and UNcertainty Analysis Methodology Implementation (TSUNAMI) modules calculate the sensitivity of keff or reactivity differences to the neutron cross-section data on an energy-dependent, nuclide-reaction-specific basis. These sensitivity data are useful for uncertainty quantification, using the comprehensive neutron cross-section-covariance data in SCALE 6. Additional modules in SCALE 6 use the sensitivity and uncertainty data to produce correlation coefficients and other relational parameters that quantify the similarity of benchmark experiments to application systems for code validation purposes. Bias and bias uncertainties are quantified using parametric trending analysis or data adjustment techniques, providing detailed assessments of sources of biases and their uncertainties and quantifying gaps in experimental data available for validation. An example application of these methods is presented for a generic burnup credit cask model.

A Statistical Sampling Method for Uncertainty Analysis with SCALE and XSUSA

A new statistical sampling sequence called Sampler has been developed for the SCALE code system. Random values for the input multigroup cross sections are determined by using the XSUSA program to sample uncertainty data provided in the SCALE covariance library. Using these samples, Sampler computes perturbed self-shielded cross sections and propagates the perturbed nuclear data through any specified SCALE analysis sequence, including those for criticality safety, lattice physics with depletion, and shielding calculations. Statistical analysis of the output distributions provides uncertainties and correlations in the desired responses, due to nuclear data uncertainties. The Sampler/XSUSA methodology is described, and example applications are shown for criticality safety and spent-fuel analysis.

Generation of a Broad-Group HTGR Library for Use with SCALE

With current and ongoing interest in high temperature gas reactors (HTGRs), the U.S. Nuclear Regulatory Commission (NRC) anticipates the need for nuclear data libraries appropriate for use in applications for modeling, assessing, and analyzing HTGR reactor physics and operating behavior. The objective of this work was to develop a broad-group library suitable for production analyses with SCALE for HTGR applications. Several interim libraries were generated from SCALE fine-group 238- and 999-group libraries, and the final broad-group library was created from Evaluated Nuclear Data File/B Version ENDF/B-VII Release 0 cross-section evaluations using new ORNL methodologies with AMPX, SCALE, and other codes. Furthermore, intermediate resonance (IR) methods were applied to the HTGR broadgroup library, and lambda factors and f-factors were incorporated into the library s nuclear data files. A new version of the SCALE BONAMI module named BONAMI-IR was developed to process the IR data in the new library and, thus, eliminate the need for the CENTRM/PMC modules for resonance selfshielding. This report documents the development of the HTGR broad-group nuclear data library and the results of test and benchmark calculations using the new library with SCALE. The 81-group library is shown to model HTGR cases with similar accuracy to the SCALE 238-group library but with significantly faster computational times due to the reduced number of energy groups and the use of BONAMI-IR instead of BONAMI/CENTRM/PMC for resonance self-shielding calculations.

A New Method for Generating Probability Tables in the Unresolved Resonance Region

Abstract A new method for constructing probability tables in the unresolved resonance region (URR) has been developed. This new methodology is an extensive modification of the single-level Breit-Wigner (SLBW) pseudo-resonance pair sequence method commonly used to generate probability tables in the URR. The new method uses a Monte Carlo process to generate many pseudo-resonance sequences by first sampling the average resonance parameter data in the URR and then converting the sampled resonance parameters to the more robust R-matrix limited (RML) format. For each sampled set of pseudo-resonance sequences, the temperature-dependent cross sections are reconstructed on a small grid around the energy of reference using the Reich-Moore formalism and the Leal-Hwang Doppler broadening methodology. The effective cross sections calculated at the energies of reference are then used to construct probability tables in the URR. The RML cross-section reconstruction algorithm has been rigorously tested for a variety of isotopes, including 16O, 19F, 35Cl, 56Fe, 63Cu, and 65Cu. The new URR method also produced normalized cross-section factor probability tables for 238U that were found to be in agreement with current standards. The modified 238U probability tables were shown to produce results in excellent agreement with several standard benchmarks, including the IEU-MET-FAST-007 (BIG TEN), IEU-MET-FAST-003, and IEU-COMP-FAST-004 benchmarks.

High-Resolution Neutron Capture and Total Cross-Section Measurements, and the Astrophysical 95Mo(n,gamma) Reaction Rate at s-process Temperatures

Abundances of Mo isotopes predicted by stellar models of the s process are, except for {sup 95}Mo, in good agreement with data from single grains of mainstream presolar SiC. Because the meteorite data seemed sound and no reasonable modification to stellar theory resulted in good agreement for {sup 95}Mo, it has been suggested that the recommended neutron capture reaction rate for this nuclide is 30% too low. Therefore, we have made a new determination of the {sup 95}Mo(n,{gamma}) reaction rate via high-resolution measurements of the neutron-capture and total cross sections of {sup 95}Mo at the Oak Ridge Electron Linear Accelerator. These data were analyzed with the R-matrix code SAMMY to obtain parameters for resonances up to E{sub n} = 10 keV. Also, a small change to our capture apparatus allowed us to employ a new technique to vastly improve resonance spin and parity assignments. These new resonance parameters, together with our data in the unresolved range, were used to calculate the {sup 95}Mo(n,{gamma}) reaction rate at s-process temperatures. We compare the currently recommended rate to our new results and discuss their astrophysical impact.

New Neutron-Induced Cross Section Measurements for Weak s Process Studies

A series of new neutron capture and transmission measurements has been undertaken at the Oak Ridge Electron Linear Accelerator (ORELA) in response to deficiencies identified in nuclear data libraries of crucial importance to the Nuclear Criticality Safety Program. New data and evaluations including covariances are required for several stable fission products as well as for materials found in mixtures with uranium. For example, chromium and nickel as constituents of stainless steel perform poorly in criticality calculations due to their relatively large neutron cross sections and substantial uncertainties in previous measurements. Therefore, new neutron-capture and total cross-section measurements are needed for 52,53Cr and 58,60Ni. These newly obtained data can be used not only to improve criticality calculations but also to serve as input parameters for the weak s-process stellar model calculations in massive stars. We will report on new experiments for these nuclides.