Hydrodynamics of two-phase flow in vertical up and down-flow across a horizontal tube bundle

Flow patterns, void fraction and friction pressure drop measurements were made for an adiabatic, vertical up-and-down, two-phase flow of air–water mixtures across a horizontal in-line, 5×20 tube bundle with a pitch-to-diameter ratio of 1.28. The flow patterns in the cross-flow zones were obtained and flow pattern maps were constructed. The data of average void fraction were less than the values predicted by a homogenous flow model and showed a strong mass velocity effect, but were well-correlated in terms of the Martinelli parameter X_tt and liquid-only Froude number Fr_LO. The two-phase friction multiplier data could be well-correlated with the Martinelli parameter.     

On the use of area-averaged void fraction and local bubble chord length entropies as two-phase flow regime indicators

In this work, the use of the area-averaged void fraction and bubble chord length entropies is introduced as flow regime indicators in two-phase flow systems. The entropy provides quantitative information about the disorder in the area-averaged void fraction or bubble chord length distributions. The CPDF (cumulative probability distribution function) of void fractions and bubble chord lengths obtained by means of impedance meters and conductivity probes are used to calculate both entropies. Entropy values for 242 flow conditions in upward two-phase flows in 25.4 and 50.8-mm pipes have been calculated. The measured conditions cover ranges from 0.13 to 5 m/s in the superficial liquid velocity j _f and ranges from 0.01 to 25 m/s in the superficial gas velocity j _g. The physical meaning of both entropies has been interpreted using the visual flow regime map information. The area-averaged void fraction and bubble chord length entropies capability as flow regime indicators have been checked with other statistical parameters and also with different input signals durations. The area-averaged void fraction and the bubble chord length entropies provide better or at least similar results than those obtained with other indicators that include more than one parameter. The entropy is capable to reduce the relevant information of the flow regimes in only one significant and useful parameter. In addition, the entropy computation time is shorter than the majority of the other indicators. The use of one parameter as input also represents faster predictions.     

The effect of magnetic field on two-phase liquid metal flow

In this paper, the basic equations of two-phase liquid metal flow in a magnetic field are derived, and specifically, two-phase liquid metal MHD flow in a rectangular channel is studied, and the expressions of velocity distribution of liquid and gas phases and the ratioK ₀ of the pressure drop in two-phase MHD flow to that in single-phase are derived. Results of calculation show that the ratioK ₀ is smaller than unity and decreases with increasing void fraction and Hartmann number because the effective electrical conductivity in the two-phase case decreases. The Project is supported by the National Natural Science Foundation of China.     

An application of the wavelet analysis technique for the objective discrimination of two-phase flow patterns

This paper presents an application of the wavelet analysis technique for two-phase flow pattern identification by using the void fraction signals obtained from a multi-channel Impedance Void Meter (IVM) in a vertical-upward air–water flow. A new method for the objective discrimination of the two-phase flow pattern has been developed to provide information regarding the local energy of void fraction signals at a given scale on the joint time–frequency diagram. The void signals are processed with Continuous Wavelet Transform (CWT) to get the local wavelet energy coefficients map on the time–frequency diagram. The effective local wavelet energy and the effective scale are then calculated. Then the criteria for flow pattern identification are, finally, obtained. A series of void fraction measurements were conducted over a wide range of air–water vertical-upward flow condition to provide an extensive database to cover several types of flow patterns. The results show that the proposed method has a high precision for characterizing different flow regimes in two-phase flow, and is considerably more promising for the online recognition of two-phase flow patterns due to the short time of data processing.     

Forced convective flow and heat transfer of upward cocurrent air–water slug flow in vertical plain and swirl tubes

This experimental study comparatively examined the two-phase flow structures, pressured drops and heat transfer performances for the cocurrent air–water slug flows in the vertical tubes with and without the spiky twisted tape insert. The two-phase flow structures in the plain and swirl tubes were imaged using the computerized high frame-rate videography with the Taylor bubble velocity measured. Superficial liquid Reynolds number (Re_L) and air-to-water mass flow ratio (AW), which were respectively in the ranges of 4000–10000 and 0.003–0.02 were selected as the controlling parameters to specify the flow condition and derive the heat transfer correlations. Tube-wise averaged void fraction and Taylor bubble velocity were well correlated by the modified drift flux models for both plain and swirl tubes at the slug flow condition. A set of selected data obtained from the plain and swirl tubes was comparatively examined to highlight the impacts of the spiky twisted tape on the air–water interfacial structure and the pressure drop and heat transfer performances. Empirical heat transfer correlations that permitted the evaluation of individual and interdependent Re_L and AW impacts on heat transfer in the developed flow regions of the plain and swirl tubes at the slug flow condition were derived.     

Heat transfer and fluid dynamics of air–water two-phase flow in micro-channels

Heat transfer, pressure drop, and void fraction were simultaneously measured for upward heated air–water non-boiling two-phase flow in 0.51 mm ID tube to investigate thermo–hydro dynamic characteristics of two-phase flow in micro-channels. At low liquid superficial velocity j_l frictional pressure drop agreed with Mishima–Hibiki’s correlation, whereas agreed with Chisholm–Laird’s correlation at relatively high j_l. Void fraction was lower than the homogeneous model and conventional empirical correlations. To interpret the decrease of void fraction with decrease of tube diameter, a relation among the void fraction, pressure gradient and tube diameter was derived. Heat transfer coefficient fairly agreed with the data for 1.03 and 2.01 mm ID tubes when j_l was relatively high. But it became lower than that for larger diameter tubes when j_l was low. Analogy between heat transfer and frictional pressure drop was proved to hold roughly for the two-phase flow in micro-channel. But satisfactory relation was not obtained under the condition of low liquid superficial velocity.     

Void fraction measurement of bubble and slug flow in a small channel using the multivision technique

Using the multivision technique, a new void fraction measurement method was developed for bubble and slug flow in a small channel. The multivision system was developed to obtain images of the two-phase flow in two perpendicular directions. The obtained images were processed—using image segmentation, image subtraction, Canny edge detection, binarization, and hole filling—to extract the phase boundaries and information about the bubble or slug parameters. With the extracted information, a new void fraction measurement model was developed and used to determine the void fraction of the two-phase flow. The proposed method was validated experimentally in horizontal and vertical channels with different inner diameters of 2.1, 2.9, and 4.0 mm. The proposed method of measuring the void fraction has better performance than the methods that use images acquired in only one direction, with a maximum absolute difference between the measured and reference values of less than 6%.     

Characteristics of developing vertical bubbly flow under normal and microgravity conditions

The axial development of the void fraction profile, interfacial area concentration and Sauter mean bubble diameter of adiabatic nitrogen-water bubbly flows in a 9 mm-diameter pipe were measured using stereo image processing under normal and microgravity conditions. The flow measurements were performed at four axial locations (axial distance from the inlet, z normalized by the pipe diameter, D, z/D = 5, 20, 40 and 60) and with various flows: superficial gas velocity of 0.00840-0.0298 m/s, and superficial liquid velocity of 0.138-0.914 m/s. The effect of gravity on radial distribution of bubbles and the axial development of two-phase flow parameters is discussed in detail based on the obtained database and visual observation. Following Serizawa-Kataoka’s phase distribution pattern criteria under normal gravity conditions, the phase distribution pattern map was developed. Similar to normal gravity two-phase flows, wall, core and intermediate void peak patterns are observed under microgravity conditions but a transition void distribution pattern is not observed in the current experimental conditions. The data obtained in the current experiment are expected to contribute to the benchmarking of CFD simulation of phase distribution pattern and interfacial area concentration in forced convective pipe flow under microgravity conditions.     

Interfacial area transport across a 90° vertical-upward elbow in air–water bubbly two-phase flow

This study develops a one-group interfacial area transport equation (IATE) for vertical-upward-to-horizontal air–water bubbly two-phase flows through a 90° elbow with a non-dimensional centerline radius of curvature of three. In order to develop the model, an extensive database is established by acquiring local two-phase flow parameters using a four-sensor conductivity probe upstream and downstream of the elbow. The data show there exist three characteristic regions in void distribution, including a bimodal-to-bimodal region, a bimodal-to-single-peaked region, and a developed horizontal flow region with void accumulated at the top of the pipe cross-section. Using the database, the preliminary dissipation length model developed by Yadav et al. (2014b) is improved by including the transition region near the exit of the elbow in addition to the dissipation region. To close the IATE model, the bubble velocity advection term and bubble interaction terms in the IATE are correlated with the parameter characterizing the “elbow-strength”. The two-phase pressure drop across the elbow is modeled using the modified Lockhart–Martinelli correlation which takes into account the minor loss effect. The closed IATE model is implemented to predict interfacial area transport in vertical-upward-to-horizontal two-phase flow. It is found that the developed model is capable of predicting interfacial area concentration with an average percent difference of less than ±6%.     

An experimental study of drag on a single tube and on a tube in an array under two-phase cross flow

In this work, the drag coefficient and the void fraction around a tube subjected to two-phase cross flow were studied for a single tube and for a tube placed in an array. The drag coefficients were determined by measuring the pressure distribution around the perimeter of the tube. Single tube drag data were taken when the tube was held both rigidly and flexibly. The test tube was made of acrylic and was 2.2 cm in diameter and 20 cm in length. In the experiments, liquid Reynolds number ranged from 430 to 21,900 for the single tube and liquid gap Reynolds number ranged from 32,900 and 61,600 for the tube placed in a triangular array. Free stream void fraction was varied from 0 to 0.4. At low Reynolds numbers, the ratio of two-phase to single-phase drag coefficient is found to be a strong function of εGr/Re2. However, at high Reynolds numbers only void fraction is the important parameter. Empirical correlations have been developed for the ratio of two-phase drag on a single tube and on a tube placed in an array.     

Liquid single- and gas—liquid two-phase flows of magnetic fluid in the presence of a high non-uniform parallel magnetic field

The effect of a non-uniform parallel high magnetic field on flow control characteristics is investigated experimentally for a magnetic fluid single-phase flow and an air—magnetic fluid two-phase flow in a vertical channel. It is found that as the magnetic field strength is increased, the friction factor of the single-phase flow increases significantly. For the two-phase flow, the friction pressure loss and the head pressure loss, which is smaller than the friction loss, are negligibly small compared with the magnetic pressure loss. In the case where air is injected 27.9d upstream from the maximum magnetic field, the air flow is blocked by the magnetic force in the entrance of the magnetic field, which leads to increases in both local void fraction and pressure drop there. In the case where air is injected 1.43d downstream from the maximum magnetic field, the air flow is accelerated, resulting in a decrease in void fraction and an increase in pressure rise. In the latter case and under the present range of experimental conditions, the magnetic pumping head reaches 0.02 MPa at the highest, and the maximum circulation flow rate reaches twice as high as non-magnetically driven flow rate.     

Experimental investigation of horizontal forced-vibration effect on air-water two-phase flow

In order to investigate the potential seismic vibrations effect on two-phase flow in an annular channel, experimental tests with air-water two-phase flow under horizontal vibrations were carried out. A low-speed eccentric-cam vibration module capable of operating at motor speed of 45–1200 rpm (f = 0.75–20 Hz) was attached to an annular channel, which was scaled down from a prototypic BWR fuel sub-channel with inner and outer diameters of 19.1 mm and 38.1 mm, respectively. The two-phase flow was operated in the ranges of 〈j_f〉 = 0.25–1.00 m/s and 〈j_g〉 = 0.03–1.46 m/s with 27 flow conditions, and the vibration amplitudes controlled by cam eccentricity (E) were designed for the range of 0.8–22.2 mm. Ring-type impedance void meters were utilized to detect the area-averaged time-averaged void fraction under stationary and vibration conditions. A systematic experimental database was built and analyzed with effective maps in terms of flow conditions (〈j_g〉-〈j_f〉) and vibration conditions (E-f and f-a), and the potential effects were expressed by regions on the maps. In the 〈j_g〉-〈j_f〉 maps, the void fraction was found to potentially decrease under vibrations in bubbly flow regime and relatively lower liquid flow conditions, which may be explained by the increase of distribution parameter. Whereas and the void fraction may increase at the region closed to bubbly-to-slug transition boundary under vibrations, which may be explained by the changes of drift velocity due to flow regime change from bubbly to slug flows. No significant change in void fraction was found in slug flow regime under the present test conditions.     

Bubbly-to-cap bubbly flow transition in a long-26 m vertical large diameter pipe at low liquid flow rate

The concurrent upward two-phase flow of air and water in a long vertical large diameter pipe with an inner diameter (D) of 200 mm and a height (z) of 26 m (z/D = 130) was investigated experimentally at low superficial liquid velocities from 0.05009 to 0.3121 m/s and the superficial gas velocities from 0.01779 to 0.5069 m/s. The resultant void fractions range from 0.03579 to 0.4059. According to the observations using a high speed video camera, the flow regimes of bubbly, developing cap bubbly and fully-developed cap bubbly flows prevailed in the flows. The developing cap bubbly flow appeared as a flow regime transition from bubbly to fully-developed cap bubble flow in the vertical large diameter pipe. The developing cap bubbly flow changes gradually and lasts for a long time period and a wide axial region in the flow direction, in contrast to a sudden transition from bubbly to slug flows in a small diameter pipe. The analysis in this study showed that the flow regime transition depends not only on the void fraction but also on the axial distance in the flow and the pipe diameter. The axial flow development brings about the transition to happen in a lower void fraction flow and the increase of pipe diameter causes the transition to happen in a higher void fraction flow. The measured void fraction showed an N-shaped axial changing manner that the void fraction increases monotonously with axial position in the bubbly flow, decreases non-monotonously with axial position in the developing cap bubbly flow, and increases monotonously again with axial position in the fully-developed cap bubbly flow. The temporary void fraction decrease phenomenon in the transition region from bubbly to cap bubbly flow can be attributed to the formation of medium to large cap bubbles and their gradual growth into the maximum size of cap bubble and/or cluster of large cap bubbles in the developing cap bubbly flow. In order to predict the N-shaped axial void fraction changing behaviors in the flow regime transition from bubbly to cap bubbly flow, the existing 12 drift flux correlation sets for large diameter pipes are reviewed and their predictabilities are studied against the present experimental data. Although some drift flux correlation sets, such as those of Clark and Flemmer (1986) and Hibiki and Ishii (2003), can predict the present experimental data with reasonable average relative deviations, no drift flux correlation set for distribution parameter and drift velocity can give a reliable prediction for the observed N-shaped axial void fraction changing behaviors in the region from bubbly to cap bubbly flow in a vertical large diameter pipe.     

Application of the ultrasonic technique and high-speed filming for the study of the structure of air–water bubbly flows

Multiphase flows are very common in industry, oftentimes involving very harsh environments and fluids. Accordingly, there is a need to determine the dispersed phase holdup using noninvasive fast responding techniques; besides, knowledge of the flow structure is essential for the assessment of the transport processes involved. The ultrasonic technique fulfills these requirements and could have the capability to provide the information required. In this paper, the potential of the ultrasonic technique for application to two-phase flows was investigated by checking acoustic attenuation data against experimental data on the void fraction and flow topology of vertical, upward, air–water bubbly flows in the zero to 15% void fraction range. The ultrasonic apparatus consisted of one emitter/receiver transducer and three other receivers at different positions along the pipe circumference; simultaneous high-speed motion pictures of the flow patterns were made at 250 and 1000 fps. The attenuation data for all sensors exhibited a systematic interrelated behavior with void fraction, thereby testifying to the capability of the ultrasonic technique to measure the dispersed phase holdup. From the motion pictures, basic gas phase structures and different flows patterns were identified that corroborated several features of the acoustic attenuation data. Finally, the acoustic wave transit time was also investigated as a function of void fraction.     

Microgravity flow regime data and analysis

To utilize the advantageous properties of two-phase flow in microgravity applications, the knowledge base of two-phase flow phenomena must be extended to include the effects of gravity. In the experiment described, data regarding the behavior of two-phase flow in a conduit under microgravity conditions (essentially zero gravity) are explored. Of particular interest, knowledge of the void fraction of the gas and liquid in a conduit is necessary to develop models for heat and mass transfer, pressure drop, and wall shear. An experiment was conducted under reduced gravity conditions to collect data by means of a capacitance void fraction sensor and high speed visual imagery. Independent parameters were varied to map the flow regime regions. These independent parameters include gas and liquid volumetric flow rates and saturation pressures. Void fraction measurements were taken at a rate of 100 Hz with six sensors at two locations along the conduit. Further, statistical parameters were developed from the void fraction measurements. Statistical parameters such as variance, signal-to-noise ratio, half height value, and linear area difference were calculated and found to have characteristics allowing flow regime identification.     

Prediction of void fraction for gas–liquid flow in horizontal,upward and downward inclined pipes using artificial neural network

In this work, the ability of artificial neural networks (ANNs) to predict void fraction of gas–liquid two–phase flow in horizontal and inclined pipes was investigated. For this purpose, an ANN model was designed and trained using a total of 301 experimental data points reported in the literature for inclination angles between –20° and +20°. Pipe inclination angle as well as superficial Reynolds number of gas (Re_sg) and liquid (Re_sl) were chosen as input parameters of different structures of multilayer perceptron (MLP) neural networks, while the corresponding void fraction was selected as their output parameter. A hyperbolic tangent sigmoid and a linear function were employed as transfer functions of hidden and output layers, respectively, and Levenberg–Marquardt back propagation algorithm was used to train the networks. By trial–and–error method, a three–layer network with 10 neurons in the hidden layer was achieved as optimal structure of the ANN which made it possible to predict the void fraction with a high accuracy. Mean absolute percent error (MAPE) of 1.81% and coefficient of determination (R²) of 0.9976 for training data and MAPE of 1.52% and R² value of 0.9948 for testing data were obtained. Also for all data, MAPE of 1.95% and R² value of 0.9972 were calculated, and 96% data were within ±5% error band. In addition, the accuracy of the proposed ANN model was compared with the predictions from 17 void fraction correlations available in the literature for different flow patterns and horizontal and inclined flows. For all cases, the proposed ANN model gave better performance than all of the studied correlations. The results confirm the very good capability of the ANNs to predict the void fractions of gas–liquid flow in inclined pipes, regardless of flow pattern. Finally, by performing interpolation using the trained network, the void fraction values for some other conditions were predicted.     

The characteristics of gas/liquid flow in large risers at high pressures

Although most of the work reported on two-phase flows are limited to small pipe diameters, two-phase flow in large risers are increasingly being encountered in the petroleum and nuclear industries. In the present work, a wire mesh sensor was employed to obtain void fraction and bubble size distribution data and visualizations of steam/water flow in a large vertical pipe (194 mm in diameter) at 46 bar. For comparison purposes, measurements were made at similar phase velocities and physical properties to a dataset for nitrogen/naphtha flow in a similar-sized riser. There exist significant differences between both sets of data. Churn-turbulent flow is observed in the present work instead of slug flow, and this differs from the intermittent and semi-annular flow patterns reported for nitrogen/naphtha data. The mean void fraction of the nitrogen/naphtha data is higher than that of the present steam/water data due to the differences in purity in the liquid phases. Furthermore, core peak distributions are observed for the present work in contrast to the flatter profiles deduced for the nitrogen/naphtha using a power law relationship.     

Bubble dynamics and interactions with a pair of micro pillars in tandem

This study investigates flow patterns and bubble dynamics of two-phase flow around two 100 μm diameter circular pillars in tandem, which were entrenched inside a horizontal micro channel. Bubble velocity, trajectory, size, and void fraction were measured using a high speed camera and analyzed using a particle tracking velocimetry method. A range of gas and liquid superficial velocities were tested, resulting in different bubbly flow patterns, which were consistent with previous studies. These flow patterns were altered as they interacted with the pillars. Depending on the relative transverse location of bubbles to the pillars, and through bubble–bubble interaction, the flow sometimes returned to its original state. It was also determined that the pillars altered both the bubble trajectory and void fraction, especially in the pillars region.     

Flow characteristics in hydraulically equilibrium two-phase flows in a vertical 2 × 3 rod bundle channel

In order to increase data on two-phase flow distribution in a multi-subchannel system, being similar to a rod bundle, experiments have been carried out using water and air at ambient pressure and temperature as the working fluids and a newly constructed 2 × 3 rod bundle channel as the test channel. The channel contained six rods in rectangular array and two-kinds of six subchannels, simulating a BWR fuel rod bundle. Experimental data on flow distribution and pressure drop along each subchannel axis were obtained in various single- and two-phase flows under a hydraulic equilibrium flow condition. From the measured pressure drop in the single-phase flow, friction factor data in each subchannel were obtained. The two-phase pressure drop data were compared with calculations by a simple, one-dimensional, one-pressure two-fluid model. In addition, Taylor bubble velocity in each subchannel in slug-churn flows was measured with a double needle contact probe. Using the bubble velocity data, we obtained a subchannel void fraction in each subchannel, and discussed a relationship of the subchannel void fractions between two different subchannels. Results of such experiments and discussions are presented in this paper.     

气液两相管流水力瞬变的数值计算及实验研究

由气液两相管流的基本方程出发,通过引入矢通量分裂,对传统的特征线差分做了较大的改进,形成了基于矢通量分裂的特征线差分解法。该法首先将控制方程组的特征值分解成正、负两部分,进而将控制方程中的矢通量雅可比矩阵分裂成两个亚矢量矩阵,对其按各自的迎风格式差分,从而建立了稳定的差分求解格式。该计算法适合于计算声速变化较大且计及液流速度的气液管流的瞬变。计算求解得到的各种不同初始空隙比的压力变化曲线、声速曲线、波速变化曲线、空隙比变化曲线及气体释放影响曲线,通过与不同初始空隙比时气液管流水力瞬变的实验结果对比分析,结果表明两者吻合较好,说明本文方法对于低空隙比的气液两相管流具有较普遍的适用性。