适用欧拉-拉格朗日方法模拟气液泡状流的气泡破碎模型

欧拉-拉格朗日方法已被广泛应用于模拟鼓泡塔等气-液反应器内的流型、气泡尺寸（或气含率）及其分布。文献中该方法主要基于临界Weber数观点来描述气泡破碎行为,且破碎后的子气泡尺寸由随机数确定。但现有实验和理论研究表明,临界Weber数约束不能体现气体密度等物性参数和泡内气体重分布对气泡破碎行为的影响。针对这些不足,提出了适用欧拉-拉格朗日框架且考虑泡内气体重分布贡献的气泡破碎机理模型,并利用开源软件OpenFOAM开发了基于新破碎模型的求解器。新模型预测结果能较好地吻合实验测量的时均轴向液速、气泡尺寸及其分布等实验数据。特别地,考虑泡内气体重分布现象的破碎机理模型成功预测了实验观测的气泡尺寸双峰分布特征。     

流体振荡器气泡碎化特性数值模拟

为了探究流体振荡器中气泡碎化特性，采用重整化群(RNG) k-ε湍流模型和流体体积函数(VOF)多相流模型研究流体振荡器中的气液两相流。通过数值模拟分析了气泡在流场中的变形和破碎行为，考察不同流动参数对气泡破碎的影响。同时进行可视化实验，用实验结果验证了模拟的准确性。结果表明,流体振荡器可作为一种新型气泡发生装置。振荡频率越高，流体振荡器中的气泡碎化效果越好；随着液速增大，气泡变形程度加剧，破碎概率增加，生成子气泡也越多；剪切应力和涡流碰撞是导致气泡破碎的主要机制；射流振荡促使气泡向振荡腔壁面偏移，气泡初始直径越大，水平偏移量越大，越容易发生破碎。     

基于介尺度PBM模型的生物反应器放大模拟及实验研究

对于通气搅拌式工业生物反应器的放大设计而言,精确预测气泡尺寸和体积传质系数非常重要,因此需要建立合适的气泡聚并和破碎模型,以保证反应器的高效操作。以5 L通气搅拌式生物反应器为对象,以气泡尺寸和体积传质系数的实验数据为基准,模拟并考察了两种聚并模型和四种破碎模型对生物反应器内流体流动行为以及传质能力的影响。结果表明,基于介尺度理论的修正聚并模型与考虑黏流剪切的破碎模型组合,所得模拟结果与实验数据吻合最好,这为大型生物反应器的桨型优化提供了模型基础。因为工业化生物发酵通常是在大型生物反应器中进行,搅拌桨型对生物反应器效能至关重要,故本研究在选定最优气泡聚并破碎模型的基础上,通过叶轮末端剪切力相等的放大原则将5 L通气搅拌式工业生物反应器放大到400 m³,同时考察了六斜叶圆盘搅拌桨、非对称式抛物线搅拌桨、布鲁马金式搅拌桨以及六直叶圆盘搅拌桨等桨型组合对气泡破碎能力和气体分散效果的影响,并通过综合对比气含率、体积传质系数等参数,得到400 m³通气搅拌式生物反应器的最优桨型组合。     

计及气泡诱导与剪切湍流的气泡破碎、湍流相间扩散及传质模型

在以鼓泡塔为代表的气液鼓泡流动中,存在着气泡诱导湍流（BIT）和剪切湍流两种湍流机制,并且二者在不同的时间、空间范围内既相互竞争又共同作用。受制于BIT动能能谱的形式和特性不够完整清晰,过去的研究中关于BIT如何对气泡破碎聚并、相间作用力、相间传热传质等相间相互作用过程产生影响的结论比较模糊。因此,本文在具有波数κ^-3特性的BIT能谱的基础上,提出了在不同工况下考虑BIT与剪切湍流共同作用的研究思路。研究结果表明,考虑两种湍流机制的气泡破碎模型和湍流相间扩散模型对BIT在整体或局部占据不同程度主导地位的情况,都能很好地捕捉气液鼓泡流动的动力学特性,为进一步准确揭示气液相间传质过程的内在机理提供了基础。     

精馏塔板相界面积的粒数衡算模型

通过对气液湍流系统中气泡的动力学行为进行分析推导出气泡破碎速率与聚并速率的表达式,在此基础上建立了描述气泡分散特性的粒数衡算模型.求解粒数衡算方程可以计算精馏塔板上气泡的粒度分布以及气液相界面积,模拟计算的结果与实验测量的数据相当一致,证明可以利用粒数衡算模型较准确地预测气液湍流系统的分散性质.     

垂直圆管内湍流泡状流的数值研究

在经典Euler/Euler型水动力模型基础上,引入考虑不同直径气泡的种群平衡方程来描述气液两相泡状流,对液相和气相分别建立了基本方程,通过对气泡的受力分析并考虑气泡之间聚合和破碎效应后给出了本构方程,建立了封闭的双流体模型并用于垂直管道湍流泡状流的三维数值模拟.模型预测值与实验数据的比较结果表明该模型能较好地模拟垂直管道湍流泡状流中的相含率分布、速度分布、湍动能分布、气泡直径分布以及气泡直径分布的演变过程.     

鼓泡塔反应器气液两相流CFD数值模拟

对圆柱形鼓泡塔反应器内的气液两相流动进行了三维瞬态数值模拟,模拟的表观气速范围为0.02～0.30 m•s^-1; 模拟采用了双流体模型,并耦合了气泡界面密度单方程模型预测气泡尺寸,该模型考虑了气泡聚并与破碎对气泡尺寸的影响。液相湍流采用考虑气相影响的修正k-ε模型,两相间的动量传输仅考虑曳力作用。模拟获得了轴向气/液相速度分布、气含率分布、湍流动能分布以及气泡表面面积密度等,对部分模拟结果与实验值进行了定量比较,结果表明模拟结果与实验结果吻合较好。     

加压鼓泡塔反应器中气液两相流CFD数值模拟

对加压气液鼓泡塔反应器内的气液两相流进行了二维数值模拟,模拟的压力为0.5～2.0 MPa,表观气速为0.120～0.312 m/s;模拟采用了Euler-Euler模型,并耦合了气泡群平衡模型（PBM）预测气泡尺寸,该模型考虑了气泡聚并与破碎对气泡的影响。液相湍流采用标准k-ε模型,两相间的作用力只考虑曳力。模拟获得了局部气含率、局部气/液相时均轴向速度及其径向分布等数据,并与实验结果进行比较。结果表明,局部气含率、局部气相速度模拟结果与实验结果吻合较好,局部液相速度径向分布特征模拟结果与文献结果相符。     

精馏塔板上气液相界面积的测量与预测

利用现代的电子光学测试仪器和计算机图像处理技术较准确地测量了精馏塔板上鼓泡液体中气泡的粒径分布和相界面积.同时，还从理论上分析了湍流液体中气泡变形与破碎的机理，提出了预测气液相界面积的多相湍流动力学模型.模拟计算结果与实验测量数据的比较证明，此模型有较高的准确性，而且形式简单，应用方便.     

液相中气泡上升行为与界面传质:实验研究与数值计算

在工业生产过程中，气泡在液相中的上升行为及气液界面的传质行为极为常见。本文针对不同条件下气泡上升过程的实验研究方法以及数值计算方法进行了总结。从实验与数值计算的角度，综述了单气泡上升过程的影响因素、多气泡上升过程聚并与破裂的现象和机理以及工业装置中气液两相流型和气泡特性，并对传质模型进行了归纳，主要关注了气侧-界面传质模型的研究现状。综述结果表明：当前对于单气泡上升行为的研究较为充分，而对于多气泡的行为机理的研究尚需深入。此外，受到研究手段的限制，进行气侧-界面传质模型研究具有一定挑战性。针对当前的相关研究进展和存在的问题，对今后气泡上升行为和传质行为的研究提出以下建议，即开展气泡聚并与破裂可控性研究，强化对气侧-界面传质过程的研究，包括泡内流体行为可视化研究和相关传质模型的建立。     

3D Single Bubble Breakup Tracking in a Stirred Tank

A detailed understanding of turbulent fluid particle breakup mechanisms is essential for the accurate modeling of gas/liquid and liquid/liquid dispersions. The design of a fully automated setup for the three‐dimensional serial examination of the single bubble breakup process in a stirred tank, ensuring high repetition rates necessary for the additionally automated statistical analysis, is described. The implementation of a three‐dimensional automatic bubble breakup tracking tool is illustrated. At last, exemplary bubble breakup trajectories that show the benefits and limitations of the developed system and method are discussed.     

A literature review of theoretical models for drop and bubble breakup in turbulent dispersions

This paper presents a literature review on mechanisms and models for the breakage of bubbles and drops (fluid particles) in turbulent dispersions. For the mechanisms, four categories are summarized, namely, turbulence fluctuation, viscous shear stress, shearing-off process and interfacial instability. The models for breakup frequency and daughter size distribution available in literature are reviewed thoroughly. The development and limitation of the existing models are studied and possible improvements are proposed.     

Experimental measurement of bubble breakup in a jet bubbling reactor

The bubble breakups in a jet bubbling reactor are captured using a high-speed camera and the velocity field is measured by particle image velocimetry. Two typical breakup patterns, jet breakup and jet-vortex breakup are observed. The breakup time interval of the jet-vortex breakup is two orders of magnitude higher than the jet breakup. The probability of the jet-vortex breakup and the jet breakup accounting in the total breakup events increases and decreases with the jet velocity and the mother bubble size, respectively. The bubble breakup region increases with the jet velocity. The bubble breakup frequency increases with the turbulent dissipation rate and the mother bubble size. The average number of daughter bubbles increases with the Weber number. An L-shaped daughter bubble size distribution is observed. Empirical correlations are established for the bubble breakup frequency, the average number of daughter bubbles and daughter bubble size distribution, and fitted well with the experimental results.     

SIMULATION OF BUBBLE BREAKUP DYNAMICS IN HOMOGENEOUS TURBULENCE

This article presents numerical simulation results for the deformation and breakup of bubbles in homogeneous turbulence under zero gravity conditions. The lattice Boltzmann method was used in the simulations. Homogeneous turbulence was generated by a random stirring force that acted on the fluid in a three-dimensional periodic box. The grid size was sufficiently small that the smallest scales of motion could be simulated for the underlying bubble-free flow. The minimum Weber number for bubble breakup was found to be about 3. Bubble breakup was stochastic, and the average time needed for breakup was much larger for small Weber numbers than for larger Weber numbers. For small Weber numbers, breakup was preceded by a long period of oscillatory behavior during which the largest linear dimension of the bubble gradually increased. For all Weber numbers, breakup was preceded by a sudden increase in the largest linear dimension of the bubble. When the Weber number exceeded the minimum value, the average surface area increased by as much as 80%.     

Coalescence efficiency of bubbles in bubble columns

Bubble size distribution was modelled by employing the population balance equation (PBE). All three bubble coalescence mechanisms (turbulence, buoyancy and laminar shear) and the main bubble breakup mechanism (breakup due to turbulent eddies) were considered in the model. Local bubble size distributions at the top and bottom of the column were obtained by solving this PBE. The results were compared with the experimental data for seven independent multiphase systems (water/air, isomax diesel/air, kerosene/air and four other liquid mixture/air) at two diverse gas velocities. The experimental adjustable constant in the coalescence efficiency function was determined by fitting the population balance to the experimental bubble size distributions. An empirical correlation was proposed for the coalescence efficiency by the dimensional analysis, which includes Reynolds and Weber numbers. © 2011 Canadian Society for Chemical Engineering     

微通道限域空间内的气泡破裂研究进展与展望

多相微反应器等微通道设备具有高效、安全等优势与广阔的应用前景,其中气泡分散相的形变、聚并、破裂等诸多流体力学行为对反应体系具有非常重要的影响,然而由于微通道的尺度特征以及多相流动非均匀性、复杂性等特点,复杂限域结构内的气泡形变与破裂机理认识仍不够充分。本文针对近年来在微尺度限域结构中的气泡形变与破裂等研究进展进行综述,首先概述了微通道多相流主要研究对象及研究方法,探讨了含有颗粒等受限狭窄空间及复杂多相条件下的传递机理研究进展,总结了气泡界面演变及不稳定破裂过程的多相相间作用机制。最后,提出研究体系构建、研究方法改进、相间作用解析以及工程化需求匹配是微尺度复杂限域空间内的气泡行为研究关键,并对下一步研究方向进行了分析与展望。     

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.     

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     

Multicompartment hydrodynamic model for slurry bubble columns

A core-annulus multicompartment two-dimensional two-bubble class model accounting for slurry recirculation and coupled with catalyst transport was developed as a part and parcel of the analysis of the behavior of slurry bubble column reactors at high gas throughputs corresponding to the churn turbulent flow regime. The model analyzed the contributions of bubble-induced turbulence closures, bubble coalescence and breakup phenomena, and catalyst axial distribution as the resultant of sedimentation, advection via liquid-solid slip, per-compartment axial dispersion and core-annulus lateral exchange of catalyst by bubble-induced turbulence. The model was also used to analyze the effects of catalyst loading, gas density and superficial velocity, and column diameter and vessel aspect ratio on the hydrodynamics of slurry bubble column reactors, namely, the per-compartment phase holdups and interstitial velocities, pressure gradient, bubble coalescence and break-up rates, and loci of velocity inversion for the gas and slurry profiles.