Robert K. Salko

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.

VERA-CS Modeling and Simulation of PWR Main Steam Line Break Core Response to DNB

The Virtual Environment for Reactor Applications core simulator (VERA-CS) being developed by the Consortium for the Advanced Simulation of Light Water Reactors (CASL) includes coupled neutronics, thermal-hydraulics, and fuel temperature components with an isotopic depletion capability. The neutronics capability employed is based on MPACT, a three-dimensional (3-D) whole core transport code. The thermal-hydraulics and fuel temperature models are provided by the COBRA-TF (CTF) subchannel code. As part of the CASL development program, the VERA-CS (MPACT/CTF) code system was applied to model and simulate reactor core response with respect to departure from nucleate boiling ratio (DNBR) at the limiting time step of a postulated pressurized water reactor (PWR) main steamline break (MSLB) event initiated at the hot zero power (HZP), either with offsite power available and the reactor coolant pumps in operation (high-flow case) or without offsite power where the reactor core is cooled through natural circulation (low-flow case). The VERA-CS simulation was based on core boundary conditions from the RETRAN-02 system transient calculations and STAR-CCM+ computational fluid dynamics (CFD) core inlet distribution calculations. The evaluation indicated that the VERA-CS code system is capable of modeling and simulating quasi-steady state reactor core response under the steamline break (SLB) accident condition, themorexa0» results are insensitive to uncertainties in the inlet flow distributions from the CFD simulations, and the high-flow case is more DNB limiting than the low-flow case.«xa0less

Assessment of the Subchannel Code CTF for Single- and Two-Phase Flows

Abstract As part of the Consortium for Advanced Simulation of Light Water Reactors, the subchannel code CTF is being used for single-phase and two-phase flow analysis under light water reactor operating conditions. Accurate determination of flow distribution, pressure drop, and void content is crucial for predicting margins to thermal crisis and ensuring more efficient plant performance. In preparation for the intended applications, CTF has been validated against data from experimental facilities comprising the General Electric (GE) 3 × 3 bundle, the boiling water reactor full-size fine-mesh bundle tests (BFBTs), the Risø tube, and the pressurized water reactor subchannel and bundle tests (PSBTs). Meanwhile, the licensed, well-recognized subchannel code VIPRE-01 was used to generate a baseline set of simulations for the targeted tests and solution parameters were compared to the CTF results. The flow split verification problem and single-phase GE 3 × 3 results are essentially in perfect agreement between the two codes. For the two-phase GE 3 × 3 cases, flow and quality discrepancies arise in the annular-mist flow regime, yet significant improvement is observed in CTF when void drift and two-phase turbulent mixing enhancement are considered. The BFBT pressure drop benchmark shows close agreement between predicted and measured results in general, although considerable overprediction by CTF is observed at relatively high void locations of the facility. This overestimation tendency is confirmed by the Risø cases. While overall statistics are satisfactory, both BFBT and PSBT bubbly-to-churn flow void contents are markedly overpredicted by CTF. The issues with two-phase closures such as turbulent mixing, interfacial and wall friction, and subcooled boiling heat transfer need to be addressed. Preliminary sensitivity studies are presented herein, but more advanced models and code stability analysis require further investigation.

Development and Testing of CTF to Support Modeling of BWR Operating Conditions

This milestone supports developing and assessing COBRA-TF (CTF) for the modeling of boiling water reactors (BWRs). This is achieved in three stages. First, a new preprocessor utility that is capable of handling BWR-specic design elements (e.g., channel boxes and large water rods) is developed. A previous milestone (L3:PHI.CTF.P12.01) led to the development of this preprocessor capability for single assembly models. This current milestone expands this utility so that it is applicable to multi-assembly BWR models that can be modeled in either serial or parallel. The second stage involves making necessary modications to CTF so that it can execute these new models. Specically, this means implementing an outer-iteration loop, specic to BWR models, that equalizes the pressure loss over all assemblies in the core (which are not connected due to the channel boxes) by adjusting inlet mass ow rate. A third stage involves assessing the standard convergence metrics that are used by CTF to determine when a simulation is steady-state. The nal stage has resulted in the implementation of new metrics in the code that give a better indication of how steady the solution is at convergence. This report summarizes these eorts and provides a demonstration of CTFs BWR-modeling capabilities. CASL-U-2016-1030-000