• Journal of Mechanical Science and Technology
• /
• v.17 no.9
• /
• pp.1380-1387
• /
• 2003

When a cold HPSI (High pressure Safety Injection) fluid associated with an overcooling transient, such as SGTR (Steam Generator Tube Rupture), MSLB (Main Steam Line Break) etc., enters the cold legs of a stagnated primary coolant loop, thermal stratification phenomena will arise due to incomplete mixing. If the stratified flow enters the downcomer of the reactor pressure vessel, severe thermal stresses are created in a radiation embrittled vessel wall by local overcooling. As general thermal-hydraulic system analysis codes cannot properly predict the thermal stratification phenomena, RG 1.154 requires that a detailed thermal-mixing analysis of PTS (pressurized Thermal Shock) evaluation be performed. Also. previous PTS studies have assumed that the thermal stratification phenomena generated in the stagnated loop side of a partially stagnated primary coolant loop are neutralized in the vessel downcomer by the strong flow from the unstagnated loop. On the basis of these reasons, this paper focuses on the development of a 3-dimensional thermal-mixing analysis model using PHOENICS code which can be applied to both partial and total loop stagnated cases. In addition, this paper verifies the fact that, for partial loop stagnated cases, the cold plume generated in the vessel downcomer due to the thermal stratification phenomena of the stagnated loop is almost neutralized by the strong flow of the unstagnated loop but is not fully eliminated.

• Journal of the Korean Society of Visualization
• /
• v.15 no.1
• /
• pp.11-18
• /
• 2017

A two-phase flow observed in a heat exchanger or nuclear power generation often has a profound effect on undesirable noise or flow characteristics. Void fraction, which refers to the ratio of gas (or liquid) to the total fluid, affects heat transfer coefficient, vibration and so forth. In other words, void fraction is one of most important parameters in two-phase flow since it contributes to comprehend the characteristics of two-phase flow. We developed a two-phase flow visualization system to measure cross-sectional and volumetric void fractions by using quick closing valves and image processing software. With this system, we could observe the plug, slug, and stratified flow patterns of two-phase flow and measure a myriad of void fractions. As a consequence of the experiment, we found that the estimated void fractions were largely coincident with the predictive values by Chisholm model.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• Nuclear Engineering and Technology
• /
• v.44 no.8
• /
• pp.871-888
• /
• 2012

This paper aims at presenting the state of the art, the recent progress, and the perspective for the future, in the modelling of two-phase flow in the horizontal legs of a PWR. All phenomena relevant for safety analysis are listed first. The selection of the modelling approach for system codes is then discussed, including the number of fluids or fields, the space and time resolution, and the use of flow regime maps. The classical two-fluid six-equation one-pressure model as it is implemented in the CATHARE code is then presented and its properties are described. It is shown that the axial effects of gravity forces may be correctly taken into account even in the case of change of the cross section area or of the pipe orientation. It is also shown that it can predict both fluvial and torrential flow with a possible hydraulic jump. Since phase stratification plays a dominant role, the Kelvin-Helmholtz instability and the stability of bubbly flow regime are discussed. A transition criterion based on a stability analysis of shallow water waves may be used to predict the Kelvin-Helmholtz instability. Recent experimental data obtained in the METERO test facility are analysed to model the transition from a bubbly to stratified flow regime. Finally, perspectives for further improvement of the modelling are drawn including dynamic modelling of turbulence and interfacial area and multi-field models.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.5
• /
• pp.381-388
• /
• 2009

In some portions of nuclear piping systems, stratification phenomena may occur due to the density difference between hot and cold stream. When the temperature difference is large, the stratified flow under diverse operating conditions can produce high thermal stress, which leads to unanticipated piping integrity issues. The objectives of this research are to examine controvertible numerical factors such as model size, grid resolution, turbulent parameters, governing equation, inflow direction and pipe wall. Parametric three-dimensional computational fluid dynamics analyses were carried out to quantify effects of these parameters on the accuracy of temperature profiles in a typical nuclear piping with complex geometries. Then, as a key finding, it was recommended to use optimized mesh of real piping with the conjugated heat transfer condition for accurate thermal stratification analyses.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.4
• /
• pp.65-72
• /
• 2011

Numerical simulation has been performed to investigate the bubbly flow in the liquid with vertical temperature gradient. The objective of this study is to establish an accurate numerical prediction program of gas-liquid two-phase flows in a vertical temperature gradient condition, whose mathematical model is formulated by the Eulerian-Lagrangian model. The present numerical results reveal the temperature mixing mechanism and the fluid dynamical characteristics induced by the bubbly flow in the liquid with stratified temperature. The effects of bubble radius, void fraction, and gas flow rate on bubble-driven convective flow are considered, too.

• Solar Energy
• /
• v.7 no.2
• /
• pp.37-53
• /
• 1987

This study has been prepared for the purpose of developing the system performance prediction method of natural circulation solar hot water system. The storage tank of the natural circulation solar hot water system equipped with flat-plate solar collector is located at higher elevation than the solar collectors. Therefor, the storage tank temperature distribution formed accordance with configuration of storage tank by flow rate of circulating fluid affect system collection efficiency. In this study measure the storage tank temperature distribution with various experimental system under real sun condition and present the theoretical prediction method of the storage tank temperature. Moreover measure the flow rate not only day-time but also night-time reverse flow rate with die injection visual flow meter. Main conclusion obtain from the present study is as follows; 1) The storage tank temperature distribution above the connecting pipe connection position is the same as that of the fully mixed tank and below the connection position is the same as that of stratified tank. 2) The system performance sensitive to the storage tank temperature distribution. Therefore detailed tank model is necessary. Average storage tank temperature can be calculate 3% and storage tank temperature profile can get less than 10% difference with this model system.

• Nuclear Engineering and Technology
• /
• v.52 no.1
• /
• pp.27-36
• /
• 2020

Natural convection and thermal stratification phenomena are found in large water pools that are being used as heat sinks for decay heat removal from the reactor core using passive heat removal systems. In this study, the two-phase (water and air) natural convection and thermal stratification phenomena with conjugate heat transfer in the rectangular enclosure were investigated numerically using ANSYS Fluent 17.2 code. The transient numerical simulations of these phenomena in the full-scale computational domain of the experimental facility were performed. Generation of water vapour bubbles around the heater rod and evaporation phenomena were included in this numerical investigation. The results of numerical simulations are in good agreement with experimental measurements. This shows that the natural convection is formed in region above the heater rod and the water is thermally stratified in the region below the heater rod. The heat from higher region and from the heater rod is transferred to the lower region via conduction. The thermal stratification disappears and the water becomes well mixed, only after the water temperature reaches the saturation temperature and boiling starts. The developed modelling approach and obtained results provide guidelines for numerical investigations of thermal-hydraulic processes in the water pools for passive residual heat removal systems or spent nuclear fuel pools considering the concreate walls of the pool and main room above the pool.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.173-173
• /
• 2012

중력류는 두 유체 사이에서 물질의 종류 또는 온도차로 인한 밀도 차이에 의해 발생되는 유체의 흐름을 말하며, 자연환경에서는 저수지에서의 중층탁류, 사막에서의 모래 폭풍, 하구에서의 염수쐐기 등이 그 예이다. 여름철에 강우가 집중되는 우리나라는 홍수에 의한 범람, 침수 등이 발생할 뿐만 아니라 댐이나 하천, 호수 등에 흙탕물이 유입되어 중간층에 정체하는 중층탁류에 의한 피해가 자주 발생하고 있다. 과거에 중층탁류는 일시적인 현상에 불과하였으나 오늘날에는 인공적으로 형성된 정체성 호수나 하천에서 장기간 체류하여 생태계의 교란과 수질 악화를 초래하고 있다. 중층탁류에 의한 피해를 저감시키기 위해서는 본래 자연 상태의 하천이나 호수와 같이 온도분포가 비선형적으로 성층화 된 상태에서 탁수가 유입되어 중층탁류가 형성되는 과정과 흐름 특성을 이해하는 것이 선행되어야 한다. 본 연구에서는 중력류 흐름 해석을 위한 실험 장치를 개발하였다. 실험 장치 중 수조의 제원은 1200 $700{\times}300mm$이며, 수조 상부에는 가열장치를 설치하였고, 수조 하부에는 냉각장치를 설치하였다. 이들 장치들을 가동하여 깊이에 따른 자연현상에서의 비선형적인 온도 단면을 재현한다. 수리 실험방법은 깊이 방향으로 온도에 따라 비선형적으로 성층화 된 유체 내로 경사가 있는 판을 따라 수조내의 성층화된 유체의 중간 온도 혹은 밀도를 가진 유체를 침입시킨다. 침입유체는 성층화된 유체 내에서 초기에는 중력에 의해 운동량을 갖고 경사면을 따라 일정한 속도로 흘러내려 간다. 이후 바닥으로 하강 한 후 부력에 의해 상승하여 동일한 밀도층의 위치에서 횡방향으로 전파하게 된다. 이 연구에서는 이와 같이 비선형적으로 성층화된 유체 내에서의 침입 중력류의 흐름 특성을 실험적으로 관측하고 분석한다.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.17 no.4
• /
• pp.382-394
• /
• 1988

Dynamic characteristics of thermally-forced stratification process in a square enclosure with a linear temperature profile at the side walls have been investigated through flow visualization experiment and numerical analysis. The experiment was performed on air with the Rayleigh numbers of order $10^5$. A particle tracer method is used for the flow visualization and to obtain a sudden linear temperature profile at the side walls copper blocks which already have a linear temperature profile are come into contact with the thin copper plates of the test section. Immediately a meridional circulation is developed and heat transfer takes place from the wall to the interior region by circulation of fluid and finally a thermal stratification is achieved. In the numerical study, QUICK scheme for convective terms, SIMPLE algorithm for pressure correction, and the implicit method for the time marching are adopted for the integration of conservation equations. Comparison of flow visualization and numerical results shows that the developing flow patterns are very similar in dynamic nature even though there is a time lag due to the inevitable time delay in setting up a linear temperature profile. For high Rayleigh numbers, the oscillatory motion is likely to take place and stratified region is extended. However, initial temperature adjustment process is much slower than that for low Rayleigh numbers.