Search

Delayed neutrons, which are described by kinetic parameters, are significant for nuclear reactor operation as they make nuclear reactors controllable. Modern lattice physics codes such as TRITON and Polaris in SCALE code system, generate kinetic parameters for a given material composition. The calculation is performed by summation of the delayed neutron data for each precursor isotope and then weighting is performed by real and/or adjoint neutron spectrum. Quantifying uncertainties in the weighted kinetic parameters is important for assembly and core calculations. Understanding uncertainty in modeling scenarios involve kinetic parameters (e.g. transients), requires propagating uncertainty in the weighted kinetic parameters due to uncertainties in the fundamental nuclear data libraries, including delayed neutron data. In this work, uncertainty analysis of the homogenized (also called weighted or macroscopic) kinetic parameters has been performed using Sampler, a module in SCALE code system, to investigate the effect of fundamental nuclear and delayed neutron data uncertainties on the weighted kinetic parameters. Two major sources of uncertainties were considered: (1) fundamental nuclear data (i.e. cross-sections, fission yield, decay data) and (2) nuclide-dependent group-wise delayed neutron data based on reported experimental measurements. In this study, a new capability was developed through SCALE code system to allow propagation of delayed neutron data uncertainties. Preliminary analysis demonstrated 7% uncertainty in β eff at Beginning of Life (BOL) and increased to 15% after fuel burnup. Decay constant groups showed lower uncertainty than delayed neutron fraction groups, both at BOL and End of Life (EOL). Delayed neutron fraction responses showed high correlation to each other which is expected to be due to the cross-section covariances reported in SCALE data libraries as well as the normalization condition of the nuclide-dependent DNF. Different sources of U-235 delayed neutron data were compared, and the results showed that the uncertainty calculated by Tuttle data was bounded by other sources. The current study can be extended to calculate kinetic parameters and their uncertainties for more advanced LWR applications.

Nuclear data are a major source of uncertainties in reactor physics calculations. The propagation of nuclear data uncertainties to important system responses is instrumental when determining appropriate safety margins in reactor safety analyses. It is also important to understand the major contributors to the observed uncertainties to make recommendations for further measurements and evaluations and aid in the understanding of the studied system. The SCALE code system allows for nuclear data uncertainty analysis based on the random sampling approach as implemented in SCALE’s Sampler sequence. Sampler was recently extended by a sensitivity analysis in terms of the calculation of two correlation-based sensitivity indices. This analysis allows for the identification of the top contributing nuclear reactions to any analyzed output uncertainty. This paper presents the sensitivity indices, along with their interpretation and limitations. It demonstrates the application in an eigenvalue and decay heat analysis for a boiling water reactor fuel assembly.

A new depletion capability has been developed in the Exnihilo radiation transport code suite. This capability enables massively parallel domain-decomposed coupling between the Shift continuous-energy Monte Carlo solver and the nuclide depletion solvers in ORIGEN to perform high-performance Monte Carlo depletion calculations. This paper describes this new depletion capability and discusses its various features, including a multi-level parallel decomposition, high-order transport-depletion coupling, and energy-integrated power renormalization. Several test problems are presented to validate the new capability against other Monte Carlo depletion codes, and the parallel performance of the new capability is analyzed.

Polaris is a 2-dimensional multigroup lattice physics capability in the SCALE code system for the analysis of light water reactor fuel designs. The goal of light water reactor lattice physics codes is to generate few-group homogenized cross sections for downstream full-core nodal diffusion calculations. Lattice physics calculations contain three primary components: the cross section processing calculation, the 2D transport calculation, and the depletion calculation. This paper summarizes the calculational methods and their implementation into Polaris, with an emphasis on implementation of the embedded self-shielding method. The accuracy of the embedded self-shielding method depends on the procedure used to generate self-shielding factors on the multigroup library. Numerical benchmarks calculations reveal that the accuracy of Polaris eigenvalue predictions is enhanced by (1) using heterogeneous unit cell models to generate the self-shielding factors on the library and (2) using self-shielding factors for within-group scattering cross section.

To more accurately represent the desired power in a constant power burnup calculation, the depletion steps of the calculation can be divided into substeps and the neutron flux renormalized on each substep to match the desired power. This paper explores how such renormalization should be performed, how large a difference it makes, and whether using renormalization affects results regarding the relative performance of different neutronics–depletion coupling schemes. When used with older coupling schemes, renormalization can provide a considerable improvement in overall accuracy. With previously published higher order coupling schemes, which are more accurate to begin with, renormalization has a much smaller effect. While renormalization narrows the differences in the accuracies of different coupling schemes, their order of accuracy is not affected.

The SCALE code system’s ability to address stochastic distributions of fuel particles within a graphite matrix has been revisited in both multigroup (MG) features and continuous-energy (CE) Monte Carlo methods. This paper describes current and emergent SCALE capabilities within the CSAS sequence to address double-heterogeneous systems and presents verification and validation studies of these methods and data. Good agreement was obtained for a high temperature gas-cooled reactor (HTGR) fuel pebble model between CSAS MG eigenvalue calculations and corresponding CE reference solutions. Code-to-code comparisons for this HTGR pebble model showed good agreement of CSAS-KENO and CSAS-Shift CE calculations and the Serpent and MCNP codes in terms of eigenvalues and reaction rate ratios. Validation studies based on two HTGR experiments resulted in good agreement between MG and CE results, as well as between experiment and calculation, although the level of agreement was significantly influenced by the applied ENDF/B nuclear data library.

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. The transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

A method for including external feed with polynomial time dependence in depletion calculations with the Chebyshev Rational Approximation Method (CRAM) is presented and the implementation of CRAM to the ORIGEN module of the SCALE suite is described. In addition to being able to handle time-dependent feed rates, the new solver also adds the capability to perform adjoint calculations. Results obtained with the new CRAM solver and the original depletion solver of ORIGEN are compared to high precision reference calculations, which shows the new solver to be orders of magnitude more accurate. Furthermore, in most cases, the new solver is up to several times faster due to not requiring similar substepping as the original one.

We present a quasidiffusion (QD) method for solving neutral particle transport problems in Cartesian XY geometry on unstructured quadrilateral meshes, including local refinement capability. Neutral particle transport problems are central to many applications including nuclear reactor design, radiation safety, astrophysics, medical imaging, radiotherapy, nuclear fuel transport/storage, shielding design, and oil well-logging. The primary development is a new discretization of the low-order QD (LOQD) equations based on cell-local finite differences. The accuracy of the LOQD equations depends on proper calculation of special non-linear QD (Eddington) factors from a transport solution. In order to completely define the new QD method, a proper discretization of the transport problem is also presented. The transport equation is discretized by a conservative method of short characteristics with a novel linear approximation of the scattering source term and monotonic, parabolic representation of the angular flux on incoming faces. Analytic and numerical tests are used to test the accuracy and spatial convergence of the non-linear method. All tests exhibit O(h^2) convergence of the scalar flux on orthogonal, random, and multi-level meshes.

Nuclear data uncertainty propagation based on stochastic sampling ( SS) is becoming more attractive while leveraging modern computer power. Two variants of the SS approach are compared in this paper. The Total Monte Carlo (TMC) method by the Nuclear Research and Consultancy Group (NRG) generates perturbed ENDF-6-formatted nuclear data by varying nuclear reaction model parameters. At Paul Scherrer Institute (PSI) the Nuclear data Uncertainty Stochastic Sampling (NUSS) system generates perturbed ACE-formatted nuclear data files by applying multigroup nuclear data covariances onto pointwise ACE-formatted nuclear data. Uncertainties of Pu-239 and U-235 from ENDF/B-VII.1, ZZ-SCALE6/COVA-44G and TENDL covariance libraries are considered in NUSS and propagated in MCNPX calculations for well-studied Jezebel and Godiva fast spectrum critical benchmarks. The corresponding uncertainty results obtained by TMC are compared with NUSS results and the deterministic Sensitivity/Uncertainty method of TSUNAMI-3D from SCALE6 package is also applied to serve as a separate verification. The discrepancies in the propagated Pu-239 and U-235 uncertainties due to method and covariance differences are discussed.