Damping of tubes due to internal two-phase flow

Two-phase internal flow is present in many piping system components. Although two-phase damping is known to be a significant constituent of the total damping, the energy dissipation mechanisms that govern two-phase damping are not well understood. In this paper, damping of three different clamped–clamped tubes subjected to two-phase air–water internal axial flow is investigated. Experimental data are reported, showing a strong dependence of two-phase damping on void fraction, flow velocity and flow regime. Data-points plotted on two-phase flow pattern maps indicate that damping is greater in a bubbly flow regime. The two-phase damping ratio reaches a maximum value at the highest void fraction before the transition to a churn flow regime. An analytical model that relates the two-phase damping ratio to the interface surface area is proposed. The model is based on rigid spherical bubbles in cubic elementary flow volumes. The analytical results are well correlated with the experiments.     

Scaling of damping induced by bubbly flow across tubes

The damping of tubes subjected to two-phase air–water bubbly cross-flow is investigated with the use of an experimental database from several authors. A new definition of damping in stagnant flow is proposed using an extrapolation of the measured values at low dimensionless flow velocities. This approach yields values of damping substantially lower than those currently defined in the literature. They are found to vary continuously with void fraction, within the bubbly flow regime. These data are used to compare several models of the equivalent viscosity of a two-phase mixture. The effect of the flow velocity is then analysed up to fluidelastic instability. It is observed that, using scaling factors based on the characteristics of the liquid phase, fluidelastic effects of bubbly flows are closely related to those known in single-phase flows.     

SOME ASPECTS OF HEAT-EXCHANGER TUBE DAMPING IN TWO-PHASE MIXTURES

Two-phase flow exists in many shell-and-tube heat exchangers. To avoid excessive flow-induced vibration, it is necessary to have information on damping. A simple experiment was undertaken to study the effect of several parameters such as void fraction, surface tension, tube frequency and confinement on damping in two-phase mixtures. The experiment consisted of a cantilevered tube immersed in a two-phase mixture generated by bubbling air through water. It is found that void fraction and flow regime are dominant parameters. Surface tension is also important. The results are presented in detail in the paper.     

Two-phase damping for internal flow: Physical mechanism and effect of excitation parameters

Two-phase flow induced-vibration is a major concern for the nuclear industry. This paper provides experimental data on two-phase damping that is crucial to predict vibration effects in steam generators. An original test section consisting of a tube subjected to internal two-phase flow was built. The tube is supported by linear bearings and compression springs allowing it to slide in the direction transverse to the flow. An excitation system provides external sinusoidal force. The frequency and magnitude of the force are controlled through extension springs. Damping is extracted from the frequency response function of the system. It is found that two-phase damping depends on flow pattern and is fairly proportional to volumetric fraction for bubbly flow. Measurements are completed by the processing of high-speed videos which allow to characterize the transverse relative motion of the gas phase with respect to the tube for bubbly flow. It is shown that the bubble drag forces play a significant role in the dissipation mechanism of two-phase damping.     

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.     

Study of two-phase flows in reduced gravity using ground based experiments

Experimental studies have been carried out to support the development of a framework of the two-fluid model along with an interfacial area transport equation applicable to reduced gravity two-phase flows. The experimental study simulates the reduced gravity condition in ground based facilities by using two immiscible liquids of similar density namely, water as the continuous phase and Therminol 59^® as the dispersed phase. We have acquired a total of eleven data sets in the bubbly flow and bubbly to slug flow transition regimes. These flow conditions have area-averaged void (volume) fractions ranging from 3 to 30% and channel Reynolds number for the continuous phase between 2,900 and 8,800. Flow visualization has been performed and a flow regime map developed which is compared with relevant bubbly to slug flow regime transition criteria. The comparison shows that the transition boundary is well predicted by the criterion based on critical void fraction. The value of the critical void fraction at transition was experimentally determined to be approximately 25%. In addition, important two-phase flow local parameters, including the void fraction, interfacial area concentration, droplet number frequency and droplet velocity, have been acquired at two axial locations using state-of-the-art multi-sensor conductivity probe. The radial profiles and axial development of the two-phase flow parameters show that the coalescence mechanism is enhanced by either increasing the continuous or dispersed phase Reynolds number. Evidence of turbulence induced particle interaction mechanism is highlighted. The data presented in this paper clearly show the marked differences in terms of bubble (droplet) size, phase distribution and phase interaction in two-phase flow between normal and reduced gravity conditions.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

Two-phase flow pattern identification in a tube bundle based on void fraction and pressure measurements,with emphasis on churn flow

Two-phase flow are frequently encountered in the industry. In particular, in steam generators of nuclear plants, water is heated so that at the top of the generator an important fraction of water flows as vapor. In this upper part, a rising co-current two phase flow transverse to the tube bundle takes place. Fluids exert significant forces on the tubes in this area which highly depend on the two-phase flow pattern. Thus, as a prerequisite, it is essential to gather information on the flow conditions associated with the different two-phase flow patterns, which can be bubbly, intermittent, or annular. Then we must analyze the potentially dangerous flow patterns. This paper presents an experimental campaign aimed at characterizing those flow patterns for a rising co-current transverse flow in a tube bundle representative of the geometry in a steam generator. A new methodology based on the understanding of key contributions to vertical two-phase flow pattern maps in tube bundles is proposed that leads to a more complete flow pattern map. Finally, the paper focuses on the churn flow, which is the flow pattern for which significant pressure fluctuations occur. For this pattern, important damages could be expected on the tubes of a steam generator. Different kinds of pressure fluctuations are observed at different frequencies depending on the flow rates and the location in the test section.     

MODELLING TWO-PHASE FLOW-EXCITED DAMPING AND FLUIDELASTIC INSTABILITY IN TUBE ARRAYS

This paper reports the results of an experimental study of the flow-induced vibration of a heat exchanger tube array subjected to two-phase cross-flow of refrigerant 11. The primary concern of the research was to develop a methodology for predicting the critical flow velocities for fluidelastic instability which better characterize the physics of two-phase flows. A new method is proposed for calculating the average fluid density and equivalent flow velocity of the two-phase fluid, using a newly developed void fraction model to account for the difference in velocity between the gas and liquid phases. Additionally, damping measurements in two-phase flow were made and compared with the data of other researchers who used a variety of modelling fluids. The results show that the two-phase damping follows a similar trend with respect to homogeneous void fraction, and when normalized, agree well with the data in the literature. The fluidelastic threshold data of several researchers who used a variety of fluids, is re-examined using the proposed void fraction model, and the results show a remarkable change in trend with flow regime. The data corresponding to the bubbly flow regime shows no significant deviation from the trend established by Connors' theory. However, the data corresponding to the intermittent flow regime show a significant decrease in stability which is nearly independent of the mass-damping parameter. It is believed that the velocity fluctuations that are inherent in the intermittent flow regime are responsible for tripping the instability, causing lower than expected stability of the bundle.     

One-group interfacial area transport equation and its sink and source terms in narrow rectangular channel

The characteristics of two-phase flow in a narrow rectangular channel are expected to be different from those in other channel geometries, because of the significant restriction of the bubble shape which, consequently, may affect the heat removal by boiling under various operating conditions. The objective of this study is to develop an interfacial area transport equation with the sink and source terms being properly modeled for the gas–liquid two-phase flow in a narrow rectangular channel. By taking into account the crushed characteristics of the bubbles a new one-group interfacial area transport equation was derived for the two-phase flow in a narrow rectangular channel. The random collisions between bubbles and the impacts of turbulent eddies with bubbles were modeled for the bubble coalescence and breakup respectively in the two-phase flow in a narrow rectangular channel. The newly-developed one-group interfacial area transport equation with the derived sink and source terms was evaluated by using the area-averaged flow parameters of vertical upwardly-moving adiabatic air–water two-phase flows measured in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm. The flow conditions of the data set covered spherical bubbly, crushed pancake bubbly, crushed cap-bubbly and crushed slug flow regimes and their superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. Good agreement with the average relative deviation of 9.98% was obtained between the predicted and measured interfacial area concentrations in this study.     

On two-phase flow patterns and transition criteria in aqueous methanol and CO2 mixtures in adiabatic,rectangular microchannels

This paper presents flow map investigations of adiabatic two-phase flow in square cross-sectioned, 200 μm deep microchannels fabricated in silicon, employing laser induced fluorescence microscopy. The influence of surface tension and nozzle geometry on the flow pattern transition was investigated using two nozzle widths (orifices of 30 μm and 50 μm, respectively) and methanol–water solutions with CO₂ as the gas phase. It was found and quantified that smaller nozzle geometries and smaller liquid surface tension promote the propagation of capillary gas bubbles at lower superficial gas and liquid velocities. Within the measurement domain of superficial gas (0.01–0.625 m/s) and liquid (0.0005–0.5000 m/s) velocities, we observed dispersed bubbly, regularly ordered bubbly, wedging, slug and annular flows, thus extending the experimental knowledge base to smaller superficial liquid velocities by almost two orders of magnitude. With the help of the flow maps presented herein, we were able to characterize the observed regularly ordered bubbly flow as the transition regime between dispersed bubbly and wedging flow. The results of the present investigation are of direct relevance to the operation of small-scale direct methanol fuel cells.     

Effect of diameter on two-phase pressure drop in narrow tubes

The effect of tube diameter on two-phase frictional pressure drop was investigated in circular tubes with inner diameters of 0.6, 1.2, 1.7, 2.6 and 3.4 mm using air and water. The gas and liquid superficial velocity ranges were 0.01-50 m/s and 0.01-3 m/s, respectively. The gas and liquid flow rates were measured and the two-phase flow pattern images were recorded using high-speed CMOS camera. Unique flow patterns were observed for smaller tube diameters. Pressure drop was measured and compared with various existing models such as homogeneous model and Lockhart-Martinelli model. It appears that the dominant effect of surface tension shrinking the flow stratification in the annular regime is important. It was found that existing models are inadequate in predicting the pressure drop for all the flow regimes visualized. Based on the analysis of present experimental frictional pressure drop data a correlation is proposed for predicting Chisholm parameter “C” in slug annular flow pattern. For all other flow regimes Chisholm’s original correlation appears to be adequate except the bubbly flow regime where homogeneous model works well. The modification results in overall mean deviation of pressure drop within 25% for all tube diameters considered. This approach of flow regime based modification of liquid gas interaction parameter appears to be the key to pressure drop prediction in narrow tubes.     

Flow regime transitions due to cavitation in the flow through an orifice

This paper presents both experimental and theoretical aspects of the flow regime transitions caused by cavitation when water is passing through an orifice. Cavitation inception marks the transition from single-phase to two-phase bubbly flow; choked cavitation marks the transition from two-phase bubbly flow to two-phase annular jet flow.

It has been found that the inception of cavitation does not necessarily require that the minimum static pressure at the vena contracta downstream of the orifice, be equal to the vapour pressure liquid. In fact, it is well above the vapour pressure at the point of inception. The cavitation number [σ = (P₃ − P_v)/(0.5 pV²); here P₃ is the downstream pressure, P_v is the vapour pressure of the liquid, ρ is the density of the liquid and V is the average liquid velocity at the orifice] at inception is independent of the liquid velocity but strongly dependent on the size of the geometry. Choked cavitation occurs when this minimum pressure approaches the vapour pressure. The cavitation number at the choked condition is a function of the ratio of the orifice diameter (d) to the pipe diameter (D) only. When super cavitation occurs, the dimensionless jet length [L/(D - d); where L is the dimensional length of the jet] can be correlated by using the cavitation number. The vaporization rate of the surface of the liquid jet in super cavitation has been evaluated based on the experiments.

Experiments have also been conducted in which air was deliberately introduced at the vena contracta to simulate the flow regime transition at choked cavitation. Correlations have been obtained to calculate the critical air flow rate required to cause the flow regime transition. By drawing an analogy with choked cavitation, where the air flow rate required to cause the transition is zero, the vapour and released gas flow rate can be predicted.     

An experimental investigation of two-phase air–water flow through a horizontal circular micro-channel

This paper is a continuation of the authors’ previous work. Two-phase air–water flow experiments are performed in a horizontal circular micro-channel. The test section is made of a fused silica tube with an inner diameter of 0.15 mm and a length of 104 mm. The flow phenomena, which are liquid/unstable annular alternating flow (LUAAF), liquid/annular alternating flow (LAAF), and annular flow, are observed and recorded by a high-speed camera mounted together with a stereozoom microscope. A flow pattern map is presented in terms of the phase superficial velocities and is compared with those of other researchers obtained from different working fluids. Image analysis is performed to determine the void fraction, which increases non-linearly with increasing volumetric quality. It is revealed that the two-phase frictional multiplier data show a dependence on flow pattern rather than mass flux. Based on the present data, a new pressure drop correlation is proposed for practical applications. According to the present study, in general the data for the two-phase air–water flow characteristics are found to comply with those of working fluids other than air–water mixtures.     

Experimental study of non-uniform inlet conditions and three-dimensional effects of vertical air–water two-phase flow in a narrow rectangular duct

To study the three-dimensional interfacial structure development in vertical two-phase flow, air–water upflow experiments were performed in a rectangular duct. Various non-uniform two-phase profiles were created by injecting air from individually controlled spargers at the duct inlet into uniformly injected water flow. A four-sensor conductivity probe was used to measure local void fraction, interfacial area concentration, bubble velocity and Sauter mean diameter at three axial locations to record the development of two-phase parameters. Experimental results showed that the lateral development across the wider dimension of the duct was significant with a non-uniform inlet profile when compared to a uniform inlet profile. It is postulated that lift, wall and turbulent forces are the major contributors to the lateral distribution of the two-phase interfacial structures making this an useful experiment for benchmarking three-dimensional two-fluid models. In examining the interfacial area, the shearing-off of group 1 bubbles (defined as the smaller spherical and distorted bubbles) from the skirt region of group 2 bubbles (defined as the bigger cap and churn bubbles), the coalescence of group 2 bubbles due to wake entrainment, and random collision are the major source and sink mechanisms of interfacial area concentration.     

Mass flow rate measurements in gas–liquid flows by means of a venturi or orifice plate coupled to a void fraction sensor

Two-phase flow measurements were carried out using a resistive void fraction meter coupled to a venturi or orifice plate. The measurement system used to estimate the liquid and gas mass flow rates was evaluated using an air–water experimental facility. Experiments included upward vertical and horizontal flow, annular, bubbly, churn and slug patterns, void fraction ranging from 2% to 85%, water flow rate up to 4000 kg/h, air flow rate up to 50 kg/h, and quality up to almost 10%. The fractional root mean square (RMS) deviation of the two-phase mass flow rate in upward vertical flow through a venturi plate is 6.8% using the correlation of Chisholm (D. Chisholm, Pressure gradients during the flow of incompressible two-phase mixtures through pipes, venturis and orifice plates, British Chemical Engineering 12 (9) (1967) 454–457). For the orifice plate, the RMS deviation of the vertical flow is 5.5% using the correlation of Zhang et al. (H.J. Zhang, W.T. Yue, Z.Y. Huang, Investigation of oil–air two-phase mass flow rate measurement using venturi and void fraction sensor, Journal of Zhejiang University Science 6A (6) (2005) 601–606). The results show that the flow direction has no significant influence on the meters in relation to the pressure drop in the experimental operation range. Quality and slip ratio analyses were also performed. The results show a mean slip ratio lower than 1.1, when bubbly and slug flow patterns are encountered for mean void fractions lower than 70%.     

Characterization of diabatic two-phase flows in microchannels: Flow parameter results for R-134a in a 0.5 mm channel

An optical measurement method for two-phase flow pattern characterization in microtubes has been utilized to determine the frequency of bubbles generated in a microevaporator, the coalescence rates of these bubbles and their length distribution as well as their mean velocity. The tests were run in a 0.5 mm glass channel using saturated R-134a at 30 °C (7.7 bar). The optical technique uses two laser diodes and photodiodes to measure these parameters and to also identify the flow regimes and their transitions. Four flow patterns (bubbly flow, slug flow, semi-annular flow and annular flow) with their transitions were detected and observed also by high speed video. It was also possible to characterize bubble coalescence rates, which were observed here to be an important phenomena controlling the flow pattern transition in microchannels. Two types of coalescence occurred depending on the presence of small bubbles or not. The two-phase flow pattern transitions observed did not compare well to a leading macroscale flow map for refrigerants nor to a microscale map for air–water flows. Time averaged cross-sectional void fractions were also calculated indirectly from the mean two-phase vapor velocities and compared reasonably well to homogeneous values.     

Flow patterns and heat transfer in vertical upward air–water flow with surfactant

Flow patterns, the pressure drag reduction and the heat transfer in a vertical upward air–water flow with the surfactant having negligible environmental impact were studied experimentally in a tube of 2.5 cm in diameter. Visual observations showed that gas bubbles in the air–water solution with surfactant are smaller in size but much larger in number than in pure air–water mixture, at the all flow regimes. The transition lines in the flow regime map for the solution of air–water mixture with surfactant of the 300 ppm concentration are mainly consistent with the experimental data obtained in clear air–water mixture. An additive of surfactant to two-phase flow reduces the total pressure drop and decrease heat transfer, especially in the churn flow regime.     

On the development of an objective flow regime indicator

A high intensity dual beam X-ray system was designed and constructed to make chordal-average void fraction measurements. This X-ray system employed a DC excited tube filament, full wave rectification and high voltage filtering to produce a stable photon source. The large photon flux produced by the X-ray system allowed analog linearization of the signal.A series of chordal-average void fraction measurements were made and used to generate probability density functions (PDF) and power spectral density (PSD) functions. The first four moments associated with these distributions were studied as possible flow regime indicators.The moments of the PDF indicated the various flow regime transitions. The moments of the PSD also show some flow regime transition information, but were sensitive to liquid phase velocity. The PDF variance, or second moment about the mean, was found to be the best indicator of flow regime. A variance of 0.04 appear to adequately discriminate between the bubbly, slug and annular flow regimes for low pressure air/water flow in a 2.54 cm I.D. vertical tube.