Mechanism on lateral distribution of small bubbles in?vertical bubbly upflows

The void fraction profile and distributing mechanism of small bubbles were investigated in detail by using the Euler?CLagrange two-way model. The continuous phase is simulated by using direct numerical simulations in Euler frame of reference. Bubble trajectories are tracked by solving the Newton equations of motion with considering forces including gravity force (buoyant force), drag force, lift force, wall lift force, pressure gradient force, virtual mass force, and inertia force in Lagrange frame of reference. Results show that the void fraction profile of small bubbles exhibits double peak values, which is caused by different interphase forces dominating over bubble movements at different bubble places, and paths of bubbles are directly related to their injection places.     

A DNS study of laminar bubbly flows in a vertical channel

Direct numerical simulations are used to examine laminar bubbly flows in vertical channels. For equal size nearly spherical bubbles the results show that at steady state the number density of bubbles in the center of the channel is always such that the fluid mixture there is in hydrostatic equilibrium. For upflow, excess bubbles are pushed to the walls, forming a bubble rich wall-layer, one bubble diameter thick. For downflow, bubbles are drawn into the channel center, leading to a wall-layer devoid of bubbles, of a thickness determined by how much the void fraction in the center of the channel must be increased to reach hydrostatic equilibrium. The void fraction profile can be predicted analytically using a very simple model and the model also gives the velocity profile for the downflow case. For the upflow, however, the velocity increase across the wall-layer must be obtained from the simulations. The slip velocity of the bubbles in the channel core and the velocity fluctuations are predicted reasonably well by results for homogeneous flows.     

Four-way coupled Eulerian–Lagrangian Direct Numerical Simulations in a vertical laminar channel flow

Direct Numerical Simulations of a laminar two-phase flow into a vertical channel are investigated. An Eulerian–Lagrangian approach allows tracking each bubble position with a four-way coupling strategy, i.e. taking into account bubble-fluid and bubble-bubble interactions. The flow configuration has been chosen to highlight the buoyancy effects due to significant values of void fraction (high numbers of bubbles); hence the bubbles collisions and wall effects are the critical parameters to ensure the dispersion of the bubble plume. The DNS approach is self-consistent and does not rely on closure relations or empirical correlations for describing the collective bubble dynamics. It is found that the DNS predicts well the behavior of the bubble plume and its back effect on the liquid phase when compared with a mixture model and experimental data. The elastic nature of collisions, the sensitivity of the mean and RMS values of velocities and void fraction to the mesh quality are explored.     

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.     

An experimental study on the effect of air bubble injection on the flow induced rotational hub

Modification of shear stress due to air bubbles injection in a rotary device was investigated experimentally. Air bubbles inject to the water flow crosses the neighbor of the hub which can rotate just by water flow shear stresses, in this device. Increasing air void fraction leads to decrease of shear stresses exerted on the hub surface until in high void fractions, the hub motion stopped as observed. Amount of skin friction decrease has been estimated by counting central hub rotations. Wall shear stress was decreased by bubble injection in all range of tested Reynolds number, changing from 50,378 to 71,238, and also by increasing air void fraction from zero to 3.06%. Skin friction reduction more than 85% was achieved in this study as maximum measured volume of air fraction injected to fluid flow while bubbles are distinct and they do not make a gas layer. Significant skin friction reduction obtained in this special case indicate that using small amount of bubble injection causes large amount of skin friction reduction in some rotary parts in the liquid phases like as water.     

Simulating air entrainment and vortex dynamics in a hydraulic jump

The air entrainment characteristics of three separate Froude number hydraulic jumps are investigated numerically using an unsteady RANS, realizable k–ε turbulence model, with a Volume of Fluid treatment for the free surface. Mean velocity profiles, average void fraction, and Sauter mean diameter compare favorably with experimental data reported in literature. In all simulations, time-averaged void fraction profiles show good agreement with experimental values in the turbulent shear layer and an accurate representation of interfacial aeration at the free surface. Sauter mean diameter is well represented in the shear layer, and free surface entrainment results indicate bubble size remains relatively unchanged throughout the depth of the jump. Several different grid resolutions are tested in the simulations. Significant improvements in void fraction and bubble size comparison are seen when the diameter to grid size ratio of the largest bubbles in the shear layer surpasses eight. A three-dimensional simulation is carried out for one Froude number jump, showing an improvement in the prediction of entrained air and bubble size compared with two-dimensional results at a substantial increase in computation time. An analysis of three-dimensional vorticity shows a complex interaction between spanwise and streamwise vortical structures and entrained air bubbles. The jump is similar to a turbulent mixing layer, constrained by the free surface, with vortex pairing and subsequent fluctuations in free surface elevation. Downstream fluctuations of the toe are associated with a roll up of the primary spanwise vortex, fluctuations of the free surface, and counter-rotating streamwise vortex pairs. The action of these flow structures is likely responsible for the improvement in three-dimensional results.     

Dispersion of Gas Bubbles in Large-Scale Homogeneous Isotropic Turbulence

An approximate analysis is given of the dispersion of gas bubbles that rise at large Reynolds number through large-scale homogeneous, isotropic turbulence, characterized by the Kraichnan energy-spectrum function. A fairly well-established equation of motion of the bubbles, originally proposed by Thomas et al. [16], is used to derive a closed set of equations for the components of the dispersion tensor of the bubbles in a manner analogous to that used by Saffman [12] for fluid particles and by Pismen and Nir [10, 11] for solid particles. The equations are then solved to obtain the diffusivities and the intensities of bubble velocity fluctuations. Analytical solutions are compared with results from simulations of the bubble motion in a Gaussian random velocity field.     

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.     

垂直壁面附近上升单气泡的弹跳动力学研究

采用高速摄影技术结合阴影法,对静止水中垂直壁面附近上升单气泡运动进行实验研究,对比气泡尺度及气泡喷嘴与壁面之间的初始无量纲距离 ($S^{\ast}$)对气泡上升运动特性的影响,分析气泡与壁面碰撞前后,壁面效应与气泡动力学机制及能量变化规律.结果表明,对于雷诺数$Re \approx 580 \sim 1100$,无量纲距离$S^{\ast } <2 \sim3$时,气泡与壁面碰撞且气泡轨迹由无约束条件下的三维螺旋转变成二维之字形周期运动;当$S^{\ast } >2 \sim3$时,壁面效应减弱,有壁面约束的气泡运动与无约束气泡运动特性趋于一致.气泡与壁面碰撞前后,壁面效应导致横向速度峰值下降为原峰值的70%,垂直速度下降50%;气泡与壁面碰撞前,通过气泡中心与壁面距离($x/R$)和修正的斯托克斯数相关式可预测垂直速度的变化规律.上升气泡与壁面碰撞过程中,气泡表面变形能量单向传输给气泡横向动能,使得可变形气泡能够保持相对恒定的弹跳运动.提出了气泡在与壁面反复弹跳时的平均阻力系数的预测模型,能够很好地描述实验数据反映出的对雷诺数${Re}$、韦伯数${We}$和奥特沃斯数${Eo}$等各无量纲参数的标度规律.      

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

An efficient method for the prediction of the motion of individual bubbles

We present a method for the simulation of two-phase flows which can be applied to problems characterised by the presence of up to several hundreds of gas bubbles. The bubble model is kept simple, requiring only six parameters to describe the shape of a single bubble. The model is coupled to a conventional time discrete finite-volume scheme for the solution of the Navier–Stokes equations by the density field which is calculated on basis of the information on the positions and the shapes of the bubbles before each time step. The motion of the bubbles is in turn calculated from an analysis of the computed flow field. Systematical errors due to simplifications are eliminated by the introduction of correction factors. For a selection of fluid dynamical problems, the results of simulations using the method are compared to experimental data. Good quantitative agreement could be found.     

On the numerical study of bubbly flow created by ventilated cavity in vertical pipe

The dispersion of bubbles into a down-liquid flow in a vertical pipe is investigated. At low flow rates, the intended design of a swarm of discrete bubbles is achieved. At high flow rates, a ventilated cavity is nonetheless formed, which is attached close to the gas sparger. Behind this ventilated cavity, three different flow regimes characterize the complex bubbly flow field downstream of the down-liquid flow: vortex region with high void fraction, transitional region and pipe flow region. In this study, a numerical model that solved the entire development of the gas–liquid flow including the extended single-phase liquid region upstream to the wall-jet and recirculating-vortex zones in order to allow a more realistic determination of the boundary conditions of the down-liquid flow was adopted. Coupling with the Eulerian–Eulerian two-fluid model to solve the respective gas and liquid phases, a population balance model was also applied to predict the bubble size distribution in the wake right below the cavity base as well as further downstream in the transitional and fully-developed pipe flow regions. The numerical model was evaluated by comparing the numerical results against the data derived from theoretical, numerical and experimental approaches. Prediction of the Sauter mean bubble diameter distributions by the population balance approach at different axial locations confirmed the dominance of breakage due to the high turbulent intensity below the ventilated cavity which led to the generation of small gas bubbles at high void fraction. Further downstream, the coalescence effect dominated leading to merging of the small bubbles to form bigger bubbles.     

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.     

Characterization of bubble mobility in channel flow with fibrous porous media walls

In composites processing, resin is introduced into a fibrous domain to cover all the empty spaces between the fibers. It is important to extract air bubbles from the domain before the resin solidifies. Failure to do so will entrap these voids in the final part, which is detrimental to its performance. Hence, there is a need to understand bubble motion in a fibrous porous domain in which the bubbles move with the resin in channels surrounded by fibrous walls. A rising bubble model is presented that consists of a single spherical void in a cylindrical axisymmetric two-phase domain of resin and air surrounded by porous media boundaries. The motion of a bubble in a channel flow with porous boundaries is modeled by replacing the walls with a slip velocity. Focus is on how the porous media permeability influences the bubble motion. A parameter called bubble mobility is defined as the ratio of bubble rise velocity to the resin free surface velocity. Results suggest that fabric permeability and fluid properties can be optimized to increase bubble mobility and ultimately lead to reduction in void content during composites processing.     

Air–water flow in a vertical pipe: experimental study of air bubbles in the vicinity of the wall

This study deals with the influence of bubbles on a vertical air–water pipe flow, for gas-lift applications. The effect of changing the bubble size is of particular interest as it has been shown to affect the pressure drop over the pipe. Local measurements on the bubbles characteristics in the wall region were performed, using standard techniques, such as high-speed video recording and optical fibre probe, and more specific techniques, such as two-phase hot film anemometry for the wall shear stress and conductivity measurement for the thickness of the liquid film at the wall. The injection of macroscopic air bubbles in a pipe flow was shown to increase the wall shear stress. Bubbles travelling close to the wall create a periodic perturbation. The injection of small bubbles amplifies this effect, because they tend to move in the wall region; hence, more bubbles are travelling close to the wall. A simple analysis based on a two-fluid set of equations emphasised the importance of the local gas fraction fluctuations on the wall shear stress.     

Simulation of void fraction profile evolution in subcooled nucleate boiling flow in a vertical annulus using a bubble-tracking approach

A three-dimensional bubble-tracking model of subcooled nucleate boiling flow in a vertical channel at low-pressure conditions is proposed with specific application to the case of boiling in an annulus with a central heating rod. Vapour is distributed in the liquid in the form of individually tracked bubbles. The overall behaviour of the liquid–vapour system results from motion, interaction, coalescence and boiling mechanisms prescribed mostly at the level of bubbles. Bubbles are nucleated at nucleation sites randomly distributed over the heated surface. After nucleation, bubbles slide on the heated wall, detach and then migrate into the lower-temperature region away from the heated surface, where they condense. The proposed model was applied to experiments on subcooled boiling from Purdue University (USA). Experimental and calculated void fraction radial profiles at different axial locations are compared.     

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.