Multi-Scale Modeling of Thermal-Fluid Phenomena Related to Loss of Forced Circulation Transient in HTGRs
Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS), an effort had been pursued to support High Temperature Gas-cooled Reactors (HTGR) technology development and its modeling and simulation needs. There was a particular need for advanced modeling & simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially the core and lower and upper plena, during safety-related transients. In this report, two main such activities were presented: 1) system level HTGR simulations using system analysis code SAM for both normal operating and accident conditions; and 2) three-dimensional (3D) computational fluid dynamics (CFD) simulations on Deteriorated Turbulent Heat Transfer (DTHT) of interest to HTGR transients, including code benchmark and closure correlation development. In this work, Modular High Temperature Gas-cooled Reactor (MHTGR) full loop simulations were performed using System Analysis Module (SAM) code. Two approaches were used to model the reactor core; (i) unit-cell or lumped parameter model, and (ii) detailed ring model, to address the different needs under different reactor conditions. Aside from the core model changes other reactor components in the primary system, such as upper/lower plenum, heat exchanger and blower etc, were kept the same. In the unit-cell core model, the lumped parameter approach was used to model the coolant channels and the heat structures using 1-D and 2-D components, respectively. The fuel assembly heat structure was modeled based on unit-cell thickness approximation. In normal and controlled operating conditions, the calculated maximum temperatures for both coolant and heat structures were within the designed limits. The unit-cell model is, however, not able to capture the radial heat conduction through fuel assemblies, which is the dominant heat transfer mechanism during accident transients, such as pressurized conduction cooldown (PCC) event. To address this modeling challenge, ring models were developed. A simplified ring model was developed first. In this simplified model, the active core was simulated with 9 circular rings consists of 6 rings for homogenized fuel heat structure and 3 rings for gas coolant. It was quickly found that this model gives rather large simulation error as it was not able to correctly capture the solid structure representative geometries. A detailed ring model was developed to address this issue. The active core was simulated with 99 circular rings consists of 66 rings for homogenized fuel heat structure and 33 rings for gas coolant. In both ring models, six additional rings were included to represent inner reflector, outer reflector, core barrel, Reactor Pressure Vessel (RPV), RPV coolant channel, and Reactive Cavity Cooling System (RCCS). The detailed ring model was applied in two different SAM simulations; (i) normal operating conditions, and (ii) accident scenarios (more specifically, PCC event), for MHTGR design. In both simulations, the SAM calculated results for coolant and heat structure temperatures were within the reactor designed limits. The newly proposed detailed ring model approach in this work has proved to be a promising method to address the challenging modeling and simulation needs in the MHTGR design. As demonstrated in this work, this model was suitable for MHTGR simulations under both normal and accident conditions. Apart from the SAM model development for gas-cooled reactors, 3D Computational Fluid Dynamics (CFD) models were developed and validated with relaminarized Deteriorated Turbulent Heat Transfer (DTHT) MIT benchmark test. Two different CFD codes, Nek5000 and STAR-CCM+, were used for this benchmark activity. Compared with MIT data, Nek5000 shows very good agreement between its simulated and experimentally measured wall temperatures. In contrast, STAR-CCM+ under-predicted wall temperature due to over predicted turbulence in the wall region. The Nek5000 code was also used to develop wall-heat transfer correlations for laminar flow in cylindrical tube. Several simulations were performed for various Reynolds flow and wall heat fluxes. Significant deviations were found when compared the predicted heat transfer coefficient and the Sieder-Tate correlation. A new set of heat transfer correlation was proposed for laminar flow which is valid for 400<Re<2000.