Bubble splitting under gas–liquid–liquid three‐phase flow in a double T‐junction microchannel

Gas–aqueous liquid–oil three‐phase flow was generated in a microchannel with a double T‐junction. Under the squeezing of the dispersed aqueous phase at the second T‐junction (T2), the splitting of bubbles generated from the first T‐junction (T1) was investigated. During the bubble splitting process, the upstream gas–oil two‐phase flow and the aqueous phase flow at T2 fluctuate in opposite phases, resulting in either independent or synchronous relationship between the instantaneous downstream and upstream bubble velocities depending on the operating conditions. Compared with two‐phase flow, the modified capillary number and the ratio of the upstream velocity to the aqueous phase velocity were introduced to predict the bubble breakup time. The critical bubble breakup length and size laws of daughter bubbles/slugs were thereby proposed. These results provide an important guideline for designing microchannel structures for a precise manipulation of gas–liquid–liquid three‐phase flow which finds potential applications among others in chemical synthesis. © 2017 American Institute of Chemical Engineers AIChE J, 63: 376–388, 2018     

Bubble breakup with permanent obstruction in an asymmetric microfluidic T‐junction

Bubble breakup with permanent obstruction in an asymmetric microfluidic T‐junction is investigated experimentally. The breakup process of bubbles can be divided into three stages: squeezing, transition, and pinch‐off stages. In the squeezing stage, the thinning of the bubble neck is mainly controlled by the velocity of the fluid flowing into the T‐junction, and the increase of the liquid viscosity can promote this process. In the transition stage, the minimum width of bubble neck decreases linearly with time. In the pinch‐off stage, the effect of the velocity of the fluid flowing into the T‐junction on the thinning of the bubble neck becomes weaker, and the increase of the liquid viscosity would delay this process. The evolution of the minimum width of the bubble neck with the remaining time before the breakup can be scaled by a power–law relationship. The bubble length has little influence on the whole breakup process of bubbles. © 2014 American Institute of Chemical Engineers AIChE J, 61: 1081–1091, 2015     

Breakup dynamics of slender bubbles in non‐newtonian fluids in microfluidic flow‐focusing devices

This study aims to investigate the breakup of slender bubbles in non‐Newtonian fluids in microfluidic flow‐focusing devices using a high‐speed camera and a microparticle image velocimetry (micro‐PIV) system. Experiments were conducted in 400‐ and 600‐μm square microchannels. The variation of the minimum width of gaseous thread with the remaining time before pinch‐off could be scaled as a power‐law relationship with an exponent less than 1/3, obtained for the pinch‐off of bubbles in Newtonian fluids. The velocity field and spatial viscosity distribution in the liquid phase around the gaseous thread were determined by micro‐PIV to understand the bubble breakup mechanism. A scaling law was proposed to describe the size of bubbles generated in these non‐Newtonian fluids at microscale. The results revealed that the rheological properties of the continuous phase affect significantly the bubble breakup in such microdevices. © 2012 American Institute of Chemical Engineers AIChE J,, 2012     

非对称Y型分岔微通道内气泡破裂与分配规律

利用高速摄像仪对气泡在非对称Y型微通道分岔口的破裂行为和分配规律进行了实验研究。采用氮气（N₂）作为分散相,含0.3%表面活性剂十二烷基硫酸钠（SDS）的蒸馏水-甘油（质量分数分别为20%、40%、50%）溶液为连续相。在分岔口处观察到了3种不同的气泡行为：无间隙的不对称破裂、有间隙的不对称破裂以及不破裂。考察了气泡破裂和不破裂行为之间的转变,并与文献进行了比较。考察了两相流率及物性对破裂气泡分配规律的影响。结果表明：破裂后两个子气泡的长度均随气相流量与气泡长度的增大而增大,随液相流量和黏度的增大而减小。随液相速度和黏度的增大,气泡破裂的不对称程度减弱。     

A unified theoretical model for breakup of bubbles and droplets in turbulent flows

Pressure has a significant effect on bubble breakup, and bubbles and droplets have very different breakup behaviors. This work aimed to propose a unified breakup model for both bubbles and droplets including the effect of pressure. A mechanism analysis was made on the internal flow through the bubble/droplet neck in the breakup process, and a mathematical model was obtained based on the Young–Laplace and Bernoulli equations. The internal flow behavior strongly depended on the pressure or gas density, and based on this mechanism, a unified breakup model was proposed for both bubbles and droplets. For the first time, this unified breakup model gave good predictions of both the effect of pressure or gas density on the bubble breakup rate and the different daughter size distributions of bubbles and droplets. The effect of the mother bubble/droplet diameter, turbulent energy dissipation rate and surface tension on the breakup rate, and daughter bubble/droplet size distribution was discussed. This bubble breakup model can be further used in a population balance model (PBM) to study the effect of pressure on the bubble size distribution and in a computational fluid dynamics‐population balance model (CFD‐PBM) coupled model to study the hydrodynamic behaviors of a bubble column at elevated pressures. © 2014 American Institute of Chemical Engineers AIChE J, 61: 1391–1403, 2015     

Three‐Dimensional Simulation of Bubble Formation Through a Microchannel T‐Junction

Three‐dimensional simulations of bubble formation in Newtonian and non‐Newtonian fluids through a microchannel T‐junction are conducted by the volume‐of‐fluid method. For Newtonian fluids, the critical capillary number Ca for the transition of the bubble breakup mechanism is dependent on the velocity ratio between the two phases and the microchannel dimension. For the power law fluid, the bubble diameter decreases and the generation frequency increases with higher viscosity parameter K and power law index n. For a Bingham fluid, the viscous force plays a more important role in microbubble formation. Due to the yield stress τ_y, a high‐viscous region is developed in the central area of the channel and bubbles deform to a flat ellipsoid shape in this region. The bubble diameter and generation frequency are almost independent of K.     

Slug bubble deformation and its influence on bubble breakup dynamics in microchannel

The deformation of moving slug bubbles and its influence on the bubble breakup dynamics in microchannel were studied. Three bubble morphologies were found in the experiment: slug, dumbbell and grenade shapes. The viscosity effect of continuous phase aggravates the velocity difference between the fluid near the wall and the bubble, resulting in that the continuous phase near the bubble head flows towards and squeezes the bubble tail, which causes the deformation of bubbles. Moreover, the experimental results show that the deformation of bubbles could significantly prolong the bubble breakup period at the downstream Y-junction. There exists the critical capillary number Ca_Cr for the asymmetric breakup of grenade bubbles, Ca_Cr increases with the rise of flow rate and viscosity of the continuous phase.     

对称分支并行微通道中气液两相流的均匀性规律

采用高速摄像系统研究了对称分支形并行微通道内气液两相流及弹状气泡均匀性规律。实验中分别采用含0.3% SDS的甘油-水溶液与氮气作为液相和气相。观察到弹状流和泡状流两种流型,作出了由两相操作条件构成的流型图及流型转变线。结果表明,气泡非均匀性主要由两微通道内流体之间的相互作用、下游通道中流体动力学的反馈作用以及通道制造误差造成。随液相黏度增大,气泡均匀性变好;在高液相流量以及低气相压力下操作,气泡尺寸分布更易达到均匀。基于压力降守恒原理和微通道内气液两相流阻力模型,构建了两通道中气泡尺寸的预测模型。     

Effect of coalescence and breakup on bubble size distributions in a two-dimensional packed bed

Bubble size distributions in a two-dimensional packed bed are investigated as a function of axial direction by using image processing techniques with a large number of bubble samples. Two inlet conditions, controlled sized bubbles and uncontrolled sized bubbles, are conducted to study the characteristic behaviors of bubbles. With the uncontrolled sized bubbles, the average bubble size corresponding to the two-dimensional bed is found, and is not affected by decrease or increase of flow rates. With the controlled sized bubbles, dominant bubble breakup and coalescence flows are separately simulated to investigate bubble breakup and coalescence rates. Unique behavior of bubble size distributions for dominant bubble breakup and coalescence has been seen, and changes in bubble size distributions along axial direction are studied with median bubble size. Near the inlet the median changes rapidly due to the dominant bubble mechanism of either coalescence or breakup, and far away from the inlet the median reaches asymptotic value due to the balance of bubble breakup and coalescence. For both dominant breakup and coalescence flows, the asymptotic values are close to the average bubble size. Therefore, the average bubble size, resulting from the balance of bubble coalescence and breakup at far downstream, is irrespective of inlet flow conditions.     

The effect of system pressure on gas‐liquid slug flow in a microchannel

Characteristics of gas‐liquid two‐phase flow under elevated pressures up to 3.0 MPa in a microchannel are investigated to provide the guidance for microreactor designs relevant to industrial application. The results indicate that a strong leakage flow through the channel corners occurs although the gas bubbles block the channel. With a simplified estimation, the leakage flow is shown to increase with an increase in pressure, leading to a bubble formation shifting from transition regime to squeezing regime. During the formation process, the two‐phase dynamic interaction at the T‐junction entrance would have a significant influence on the flow in the main channel as the moving velocity of generated bubbles varies periodically with the formation cycle. Other characteristics such as bubble formation frequency, bubble and slug lengths, bubble velocities, gas hold‐up, and the specific surface area are also discussed under different system pressures. © 2013 American Institute of Chemical Engineers AIChE J, 60: 1132–1142, 2014     

On the leakage flow around gas bubbles in slug flow in a microchannel

The leakage flow is that liquid does not push gas bubbles and leaks through the channel corners. This leakage flow was confirmed by tracking particles moving in the liquid film with a double light path method and was quantified by tracking the gas–liquid interface movement. The results show that leakage flow varies during bubble formation process. The average net leakage flow Q_net‐leak in a bubble formation cycle at T‐junction can be as large as 62.4% of the feeding liquid flow rate, depending on the liquid properties. Q_net‐leak for regular flow at main channel is much smaller, ranging from about 0 to 30% of the feeding liquid flow rate. The difference between the two leakage flows would lead to an increase in liquid slug length after generation. Finally, the effects of parameters such as phase flow rates, surface tension, and viscosity were investigated. © 2015 American Institute of Chemical Engineers AIChE J, 61: 3964–3972, 2015     

A visualized investigation of bubble breakup in a swirl-venturi bubble generator

The dynamics and breakup of bubbles in swirl-venturi bubble generator (SVBG) are explored in this work. The three-dimensional movement process and breakup phenomena of bubbles are captured by one high-speed camera system with two cameras while the distribution of swirling flow field is recorded through Particle Image Velocimetry technology. It is revealed that bubbles have two motion trajectories, which are deeply related to bubble breakup. One trajectory is that mother bubble moves upward in an axial direction of the SVBG to the diverging section, and the other trajectory is that mother bubble rotates obliquely upward to another side-wall along the radial direction. Meanwhile, binary breakup, shear-off-induced breakup, static erosive breakup, and dynamic erosive breakup are observed. For relatively high liquid Reynolds number, vortex flow regions are extended and the bubble size is reduced. Furthermore, it is worth noting that the number of microbubbles increases significantly for intensive swirling flow.     

Bubble mechanisms and characteristics at pore scale in a packed-bed reactor

An image processing technique was used to study dominant bubble mechanisms in a two-dimensional packed-bed at pore level under the bubbly flow regime. Bubble breakup and coalescence were identified as dominant mechanisms using a large number of image samples. Two types of coalescence mechanisms were identified that occur due to compression and deceleration associated with the bubbles and three breakup mechanisms were identified that are result of liquid shear force, bubble acceleration, and bubble impact. Data on various two-phase parameters, such as local void fraction, bubble velocity, size, number, and shape were obtained from the images. Results indicated that when a flow regime changed from bubbly to either trickling or pulsing flow, the number of average sized bubbles significantly decreased and the shape of the majority of the bubbles was no longer spherical. Although a mean bubble velocity of all sized bubbles was uniform for given gas and liquid superficial velocities, individual bubble velocities were quite different depending on the bubble location in the pore. The present bubble size distributions were compared with previous studies and the results on bubble size are in general agreement.     

CFD-PBM耦合模型模拟气液鼓泡床的通用性研究

通过对不同操作压力和不同液体性质气液鼓泡床的模拟值与实验数据进行对比，从而验证CFD-PBM耦合模型的通用性。结果表明，CFD-PBM耦合模型在加入了气泡破碎修正因子后，可以很好地预测压力对鼓泡床流体力学行为的影响趋势，当压力升高时，气含率显著升高。不同液体黏度和表面张力条件下CFD-PBM耦合模型的模拟结果与实验结果均吻合较好。随液体黏度增大，气泡破碎速率减小，气泡尺寸分布变宽，曳力显著下降，气含率随之降低。随表面张力减小，气泡破碎速率增大，气泡变小，气含率升高。CFD-PBM耦合模型具有很好的通用性，原因在于考虑了压力、液体黏度和表面张力对气泡聚并、破碎和气液相间作用力的影响。     

Generation of micromonodispersed droplets and bubbles in the capillary embedded T‐junction microfluidic devices

This work focuses on the dispersion of micromonodispersed droplets and bubbles in the capillary embedded T‐junction microfluidic devices. The effects of the microchannel structure, operating conditions, and physical properties on the dispersion rules were carefully investigated. It was found that the extended capillary could greatly affect the dispersion rules, which was favorable for reducing the dispersed size. The dispersed size was mainly dominated by the Ca number, and the effects of dispersed phase flow rate and viscosity ratio of the two phases were also very important. The dispersion mechanism and size rules in the capillary embedded microfluidic devices were discussed seriously by comparing the similarities and differences of the liquid/liquid and gas/liquid dispersion processes. © 2010 American Institute of Chemical Engineers AIChE J, 2011     

Regime transition in viscous and pseudo viscous systems: A comparative study

A comprehensive quantitative study on the effect of liquid viscosity (1 ≤ µ_L ≤ 1149 mPa‐s) on the local flow phenomena of the gas phase in a small diameter bubble column is performed using ultrafast electron beam X‐ray tomography. The internal dynamic flow structure and the bubble size distribution shows a dual role of the liquid viscosity on the hydrodynamics. Further, the effect of solid concentration (C_s = 0.05, 0.20) on the local flow behavior of the gas phase is studied for the pseudo slurry viscosities similar to the liquid viscosities of the gas–liquid systems. The effects of liquid and pseudo slurry viscosities on flow structure, bubble size distribution, and gas phase distribution are compared. The bubble coalescence is significantly enhanced with the addition of particles as compared to the system without particles for apparently same viscosity. The superficial gas velocity at which transition from homogeneous bubbly to slug flow regime occurs is initiated by the addition of particles as compared to the particle free system for apparently same viscosity. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3079–3090, 2014     

Dynamics of bubble breakup in a microfluidic T-junction divergence

The aim of this work is to investigate experimentally the bubble breakup in a microfluidic T-junction divergence using a high-speed digital camera and a micro-Particle Image Velocimetry (micro-PIV) system. The breakup and non-breakup of N₂ bubbles in glycerol–water mixtures with several concentrations of sodium dodecyl sulphate (SDS) as surfactant were studied with capillary number ranging from 0.001 to 0.1. The cross section of PMMA square microchannel is 400 μm wide and 400 μm deep. Four various flow patterns were observed at the T-junction by changing gas and liquid flow rates. The dynamics of three various types of symmetric breakup of bubbles were investigated. The symmetric breakup of bubbles type I is mainly controlled by the augmented pressure in liquid phase. The symmetric breakup of bubbles type II is controlled by both the increased pressure and viscous forces. In the symmetric breakup of bubbles type III, a scaling law for the minimum bubble neck and the remaining time during bubble breaking process were found. The transitions between breakup and non-breakup of bubbles were investigated, and a power–law relationship between bubble extension and capillary number was proposed to predict the transitions between adjacent regimes. Our experimental results reveal that the bubble breakup in a microfluidic T-junction divergence is similar to the droplet behaviours in such a device ( [Jullien et al., 2009] , [Leshansky and Pismen, 2009] and [Link et al., 2004] ).     

Study of Air‐Liquid Flow Patterns in Hydrocyclone Enhanced by Air Bubbles

In order to improve the oil‐water separation efficiency of a hydrocyclone, a new process utilizing air bubbles has been developed to enhance separation performance. Using the two‐component phase Doppler particle analyzer (PDPA) technique, the velocities of two phases, air and liquid, and air bubble diameter were measured in a hydrocyclone. The air‐liquid mixing pump can produce 15 to 60 μm‐diameter air bubbles in water. There is an optimum air‐liquid ratio for oil‐water separation of a hydrocyclone enhanced by air bubbles. An air core occurs in the hydrocyclone when the air‐liquid ratio is more than 1 %. The velocities of air bubbles have a similar flow pattern to the water phase. The axial and tangential velocity differences of the air bubbles at different air‐liquid ratio are greater near the wall and near the core of the hydrocyclone. The measured results show that the size distribution of the air bubbles produced by the air‐liquid mixing pump is beneficial to the process where air bubbles capture oil droplets in the hydrocyclone. These studies are helpful to understand the separation mechanism of a hydrocyclone enhanced by air bubbles.     

Bubble size in a forced circulation loop reactor

BACKGROUND: The bubble size distribution in gas‐liquid reactors influences gas holdup, residence time distribution, and gas‐liquid interfacial area for mass transfer. This work reports on the effects of independently varied gas and liquid flow rates on steady‐state bubble size distributions in a new design of forced circulation loop reactor operated with an air–water system. The reactor consisted of a cylindrical vessel (～26 L nominal volume, gas‐free aspect ratio ≈ 6, downcomer‐to‐riser cross‐sectional area ratio of 0.493) with a concentric draft tube and an annular riser zone. Both gas and liquid were in forced flow through a sparger that had been designed for minimizing the bubble size. RESULTS: Photographically measured bubble size distributions in the riser zone could be approximated as normal distributions for the combinations of gas and liquid flow rates used. This contrasted with other kinds of size distributions (e.g. bimodal, Gaussian) that have been reported for other types of gas‐liquid reactors. Most of the bubbles were in the 3 to 5 mm diameter range. At any fixed low value of aeration rate (≤1.8 × 10^?4 m³s^?1), increase in the liquid flow rate caused earlier detachment of bubbles from the sparger holes to reduce the Sauter mean bubble size in the riser region. CONCLUSION: Unlike in conventional bubble columns where bimodal and Gaussian bubble size distributions have been reported, a normal bubble size distribution is attained in forced circulation loop reactors with an air–water system over the entire range of operation. Copyright © 2007 Society of Chemical Industry     

Theoretical prediction of flow regime transition in bubble columns by the population balance model

Theoretical prediction of flow regime transition in bubble columns was studied based on the bubble size distribution by the population balance model (PBM). Models for bubble coalescence and breakup due to different mechanisms, including coalescence due to turbulent eddies, coalescence due to different bubble rise velocities, coalescence due to bubble wake entrainment, breakup due to eddy collision and breakup due to large bubble instability, were proposed. Simulation results showed that at relatively low superficial gas velocities, bubble coalescence and breakup were relatively weak and the bubble size was small and had a narrow distribution; with an increase in the superficial gas velocity, large bubbles began to form due to bubble coalescence, resulting in a much wider bubble size distribution. The regime transition was predicted to occur when the volume fraction of small bubbles sharply decreased. The predicted transition superficial gas velocity was about 4 cm/s for the air-water system, in accordance with the values obtained from experimental approaches.