Kang Seog Kim

High-fidelity light water reactor analysis with the numerical nuclear reactor

Abstract The Numerical Nuclear Reactor (NNR) was developed to provide a high-fidelity tool for light water reactor analysis based on first-principles models. High fidelity is accomplished by integrating full physics, highly refined solution modules for the coupled neutronic and thermal-hydraulic phenomena. Each solution module employs methods and models that are formulated faithfully to the first principles governing the physics, real geometry, and constituents. Specifically, the critical analysis elements that are incorporated in the coupled code capability are a direct whole-core neutron transport solution and an ultra-fine-mesh computational fluid dynamics/heat transfer solution, each obtained with explicit (sub-fuel-pin-cell level) heterogeneous representations of the components of the core. The considerable computational resources required for such highly refined modeling are addressed by using massively parallel computers, which together with the coupled codes constitute the NNR. To establish confidence in the NNR methodology, verification and validation of the solution modules have been performed and are continuing for both the neutronic module and the thermal-hydraulic module for single-phase and two-phase boiling conditions under prototypical pressurized water reactor and boiling water reactor conditions. This paper describes the features of the NNR and validation of each module and provides the results of several coupled code calculations.

A Full-Core Resonance Self-Shielding Method Using a Continuous-Energy Quasi–One-Dimensional Slowing-Down Solution that Accounts for Temperature-Dependent Fuel Subregions and Resonance Interference

Abstract A correction-based resonance self-shielding method is developed that allows annular subdivision of the fuel rod. The method performs the conventional iteration of the embedded self-shielding method (ESSM) without subdivision of the fuel to capture the interpin shielding effect. The resultant self-shielded cross sections are modified by correction factors incorporating the intrapin effects of radial variation of the shielded cross section, radial temperature distribution, and resonance interference. A quasi–one-dimensional slowing-down equation is developed to calculate such correction factors. The method is implemented in the DeCART code and compared with the conventional ESSM and subgroup method with benchmark MCNP results. The new method yields substantially improved results for both spatially dependent reaction rates and eigenvalues for typical pressurized water reactor pin cell cases with uniform and nonuniform fuel temperature profiles. The new method is also proved effective in treating assembly heterogeneity and complex material composition such as mixed oxide fuel, where resonance interference is much more intense.

VERA Core Simulator methodology for pressurized water reactor cycle depletion

This paper describes the methodology developed and implemented in the Virtual Environment for Reactor Applications Core Simulator (VERA-CS) to perform high-fidelity, pressurized water reactor (PWR), multicycle, core physics calculations. Depletion of the core with pin-resolved power and nuclide detail is a significant advance in the state of the art for reactor analysis, providing the level of detail necessary to address the problems of the U.S. Department of Energy Nuclear Reactor Simulation Hub, the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS has three main components: the neutronics solver MPACT, the thermal-hydraulic (T-H) solver COBRA-TF (CTF), and the nuclide transmutation solver ORIGEN. This paper focuses on MPACT and provides an overview of the resonance self-shielding methods, macroscopic-cross-section calculation, two-dimensional/one-dimensional (2-D/1-D) transport, nuclide depletion, T-H feedback, and other supporting methods representing a minimal set of the capabilities needed to simulate high-fidelity models of a commercial nuclear reactor. Results are presented from the simulation of a model of the first cycle of Watts Bar Unit 1. The simulation is within 16u2009parts per million boron (ppmB) reactivity for all state points compared to cycle measurements, with an average reactivity bias of <5u2009ppmB for the entire cycle. Comparisons to cycle 1 flux map data are also provided, and the average 2-D root-mean-square (rms) error during cycle 1 is 1.07%. To demonstrate the multicycle capability, a state point at beginning of cycle (BOC) 2 was also simulated and compared to plant data. The comparison of the cycle 2 BOC state has a reactivity difference of +3u2009ppmB from measurement, and the 2-D rms of the comparison in the flux maps is 1.77%. These results provide confidence in VERA-CSs capability to perform high-fidelity calculations for practical PWR reactor problems.

Subgroup Benchmark Calculations for the Intra-Pellet Nonuniform Temperature Cases

A benchmark suite has been developed by Seoul National University (SNU) for intrapellet nonuniform temperature distribution cases based on the practical temperature profiles according to the thermal power levels. Though a new subgroup capability for nonuniform temperature distribution was implemented in MPACT, no validation calculation has been performed for the new capability. This study focuses on bench-marking the new capability through a code-to-code comparison. Two continuous-energy Monte Carlo codes, McCARD and CE-KENO, are engaged in obtaining reference solutions, and the MPACT results are compared to the SNU nTRACER using a similar cross section library and subgroup method to obtain self-shielded cross sections.