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.     

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.     

Measured Reynolds stress distributions and energy budgets of a fully pulsed round air jet

 Most particle-tracking velocimetry (PTV) algorithms are not suitable for calculating the velocity vectors of a fluid flow subjected to strong deformation, because these algorithms deal only with flows due to translation. Accordingly, it is necessary to develop a novel algorithm applicable to flows subjected to strong deformations such as rotation, shear, expansion and compression. This paper proposes a novel particle tracking algorithm using the velocity gradient tensor (VGT) which can deal with strong deformations and demonstrates that this algorithm is applicable to some basic fluid motions (rigidly rotating flow, Couette flow, and expansion flow). Furthermore, the performance of this algorithm is compared with the binary image cross-correlation method (BICC), the four-consecutive-time-step particle tracking method (4-PTV), and the spring model particle tracking algorithm (SPG) using simulations and experimental data. As a result, it is shown that this novel algorithm is useful and applicable for the highly accurate measurement and analysis of fluid flows subjected to strong deformations. Received: 9 February 1999/Accepted: 22 November 1999     

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.     

Investigation of turbulence modification in a non-reactive two-phase flow

In a two-phase flow the influence of a dispersed phase on the turbulence properties of a continuous phase, known as turbulence modification, is investigated. An experimental approach is discussed that is suitable for studying the decay of grid-generated turbulence in a vertically orientated wind tunnel with a cross-section large enough to avoid influences from walls. Phase Doppler anemometry is used to characterize both single and two-phase flow by measuring mean axial and radial velocity components, velocity fluctuations, turbulent kinetic energy, and integral time scales. By direct comparison of results from single- and two-phase flows, the feedback of the dispersed phase on the continuous phase can be isolated. The data is used to deduce a source term for particle-induced turbulence production appropriate for a numerical simulation of the flow, based on the Reynolds-averaged Navier-Stokes equations. Although of special importance for a detailed understanding of turbulent two-phase combustion, additional complexity introduced by evaporation and chemical reactions is avoided by using glass beads as dispersed phase.     

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.     

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.     

Application of photobleaching molecular tagging velocimetry to turbulent bubbly flow in a square duct

To evaluate turbulence energy budget in bubbly flows, an image processing method in a photobleaching molecular tagging velocimetry is improved for accurate evaluation of velocity gradients. Turbulence properties in single-phase and two-phase dilute-bubbly flows in a square duct are measured using the improved method. As a result, the following conclusions are obtained: (1) The axial velocity and axial turbulent intensity measured by the present method agree well with those measured by laser Doppler velocimetry not only for the single-phase flow but also for the dilute-bubbly flow. (2) The present method can measure velocity components and velocity gradients in the vicinity of the wall, and therefore the present method is of great use in understanding the mechanism of turbulence generation and dissipation near the wall. (3) The present method can provide detailed information on turbulence structure such as turbulence kinetic energy budget. (4) Bubbles tend to increase not only the turbulence production but also the turbulence dissipation.     

Theoretical and numerical studies on the flow multiplicity phenomenon for gas–solids two-phase flows in CFB risers

The dependence of the fully-developed flow profiles on the inlet flow conditions for gas–solids two-phase flows, i.e. the flow multiplicity phenomenon, in circulating fluidized bed (CFB) risers was proposed and discussed in this article. The flow multiplicity phenomenon for gas–solids two-phase flows was first proved mathematically based on the conservation equations of mass and momentum. Then the CFD model using Eulerian–Eulerian approach with k–ε turbulence model for each phase was further adopted to analyze the details of this flow multiplicity phenomenon. It is theoretically and numerically revealed that for gas–solids two-phase flows, the flow profiles in the fully-developed region are always dominated by the flow profiles at the inlet. The solids concentration profile is closely coupled with the velocity profile, and the inlet solids concentration and velocity profiles can largely influence the fully-developed concentration and velocity profiles.     

Benchmarking of the one-dimensional one-group interfacial area transport equation for reduced-gravity bubbly flows

In order to properly design and safely operate two-phase flow systems, especially those deployed on future space missions, it is necessary to have accurate predictive capabilities. The application of a novel predictive method, the interfacial area transport equation (IATE), to dynamically predict the change of interfacial area concentration for reduced-gravity two-phase flows is described in this paper. Fluid particle interaction mechanisms such as coalescence and breakup that are present in reduced-gravity two-phase flows have been studied experimentally as reported in a previous paper by the current authors [Vasavada et al., 2007]. These mechanisms represent the source and sink terms in the IATE and their mechanistic models are benchmarked using experimental data obtained in a 25 mm inner diameter ground-based test section wherein reduced-gravity conditions were simulated. The comparison of the predictions from the model against experimental data shows good agreement. It has been found that, in contrast to the hypothesis extended in the literature, the wake entrainment based coalescence mechanism is present in reduced-gravity two-phase flows and in some cases is more important than coalescence due to random collision. Physics based arguments are extended to support this conclusion.     

Gas-Dynamic Processes in Two-Phase Flows in MHD Generators

A plane problem of a two-phase monodisperse flow of combustion products of plasma-forming composite solid propellants in the duct of a Faraday's MHD generator with continuous electrodes, including an accelerating nozzle, MHD channel, and diffuser, is considered. An algorithm based on the pseudo-transient method is developed to solve the system of equations describing the two-phase flow. Gas-dynamic processes in the channels of the Pamir-1 setup are numerically studied. It is shown that shock-free deceleration of a supersonic flow to velocities close to the equilibrium velocity of sound in a two-phase mixture and significantly lower than the velocity of sound in the gas is possible in two-phase flows.     

Measurement and simulation of the two-phase velocity correlation in sudden-expansion gas-particle flows

In this paper the present authors measured the gas-particle two-phase velocity correlation in sudden expansion gas-particle flows with a phase Doppler particle anemometer(PDPA) and simulated the system behavior by using both a Reynolds-averaged Navier-Stokes(RANS) model and a large-eddy simulation(LES).The results of the measurements yield the axial and radial time-averaged velocities as well as the fluctuation velocities of gas and three particle-size groups(30 μ m,50 μ m,and 95 μ m) and the gas-particle velocity correlation for 30 μ m and 50 μ m particles.From the measurements,theoretical analysis,and simulation,it is found that the two-phase velocity correlation of sudden-expansion flows,like that of jet flows,is less than the gas and particle Reynolds stresses.What distinguishes the two-phase velocity correlations of sudden-expansion flow from those of jet and channel flows is the absence of a clear relationship between the two-phase velocity correlation and particle size in sudden-expansion flows.The measurements,theoretical analysis,and numerical simulation all lead to the above-stated conclusions.Quantitatively,the results of the LES are better than those of the RANS model.     

Turbulent cascade modeling of single and bubbly two-phase turbulent flows

The analysis of turbulent two-phase flows requires closure models in order to perform reliable computational multiphase fluid dynamics (CMFD) analyses. A spectral turbulence cascade-transport model, which tracks the evolution of the turbulent kinetic energy from large to small liquid eddies, has been developed for the analysis of the homogeneous decay of isotropic single and bubbly two-phase turbulence. This model has been validated for the decay of homogeneous, isotropic single and two-phase bubbly flow turbulence for data having a 5 mm mean bubble diameter. The Reynolds number of the data based on bubble diameter and relative velocity is approximately 1400.     

Interfacial area concentration in gas–liquid bubbly to churn flow regimes in large diameter pipes

This study performed a survey on existing correlations for interfacial area concentration (IAC) prediction and collected an IAC experimental database of two-phase flows taken under various flow conditions in large diameter pipes. Although some of these existing correlations were developed by partly using the IAC databases taken in the low-void-fraction two-phase flows in large diameter pipes, no correlation can satisfactorily predict the IAC in the two-phase flows changing from bubbly, cap bubbly to churn flow in the collected database of large diameter pipes. So this study presented a systematic way to predict the IAC for the bubbly-to-churn flows in large diameter pipes by categorizing bubbles into two groups (group 1: spherical or distorted bubble, group 2: cap bubble). A correlation was developed to predict the group 1 void fraction by using the void fraction for all bubble. The group 1 bubble IAC and bubble diameter were modeled by using the key parameters such as group 1 void fraction and bubble Reynolds number based on the analysis of Hibiki and Ishii (2001, 2002) using one-dimensional bubble number density and interfacial area transport equations. The correlations of IAC and bubble diameter for group 2 cap bubbles were developed by taking into account the characteristics of the representative bubbles among the group 2 bubbles and the comparison between a newly-derived drift velocity correlation for large diameter pipes and the existing drift velocity correlation of Kataoka and Ishii (1987) for large diameter pipes. The predictions from the newly-developed two-group IAC correlation were compared with the collected experimental data in gas–liquid bubbly to churn flow regimes in large diameter pipes and their mean absolute relative deviations were obtained to be 28.1%, 54.4% and 29.6% for group 1, group 2 and all bubbles respectively.     

Velocity measurements of flow through a step stenosis using Magnetic Resonance Imaging

 Magnetic resonance imaging (MRI) is a versatile noninvasive tool for achieving full-field quantitative visualization of complex fluid flows. The MRI signal results from the interaction of radio-frequency (RF) pulses with nuclear spins exposed to a strong static magnetic field. The two main techniques of MRI velocimetry are time-of-flight and phase contrast techniques. Time-of- flight techniques involve tagging and tracking a material volume of fluid, whereas phase contrast techniques use magnetic field gradients to encode velocity information into the phase of the MRI signal. In this study, both techniques are used to probe the pressure-driven steady flow of water in a pipe with a step stenosis. The velocity measurements were then compared with computational results obtained using the FIDAP software package. The experiments show that the phase contrast method gives more accurate results, with 90% of the measurements within 10% of the local computational fluid dynamics (CFD) velocity predictions at Re = 100 and 94% of the measurements within 10% of the local CFD predictions at Re = 258. Although the time-of-flight experiments were not as accurate, they provide a good qualitative image of the flow field. Sources of the discrepancies between the MRI data and the CFD results are also discussed, including acceleration and spin flow-through artifacts. Received: 7 April 1999/Accepted: 20 December 1999     

Microgravity experiments of two-phase flow patterns aboard MIR space station

A first experimental study on two-phase flow patterns at a long-term, steady microgravity condition was conducted on board the Russian Space Station “MIR” in August 1999. Carbogal and air are used as the liquid and the gas phase, respectively. Bubble, slug, slug-annular transitional, and annular flows are observed. A new region of annular flow with lower liquid superficial velocity is discovered, and the region of the slug-annular transitional flow is wider than that observed by experiments on board the parabolic aircraft. The main patterns are bubble, slug-annular transitional and annular flows based on the experiments on board MIR space station. Some influences on the two-phase flow patterns in the present experiments are discussed. The project supported by the National Natural Science Foundation of China (19789201), the Ministry of Science and Technology of China (95-Yu-34), and the Post-doctoral Science Foundation of China     

Flow regime identification in a two-phase flow using wavelet transform

This study addresses the problem of the automatic flow regime identification in two-phase flows in pipes. A novel wavelet transform-based approach is proposed and validated using time series of differential pressure fluctuations. The experimental data on the differential pressure measured in a vertically installed Venturi meter for air-water flow were analyzed and found to be appropriate for flow regime identification. The wavelet spectrum of the measured signal is shown to characterize the flow patterns completely, and the vector of the wavelet variances is proposed as the characteristic vector for use in an on-line flow regime identification system.     

Improved subgrid scale model for dense turbulent solid-liquid two-phase flows

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules. Assuming that the solid-phase velocity distributions obey the Maxwell equations, the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow, the solid-particle‘s governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations. Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale (SGS) model, in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor, is proposed to model the two-phase governing equations by applying dimension analyses. Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls, the velocity and pressure fields, and the volumetric concentration are calculated. The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.     

Separated two-phase flow in a nozzle

This work was performed to extend and further test the method of handling separated two-phase flow by studying each phase separately and, particularly, by placing emphasis on the study of the gas phase with interface transport expressions showing the influence of the liquid phase on it. A one-dimensional flow model for accelerating flows was used in conjunction with experimental data to obtain the pressure distribution and velocity distribution in a converging nozzle for several values of flow quality and nozzle inlet stagnation pressure. The results tend to support the use of the model (which includes the assumption that the gas is in critical flow when the two-phase mixture is in critical flow) and give some insight regarding the nature of the liquid distribution near the nozzle throat.     

Influence of the wall boundary layer thickness on a gas jet injected into a liquid crossflow

The effects of the wall boundary layer thickness on the development of an axisymmetric gas jet injected into a confined vertical water flow were investigated. The variations in the wall boundary layer were made by using suction at the wall through rectangular profiled slots. The water velocity around the two-phase jet was studied for several boundary layer thickness values by laser Doppler velocimetry. The gas jet outline was extracted by image processing applied on visualisations for a wide variety of water, gas and suction conditions. These comparisons showed that the boundary layer has no influence on jet development. The data showed that the interactions between flows near the injection do not develop downstream, which accounts for the absence of the two classical contrarotating vortices in the medium field of the jet. The influence of the gas pocket, and then of buoyancy, is predominant over other phenomena. Received: 15 December 1999 / Accepted: 29 August 2000