Development of a robust image processing technique for bubbly flow measurement in a narrow rectangular channel

This paper presents a robust image processing technique for bubbly flow measurement over a wide range of void fractions. The proposed algorithm combines geometrical, optical and topological information recorded with high speed cameras to separate and reconstruct the overlapping bubbles. The common difficulties such as overlapping, irregular bubble shape, surface deformation and large clustering in digital image processing are solved by combining different information based on a preset decision table and flow chart. Test with synthetic bubble images is performed to evaluate the reliability of the algorithm and quantify the uncertainty of the data. The result shows that the proposed algorithm can accurately measure bubbly flows with void fraction up to 18% for large bubbles. Four runs of bubbly flow images in a 30 mm  ×  10 mm rectangular channel are then recorded by three high speed cameras. The area-averaged void fraction of these test runs range from 2.4% to 9.1%. The axial and lateral distributions of bubble number density are obtained by the present algorithm for studying the characteristics of these flows.     

Vector post-processing algorithm for phase discrimination of two-phase PIV

A simple phase separation method using vector post-processing techniques is evaluated to measure velocity fields in a bubble plume. To provide for validation, fluorescent seeding is used, and two sets of synoptic images are obtained: mixed-phase images containing bubbles and fluorescent particles, and fluid-phase images containing only fluorescent particles. A third dataset is derived by applying a digital mask to remove bubbles from the mixed-phase images. All datasets are processed using cross-correlation particle image velocimetry (PIV). The resulting vector maps for the raw, mixed-phase data contain both bubble and continuous-phase velocity vectors. To separate the phases, a vector post-processing algorithm applies a maximum velocity threshold for the continuous-phase velocities coupled with the vector median filter to identify remaining bubble-velocity vectors and remove them from the mixed-phase velocity field. To validate the phase separation algorithm, the post-processed fluid-phase vectors are compared to PIV results obtained from both the optically separated and digitally masked data. The comparison among these methods shows that the post-processed mixed-phase data have small errors in regions near some bubbles, but for dilute environmental flows (low void fraction and slip velocity approximately equal to the entrained fluid velocity), the algorithm predicts well both instantaneous and time average statistical quantities. The method is reliable for flows having 10% or less of the field of view occupied by bubbles. The resulting instantaneous data provide information on plume wandering and eddy-size distributions within the bubble plume. By comparison among the datasets, it is shown that the patchiness of the vector-post processed and image masked data limit the diameter of identifiable eddy structures to the average distance between bubbles in the image, and that both datasets give identical probability density functions of eddy size. The optically filtered data have better data coverage and predict a greater probability of larger eddies as compared to the other two datasets.     

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

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

Computational study of steady streaming from oscillating microbubbles with uniform and wavy wall motions

Steady streaming flow fields of a 5 μm bubble oscillating with uniform radial wall motion and a 500 μm bubble oscillating with wavy wall motion were simulated using a computational fluid dynamics method that incorporated fluid–structure interactions. The steady streaming flow fields for both bubbles were calculated, and they exhibited upward jet flow with two symmetrical counter-rotating vortices. The maximum streaming velocity ranged from a few to tens of millimeters per second. The simulated flow fields were compared with the theory and experimental measurements using particle image velocimetry. The simulation results agreed well with the theoretical and experimental data. Therefore, the proposed computational method would provide a useful tool to predict steady streaming flow fields of oscillating bubbles.     

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.     

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.     

Effect of bubbles on turbulent kinetic energy transport in downward flow measured by time-resolved PTV

A time-resolved particle tracking velocimetry (PTV) system and a shape projection imaging system were used to investigate the turbulence modifications by bubbles in a downward bubbly flow. Two bubble sizes and three mean void fractions were tested at a Reynolds number of about 20,000. The strong modifications in the mean velocity, turbulent kinetic energy (TKE) budget, and velocity spectra are observed in the central region of the pipe that has a high local void fraction. In particular, kinetic energy production decreased, whereas the TKE dissipation rate increased. This suggests that the transfer of energy due to bubbles has a very large effect on the TKE budget. Moreover, velocity spectra reveal that the presence of bubbles modifies the length scales of turbulent eddies, which contain, transfer, and dissipate energy.     

Comprehensive experimental investigation of the hydrodynamics of large-scale, 3D,oscillating bubble plumes

An extensive study of the most important hydrodynamic characteristics of fairly large-scale bubble plumes was conducted using several measurement techniques and a variety of tools to analyze the data. Particle image velocimetry (PIV), double-tip optical probes (OP) and photographic techniques were extensively applied to measure bubble and liquid velocities, void-fraction and bubble sizes. PIV measurements in a vertical plane crossing the centre of the injector provided the instantaneous velocity fields for both phases, as well as hydrodynamic parameters, such as the movement of the axis of the plume and its instantaneous shape. Statistical studies were performed using image processing to determine the distribution of the apparent instantaneous plume diameter and centreline position. An important finding was that there is little instantaneous spreading of the bubble plume core; the spreading of the time-averaged plume width (as measured from the time-averaged void-fraction and time-averaged liquid velocity fields) is largely due to plume meandering and oscillations. The liquid-phase stress tensor distributions obtained from the instantaneous velocity data indicate that, for the continuous phase, these stresses scale linearly with the local void-fraction in the range of 0.5% < α < 2.5%. The bubbles were found to be ellipsoidal, with shape factor e ≈ 0.5.     

Local gas and liquid phase velocity measurement in a miniature stirred vessel using PIV combined with a new image processing algorithm

A method which combines standard two-dimensional particle image velocimetry (PIV) with a new image processing algorithm has been developed to measure the average local gas bubble velocities, as well as the local velocities of the liquid phase, within small stirred vessel reactors. The technique was applied to measurements in a gas–liquid high throughput experimentation (HTE) vessel of 45 mm diameter, but it is equally suited to measurements in larger scale reactors. For the measurement of liquid velocities, 3 μm latex seeding particles were used. For gas velocity measurements, a separate experiment was conducted which involved doping the liquid phase with fluorescent Rhodamine dye to allow the gas–liquid interfaces to be identified. The analysis of raw PIV images enabled the detection of bubbles within the laser plane, their differentiation from obscuring bubbles in front of the laser plane, and their use in lieu of tracer particles for gas velocity analysis using cross-correlation methods. The accuracy of the technique was verified by measuring the velocity of a bubble rising in a vertical glass column. The new method enabled detailed velocity fields of both phases to be obtained in an air–water system. The overall flow patterns obtained showed a good qualitative agreement with previous work in large scale vessels. The downward liquid velocities above the impeller were greatly reduced by the addition of the gas, and significant differences between the flow patterns of the two-phases were observed.     

Bubble and liquid turbulence characteristics of bubbly flow in a large diameter vertical pipe

The bubble and liquid turbulence characteristics of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. The bubble characteristics were measured using a dual optical probe, while the liquid-phase turbulence was measured using hot-film anemometry. Measurements were performed at six liquid superficial velocities in the range of 0.2–0.68 m/s and gas superficial velocity from 0.005 to 0.18 m/s, corresponding to an area average void fraction from 1.2% to 15.4%. At low void fraction flow, the radial void fraction distribution showed a wall peak which changed to a core peak profile as the void fraction was increased. The liquid average velocity and the turbulence intensities were less uniform in the core region of the pipe as the void fraction profile changed from a wall to a core peak. In general, there is an increase in the turbulence intensities when the bubbles are introduced into the flow. However, a turbulence suppression was observed close to the wall at high liquid superficial velocities for low void fractions up to about 1.6%. The net radial interfacial force on the bubbles was estimated from the momentum equations using the measured profiles. The radial migration of the bubbles in the core region of the pipe, which determines the shape of the void profile, was related to the balance between the turbulent dispersion and the lift forces. The ratio between these forces was characterized by a dimensionless group that includes the area averaged Eötvös number, slip ratio, and the ratio between the apparent added kinetic energy to the actual kinetic energy of the liquid. A non-dimensional map based on this dimensionless group and the force ratio is proposed to distinguish the conditions under which a wall or core peak void profile occurs in bubbly flows.     

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.     

Energy spectra and bubble velocity distributions in pseudo-turbulence: Numerical simulations vs. experiments

Direct numerical simulations (DNS) are performed to study the behavior of a swarm of rising air bubbles in water, employing the front tracking method, which allows to handle finite-size bubbles. The swarms consist of monodisperse deformable 4 mm bubbles with a gas fraction of 5% and 15%. This paper focuses on the comparison of the liquid energy spectra and bubble velocity probability density functions (PDFs) with experimental data obtained by phase-sensitive constant-temperature anemometry (CTA) and three-dimensional particle tracking velocimetry (PTV), respectively.     

Flow structure of microbubble-laden turbulent channel flow measured by PIV combined with the shadow image technique

The turbulence structure of a horizontal channel flow with microbubbles is experimentally investigated using combined particle image velocimetry (PIV) in order to clarify the mechanism of drag reduction caused by microbubbles. A new system which simultaneously measures the liquid phase and the dispersed bubbles is proposed, based on a combination of particle tracking velocimetry (PTV), laser-induced fluorescence (LIF) and the shadow image technique (SIT). To accurately obtain the velocity of the liquid phase, tracer particles which overlap with the bubble shadow images are almost entirely eliminated in the post-processing. Finally, the turbulence characteristics of the flow field are presented, including measurements for both phases, and the bubble effect on the turbulence is quantified.     

Effect of bubbles on the turbulence near the exit of a liquid jet

The objective of this paper is to examine the effect of bubbles on the turbulence levels of a water jet. Simultaneous measurements of the axial and radial velocity components were taken in a bubbly jet with a Laser Doppler Velocimeter (LDV) and then compared to the velocities of a single phase jet at the same liquid flow rate. Mean bubble diameters ranged from 0.6 to 2 mm and the void fractions were up to about 20%. The liquid Reynolds numbers were from 5,000 to 10,000 approximately. The measurements extended to from an axial distance of 4–12 cm. It was observed that bubbles did not affect significantly the average velocity profiles in the jet. However bubbles increased the turbulence intensities in the core of the jet near the jet exit. The increase in turbulence intensities was more pronounced at lower Reynolds numbers and at higher void fractions.     

Study of periodically excited bubbly jets by PIV and double optical sensors

Interactions between large coherent structures and bubbles in two-phase flow can be systematically observed in a periodically excited bubbly jet. Controlled excitation at fixed frequency causes large eddy structures to develop at regular intervals. Thus, interactions between large vortices and bubbles can be studied with PIV and double optical sensors (DOS) using phase-averaging techniques. A number of results on the time and space dependence of velocities and void fractions are presented revealing physical interactions between the liquid flow field and bubble movement as well as feedbacks from bubble agglomeration on the development of flow structures. A clear indication of bubble trapping inside the vortex ring is the generation of a bubble ring that travels with the same velocity as the vortex ring. The DOS results indicate clustering of the bubbles in coherent vortex structures, with a periodic variation of void fraction during the excitation period.     

Effect of gas expansion on the velocity of a Taylor bubble: PIV measurements

The effect of gas expansion on the velocity of a Taylor bubble was studied experimentally. The velocity field in the liquid ahead of a Taylor bubble was measured by particle image velocimetry (PIV), and the bubble velocity was measured with two pairs of laser diodes and photocells. The experiments were done in a 7.0 m long vertical tube with a 32 mm internal diameter. Solutions of carboxymethylcellulose (CMC) polymer with weight percentages between 0.01% and 0.1% were used. The expansion of slug gas induces an increase in the bubble velocity and a corresponding displacement of the liquid ahead of the bubble. The velocity of the bubble increases by an amount equal to the maximum velocity in the liquid displaced. For the solutions studied, the induced velocity profile was parabolic and the bubble velocity increase was equal to the liquid velocity at the tube axis, i.e., twice the mean velocity in the liquid displaced. The corrected velocity obtained by subtracting the velocity increase from the value of the bubble velocity is independent of the bubble length.     

Green water void fraction due to breaking wave impinging and overtopping

The present study uses laboratory measurements to investigate the void fraction of an overtopping flow on a structure. The overtopping flow, also called green water, was generated by the impingement of a plunging breaking wave on the structure following the Froude similarity of an extreme hurricane wave and a simplified offshore structure. The flow is multi-phased and turbulent with significant aeration. A fiber optic reflectometer (FOR) and bubble image velocimetry (BIV) were employed to measure the void fraction and velocity in the flow, respectively, and to determine the water level on the deck. Mean properties of void fraction and velocity were obtained by ensemble-averaging and time-averaging the repeated instantaneous measurements. The temporal and spatial distributions of void fraction reveal that the flow is very highly aerated near the front of green water and has relatively low aeration near the deck surface. The mean void fraction and velocity distributions were also depth-averaged for simplicity and potential use in engineering applications. Using the measured data, similarity profiles for depth-averaged void fraction, depth-averaged velocity, and water level were found. The study suggests that using only the velocity data is insufficient if the flow momentum or the flow rate is to be determined. The accuracy of the void fraction measurements was validated by comparing the directly measured water volume of the overtopping flow with the calculated water volume based on the measured velocity and void fraction.     

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

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

Laser-doppler measurement of the velocity and diameter of bubbles using the triple-peak technique

The velocity and diameter of nitrogen bubbles rising in quiescent water were simultaneously measured by a dual-beam laser-Doppler anemometer. The range of bubble diameters was 0.8–1.8 mm. The triple-peak technique was used to obtain the bubble diameter from the anemometer signal. Using a single stream of bubbles, the horizontal area of the bubble-detection region was proportional to the square of the bubble diameter. Using additional streams of bubbles to generate an optical disturbance, the probability of detecting a valid signal decreased linearly with increasing bubble flowrate of the additional streams and was independent of their location and number. The maximum void fraction was 0.01.     

Experimental study on turbulent natural convection heat transfer in water with sub-millimeter-bubble injection

Using thermocouples and a particle tracking velocimetry technique, temperature and velocity measurements are conducted to investigate flow and heat transfer characteristics of turbulent natural convection from a vertical heated plate in water with sub-millimeter-bubble injection. Hydrogen-bubbles generated by the electrolysis of water are used as the sub-millimeter-bubbles. In the turbulent region, the heat transfer deterioration occurs for a bubble flow rate Q = 33 mm³/s, while the heat transfer enhancement occurs for Q = 56 mm³/s. Temperature and velocity measurements suggest that the former is caused by a delay of the transition due to the bubble-induced upward flow. On the other hand, the latter is mainly due to two factors: one is the enhancement of the rotation of eddies in the outer layer, and the other is the increase in the gradient of the streamwise liquid velocity at the heated wall. These are caused by bubbles, which are located in the inner layer, rising at high speed.