Engineering Papers⌕ Search

Engineering topics

Clarno, Kevin T.

Publications and source records attributed to Clarno, Kevin T..

High-level fuel fabrication facility designs from discrete-event simulation

Like other industrial processes, the production of metallic nuclear fuels (MNF) requires that fabrication facilities be able to reliably meet production demands, operate efficiently, and adhere to federal and local safety regulations. In turn, the set of design variables employed by such a facility, such as operations policies, infrastructure and machinery purchased, and the type and number of staff hired, directly impact a facility’s ability to satisfy these goals. Therefore, facility designers must carefully determine which set of design variable values optimally satisfies these constraints. In this paper, we explore how values for these high-level design variables, namely hiring requirements, can be determined in the context of nuclear fuel manufacturing through the coupling of physics-based and discrete-event simulation technologies. Using the Versatile Test Reactor (VTR) program as a case study, we demonstrate how SCALE and MCNP nuclear physics model outputs can be integrated into ExtendSim discrete-event simulation (DES) models of the fuel fabrication process to determine the optimal number of staff hired to ensure fuel production goals are met, operations comply with effective dose limit regulations, and overall project costs are reduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermochemically-informed mass transport model for interdiffusion of U and Zr in irradiated U-Pu-Zr fuel with fission products

A new formulation for redistribution of constituents in metallic U-Pu-Zr nuclear fuel is presented that can incorporate the contributions from fission products. The formulation is based on the thermodynamic driving forces derived from the generalized chemical potential that includes effects of composition and temperature. As a result, the redistribution model can readily account for the composition changes due to the generation of fission products while using only a limited set of transport coefficients. The thermodynamic model for the metallic fuel and thermochemistry solver Thermochimica were coupled with the nuclear fuel performance code BISON to implement the redistribution model. The simulations reproduce an experimentally observed Zr-depletion zone in the mid-radius region of an irradiated fuel slug. The generation of fission products during burnup is shown to have a stabilizing effect on the fuel chemistry, slowing the rate of U-Zr interdiffusion and reducing the size of the Zr-depleted zone. Here, the overall objective of the work is to advance predictive capabilities of fission product behaviour in U-Zr metallic fuel in the context of fuel performance and safety.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and simulation of oxygen transport in high burnup LWR fuel

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. Here, 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.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of the multigroup cross section library for the CASL neutronics simulator MPACT: Method and procedure

The MPACT neutronics module of the Consortium for Advanced Simulation of Light Water Reactors (CASL) core simulator is a 3D whole core transport code being developed for the CASL toolset known as the Virtual Environment for Reactor Analysis (VERA). Key characteristics of the MPACT code include a subgroup method for resonance self-shielding and a whole-core transport solver with a 2D/1D synthesis method. The AMPX and SCALE code packages have been significantly improved to include intermediate resonance parameters, as well as new modules to enhance the accuracy of resonance data and scattering matrices for the Bondarenko approach and transport cross sections of hydrogen. In addition, CASL XSTools have been developed to generate the MPACT multigroup cross section library. A simple super-homogenization method for 238U has been developed to resolve reaction rate discrepancy issues caused by angle-dependent total cross sections, poor scattering matrices, and a poor resonance interference model at coarse group structure. Here, due to the limitation of computing capacity that can occur even in high performance computing, a new 51-group structure has been developed for efficient simulation that is applicable for both pressurized water reactor (PWR) and boiling water reactor (BWR) simulations. A new ENDF/B-VII.1 MPACT 51-group library has been developed using the AMPX/SCALE and XSTools code packages, and it has been verified and validated by performing various benchmark calculations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Implementation of a Spacer Grid Rod Thermal-Hydraulic Reconstruction (ROTHCON) Capability into the Thermal-Hydraulic Subchannel Code CTF

The Consortium for Advanced Simulation of Light Water Reactors is developing a core simulator capability known as the Virtual Environment for Reactor Applications (VERA) to address nuclear industry challenge problems such as crud-induced power shift (CIPS). The CTF thermal-hydraulic (T/H) subchannel code provides thermal feedback in the coupled neutronics, T/H, crud chemistry simulation that VERA performs. It has been discovered that the coarse meshing approach used by CTF (in which fuel rods are discretized into four azimuthal segments) can be a source of error in predicting crud growth and boron distribution in VERA CIPS calculations. Spacer grid effects lead to complex rod-to-fluid heat transfer behavior that, when not resolved, can lead to error in the prediction of crud growth and boron deposition. A higher-fidelity computational fluid dynamics approach can be used instead of CTF, but this leads to excessive simulation times. This paper presents an approach for using high-fidelity computational fluid dynamics data to create shape functions that are used in CTF to reconstruct rod surface heat transfer behavior as a function of spacer grid geometry. The approach is demonstrated for a 5 × 5 rod bundle facility with five mixing vane grids under a range of operating conditions encountered in nominal pressurized water reactor conditions. Furthermore, it is demonstrated that the grid heat transfer maps are successful at introducing a higher-fidelity heat transfer modeling capability into CTF.

CFD-informed↗

Depletion, Chemical Reaction and Transport in High Burnup Nuclear Fuel

We have developed a formulation and computational model for oxygen transport in Light Water Reactor (LWR) uranium dioxide fuel. The overall model couples the burnup simulation isotopic composition with a thermochemistry model of the fuel phase, and oxygen transport using the driving forces from the thermodynamic calculations. The diffused oxygen is accounted for in the thermochemistry model in order to establish consistent material and thermodynamic conditions in the fuel undergoing burnup. The model has been implemented in the nuclear fuel performance code Bison. The formulation and the model have been successfully demonstrated on fuel burnup data from the open literature. The developed capability enables consideration of complex chemical composition of irradiated fuels beyond available burnup models in Bison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ONRAMP, ORNL Nuclear Resources Analysis and Modeling Portfolio

Many teams from US and international organizations, government, industry, research organizations, and academic institutions have established near-term technical exchanges and long-term partnerships with ORNL, realizing previously unforeseen innovations and efficiencies. Through ONRAMP, we invite you to learn more about the rich resources available at ORNL that can enhance your current approaches. Learn what is possible using customized, advanced techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new fuel modeling capability, CTFFuel, with a case study on the fuel thermal conductivity degradation

Here, a new fuel modeling capability, CTFFuel, is developed from the subchannel code, CTF. This code is a standalone interface to the CTF fuel rod models, allowing for fuel rod simulations to be run independently from the fluid. This paper provides an overview of the code with a case study on the thermal conductivity degradation of LWR fuels to demonstrate its capabilities. The modeling of fuel thermal conductivity degradation in the code is improved through the addition of new modeling options to account for the irradiation effects via globally defined parameters. After the initial implementation, a variety of order-of-accuracy tests and code comparisons are performed to test software quality. A controlled analysis is allowed by CTFFuel to verify the numerical scheme of CTF’s conduction solution and to benchmark its fuel temperature predictions against FRAPCON-4.0’s. Overall, the software quality and verification procedure ensures that the new model is coded correctly, that it properly interacts with the rest of the code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Neutron Capture Energies for Flux Normalization and Approximate Model for Gamma-Smeared Power

The Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) neutronics simulator MPACT has used a single recoverable fission energy for each fissionable nuclide assuming that all recoverable energies come only from fission reaction, for which capture energy is merged with fission energy. This approach includes approximations and requires improvement by separating capture energy from the merged effective recoverable energy. This report documents the procedure to generate recoverable neutron capture energies and the development of a program called CapKappa to generate capture energies. Recoverable neutron capture energies have been generated by using CapKappa with the evaluated nuclear data file (ENDF)/B-7.0 and 7.1 cross section and decay libraries. The new capture kappas were compared to the current SCALE-6.2 and the CASMO-5 capture kappas. These new capture kappas have been incorporated into the Simplified AMPX 51- and 252-group libraries, and they can be used for the AMPX multigroup (MG) libraries and the SCALE code package. The CASL VERA neutronics simulator MPACT does not include a gamma transport capability, which limits it to explicitly estimating local energy deposition from fission, neutron, and gamma slowing down and capture. Since the mean free path of gamma rays is typically much longer than that for the neutron, and the total gamma energy is about 10% to the total energy, the gamma-smeared power distribution is different from the fission power distribution. Explicit local energy deposition through neutron and gamma transport calculation is significantly important in multi-physics whole core simulation with thermal-hydraulic feedback. Therefore, the gamma transport capability should be incorporated into the CASL neutronics simulator MPACT. However, this task will be timeconsuming in developing the neutron induced gamma production and gamma cross section libraries. This study is to investigate an approximate model to estimate gammasmeared power distribution without performing any gamma transport calculation. A simple approximate gamma smearing model has been investigated based on the facts that pinwise gamma energy depositions are almost flat over a fuel assembly, and assembly-wise gamma energy deposition is proportional to kappa-fission energy deposition. The approximate gamma smearing model works well for single assembly cases, and can partly improve the gamma smeared power distribution for the whole core model. Although the power distributions can be improved by the approximate gamma smearing model, still there is an issue to explicitly obtain local energy deposition. A new simple approach or gamma transport/diffusion capability may need to be incorporated into MPACT to estimate local energy deposition for more robust multi-physics simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment of the analysis capability for core-wide PWR pellet-clad interaction screening of Watts Bar Unit 1

This report illustrates the application of the Virtual Environment for Reactor Applications simulation environment to pressurized water reactor core-wide pellet-clad interaction screening, with emphasis on the multi-dimensional analysis capabilities of the BISON fuel performance code. The Virtual Environment for Reactor Applications codes utilized in this effort include the BISON fuel performance code and the VERA-CS sub-element, which combines the MPACT neutron transport code and the COBRA-TF thermal-hydraulics code. The target reactor for this work is Watts Bar Unit 1, cycle 6 and cycle 7, selected on the basis that pellet-clad interaction failures had occurred during power start up in both cycles, despite the moderate power levels. These instances of failure at Watts Bar Unit 1 offered an interesting test of the code’s veracity and a particularly challenging condition for this assessment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Verification and Validation of the ENDF/B-VII.1 v4.3m1 MPACT 51-Group Cross Section Library

The MPACT neutronics module of the Consortium for Advanced Simulation of Light Water Reactors (CASL) core simulator is a 3-D whole core transport code being developed for the Virtual Environment for Reactor Analysis (VERA) CASL toolset. Key characteristics of the MPACT code include (1) a subgroup method for resonance selfshielding and (2) a whole-core transport solver with a 2-D/1-D synthesis method. The MPACT code requires a cross section library to support all of its core simulation capabilities; a cross section library would be MPACT’s most influencing component for simulation accuracy. Due to the limitation of computing capacity even in high performance computing, the cross section library must be developed to enhance computational efficiency in memory and computing time without losing accuracy and generality. This requires a very coarse energy group structure with about 50 groups. Multigroup (MG) neutron cross section libraries for the CASL VERA-CS neutronics simulator MPACT have been created by using the AMPX/SCALE code package developed at Oak Ridge National Laboratory (ORNL). A new 51-group structure was developed for efficient simulation to be applicable for both pressurized water reactor (PWR) and boiling water reactor (BWR) simulations. New v4.3m1 MPACT 51-group library with the ENDF/B- VII.1 nuclear data were developed for MPACT. Since SCALE KENO has some issues in the probability table and has a low cutoff energy for the free gas model, this issue has been corrected. This study focused on the development of the ENDF/B-VII.1 v4.3m1 MPACT MG library, including verification and validation. This document includes the library generation methodology and procedure, the verification procedure, and benchmark results compared to continuous energy Monte Carlo results using newly developed cross section libraries.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of Ultra-Fine Multigroup Cross Section Library of the AMPX/SCALE Code Packages

The Consortium for Advanced Simulation of Light Water Reactors Virtual Environment for Reactor Applications (VERA) neutronic simulator MPACT is being developed by Oak Ridge National Laboratory and the University of Michigan for various reactor applications. The MPACT and simplified MPACT 51- and 252-group cross section libraries have been developed for the MPACT neutron transport calculations by using the AMPX and Standardized Computer Analyses for Licensing Evaluations (SCALE) code packages developed at Oak Ridge National Laboratory. It has been noted that the conventional AMPX/SCALE procedure has limited applications for fast-spectrum systems such as boiling water reactor (BWR) fuels with very high void fractions and fast reactor fuels because of its poor accuracy in unresolved and fast energy regions. This lack of accuracy can introduce additional error sources to MPACT calculations, which is already limited by the Bondarenko approach for resolved resonance self-shielding calculation. To enhance the prediction accuracy of MPACT for fast-spectrum reactor analyses, the accuracy of the AMPX/SCALE code packages should be improved first. The purpose of this study is to identify the major problems of the AMPX/SCALE procedure in generating fast-spectrum cross sections and to devise ways to improve the accuracy. For this, various benchmark problems including a typical pressurized water reactor fuel, BWR fuels with various void fractions, and several fast reactor fuels were analyzed using the AMPX 252-group libraries. Isotopic reaction rates were determined by SCALE multigroup (MG) calculations and compared with continuous energy (CE) Monte Carlo calculation results. This reaction rate analysis revealed three main contributors to the observed differences in reactivity and reaction rates: (1) the limitation of the Bondarenko approach in coarse energy group structure, (2) the normalization issue of probability tables, and (3) neglect of the self-shielding effect of resonance-like cross sections at high energy range such as (n,p) cross section of Cl35. The first error source can be eliminated by an ultra-fine group (UFG) structure in which the broad scattering resonances of intermediate-weight nuclides can be represented accurately by a piecewise constant function. A UFG AMPX library was generated with modified probability tables and tested against various benchmark problems. The reactivity and reaction rates determined with the new UFG AMPX library agreed very well with respect to Monte Carlo Neutral Particle (MCNP) results. To enhance the lattice calculation accuracy without significantly increasing the computational time, performing the UFG lattice calculation in two steps was proposed. In the first step, a UFG slowing-down calculation is performed for the corresponding homogenized composition, and UFG cross sections are collapsed into an intermediate group structure. In the second step, the lattice calculation is performed for the intermediate group level using the condensed group cross sections. A preliminary test showed that the condensed library reproduces the results obtained with the UFG cross section library. This result suggests that the proposed two-step lattice calculation approach is a promising option to enhance the applicability of the AMPX/SCALE system to fast system analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Review of Bison Usage in VERA

Bison is being used in VERA in a variety of ways; this milestone will document an independent review of the current usage of Bison and provide guidance that will improve the accuracy, performance, and consistency in the ways that Bison is used. This task will entail running a suite of small, single and multi-cycle problems with VERA-CS, followed by Bison, and Tiamat (inline) and evaluating the usage. It will also entail performing several detailed ramp to full power solutions to compare the one-way coupled solver with the fully-coupled Tiamat. This will include at least one iteration with the PHI team to incorporate some of the feedback and improve the usage. This work will also be completed in conjunction with an FMC task to evaluate the ability of Bison to model load-follow in a PWR

97 MATHEMATICS AND COMPUTING↗

Demonstration of Coupled Tiamat Single Assembly Calculations

This report corresponds to milestone L3:PHI.PCI.P15.03, which was originally intended to investigate the time discretization approaches with the newly developed fully coupled Tiamat capability, targeting single assembly problems.

97 MATHEMATICS AND COMPUTING↗

Software Quality Assurance and Verification for the MPACT Library Generation Process

This report fulfills the requirements for the Consortium for the Advanced Simulation of Light-Water Reactors (CASL) milestone L2:RTM.P14.02, “SQA and Verification for MPACT Library Generation,” by documenting the current status of the software quality, verification, and acceptance testing of nuclear data libraries for MPACT. It provides a brief overview of the library generation process, from general-purpose evaluated nuclear data files (ENDF/B) to a problem-dependent cross section library for modeling of light-water reactors (LWRs). The software quality assurance (SQA) programs associated with each of the software used to generate the nuclear data libraries are discussed; specific tests within the SCALE/AMPX and VERA/XSTools repositories are described. The methods and associated tests to verify the quality of the library during the generation process are described in detail. The library generation process has been automated to a degree to (1) ensure that it can be run without user intervention and (2) to ensure that the library can be reproduced. Finally, the acceptance testing process that will be performed by representatives from the Radiation Transport Methods (RTM) Focus Area prior to the production library’s release is described in detail.

97 MATHEMATICS AND COMPUTING↗

Development of the V4.2m5 and V5.0m0 Multigroup Cross Section Libraries for MPACT for PWR and BWR

The MPACT neutronics module of the Consortium for Advanced Simulation of Light Water Reactors (CASL) core simulator is a 3-D whole core transport code being developed for the CASL toolset, Virtual Environment for Reactor Analysis (VERA). Key characteristics of the MPACT code include (1) a subgroup method for resonance selfshielding and (2) a whole-core transport solver with a 2-D/1-D synthesis method. The MPACT code requires a cross section library to support all the MPACT core simulation capabilities which would be the most influencing component for simulation accuracy.

97 MATHEMATICS AND COMPUTING↗

MPACT Subgroup Self-Shielding Efficiency Improvements

Recent developments to improve the efficiency of the MOC solvers in MPACT have yielded effective kernels that loop over several energy groups at once, rather that looping over one group at a time. These kernels have produced roughly a 2x speedup on the MOC sweeping time during eigenvalue calculation. However, the self-shielding subgroup calculation had not been reevaluated to take advantage of these new kernels, which typically requires substantial solve time. The improvements covered in this report start by integrating the multigroup kernel concepts into the subgroup calculation, which are then used as the basis for further extensions. The next improvement that is covered is what is currently being termed as “Lumped Parameter MOC”. Because the subgroup calculation is a purely fixed source problem and multiple sweeps are performed only to update the boundary angular fluxes, the sweep procedure can be condensed to allow for the instantaneous propagation of the flux across a spatial domain, without the need to sweep along all segments in a ray. Once the boundary angular fluxes are considered to be converged, an additional sweep that will tally the scalar flux is completed. The last improvement that is investigated is the possible reduction of the number of azimuthal angles per octant in the shielding sweep. Typically 16 azimuthal angles per octant are used for self-shielding and eigenvalue calculations, but it is possible that the self-shielding sweeps are less sensitive to the number of angles than the full eigenvalue calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗