Stratification and Heat Transfer of Molten Corium Pool for In-Vessel Retention

Author(s):  
Qingan Xiang ◽  
Jian Deng ◽  
Dahuan Zhu ◽  
Xiaoli Wu ◽  
Jinsheng Bi ◽  
...  

Abstract In-vessel retention (IVR) consists in cooling the molten corium contained in the lower head of reactor vessel by natural convection and reactor cavity flooding. The general approach which is used to study IVR problems is a “bounding” approach which consists in assuming a specified corium pool stratification in the vessel and then demonstrating that the vessel can cope with the resulting thermal and mechanical loads. Thermal loading on the vessel is controlled by the convective heat transfer inside the molten corium pool. Traditionally molten corium pool in the lower head was expected to stratify into two-layer with the dense oxide pool at the bottom and the light metal layer on the top. Based on the MASCA experiments, the increased density of the metal layer attributed to a transfer of uranium metal leads to inverse stratification with a heavy metal layer relocating below the oxide pool. This behavior can be explained by physicochemical interaction between the oxidic and metallic phases of the corium pool. Therefore, a methodology which couples physicochemical effects and thermal hydraulics has been developed to address the IVR issue. The main purpose of this paper is to present this methodology and its application for calculate stratification probability of two-layer and three-layer configuration, analyze the safety margin of IVR in two-layer and three-layer configuration, and evaluate the lower head heat thermal failure probability.

Author(s):  
Walter Villanueva ◽  
Chi-Thanh Tran ◽  
Pavel Kudinov

An in-vessel stage of a severe core melt accident in a Nordic type Boiling Water Reactor (BWR) is considered wherein a decay-heated pool of corium melt inflicts thermal and mechanical loads on the lower-head vessel wall. This process induces creep leading to a mechanical failure of the reactor vessel wall. The focus of this study is to investigate the effect of Control Rod Guide Tube (CRGT) and top cooling on the modes of global vessel failure of the lower head. A coupled thermo-mechanical creep analysis of the lower head is performed and cases with and without CRGT and top cooling are compared. The debris bed heat-up, re-melting, melt pool formation, and heat transfer are calculated using the Phase-change Effective Convectivity Model and transient heat transfer characteristics are provided for thermo-mechanical strength calculations. The creep analysis is performed with the modified time hardening creep model and both thermal and integral mechanical loads on the reactor vessel wall are taken into account. Known material properties of the reactor vessel as a function of temperature, including the creep curves, are used as an input data for the creep analysis. It is found that a global vessel failure is imminent regardless of activation of CRGT and top cooling. However, if CRGT and top cooling is activated, the mode and timing of failure is different compared to the case with no CRGT and top cooling. More specifically, with CRGT and top cooling, there are two modes of global vessel failure depending on the size of the melt pool: (a) ‘ballooning’ of the vessel bottom for smaller pools, and (b) ‘localized creep’ concentrated within the vicinity of the top surface of the melt pool for larger pools. Without CRGT and top cooling, only a ballooning mode of global vessel failure is observed. Furthermore, a considerable delay (about 1.4 h) on the global vessel failure is observed for the roughly 30-ton debris case if CRGT and top cooling is implemented. For a much larger pool (roughly 200-ton debris), no significant delay on the global vessel failure is observed when CRGT and top cooling is implemented, however, the liquid melt fraction and melt superheat are considerably higher in non-cooling case.


Author(s):  
Jun Yeong Jung ◽  
Yong Hoon Jeong

In-Vessel Retention by External Reactor Vessel Cooling (IVR-ERVC) is method of removing the decay heat by cooling reactor vessel after corium relocation, and is also one of severe accident management strategies. Estimating heat transfer coefficients (HTCs) is important to evaluate heat transfer capability of the ERVC. In this study, the HTCs of outer wall of reactor vessel lower head were experimentally measured under the IVR-ERVC situation of Large Loss of Coolant Accident (LLOCA) condition. Experimental equipment was designed to simulate flow boiling condition of ERVC natural circulation, and based on APR+ design. This study focused on effects of real reactor vessel geometry (2.5 m of radius curvature) and material (SA508) for the HTCs. Curved rectangular water channel (test section) was design to simulate water channel which is between the reactor vessel lower head outer wall and thermal insulator. Radius curvature, length, width and gap size of the test section were respectively 2.5 m, 1 m, 0.07 m and 0.15 m. Two connection parts were connected at inlet and outlet of the test section to maintain fluid flow condition, and its cross section geometry was same with one of test section. To simulate vessel lower head outer wall, thin SA508 plate was used as main heater, and test section supported the main heater. Thickness, width, length and radius curvature of the main heater were 1.2 mm, 0.07 m, 1 m and 2.5 m respectively. The main heater was heated by DC rectifier, and applied heat flux was under CHF value. The test section was changed for each experiment. The HTCs of whole reactor vessel lower head (bottom: 0 ° and top: 90 °) were measured by inclining the test section, and experiments were conducted at four angular ranges; 0–22.5, 22.5–45, 45–67.5 and 67.5–90 °. DI water was used as working fluid in this experiment, and all experiments were conducted at 400 kg/m2s of constant mass flux with atmospheric pressure. The working fluid temperatures were measured at two point of water loop by K-type thermocouple. The main heater surface temperatures were measured by IR camera. The main heater was coated by carbon spray to make uniform surface emissivity, and the IR camera emissivity calibration was also conducted with the coated main heater. The HTCs were calculated by measured main heater surface temperature. In this research, the HTC results of 10, 30, 60 and 90 ° inclination angle were presented, and were plotted with wall super heat.


Author(s):  
D. L. Knudson ◽  
J. L. Rempe

Molten core materials may relocate to the lower head of a reactor vessel in the latter stages of a severe accident. Under such circumstances, in-vessel retention (IVR) of the molten materials is a vital step in mitigating potential severe accident consequences. Whether IVR occurs depends on the interactions of a number of complex processes including heat transfer inside the accumulated molten pool, heat transfer from the molten pool to the reactor vessel (and to overlying fluids), and heat transfer from exterior vessel surfaces. SCDAP/RELAP5-3D© has been developed at the Idaho National Engineering and Environmental Laboratory to facilitate simulation of the processes affecting the potential for IVR, as well as processes involved in a wide variety of other reactor transients. In this paper, current capabilities of SCDAP/RELAP5-3D© relative to IVR modeling are described and results from typical applications are provided. In addition, anticipated developments to enhance IVR simulation with SCDAP/RELAP5-3D© are outlined.


2018 ◽  
Vol 111 ◽  
pp. 293-302 ◽  
Author(s):  
Luteng Zhang ◽  
Simin Luo ◽  
Yapei Zhang ◽  
Wenxi Tian ◽  
G.H. Su ◽  
...  

Author(s):  
Michel Arnal ◽  
Christian Precht ◽  
Thomas Sprunk ◽  
Tobias Danninger ◽  
John Stokes

The present paper outlines a practical methodology for improved virtual prototyping, using as an example, the recently re-engineered, internally-cooled 1st stage blade of a 40 MW industrial gas turbine. Using the full 3-D CAD model of the blade, a CFD simulation that includes the hot gas flow around the blade, conjugate heat transfer from the fluid to the solid at the blade surface, heat conduction through the solid, and the coolant flow in the plenum is performed. The pressure losses through and heat transfer to the cooling channels inside the airfoil are captured with a 1-D code and the 1-D results are linked to the three-dimensional CFD analysis. The resultant three-dimensional temperature distribution through the blade provides the required thermal loading for the subsequent structural finite element analysis. The results of this analysis include the thermo-mechanical stress distribution, which is the basis for blade life assessment.


2005 ◽  
Author(s):  
H. K. Cho ◽  
D. U. Seo ◽  
M. O. Kim ◽  
G. C. Park

In the HTGR (High Temperature Gas Cooled Reactor), the Reactor Cavity Cooling System (RCCS) is equipped to remove the heat transferred from the reactor vessel to the structure of the containment. The function of the RCCS is to dissipate the heat from the reactor cavity during normal operation including shutdown. The system also removes the decay heat during the loss of forced convection (LOFC) accident. A new concept of the water pool type RCCS was proposed at Seoul National University. The system mainly consists of two parts, water pool located between the containment and reactor vessel and five trains of air cooling system installed in the water pool. In normal operations, the heat loss from the reactor vessel is transferred into the water pool via cavity and it is removed by the forced convection of air flowing through the cooling pipes. During the LOFC accident, the after heat is passively removed by the water tank without the forced convection of air and the RCCS water pool is designed to provide sufficient passive cooling capacity of the after heat removal for three days. In the present study, experiments and numerical calculations using CFX5.7 for the water pool and cooling pipe were performed to investigate the heat transfer characteristics and evaluate the heat transfer coefficient model of the MARS-GCR (Multi-dimensional Analysis of Reactor Safety for Gas Cooled Reactor Analysis) which was developed for the safety analysis of the gas cooled reactor. From the results of the experiments and CFX calculations, heat transfer coefficients inside the cooling pipe were calculated and those were used for the assessment for the heat transfer coefficient model of the MARS-GCR.


1981 ◽  
Vol 40 (3) ◽  
pp. 235-238
Author(s):  
A. M. Arkharov ◽  
A. I. Ageev ◽  
V. I. Pryanichnikov ◽  
N. B. Rubin

Author(s):  
Yan Dapeng ◽  
Ying Luo

Metallic insulation is commonly used in reactor vessel because of its resistance to radiation and corrosion. Since the main mode of heat loss of reactor vessel is thermal radiation, the ability to prevent radiation heat transfer is important for metallic insulation. But the thermal conductivity of metallic insulation is difficult to calculate owing to their complex geometry. This article uses FLUENT 14.0 to obtain the important parameter “view factor”, and then develops a computational model of effective conductivity of metallic insulation. Heat transfer test of metallic insulation was done, and the numerical simulation of metallic insulation was also performed. Based on results of test and simulation, the computational model is modified. The modified model can fit the test result better. Based on the modified model, the effective conductivity of metallic insulation increases with the increase of temperature of hot side and cold side, among which the temperature of hot side influences more. And when the temperature is high, the effective conductivity increases much faster.


1999 ◽  
Author(s):  
Larry James Siefken ◽  
Eric Wesley Coryell ◽  
Seungho Paik ◽  
Han Hsiung Kuo
Keyword(s):  

Author(s):  
Ewelina Pazera

This work is about a heat transfer phenomenon in relation to the periodically laminated composite. The specific type of thermal loading, analyzed in this paper, require formulation of Robin boundary conditions. To consider a layered structure of analyzed composite, the tolerance averaging technique is used. This method allows to take into account a thickness of the layers and obtain the equations with continuous coefficients. To solve these equations, the finite difference method is used, because an analytical solution is not available in this case (in contrast to the analogous issue in relation to a homogeneous layer).


Sign in / Sign up

Export Citation Format

Share Document