Drift-flux modeling of void fraction for boiling two-phase flow in a tight-lattice rod bundle

Author(s):  
Takashi Hibiki ◽  
Shuichiro Miwa ◽  
Kenichi Katono
Author(s):  
Hiroyuki Yoshida ◽  
Takeharu Misawa ◽  
Kazuyuki Takase

Two-fluid model can simulate two phase flow less computational cost than inter-face tracking method and particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D, which adopts boundary fitted coordinate system in order to simulate complex shape channel flow. In this paper, boiling two-phase flow analysis in a tight lattice rod bundle is performed by ACE-3D code. The parallel computation using 126CPUs is applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. At height z = 0.5 m, void fraction in the gap region is higher in comparison with that in center region of the subchannel. However, at height of z = 1.1m, higher void fraction distribution exists in center region of the subchannel in comparison with the gap region. The tendency of void fraction to concentrate in the gap region at vicinity of boiling starting point, and to move into subchannel as water goes through rod bundle, is qualitatively agreement with the measurement results by neutron radiography. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight lattice rod bundle with no lift force model (neglecting lift force acting on bubbles) is also performed. From the comparison of numerical results, it is concluded that the effects of lift force model are not so large on overall void fraction distribution in tight lattice rod bundle. However, higher void fraction distribution in center region of the subchannel was not observed in this simulation. It is concluded that the lift force model is important for local void fraction distribution in rod bundles.


2018 ◽  
Vol 115 ◽  
pp. 480-486 ◽  
Author(s):  
Bin Yu ◽  
Wenxiong Zhou ◽  
Liangming Pan ◽  
Hang Liu ◽  
Quanyao Ren ◽  
...  

2019 ◽  
Vol 109 ◽  
pp. 109881 ◽  
Author(s):  
Miao Gui ◽  
Zhaohui Liu ◽  
Bo Liao ◽  
Teng Wang ◽  
Yong Wang ◽  
...  

Author(s):  
Xu Han ◽  
Xiuzhong Shen ◽  
Toshihiro Yamamoto ◽  
Ken Nakajima ◽  
Takashi Hibiki

2012 ◽  
Vol 40 ◽  
pp. 166-177 ◽  
Author(s):  
Shao-Wen Chen ◽  
Yang Liu ◽  
Takashi Hibiki ◽  
Mamoru Ishii ◽  
Yoshitaka Yoshida ◽  
...  

Author(s):  
Quanyao Ren ◽  
Liangming Pan ◽  
Wenxiong Zhou ◽  
Tingpu Ye ◽  
Hang Liu ◽  
...  

In order to simulate the transfer of mass, momentum and energy in the gas-liquid two-phase flow system, tremendous work focused on the phenomenon, mechanisms and models for two-phase flow in different channels, such as circular pipe, rectangular channel, rod bundle and annulus. Drift-flux model is one of the widely used models for its simplicity and good accuracy, especially for the reactor safety analysis codes (RELAP5 and TRAC et al.) and sub-channel analysis code (COBRA, SILFEED and NASCA et al.). Most of the adopted drift-flux models in these codes were developed based on the void fraction measured in pipe and annulus, which were different with the actual nuclear reactor. Although some drift-flux models were developed for rod bundles, they were based on the void fraction on the whole cross-section not in subchannel in rod bundles due to the lack of effective measuring methods. A novel sub-channel impedance void meter (SCIVM) has been developed to measure the void fraction in sub-channel of 5 × 5 rod bundles, which is adopted to evaluate these existing drift-flux models for rod bundles. By comparison, the values of drift-flux parameters have large differences among different correlations, which are suggested to be reconsidered. Based on the experimental data and physical laws, Lellouche-Zolotar and Chexal-Lellouche correlations show a better performance for drift velocity. If the predicting error of void fraction is the only concerned parameter, Chen-Liu, Ishizuka-Inoue and Chexal-Lellouche correlations are recommended for averaged relative error less than 30%. More experiments are suggested to focus on the distribution parameter and drift velocity through their definition.


Author(s):  
Y. S. Lim ◽  
Simon C. M. Yu

Single phase and two phase flow characteristics in micro-sized glass tubes with i.d. (inner diameter) of 300 and 500 μm have been examined experimentally. Single phase pressure drop measurements are found generally in good agreement with Poiseulle flow theory. Transitional flow is found to start earlier at Reynolds number about 1600 as compared to the onset of transitional flow at Reynolds number of 2300 for macro-scale tubes. In addition, these glass tubes are employed for the investigation of adiabatic two phase flow characteristic by introducing gas phase via a stainless steel tube inserted at the center of the glass tube. Real time flow visualization obtained under the same flow condition are analyzed by both cross sectional void fraction (one dimensional drift flux model) and volumetric void fraction (image processing method). The analysis shows that the void fraction estimated by drift flux model (DFM) agrees with homogeneous correlation (α = β) and Armand correlation (α = 0.833β). However image processing method seems to reveal that the slip ratio for the two phase flow is more significant and that the void fraction results are clustering between slip ratio of 3 and 7. Additionally, two phase frictional pressure losses are compared with the convention correlation for macro-sized tube (Lockhart-Martinelli model). It is found that measurements of the two phase frictional pressure drop can serve as a flow map to predict the flow patterns when the flow in the channel is not transparent.


2021 ◽  
Vol 39 (4) ◽  
pp. 1343-1350
Author(s):  
Tat Thang Nguyen

The drift-flux model is widely used in study, calculation and design of two-phase flow. It is a highly efficient model that requires little computation resources. In the model, accurate calculation of the distribution parameter C0 and the drift velocity Vgj is a critically important factor. The calculation requires simultaneously measured data of phase velocity and void fraction distributions or profiles. By using currently widely used methods for two-phase flow measurement, satisfying the requirement is highly difficult. This paper presents novel results of simultaneous measurement of the phase velocity and void fraction profiles in a vertical round tube of 50 mm inner diameter. A combination measurement method has been developed. It comprises the multiwave Ultrasonic Velocity Profile (multiwave UVP) method and the Wire Mesh Tomography (WMT). Based on the measured data, C0 and Vgj have been calculated. They have been compared with those of the published experimental data and correlations. Analyses of the measured data have been carried out. For the first time, the analysis results reveal the variation of C0 and Vgj in the measured flow conditions. More importantly, the data obtained are also useful for the development and validation of the computational codes for two-phase flow.


Author(s):  
Pei-Syuan Ruan ◽  
Shao-Wen Chen ◽  
Min-Song Lin ◽  
Jin-Der Lee ◽  
Jong-Rong Wang

Abstract This paper presents the experimental results and analyses of the structure velocity of air-water two-phase flow in a 3 × 3 rod bundle channel. A total of 56 flow conditions were tested and investigated for rod-gap, sub-channel, rod-wall and global regions of rod bundle geometry. The experimental tests were carried out under bubbly and cap-bubbly flow regimes with superficial gas and liquid velocities of 0–1 m/s and 1–1.7 m/s, respectively. The conductivity probes were set at different heights to measure the global and local void fractions. The structure velocity of air-water two-phase flow is the average bubble velocity calculated by the method in this study. The structure velocity were determined by utilizing the cross-correlation technique to analyze the time lags of the bubbles passing through the conductivity probes. The results of this study indicated that the structure velocity may increase with increasing superficial gas and liquid velocities. In low superficial gas velocity region, the structure velocity may first slightly increase and follow by a sudden jump which appear in most regions. After the sudden jump, the structure velocity may keep increasing mildly. The present structure velocities have been compared with the area-averaged gas velocities predicted by the drift flux model, and it appears that most structure velocities show a good agreement with the averaged gas velocities from the drift flux model after the jump.


Sign in / Sign up

Export Citation Format

Share Document