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.     

Updated results on hydrodynamic mass and damping estimations in tube bundles under two-phase crossflow

Flow-Induced Vibration (FIV) is the most critical dynamic issue in the design of shell-and-tube heat exchangers. This fluid-structure phenomenon may generate high amplitude vibration of tubes or structural parts, which leads to fretting wear between the tubes and supports, noise or even fatigue failure of internal components. The study of this phenomenon is more challenging if considered that two-phase crossflow exists in many shell-and-tube heat exchangers. In this framework, the analysis of the influence of void fraction and flow patterns on FIV is of particular interest. In fact, void fraction and flow patterns do affect the dynamic parameters involved in tube vibration and, hence, the current vibration mechanism. However, in spite of the importance of devices subjected to two-phase flow, FIV under these conditions have not been entirely understood. In this paper, the results of an extensive experimental campaign, aiming at validating the flow pattern maps found in open literature, are presented. For this purpose, a normal triangular (transversal pitch per diameter ratio of 1.26) tube bundle subjected to two-phase air - water vertical upward crossflow is used. Structural sensors are used to measure the tube dynamic responses and estimate parameters such as hydrodynamic mass and damping ratios, which are strongly dependent on flow conditions. Theoretical models and data previously published are compared with the present experimental results, showing good agreement.     

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.     

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.     

Fluidelastic instability study in a rotated triangular tube array subject to two-phase cross-flow. Part I: Fluid force measurements and time delay extraction

Fluidelastic instability is a key issue in steam generator tube bundles subjected to cross-flow. The extension to two-phase flow of the existing theoretical models, developed and tested mostly for single phase flow, is investigated in this paper. The time delay is one of the key parameter for modeling fluidelastic instability, especially the damping controlled mechanism. The direct measurement of the time delay between the tube motion and the fluid force faces certain difficulties in two-phase flow since the high turbulence due to the interaction of the two components of the flow may increase the randomness of the measured force. To overcome this difficulty, an innovative method for extracting the time delay inherent to the quasi-steady model for fluidelastic instability is proposed in this study.Firstly, experimental measurements of unsteady and quasi-static fluid forces (in the lift direction) acting on a tube subjected to air–water two-phase flow were conducted. The unsteady fluid forces were measured by exciting the tube using a linear motor. These forces were measured for a wide range of void fractions, flow velocities and excitation frequencies. The experimental results showed that the unsteady fluid forces could be represented as single valued function of the reduced flow velocity. It was also found that for a given frequency, the unsteady fluid force phase was weakly dependent on the void fraction for the range of flow velocities considered.The time delay was determined by equating the unsteady fluid forces with the quasi-steady forces. The results given by this innovative method of measuring the time delay in two-phase flow were consistent with theoretical expectations. The time delay could be expressed as a linear function of the convection time and the time delay parameter was determined for void fractions ranging from 60% to 90%.     

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.     

垂直向上气液两相流中两相斯托拉赫数的研究

试验研究了三角形、T形两种形状4种规格的物体,在垂直上升气液两相流中,发生气液两相涡街时,气液两相斯托拉赫数的变化规律,在测得大量数据的基础上,得出了发生气液两相涡街时,气液两相斯托拉赫数的通用关系式,研究表明,气液两相斯托拉赫数在两相工况下为一变数,其值与来流截面含气率、涡街发生体形状和特征尺寸、来流方向等因素有关,应用此关系式,根据测得的两相涡街频率可将涡街发生体作为测量两相流流量与组分的测量元件。     

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.     

Analysis of frictionless two-phase flow in tubes under centrifugal acceleration and minimum heating for choked flow

To determine the void fraction in a tube of a rotating heat exchanger, an analytical investigation was undertaken to model frictionless two-phase flow boiling. Steady, one-dimensional separated two-phase conservation equations in differential form, were first applied to a stationary system. The equations were integrated between the inlet and exit of the flow channel to yield three coupled algebraic equations. The algebraic equations were then modified to represent rotating systems. To obtain closure, the velocity ratio, mass quality and void fraction are defined as a function of pressure.

A numerical technique was used to solve the equations. Sample results are presented in a graph of mass quality versus void fraction. The graph demonstrates that a minimum heat input must be exceeded to change from a single-phase flow to saturated two-phase flow boiling. Also, the void fraction was found to increase for increasing heat input, decreasing mass flow rate, increasing inlet mass quality and decreasing pressure difference between the inlet and exit.     

The analysis of two-dimensional two-phase flow in horizontal heated tube bundles using drift flux model

This paper presents the experimental study and numerical simulation of two-dimensional two-phase flow in horizontal heated tube bundles. In the experiments, two advanced measuring systems with a single-fibre optical probe and a tri-fibre-optical-probe were developed to measure respectively the local void fraction and vapor bubble velocities among the heated tube bundles. In accordance with the internal circulation characteristics of two-phase flow in the tube bundles, a mathematical model of two-dimensional two-phase low Reynolds number turbulent flow based on the modified drift flux model and the numerical simulation method to analyze the two-phase flow structures have been developed. The modified drift flux model in which both the acceleration by gravity and the acceleration of the average volumetric flow are taken into account for the calculation of the drift velocities enables its application to the analysis of multi-dimensional two-phase flow. In the analysis the distributions of the vapor-phase velocity, liquid-phase velocity and void fraction were numerically obtained by using the modified drift flux model and conventional drift flux model respectively and compared with the experimental results. The numerical analysis results by using the modified drift flux model agree reasonably well with the experimental investigation. It is confirmed that the modified drift flux model has the capability of correctly simulating the two-dimensional two-phase flow. Received on 3 September 1998     

Propagation of a pressure wave in two-phase flow with very high void fraction

An experimental and theoretical study of a finite amplitude pressure wave propagating through a two-phase media of about 0.9999–0.99999 void fraction is performed. This two-phase media consists of many parallel liquid films in a gas. The films are perpendicular to the wave propagation direction and result in a two-phase fluid of extremely high void fraction. Experiments are done in a vertical shock tube and show that the shock wave is broken down into an initial sharply rising wave and a second gradually rising wave. The velocity of the first wave agrees well with the theoretical prediction assuming an adiabatic thermal equilibrium change, which approaches the gas sonic velocity in the two-phase flow in the low void fraction region. The second wave is caused by the complex reflection and destruction of the waves.     

Void fraction and pressure drop during external upward two-phase cross flow in tube bundles – part II: Predictive methods

The present paper is the Part II of a broad study concerning void fraction and pressure drop for air-water upward external flow across tube bundles. In the Part I, the experimental facility and the data regression procedures were described and the experimental results are presented and discussed. Initially, Part II presents a literature review concerning void fraction and pressure drop predictive methods available in the open literature for two-phase upward flow across tube bundles. Next, the methods from literature are compared among them and with the database presented in paper Part I. Significant discrepancies are observed among the predictive methods, and deviations as high as two orders of magnitude are verified among the predicted values of pressure drop. Then, a new void fraction predictive method is proposed based on the experimental results and on the minimum kinetic energy principle. This method provides satisfactory predictions of the results described in paper Part I and also of independent data from the literature. A new predictive method for frictional pressure drop during two-phase flow based on two-phase multiplier is also proposed. This method predicted 94% of the experimental data obtained in the present study within an error margin of ± 30%, and also provides accurate predictions of independent results for triangular tube bundles gathered in the open literature.     

AN EXPERIMENTAL STUDY ON THE FLUIDELASTIC FORCES ACTING ON A SQUARE TUBE BUNDLE IN TWO-PHASE CROSS-FLOW

A tube in a square tube bundle of P/D=1·42 was oscillated in the lift direction in air–water two-phase cross-flow, and fluidelastic forces acting on the oscillated tube were measured. First, the tube amplitude was fixed to 3 mm (=0·136 D), and added mass, damping, and stiffness coefficients were obtained as a function of two-phase mixture characteristics such as nondimensional gap velocity and void fraction. When reference mixture density and velocity were estimated, the drift–flux model, in which the relative velocity between the gas and liquid phases was estimated, generated better results than the homogeneous model. The added mass coefficient was obtained from quiescent two-phase flow as a function of void fraction. Using the added mass coefficient, the added stiffness coefficient converged to zero with decreasing nondimensional gap velocity. This overcame the contradiction in the added stiffness estimation without added mass, in which the added stiffness coefficient did not converge to zero with decreasing nondimensional gap velocity. Next, the effects of the vibration amplitude on the fluidelastic force coefficients were considered. When the tube amplitude was 3 mm (=0·136 D) or less, the equivalent added stiffness and damping coefficients were almost constant and nonlinearity was small. This showed the validity of the fluidelastic force coefficients obtained based on the data of amplitude of 3 mm. The linearity did not exist when the tube displacement amplitude was 4·5 mm (=0·205 D) or more; a remarkable nonlinearity appeared in the equivalent added damping coefficient. A method to estimate the limit-cycle amplitude of the fluidelastic vibration was proposed when only one tube in the tube bundle was able to vibrate in the lift direction. The amplitude could be obtained from the amplitude at which the equivalent added damping coefficient changed from negative to positive with increase in the tube amplitude.     

Interfacial measurements in stratified types of flow. Part I: New optical measurement technique and dry angle measurements

The development of a new non-intrusive computerized image analysis and optical observation method for accurately detecting the vapor–liquid interface in stratified two-phase flows is presented. This technique is applied to a round horizontal sight glass tube using a monochromatic laser sheet for observing two-phase flow patterns and for measuring cross-sectional dry angles and void fractions in these types of flow. The cross-sectional image observed externally is distorted by refraction and is thus reconstructed by computer. From this image, the shape of the vapor–liquid interface is detected and the dry angle and void fraction are accurately determinable over a wide range of conditions for a glass tube of 13.6 mm diameter. Results for dry angle are reported here while the test facility and void fraction measurements are presented in Part II. R-22 and R-410A are the test fluids. Dry angles were quite close to values predicted for stratified flow and much larger than comparable values for air–water flows.     

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.     

Flow pattern and flow characteristics for counter-current two-phase flow in a vertical round tube with wire-coil inserts

Flow pattern, void fraction and slug rise velocity on counter-current two-phase flow in a vertical round tube with wire-coil inserts are experimentally studied. Flow pattern and slug rise velocity are measured visually with a video camera. The void fraction is measured by the quick-closing valve method. Four kinds of coils with different coil pitches and coil diameters are used as inserts. The presence of wire-coil inserts induces disturbance into gas and liquid flows so that the shape and motion of gas slug or bubbles in a wire-coil inserted tube are quite different from those observed in a smooth tube without insert. The bubbly flow occurs in the low gas superficial velocity region in the wire-coil inserted tube, while the slug or churn/annular flow only appears in the smooth tube without insert over the all test range. The measured slug rise velocity in the wire-coil inserted tube is higher than that in the smooth tube. With modified mean flow velocity calculated with core area, the slug rise velocity in wire-coil tube inserted is in good agreement with Nicklin's correlation. The void fraction in a wire-coil inserted tube is lower than that in a smooth tube in the range of high gas superficial velocities. By introducing a simple assumption on considering the effective flowing area, the measured void fractions in a wire-coil inserted tube are in relatively good agreement with the predicted result based on the drift flux model proposed by others with the correlation for slug rise velocity given by others when the coil pitch is dense.     

Two-phase damping and interface surface area in tubes with vertical internal flow

Two-phase flow is common in the nuclear industry. It is a potential source of vibration in piping systems. In this paper, two-phase damping in the bubbly flow regime is related to the interface surface area and, therefore, to flow configuration. Experiments were performed with a vertical tube clamped at both ends. First, gas bubbles of controlled geometry were simulated with glass spheres let to settle in stagnant water. Second, air was injected in stagnant alcohol to generate a uniform and measurable bubble flow. In both cases, the two-phase damping ratio is correlated to the number of bubbles (or spheres). Two-phase damping is directly related to the interface surface area, based on a spherical bubble model. Further experiments were carried out on tubes with internal two-phase air–water flows. A strong dependence of two-phase damping on flow parameters in the bubbly flow regime is observed. A series of photographs attests to the fact that two-phase damping in bubbly flow increases for a larger number of bubbles, and for smaller bubbles. It is highest immediately prior to the transition from bubbly flow to slug or churn flow regimes. Beyond the transition, damping decreases. It is also shown that two-phase damping increases with the tube diameter.     

Propagation of pressure waves in two-phase flow

We deal with a pressure wave of finite amplitude propagating in a gas and liquid medium or in the fluid in an elastic tube. We study the effects of pipe elasticity on the propagation velocity of the pressure wave. Pressure waves of finite amplitude progressing in the two-phase flow are treated considering the void fraction change due to pressure rise. The propagation velocity of the two-phase shock wave is also investigated, and the behavior of the reflection of the pressure wave at the rigid wall is analyzed and compared to that in a pure gas or liquid. The results are compared to experimental data of a pressure wave propagating in the two-phase flow in a vertical shock tube.     

DYNAMICS OF AN IN-LINE TUBE ARRAY SUBJECTED TO STEAM–WATER CROSS-FLOW. PART III: FLUIDELASTIC INSTABILITY TESTS AND COMPARISON WITH THEORY

This paper presents the results of comprehensive flow-induced vibration tests conducted on an in-line array in steam–water two-phase flow. The responses of three essentially isolated flexible cylinders at different depths within the array were simultaneously measured. The main test parameters were, ambient pressure (and saturation temperature), in the range 0·5–5·8 MPa, void fraction, 0·70–0·96, and phase flow velocity. Tests reported here were conducted simultaneously with the damping tests reported in Part I of this study. At the highest pressures (3·0 and 5·8 MPa), strong instabilities, in homogeneous flow akin to single-phase flow occurred. The test tube located in the central region of the array was the most susceptible to instability. This was attributed partly to reduced two-phase damping deep in the array, while differences in local fluid forces at different locations in the array are not ruled out. The flow at 0·5 MPa was a nonhomogeneous intermittent slug-type flow. Strong turbulence excitation obscured clear fluidelastic instability; intermittent instability was, however, ascertained. Stability boundary calculations were done using unsteady fluid forces presented in Part II of this series of papers. Results for the case of P=5·8 MPa show good agreement with the measured instability boundary.     

水平管气液两相间歇流含气率研究

对大直径水平管内气液二相流动进行了实验，实验结果表明大直径水平管与小直径水平管具有不同的气液二相流动特征。分析和提出了适用于大直径水平管内间歇流液弹含气率模型，其计算值与实验值非常接近。