Search NASA⌕ Search

DOE OSTI · 1997830

SCALE 6.2 Lattice Physics Performance Assessment

Abstract

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Mertyurek, Ugur, Jessee, Matthew A., Betzler, Benjamin R., Bowman, Stephen M.. 2023-03-01. SCALE 6.2 Lattice Physics Performance Assessment. https://doi.org/10.2172/1997830

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Design-informed neutronics assessment of long-lived fission product transmutation in a tokamak fusion reactor blanket

This study presents a neutronics-based assessment of the feasibility and viability of transmuting six major long-lived fission products (LLFPs) from light-water reactors, namely 99 Tc, 129 I, 79 Se, 93 Zr, 126 Sn, and 135 Cs, within the blanket region of a tokamak fusion reactor, using the MIT ARC design as a concrete fusion configuration. Monte Carlo neutronics simulations were performed to evaluate LLFP transmutation and to compare the results with a reference boiling water reactor (BWR). The results indicate that transmutation of all six LLFPs is neutronics-feasible in fusion reactors, with transmutation half-lives significantly shorter than their natural decay half-lives. For elemental targets, transmutation of 135 Cs, 126 Sn, and 93 Zr was found potentially viable, as the net mass transmuted exceeded that achievable in the reference BWR under identical target volume and irradiation time. When isotopically separated targets were considered, transmutation of 126 Sn and 93 Zr appeared potentially viable. A parametric study demonstrated that plasma geometry modifications can enhance local neutron flux, increasing the transmuted 93 Zr mass by approximately 33% and reducing the transmutation half-life from approximately 240 years to 180 years. Repositioning the target and adjusting material layer thickness reduced the transmutation half-life of 93 Zr to 67 years and increased the net mass transmuted by a factor of 50. Furthermore, these results demonstrate that fusion reactors can enable LLFP transmutation beyond the practical limits of thermal fission reactors and highlight the critical role of reactor and blanket design optimization. Engineering and fuel-cycle considerations required for deployment are beyond the scope of this neutronics-focused study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗