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.     

Effect of bubble size on void fraction fluctuations in dispersed bubble flows

It is known that bubble size affects seriously the average void fraction in bubbly flows where buoyant velocities vary considerably with bubble size. On the contrary, there is no systematic literature report about bubble size effects on the intensity and frequency of void fraction fluctuations around the average void fraction. This work aims to provide such information. An electrical impedance technique is employed along with non-intrusive ring electrodes to register void fraction fluctuations down to 10⁻⁵. Bubble size fluctuations are estimated from high resolution optical images. Experiments are conducted in co-current upward dispersed bubble flow inside a 21 mm tube with average bubble size between ∼50 and ∼700 μm. Water and blood simulant are used as test liquids with velocity from ∼3 to ∼30 cm s⁻¹. The above resemble conditions of Decompression Sickness (DCS) in the bloodstream of human vena cava. It is found that the intensity and frequency of void fraction fluctuations vary appreciably with bubble size at constant gas and liquid flow rates. Moreover, these variations are not random but scale with bubble size. As a first step to quantify this effect, an empirical expression is derived that relates average bubble size to the ratio standard deviation/average value of void fraction.     

Bubble–wall interaction and two-phase flow parameters on a full-scale boat boundary layer

Full scale bubbly flow experiments were performed on a 6 m flat bottom survey boat, measuring the void fraction, bubble velocity and size distributions as the bubbles naturally entrained at the bow of the boat interact with the boat’s boundary layer. Double-tip sapphire optical probes capable of measuring bubbles down to 50 μm in diameter were specifically designed and built for this experiment. The probes were positioned under the hull at the bow near the bubble entrainment region and at the stern at the exit of the bottom flat plate. Motorized positioners were used to vary the probe distance to the wall from 0 to 50 mm. The experiments were performed in fresh water (Coralville Lake, IA) and salt water (Panama City Beach, FL), at varying velocities with most data analysis performed at 10, 14 and 18 knots. The results indicate that the bubbles interact significantly with the boundary layer. At low velocity in fresh water, bubble accumulation under the hull and coalescence are evident by the presence of large bubbles at the stern. At high speeds bubble breakup dominates and very small bubbles are produced near the wall. It is also observed that salt water inhibits coalescence, even at low boat speeds. The void fraction increases with speed beyond 10 knots and peaks near the wall. Bubble velocities show slip with the wall at all speeds and exhibit large RMS fluctuations, increasing near the wall.     

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.     

Planar fluorescence for round bubble imaging and its application for the study of an axisymmetric two-phase jet

A correlation-based processing algorithm for bubble identification by a planar fluorescence for bubble imaging (PFBI) technique is presented in this paper. The algorithm includes procedures to identify bubble positions and sizes, as well as to track bubbles and correct bubble displacement vectors. Moreover, several schemes for calculation time optimisation were realised to achieve a reliable calculation time. The developed algorithm identifies and tracks overlapping bubble images or images with non-uniform intensity distributions. The employed correlation and iterative passing approach provides sub-pixel accuracy of bubble displacement estimation. In addition, the presented algorithm for bubble ring detection can be easily applied to shadow photography images of bubbles, after the application of a derivative filter. The PFBI technique, combined with the particle image velocimetry and particle tracking velocimetry algorithms, was applied for the experimental study of bubbly free jet two-phase flows at Re = 12,000. Four cases of volumetric gas content in the jet core were studied: 0, 1.2, 2.4 and 4.2%, with the same mean bubble diameter—0.85 mm. The developed technique measures two-dimensional distributions of instantaneous void fractions, as well as both gaseous and liquid-phase velocities. Consequently, the mean void fraction and velocity fields and a set of second-order statistical moments were obtained, including correlations of void fraction and velocity pulsations. It was shown that the increase in volumetric gas content leads to the suppression of liquid-phase velocity fluctuations in the jet mixing layer.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Vertical stability of bubble chain: Multiscale approach

Linear stability is investigated of a uniform chain of equal spherical gas bubbles rising vertically in unbounded stagnant liquid at Reynolds number Re = 50–200 and bubble spacing s > 2.6 bubble radii. The equilibrium bubble positions are questioned for their stability with respect to small displacements in the vertical direction, parallel to the chain motion. The transverse displacements are not considered, and the chain is assumed to be laterally stable. The bubbles are subjected to three kinds of forces: buoyant, viscous, inviscid. The viscous and inviscid forces have both pairwise (local) and distant (nonlocal) components. The pairwise forces are expressed by the leading-order formulas known from the literature. The distant forces are expressed as a linear superposition of the pairwise forces taken over several farther neighbours. The stability problem is addressed on three different length scales corresponding to: discrete chain (microscale), continuous chain (mesoscale), bubbly chain flow (macroscale). The relevant governing equations are derived for each scale. The microscale equations are a set of ODE’s, the Newton force laws for the individual discrete bubbles. The mesoscale equation is a PDE for bubbles continuously distributed along a line, obtained by taking the continuum limit of the microscale equations. The macroscale equations are two PDEs, the mass and momentum conservation equations, for an ensemble of noninteracting mesoscale chains rising in parallel. This transparent two-step process (micro → meso → macro) is an alternative to the usual one-step averaging, in obtaining the macroscale equations from microscale information. Here, the scale-up methodology is demonstrated for 1D motion of bubbles, but it can be used for behaviour of 2D and 3D lattices of bubbles, drops, and solids.It is found that the uniform equilibrium spacing results from a balance between the attractive and repulsive forces. On all three length scales, the equilibrium is stabilized by the viscous drag force, and destabilized by the viscous shielding force (shielding instability). The inviscid forces are stability neutral and generate conservative oscillations and concentration waves. The stability region in the parameter plane s − Re is determined for each length scale. The stable region is relatively small on the microscale, larger on the mesoscale, and shrinks to zero on the macroscale where the bubbly chain flow is inherently unstable.The shielding instability is expected to occur typically in intermediate Re flows where the vertical bubble interactions dominate over the horizontal interactions. This new kind of instability is studied here in a great detail, likely for the first time. Its relation to the elasticity properties of bubbly suspension on different length scales is discussed too. The shielding force takes the form of a negative bulk modulus of elasticity of the bubbly mixture.     

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%.     

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.     

Enhancement of channel wall vibration due to acoustic excitation of an internal bubbly flow

The effect of an internal turbulent bubbly flow on vibrations of a channel wall is investigated experimentally and theoretically. Our objective is to determine the spectrum and attenuation rate of sound propagating through a bubbly liquid flow in a channel, and connect these features with the vibrations of the channel wall and associated pressure fluctuations. Vibrations of an isolated channel wall and associated wall pressure fluctuations are measured using several accelerometers and pressure transducers at various gas void fractions and characteristic bubble diameters. A waveguide-theory-based model, consisting of a solution to the three-dimensional Helmholtz equation in an infinitely long channel with the effective physical properties of a bubbly liquid is developed to predict the spectral frequencies of the wall vibrations and pressure fluctuations, the corresponding attenuation coefficients and propagation phase speeds. Results show that the presence of bubbles substantially enhances the power spectral density of the channel wall vibrations and wall pressure fluctuations in the 250–1200 Hz range by up to 27 and 26 dB, respectively, and increases their overall rms values by up to 14.1 and 12.7 times, respectively. In the same frequency range, both vibrations and spectral frequencies increase substantially with increasing void fraction and slightly with increasing bubble diameter. Several weaker spectral peaks above that range are also observed. Trends of the frequency and attenuation coefficients of spectral peaks, as well as the phase velocities are well predicted by the model. This agreement confirms that the origin of enhanced vibrations and pressure fluctuations is the excitation of streamwise propagating pressure waves, which are created by the initial acoustic energy generated during bubble formation.     

Flow pattern characterization for R-245fa in minichannels: Optical measurement technique and experimental results

An optical measurement method using image processing for two-phase flow pattern characterization in minichannel is developed. The bubble frequency, the percentage of small bubbles as well as their velocity are measured. A high-speed high-definition video camera is used to measure these parameters and to identify the flow regimes and their transitions. The tests are performed in a 3.0 mm glass channel using saturated R-245fa at 60 °C (4.6 bar). The mass velocity is ranging from 100 to 1500 kg/m² s, the heat flux is varying from 10 to 90 kW/m² and the inlet vapor quality from 0 to 1. Four flow patterns (bubbly flow, bubbly–slug flow, slug flow and annular flow) are recognized. The comparison between the present experimental intermittent/annular transition lines and five transition lines from macroscale and microscale flow pattern maps available in the literature is presented. Finally, the influence of the flow pattern on the heat transfer coefficient is highlighted.     

Evaluation of flow patterns and elongated bubble characteristics during the flow boiling of halocarbon refrigerants in a micro-scale channel

In the present study, quasi-diabatic two-phase flow pattern visualizations and measurements of elongated bubble velocity, frequency and length were performed. The tests were run for R134a and R245fa evaporating in a stainless steel tube with diameter of 2.32 mm, mass velocities ranging from 50 to 600 kg/m² s and saturation temperatures of 22 °C, 31 °C and 41 °C. The tube was heated by applying a direct DC current to its surface. Images from a high-speed video-camera (8000 frames/s) obtained through a transparent tube just downstream the heated sections were used to identify the following flow patterns: bubbly, elongated bubbles, churn and annular flows. The visualized flow patterns were compared against the predictions provided by Barnea et al. (1983) [1], Felcar et al. (2007) [10], Revellin and Thome (2007) [3] and Ong and Thome (2009) [11]. From this comparison, it was found that the methods proposed by Felcar et al. (2007) [10] and Ong and Thome (2009) [1] predicted relatively well the present database. Additionally, elongated bubble velocities, frequencies and lengths were determined based on the analysis of high-speed videos. Results suggested that the elongated bubble velocity depends on mass velocity, vapor quality and saturation temperature. The bubble velocity increases with increasing mass velocity and vapor quality and decreases with increasing saturation temperature. Additionally, bubble velocity was correlated as linear functions of the two-phase superficial velocity.     

An advanced LIF-PLV system for analysing the hydrodynamics in a laboratory bubble column at higher void fractions

Bubble columns are widely used in the chemical industry and biotechnology. Flow and turbulence in such an apparatus are induced by the bubble rise, and the bubble behaviour is strongly affected by swarm effects (i.e. the interaction between bubbles). For analysing the bubble swarm behaviour and simultaneously evaluating the flow structure and bubble-induced turbulence, a bubble column of 140 mm diameter and a height of 650 mm or 1,400 mm (initial water level) were considered. The bubble column was aerated with relatively fine bubbles having a mean size between about 0.5 and 4.0 mm. The gas hold-up was varied in the range between 0.5 and 19%. A two-phase pulsed-light velocimetry (PLV) system was developed to evaluate instantaneous flow fields of both rising bubbles and the continuous phase. The measurement of the liquid velocities in the bubble swarm was achieved by adding fluorescing seed particles. Images of bubbles and fluorescing tracer particles were acquired by two CCD cameras. Hence, the images from tracers and bubbles were easily separated by optical interference filters with a bandwidth corresponding to the emitting wavelength of the fluorescing tracer particles and the wavelength of the applied Nd-YAG pulsed laser, respectively. To improve the phase separation of the system, the CCD cameras were additionally placed in a non-perpendicular arrangement with respect to the light sheet. The acquired images were evaluated with the minimum-quadratic-difference algorithm. The potential of this technique for the analysis of bubbly flows with higher void fraction was explored. In order to obtain averaged velocity maps of bubble and fluid within the entire column, about 1,000 image pairs were recorded and evaluated for each phase. In addition, turbulence intensities of the fluid were deduced from the measurements. The turbulence properties were used to characterise bubble-induced turbulence for various bubble mean diameters and gas hold-ups. Moreover, the determination of the average bubble slip velocity within the bubble swarm was possible.     

Three-dimensional measurement of void fraction in a bubble plume using statistic stereoscopic image processing

A method for measuring the instantaneous 3-D bubble distribution in bubbly flows is proposed by making use of stereoscopic image processing. The method applies to the case of higher rather than dilute void fractions, because it is based on the following two statistic-based procedures. The first process is to obtain the 2-D distribution of local void fraction from visualized images, where a correlation between the bubble shadow fraction and the projection void fraction is derived by introducing the new idea of overlapping probability of the bubble-images under various conditions. The second process is to reconstruct the 3-D distribution of the local void fraction. In this step, a logical method is proposed in order to obtain the bubble-existence probability function which simultaneously satisfies the two distributions of the two projection void fractions. The present method is applied for instantaneous measurement of the 3-D structure of a bubble plume. The measured results show its particular structure, such as the hollow distribution of the bubbles, and the local spiral structure, which cannot be captured by point-measurement experiments. Received: 5 January 1999/Accepted: 3 February 2000     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Dynamics of a Taylor bubble in steady and pulsatile co-current flow of Newtonian and shear-thinning liquids in a vertical tube

A computational analysis is carried out to ascertain the effects of steady and pulsatile co-current flow, on the dynamics of an air bubble rising in a vertical tube containing water or a solution of Carboxymethylcellulose (CMC) in water. The mass fraction (m_f) of CMC in the solution is varied in the range 0.1% ⩽ m_f ⩽ 1% to accommodate zero-shear dynamic viscosities in the range 0.009–2.99 Pa-s. It was found that the transient and time-averaged velocities of Taylor bubbles are independent of the bubble size under both steady as well as pulsatile co-current flows. The lengths of the Taylor bubbles under the Newtonian conditions are found to be consistently greater than the corresponding shear-thinning non-Newtonian conditions for any given zero-shear dynamic viscosity of the liquid. In contrast to observations in stagnant liquid columns, an increase in the dynamic viscosity of the liquid (under Newtonian conditions) results in a concomitant increase in the bubble velocity, for any given co-current liquid velocity. In shear-thinning liquids, the change in the bubble velocity with an increase in m_f is found to be relatively greater at higher co-current liquid velocities. During pulsatile shear-thinning flows, distinct ripples are observed to occur on the bubble surface at higher values of m_f, the locations of which remain stationary with reference to the tube for any given pulsatile flow frequency, while the bubble propagated upwards. In such a pulsatile shear-thinning flow, a localised increase in dynamic viscosity is accompanied near each ripple, which results in a localised re-circulation region inside the bubble, unlike a single re-circulation region that occurs in Newtonian liquids, or shear-thinning liquids with low values of m_f. It is also seen that as compared to frequency, the amplitude of pulsatile flow has a greater influence on the oscillating characteristics of the rising Taylor bubble. The amplitude of oscillation in the bubble velocity increases with an increase in the CMC mass fraction, for any given value of pulsatile flow amplitude.     

Air–water flow characteristics in high-velocity free-surface flows with 50% void fraction

High-velocity free-surface flows are complex two-phase flows and limited information is available about the interactions between air and water for void fractions of about 50%. Herein a detailed experimental study was conducted in the intermediate flow region (C ∼ 50%) on a stepped spillway and the microscopic air–water flow characteristics were investigated. The results showed differences in water and droplet chord times with comparatively larger number of air chord times (0–2 ms), and larger number of water chord times (2–6 ms). A monotonic decrease of particle chord modes was observed with increasing bubble count rates. Several characteristic time scales were identified based upon inter-particle arrival time analyses of characteristic chord time classes as well as spectral analyses of the instantaneous void fraction signal. Chord times of 3–5 ms appeared to be characteristic time scales of the intermediate flow region having similar time scales compared to the local correlation and integral turbulent time scales and to time scales associated with bubble break-up and turbulent velocity fluctuations. A further characteristic time scale of 100 ms was identified in a frequency analysis of instantaneous void fraction. This time scale was of the same order of magnitude as free-surface auto-correlation time scales suggesting that the air–water flow structure was affected by the free-surface fluctuations.     

Void fraction and pressure drop during external upward two-phase crossflow in tube bundles – part I: Experimental investigation

The present paper is the part I of a broad study concerning void fraction and pressure drop for air-water upward external flow across tube bundles. Experimental results were obtained for liquid and gas superficial velocities ranging from 0.02 to 1.50 m/s and 0.20 to 10.00 m/s, respectively. Void fraction measurements were performed for bubbly flow using a capacitive probe. The test section consisted of a triangular tube bundle counting with 19 mm OD tube and transverse pitch of 24 mm. Initially, the paper describes the test facility and the data regression and experimental procedures. Then, the pressure drop and void fraction measurements are validated based on tests for single-phase flow and quiescent liquid conditions, respectively. Finally, the experimental data are presented and analyzed. In the second part of this study (Part II), a literature review on predictive methods for void fraction and pressure drop is presented. Additionally, these methods are compared with the database presented in Part I and new predictive methods for void fraction and frictional pressure drop are proposed.     

Experimental and computational studies of shock wave-to-bubbly water momentum transfer

Momentum transfer from shock waves (SWs) of various intensity (from 0.05 MPa to 0.5 MPa in amplitude) to water containing air bubbles 2.5 to 4 mm of mean diameter is studied both experimentally and by means of numerical simulation. Experiments are performed in a vertical shock tube of a 50 × 100 mm² rectangular cross section consisting of a 495-mm long high-pressure section (HPS), 495-mm long low-pressure section (LPS), and 990 mm long test section (TS) equipped with an air bubbler and filled with water. Experiments have shown that as the initial gas volume fraction in water increases from 0 to 0.3 the momentum imparted in bubbly water by SWs increases monotonically, gradually levelling off at an air volume fraction of about 0.30. The experimental data are confirmed by two-dimensional (2D) simulation of SW propagation in bubbly water in terms of the SW velocity versus the air content, pressure profiles, as well as liquid and gas velocity behind the shock front.