Analysis of Turbulence Structure and Void Fraction Distribution in Gas-Liquid Two-Phase Flow Under Bubbly and Slug Flow Regime

Author(s):  
Isao Kataoka ◽  
Kenji Yoshida ◽  
Tsutomu Ikeno ◽  
Tatsuya Sasakawa ◽  
Koichi Kondo

Accurate analyses of turbulence structure and void fraction distribution are quite important in designing and safety evaluation of various industrial equipments using gas-liquid two-phase flow such as nuclear reactor, etc. Using turbulence model of two-phase flow and models of bubble behaviors in bubble flow and slug flow, systematic analyses of distributions of void fraction, averaged velocity and turbulent velocity were carried out and compared with experimental data. In bubbly flow, diffusion of bubble and lift force are dominant in determining void fraction distribution. On the other hand, in slug flow, large scale turbulence eddies which convey bubbles into the center of flow passage are important in determining void fraction distribution. In turbulence model, one equation turbulence model is used with turbulence generation and turbulence dissipation due to bubbles. Mixing length due to bubble is also modeled. Using these bubble behavior models and turbulence models, systematic predictions were carried out for void distributions and turbulence distributions for wide range of flow conditions of two phase flow including bubbly and slug flow. The results of predictions were compared with experimental data in round straight tube with successful agreement. In particular, concave void distributions in bubbly flow and convex distribution in slug flow were well predicted based on the present model.

Author(s):  
Hiroyuki Yoshida ◽  
Taku Nagatake ◽  
Kazuyuki Takase ◽  
Akiko Kaneko ◽  
Hideaki Monji ◽  
...  

An earthquake is one of the most serious phenomena to consider for the safety of a nuclear reactor in Japan. Therefore, structural safety of nuclear reactors has been studied and nuclear reactors were contracting with structural safety for a big earthquake. However, it is not enough for safety operation of nuclear reactors because thermal-fluid safety is not confirmed under the earthquake. For instance, behavior of gas-liquid two-phase flow is unknown in seismic conditions. Especially, fluctuation of void fraction is an important factor for the safety operation of the nuclear reactor. In previous work, fluctuation of void faction in bubbly flow was studied experimentally and theoretically to investigate the stability of the bubbly flow. In such studies, flow rate or void fraction fluctuations were given to the steady bubbly flow. In case of the earthquake, the fluctuation is not only the flow rate, but also a body force on the two-phase flow and shear force through the pipe wall. Interactions of gas and liquid through their interface also act on the behavior of the two-phase flow. The fluctuation of the void fraction is not clear for such complicated situation during the earthquake. Therefore, the behavior of gas-liquid two-phase flow is investigated experimentally and numerically in a series of studies. In this study, to develop the predictive technology of two-phase flow dynamics under earthquake acceleration, a detailed two-phase flow simulation code with an advanced interface tracking method TPFIT (Two-Phase Flow simulation code with Interface Tracking) was expanded to two-phase flow simulation in seismic conditions. In a previous study, we performed a numerical simulation of a two-phase bubbly flow in a horizontal pipe and a vertical bubble motion in a water tank in seismic conditions. And it was confirmed that the modified TPFIT can be applicable to the bubbly flow in seismic conditions. In this paper, the two-phase bubbly flow in a simulated single-subchannel excited by oscillation acceleration was simulated by using the expanded TPFIT. A calculation domain used in this simulation was a simplified subchannel in a BWR core. And time-series of void fraction distributions were evaluated based on predicted bubble distributions. When no oscillation acceleration was added, void fraction concentrated in a region near the wall. When oscillation acceleration was added, void fraction distribution was changed by time. And coalesces of bubbles occurred in the numerical simulation, and bubbles with relatively large diameter were observed. In the results, complicated void fraction distribution was observed, because the response of void fraction distribution on the oscillation acceleration was dependent on not only imposed acceleration, but also the bubble diameter.


Author(s):  
Kenichi Katono ◽  
Jun Nukaga ◽  
Takuji Nagayoshi ◽  
Kenichi Yasuda

We have been developing a void fraction distribution measurement technique using the three-dimensional (3D) time-averaged X-ray CT (computed tomography) system to understand two-phase flow behavior inside a fuel assembly for BWR (boiling water reactor) thermal hydraulic conditions of 7.2 MPa and 288 °C. Unlike CT images of a normal standstill object, we can obtain 3D CT images that are reconstructed from time-averaged X-ray projection data of the intermittent two-phase flow. We measured the 3D void fraction distribution in a vertical square (5 × 5) rod array that simulated a BWR fuel assembly in the air-water test. From the 3D time-averaged CT images, we confirmed that the void fraction at the center part of the channel box was higher than that near the channel box wall, and the local void fraction at the central region of a subchannel was higher than that at the gap region of the subchannel. A comparison of the volume-averaged void fractions evaluated by the developed X-ray CT system with those evaluated by a differential pressure transducer in a void fraction range from 0.05 to 0.40 showed satisfactory agreement within a difference of 0.03.


Sign in / Sign up

Export Citation Format

Share Document