Wenxi Tian

Comparison of Hydrogen Generation Rate between CORA-13 Test and MELCOR Simulation: Clad Solid-Phase Oxidation Models Using Self-Developed Code MYCOAC

Abstract To better understand the MELCOR oxidation and degradation models, past work compared the MELCOR model to a CORA experiment (CORA Test 13). These MELCOR analyses specifically focused on fuel bundle heatup and clad oxidation when compared to CORA test data. The comparison allowed the authors to investigate differences between hydrogen generation data and simulation results. Several potential reasons were considered for hydrogen generation rate differences, including MELCOR input power, heat transfer modeling, the clad solid-phase oxidation model, and the gaseous steam diffusion model. This work focuses on the possible uncertainty in the clad oxidation models used in MELCOR. First, the MELCOR nodalization approach for the CORA test was reviewed. Then, the temperature history and spatial variation were examined. One main focus was to consider other clad solid-phase oxidation models to compare the MELCOR models. This was accomplished by developing a separate model, MYCOAC, using MELCOR temperature predictions as input. Finally, the mass transfer resistance of steam diffusion to the clad surface was examined and found to be a small effect. While the Baker-Just solid-phase oxidation model showed better agreement with CORA data at low temperatures, the conclusion in this paper is that the oxidation models are not the major source of uncertainty in hydrogen generation rate differences. Future work will focus on heat transfer modeling of the CORA test.

Development of a Thermal-Hydraulic Analysis Code and Transient Analysis for a FHTR

The Fluoride-salt-cooled High-temperature Reactor (FHR) is an advanced reactor concept that uses high temperature TRISO fuel with a low-pressure liquid salt coolant. Design of Fluoride-salt-cooled High-temperature Test Reactor (FHTR) is a key step in the development of the FHR technology and is currently in progress both in China and the United States. An FHTR based on pebble bed core design with coolant temperature 600–700 °C is being planned for construction by the Chinese Academy of Sciences (CAS)’s Thorium Molten Salt Reactor (TMSR) Research Center, Shanghai Institute of Applied Physics (SINAP). This paper provides preliminary thermal hydraulic transient analyses of an FHTR using SINAP’s pebble core design as a reference case. A point kinetic model is calculated by developing a microcomputer code coupling with a simplified porous medium heat transfer model in the core. The founded models and developed code are applied to analyze the safety characteristics of the FHTR by simulating basic transient conditions including the unprotected loss of flow, unprotected overcooling, and unprotected transient overpower accidents. The results show that the SINAP’s pebble core design is an inherently safe reactor design.© 2014 ASME

Experimental Investigation of Air-Water CCFL in the Pressurizer Surge Line of AP1000

Abstract Countercurrent flow limitation (CCFL) may occur under certain flow conditions in the surge line, restricting the draining of water from the pressurizer and thus affecting the coolant inventory and water level in the reactor pressure vessel (RPV). The complexity of the AP1000 pressurizer surge line structure makes predicting CCFL fairly difficult, and there are still not enough CCFL studies on this complex structure. Based on an extensive literature survey, the authors of this paper are particularly aware of the need for improved CCFL models for the pressurizer surge line of AP1000. To investigate the CCFL phenomenon in the surge line assembly fixture of AP1000, a whole-visual test model of the surge line is designed with a scaling ratio of 1:4, and a test loop is established to carry out visualization experiments with air-water countercurrent flow (CCF). The whole-visual test section made of acrylic material is composed of a pressurizer simulator, a surge line tube, a hot leg T-type tube, and an RPV simulator. The air-water CCF experiments are conducted at atmospheric pressure and room temperature with the pressurizer simulator water level varying from 150 to 900 mm. The visual CCF experimental processes and CCFL phenomena are filmed by a high-speed camera and analyzed in detail. The pressure drops at different CCFL locations are measured and evaluated to explore the relationships between the CCFL characteristics and flow patterns in the surge line. The development process of the CCFL is defined as the CCFL region, which can be divided into different regions according of the changes in water mass flow and CCF flow behavior. The CCFL data are analyzed and compared using the air and water superficial velocities to study the effects of hysteresis and water level. Small discrepancies are found between the data of different water levels, reflecting the small but not-negligible influence of the upper tank water level. Empirical models for the CCFL in the surge line assembly fixture are explored preliminarily using Kutateladze-type correlation and Froude-Ohnesorge correlation. Deficiencies still exist in the present semiempirical models, inspiring a more in-depth study on the empirical models for CCFL in the surge line assembly fixture that considers the complex two-phase flow behaviors in the upper tank and near the joint between the upper tank and surge line tube. The present CCFL data are compared broadly and in detail with groups of CCFL data of similar former experiments to demonstrate the applicability of the present air-water CCFL data to the development of a CCFL prediction model for the prototype large-diameter surge line assembly fixture of the AP600/AP1000. We will perform much more experimental and theoretical work to study the detailed mechanism of these special phenomena and to develop a more applicable CCFL model for the geometry and conditions of the prototype large-diameter surge line assembly fixture.

3D Flow Field Analysis for a Traditional II+ PWR Containment Under Normal Condition

Flow field analysis is a foundation to many thermal-dynamic phenomena in the nuclear containment. There are several ventilation systems under normal condition to assure the proper environment for staff and facilities, and the two main ventilation systems are the Reactor Pit Ventilation system (EVC) and the Containment Continuous Ventilation system (EVR). The fans of the ventilation systems recycle the air in the containment to various rooms, and cooling coils are cooled by the Nuclear Island Chilled Water system (DEG). In this paper, we simulated the 3D flow field in a Chinese traditional generation II+ PWR reactor containment under normal conditions with a commercial CFD software ANSYS FLUENT.According the actual geometry data of the reactor plant, we built the geometry model, including the EVC and EVR system, and they were checked by plant experts to assure authenticity. Proper maximum mesh sizes were set for different parts at the ICEM CFD, and the grid number was about ten millions. We used the fan model in the ANSYS FLUENT to simulate fans in the ventilations. The calculated values of flow rate in ventilation systems were in good agreement with the design values.Copyright

The Research on Core Melting Process: Oxidation

Beyond-design basis accidents in the AP1000 may result in reactor core melting and are therefore termed core melt accidents. The aim of this work is to develop a code to calculate and analyze the oxidation of a single fuel rod with total failures of engineered safeguard systems under a certain beyond-design basis accident such as a gigantic earthquake which can result in station blackout and then total loss of coolant flow. Using the code, the responses of the most dangerous fuel rod in the AP1000 were calculated under the accident. A discussion involving fuel pellets melting, cladding rupture and oxidation, and hydrogen production then was carried out, focused on DNBR during coolant pump coastdown, the cladding intactness under different flow rates in natural circulation, and the delay effect on cladding rupture due to cladding oxidation. By the analysis of calculated results, several suggestions on guaranteeing the security of fuel rods were provided.Copyright

Transient Safety Analysis of a Transportable Fluoride-Salt-Cooled High-Temperature Reactor Using RELAP5-3D

Abstract A transportable fluoride-salt-cooled high-temperature reactor (TFHR) design with 20-MW(thermal) rated power and 18-month fuel cycle is proposed for off-grid applications. One of the design goals of the compact reactor core is potential transport by truck, rail, or air. Full-core thermal-hydraulic analyses and improvements using three-dimensional computational fluid dynamics (CFD) were performed previously to demonstrate the feasibility of a TFHR design at a nominal power of 20 MW(thermal). In this paper, the best-estimate system code Reactor Excursion Leak Analysis Program (RELAP5-3D) is adopted to study the transient behavior of this TFHR design and the safety characteristics of the primary loop system during accident conditions. The modeling results of the steady state were verified using CFD results with consideration of radial heat conduction between heat transfer unit cells. Four most challenging accidents of anticipated transient without scram were analyzed, as well as parametric studies of some key factors. These accidents include unprotected reactivity insertion accident (URIA), unprotected loss of heat sink (ULOHS), unprotected loss of flow (ULOF), and a combination accident of ULOF and ULOHS. The results indicate that transient temperature limits are not exceeded during the most severe accidents. They indicate satisfactory transient performance of the TFHR design. The transient temperature limit of structure material Hastelloy N, based on embrittlement phenomena, poses the most limiting constraint due to the small temperature margin of about 20 K in the accident combination of ULOF and ULOHS. Overall, TFHR is a sound reactor design from a thermal-hydraulic viewpoint.

Analysis of the Particulate Debris Bed Quenching During Top and Bottom Flood

The quenching characteristics of particulate debris bed during bottom and top flood is analyzed in this paper. The top flood model is formulated by dividing the quenching process into downward frontal period and upward frontal period, which are controlled by the counter-current flow limitation (CCFL) condition and effects of the incoming coolant subcooling and steam cooling in dry channels during quenching process. The bottom flood model is based on porous media theory under the assumption that the height of the two phase region is negligible and the particulate debris bed is divided into single phase liquid and single phase vapor region. The results calculated by these models were compared with the experimental data. The influences of porosity, initial debris temperature and other parameters on both the top and bottom quenching process were studied in this paper. During the top flood, the quenching velocity increased with the increase of the porosity and the decrease of the initial debris temperature. The porosity and initial debris temperature had a larger influence on quenching velocity compared with other parameters, such as initial coolant temperature and coolant flow rate. During the bottom flood, the quenching velocity also increased with the increase of the porosity and the decrease of the initial debris temperature. However, the coolant flow rate had a large influence on the quenching velocity unlike that during the top flood. Quenching from bottom may be superior to the quenching from top. The results can be expected to be useful to evaluate the quenching process of the particulate debris bed.Copyright

Loss-of-Flow-Accidents (LOFA) Study for 100 MW IPWR

The inherent system safety of the 100 MW integral pressurized water reactor (IPWR) can be improved by placing the core, the efficient once-through steam generators and the main coolant pumps in the reactor pressure vessel, and omitting some large pipes and valves in the primary coolant system which can prevent the occurrence of large break loss of coolant accident and reduce the possibility of core melt accident. The application of the passive safety systems simplifies the structures of IPWR and improves the economy of the reactor. In case of accidents, the primary coolant system establishes natural circulation to take the core decay heat away by passive safety systems using gravity and other natural driving forces, thereby enhancing the safety and reliability of the system IPWR. It’s of great significance to establish reasonable and correctable models, including the primary coolant system model, the second loop model and passive core cooling system model, to study thermal-hydraulic phenomena under steady state, transient state and accident conditions.Based on transient safety analysis program RELAP5/MOD3.4, 100 MW IPWR system was simulated. A series of models of reactor coolant system and passive safety systems were established. The main system models are composed of primary coolant system model, part of second loop model, passive safety injection system model and passive residual heat removal system model. The primary coolant system model includes core, lower plenum, downcomer, region of steam generators, upper plenum, riser, pressurizer, and surge line; the second loop model includes the main feed water line, the steam line, and steam generator tubes; passive safety injection system model includes core makeup tank, accumulator, automatic depressurization system, direct vessel injection line; and passive residual heat removal system model includes passive residual heat removal heat exchanger in containment refueling water storage tank. Based on the established models, the steady state was debugged with the RELAP5 input card.Steady state calculation was performed and the results agree well with designed values, which verifies the validity of the model and the input card. Using the steady state results as initial conditions, transient calculation was performed. Typical accidents (loss of main water accident) were calculated, which were relieved by auxiliary feedwater system (AFWS) and passive residual heat removal system (PRHR SYSTEM). The results, obtained from AFWS and PRHR SYSTEM, were contrasted and process of accident and thermal-hydraulic phenomena were analyzed according to transient calculation results. The transient calculation results showed that the integral PWR system and the passive safety system model can provide a reference for IPWR transient safety analysis.Copyright

Sub-Channel Analysis of Pb-Bi-Cooled Reactor With Modified COBRA-EN

Lead-bismuth eutectic (LBE) cooled fast reactor, one of the six types of reactors in Gen-IV, has very good inherent safety and significant advantages in reducing and burning nuclear wastes, enhancing economy. Also LBE cooled accelerator driven system (ADS) has been a very innovative and potential waste burner. COBRA-EN is a mature, stable and widely-used sub-channel analysis code for light water cooled reactor but it couldn’t be applied in Pb-Bi-cooled reactor directly. Some modifications were made for COBRA-EN in the present work, then the code was named COBRA-PB and was suitable for the sub-channel analysis of Pb-Bi-cooled reactor. The modified code was verified and validated with CFX and experimental results. There was a good agreement between the two results. Then sub-channel analysis of Pb-Bi-cooled reactor was done with the modified code.Copyright

Theoretical Research on Flow Instability in Parallel Channels Under Motion Conditions

In motion conditions, in addition to gravitational acceleration, a new acceleration was developed and it was added to the thermal hydraulics characteristics in flow channels. The additional acceleration leads to the different thermal hydraulic characteristics and will trigger the flow oscillation and even flow instability in parallel channels. In order to study the effect of the additional acceleration on the flow oscillation, the corresponding physical models are established in this work. Through the deduction of the mathematical model, the code for flow instability under motion conditions with Gear algorithm is developed. The flow oscillation curves, critical power, marginal stability boundary (MSB) are obtained. After comparison and analysis, it is found that some motion conditions lead to flow periodic oscillation. Different flow passage position results in different oscillation amplitudes. The marginal stability boundaries (MSB) under different motion conditions fit well, that is, the effect of motion conditions on MSB is small. Number of channels has little effect; however, channel arrangement influences the flow in every channel. These conclusions are of great significance in marine reactor design.Copyright