Analytical Evaluation of Debris Cooling and Spreading Behaviors at Molten Core in Severe Accident

Author(s):  
Akihiro Kobayashi ◽  
Shuichiro Miwa ◽  
Michitsugu Mori

On March 11, 2011, severe accident occurred at Fukushima Daiichi Nuclear Power Plant, and Units 1 to 3 of the plant have led to core melt. That is to say, melted fuel rods and core internals fell to the bottom of the Reactor Pressure Vessel (RPV). It is also believed that molten core has leaked into the reactor containment vessel. In order to plan for a safe molten core removal from the reactor, it is important to estimate the conditions of molten core by conducting analysis. Particular importance of the analysis is to understand the mechanisms of molten core spreading-cooling processes. However, sufficient understanding of this process has not been obtained yet. The main purpose of this study is to evaluate molten metal spreading-cooling phenomena and subsequently, estimate the conditions of the molten metal. In order to achieve the purpose, the Computational Fluid Dynamics (CFD) for thermal fluid analysis, STAR-CCM+ was utilized. In the simulation of the unsteady two-phase flow, the volume of fluid model was applied for the spreading and interfacial surface formation of molten metal with the surrounding air. The key parameter for the molten metal spreading is the temperature dependent viscosity of molten metal. To assess the validity of this model, the analysis of the VULCANO VE-U7, molten metal spreading experiment, has been compared with simulation results.

2021 ◽  
Author(s):  
Chandan Kumawat ◽  
Bhupendra Kumar Sharma ◽  
Khalid Saad Mekheimer

Abstract A two-phase blood flow model is considered to analyze the fluid flow and heat transfer in a curved tube with time-variant stenosis. In both core and plasma regions, the variable viscosity model ( Hematocrit and non linear temperature-dependent, respectively) is considered. A toroidal coordinate system is considered to describe the governing equations. The perturbation technique in terms of perturbation parameter ε is used to obtain the temperature profile of blood flow. In order to find the velocity, wall shear stress and impedance profiles, a second-order finite difference method is employed with the accuracy of 10−6 in the each iteration. Under the conditions of fully-developed flow and mild stenosis, the significance of various physical parameters on the blood velocity, temperature, wall shear stress (WSS) and impedance are investigated with the help of graphs. A validation of our results has been presented and comparison has been made with the previously published work and present study, and it revels the good agreement with published work. The present mathematical study suggested that arterial curvature increase the fear of deposition of plaque (atherosclerosis), while, the use of thermal radiation in heat therapies lowers this risk. The positive add in the value of λ1 causes to increase in plasma viscosity; as a result, blood flow velocity in the stenosed artery decreases due to the assumption of temperature-dependent viscosity of the plasma region. Clinical researchers and biologists can adopt the present mathematical study to lower the risk of lipid deposition, predict cardiovascular disease risk and current state of disease by understanding the symptomatic spectrum, and then diagnose patients based on the risk.


Author(s):  
Miki Saito ◽  
Taizo Kanai ◽  
Satoshi Nishimura ◽  
Yoshihisa Nishi

Abstract Understanding the mechanism of fission product (FP) removal by pool scrubbing is essential for improving the prediction accuracy of FP emissions concerning severe accident (SA) in a nuclear power plant. Since FP migrates from a gas-phase to a liquid-phase via a gas-liquid interface, the FP removal efficiency by pool scrubbing is largely affected by the flow regime of gas-liquid two-phase flow. In order to gain a deeper understanding of the influence of gas properties on flow regimes, experiments were performed by injecting helium (He) and nitrogen (N2) gas mixtures of several volumetric ratios through a pool of stagnant water. The result suggests clear effects of gas compositions on gas-liquid two-phase flow, where both void and holdup fractions were found to increase with N2 fraction in the supplied gas. The results were compared with previous studies, and a detailed analysis of bubble characteristics for different compositions of gases was performed using a wire-mesh sensor (WMS). This paper also illustrates further research aspects needed to discuss the effect of its results on FP removal efficiency in a SA, and to acquire comprehensive physics behind such gas property influences on two-phase flow.


Author(s):  
Jiwon Mun ◽  
Jaehyung Ju ◽  
Byoung-Gwan Yun ◽  
Byung-Moon Chang ◽  
Doo-Man Kim

Investment casting processes are influenced by a variety of parameters. Many researches considering viscosity as a constant have been conducted up to this point. In particular, however, viscosity with temperature change has not been much accounted for solidification and heat transfer simulation of molten metal in the investment casting process. In addition, analysis of behavior of metal flow as well as air gap problems for complex network structures have not been investigated much. The aim of this study is to build transient metal flow and velocity profile models considering temperature dependent viscosity in investment casting processes of cellular structures. In this study, a Computational Fluid Dynamics (CFD) modeling tool was used for metal flow and velocity profile in investment casting processing using User Defined Function (UDF) for temperature dependent viscosity. The results of the metal flow and velocity profile inside of the simple cylindrical geometry are represented. It is shown that for the validation of the numerical simulation, the velocity profile between analytical and numerical approaches showed very good agreement. Analytical approaches showed that velocity was reduced with the increase in viscosity, which is applied as a function of temperature. In particular, rapid decreasing in velocity was shown from under the melting temperature of the molten metal. There was no movement on metal flow at the room temperature. Numerical approaches showed that the liquid metal began to be solidified from the wall surface inside of the mold. For the same simulation time, it was shown that the metal flow in a cylinder that has 1mm diameter showed better fluidity rather than that of the cylinder that has 2mm diameter due to the increase in adhesion between liquid metal and the surface of the mold and surface tension between molten metal and air. The effective diameter by solidification is decreased with the time change.


Author(s):  
Chue Cui Ting ◽  
Afiq Mohd Laziz ◽  
Khoa Dang Dang Bui ◽  
Ngoc Thi Nhu Nguyen ◽  
Pha Ngoc Bui ◽  
...  

Microfluidic systems undergo rapid expansion of its application in different industries over the few decades as its surface tension-dominated property provides better mixing and improves mass transfer between two immiscible liquids. Synthesis of biodiesel via transesterification of vegetable oil and methanol in microfluidic systems by droplet flow requires separation of the products after the reaction occurred. The separation technique for multiphase fluid flow in the microfluidic system is different from the macro-system, as the gravitational force is overtaken by surface force. To understand these phenomena completely, a study on the hydrodynamic characteristics of two-phase oil-methanol system in microchannel was carried out. A multiphase Volume of Fluid model was developed to predict the fluid flow in the microchannel. An inline separator design was proposed along with its variable to obtain effective separation for the oil-methanol system. The separation performance was evaluated based on the amount of oil recovered and its purity. The capability of the developed model has been validated through a comparison of simulation results with published experiment. It was predicted that the purity of recovered oil was increased by more than 46% when the design with side openings arranged at both sides of the microchannel. The highest percentage recovery of oil from the mixture was simulated at 91.3% by adding the number of side openings to ensure the maximum recovery. The oil that was separated by the inline separator was predicted to be at 100% purity, which indicates that no methanol contamination throughout the separation process. The purity of the separated product can be increased by manipulating the pressure drop across the side openings. Hence, it can be concluded that the separation in a large diameter microchannel system is possible and methodology can be tuned to achieve the separation goal. Finally, the simulation results showed that the present volume of fluid model had a good agreement with the published experiment.


2015 ◽  
Vol 2015 ◽  
pp. 1-7
Author(s):  
Dao-gang Lu ◽  
Fan Zhang ◽  
Dan-ting Sui ◽  
Xue-zhang Xi ◽  
Lei-bo Yu

Chinese large-capacity advanced PWR under construction in China is a new and indispensable reactor type in the developing process of NPP fields. At the same time of NPP construction, accident sequences prediction and operators training are in progress. Since there are some possible events such as feedwater pumps trip in secondary circuit may lead to severe accident in NPP, training simulators and engineering simulators of CI are necessary. And, with an increasing proportion of nuclear power in China, NPP will participate in regulating peak load in power network, which requires accuracy calculation and control of secondary circuit. In order to achieve real-time and full scope simulation in the power change transient and accident scenarios, RELAP5/MOD 3.4 code has been adopted to model the secondary circuit for its advantage of high calculation accuracy. This paper describes the model of steady state and turbine load transient from 100% to 40% of secondary circuit using RELAP5 and provides a reasonable equivalent method to solve the calculation divergence problem caused by dramatic two-phase condition change while guaranteeing the heat transfer efficiency. The validation of the parameters shows that all the errors between the calculation values and design values are reasonable and acceptable.


Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1779
Author(s):  
Piotr Siwicki ◽  
Marcin Krukowski ◽  
Jan Studziński ◽  
Bartosz Szeląg ◽  
Rafał Wojciechowski

The paper analyzes the possibility of using the CFD (Computational Fluid Dynamics) method to predict the amount of sewage remaining in siphons after a full air blast of the pressure sewer system. For this purpose, the results from measurements carried out on a laboratory installation were compared with the results obtained from modelling using a spatial model (3D) and a plain model (2D) of the installation. To determine these models, the structure of the VOF (Volume of Fluid) model was used in the CFD method. The simulation calculations carried out make it possible to state that the use of the plain model with the development of the installation modelled in the plan does not result in significant deterioration of the obtained results. The possibility of using 2D models for modelling pumped sewer systems allows for a significant shortening of the calculation time, which, in practice, results in the possibility of modelling much larger and longer installations than is possible with 3D models.


Author(s):  
A. M. Razmi ◽  
B. Firoozabadi

In the present study, the presence of a baffle and its effect on the hydrodynamics of the flow in a primary settling tank has been investigated experimentally by ADV (Acoustic Doppler Velocimeter). On the other hand, the characteristics of this flow field were simulated by an unsteady two-phase finite volume method, and VOF (Volume of Fluid) model; and results were evaluated by the experimental data. The numerical calculation performed by using k–ε RNG model agrees well with experiments. It depicts the ability of this method in predicting the velocity profile and flow structure. In addition, the optimum position of the baffle to achieve the best performance of the tank was determined by applying the above mentioned numerical scheme.


Author(s):  
Alexandre Zanchetti ◽  
Mickael Hassanaly ◽  
Hervé Cordier ◽  
Antonio Sanna ◽  
Namane Mechitoua ◽  
...  

The Fukushima accident reminded us of the possible consequences in terms of radiological release that can result from a hydrogen explosion in a nuclear power plant, and, specifically, within the containment of a water cooled reactor building. Some mitigation means against hydrogen hazards exist but performance improvements in numerical tools simulating thermal-hydraulic flows and hydrogen combustion are necessary to allow realistic assessments of severe accident consequences in the containment. In this context, EDF works on CFD simulation of hydrogen distribution in penalized conditions. After dealing with cases for which the water spray system was assumed to be unavailable, and so treated with single-phase CFD code [1] [2], the present paper content is now about simulation and analysis of the local hydrogen concentration in the case of a severe accident for which the water spray system is available. Numerical developments of a multi-phase CFD code (Neptune_CFD) and code validation lead to consistent simulations. The numerical simulation performed by EDF confirms the favorable safety impact of water spray on pressure and temperature for a LOCA scenario occurring on a 1300 MWe Pressurized Water Reactor. Nevertheless, CFD results show that the activation of the spray system before hydrogen injection gives greater hydrogen concentration. So, in the future, to better assess hydrogen risk, EDF will perform computations at CFD taking into account the interaction between combustion and water sprays.


Author(s):  
Milan Amižić ◽  
Estelle Guyez ◽  
Jean-Marie Seiler

In the frame of severe accident research for the second and the third generation of nuclear power plants, some aspects of the concrete cavity ablation during the molten corium–concrete interaction are still remaining issues. The determination of heat transfer along the interfacial region between the molten corium pool and the ablating basemat concrete is crucial for the assessment of concrete ablation progression and eventually the basemat melt-through. For the purpose of experimental investigation of thermal-hydraulics inside a liquid pool agitated by gas bubbles, the CLARA project has been launched jointly by CEA, EDF, IRSN, GDF-Suez and SARNET. The CLARA experiments are performed using simulant materials and they reveal the influence of superficial gas velocity, liquid viscosity and pool geometry on the heat transfer coefficient between the internally heated liquid pool and vertical and horizontal pool walls maintained at uniform temperature. The first test campaign has been conducted with the smallest pool configuration (50 cm × 25 cm × 25 cm). The tests have been performed with liquids covering a wide range of dynamic viscosity from approximately 1 mPa s to 10000 mPa s. This paper presents some preliminary conclusions deduced from the experiments which involve a liquid pool with the gas injection only from the bottom plate. A comparison with existing models for the assessment of heat transfer has also been carried out.


Sign in / Sign up

Export Citation Format

Share Document