双气泡相群平衡模型模拟鼓泡塔气液两相流

基于对鼓泡塔气泡行为的现有认识,把气泡分成大、小气泡. 首次建立了完整的双气泡相群平衡模型(TBPBM),以预测气泡尺寸,并耦合TBPBM与CFD双流体模型对D=440 mm鼓泡塔进行数值模拟,获得了气泡尺寸体积概率分布、时均气含率与液相速度径向分布、大小气泡相尺寸分布,对部分模拟结果与实验值及文献模拟结果进行了比较. 结果表明,TBPBM-CFD模型预测的时均气含率和液相速度分布与实验结果吻合最好,较SBPBM、平均气泡尺寸模型的模拟结果有明显改善. 与实验值相比,TBPBM模型的整体气含率模拟误差为5.7%,而SBPBM模型和平均气泡尺寸模型的误差分别为27.2%和17.3%.     

加温加压下CFD-PBM耦合模型空气-水两相流数值模拟研究

计算流体力学与群体平衡模型（CFD-PBM）结合可有效地模拟鼓泡塔内流体行为,较准确地预测流场、相含率以及局部气泡尺寸分布。以直径100 mm、高1.3 m的加温加压鼓泡塔为模拟对象,在系统压力为1 MPa、表观气速为0.08~0.24 m/s、温度为30~160℃条件下系统地考察了空气-水体系的表观气速、温度以及固含率对平均气含率、大小气泡气含率、气泡直径和气泡尺寸分布等参数的影响。结果表明,平均气含率的模拟结果和实验值在10%的误差范围内吻合较好;温度的变化主要影响了塔内气泡的聚并和破碎,并用聚并破碎的机理解释了温度对其流体行为的影响。     

不同分布器对鼓泡塔气液两相流影响的CFD模拟

对不同分布器的鼓泡塔反应器内的气液两相流体动力学行为进行了三维瞬态数值模拟.模型采用欧拉-欧拉双流体模型以及RNG κ-ε湍流模型,研究了不同分布器设计对鼓泡塔反应器气液两相流的影响.模拟得到了鼓泡塔内整体气含率、液相速度矢量分布、时均气含率以及时均轴向液速的径向分布等结果,并对部分模拟结果与文献实验结果进行了比较,其结果吻合得较好,证实了随着分布器开孔越均匀则有利于加强流形的非对称性,从而加剧了气液两相径向混合,但对整体气含率的影响不大.     

鼓泡塔内空气-醋酸体系流体力学参数的CFD-PBM耦合模型数值模拟

通过二维和三维CFD-PBM耦合模型对空气-醋酸体系中流体力学参数进行数值模拟,采用表面张力修正曳力模型与聚并模型,考察了醋酸浓度对鼓泡塔内气含率、气泡大小分布及轴向液速等参数的影响,与差压法、光纤探针和电阻层析成像技术（ERT）测量的实验数据进行了对比,并讨论分析了气含率和气泡直径等流体力学参数的模拟结果。结果表明,醋酸浓度在70%~80%（质量分数）范围内平均气含率存在最大值,且平均气含率的预测值在±10%误差内,三维模拟结果和ERT实验值吻合较好,说明修正后的模型在不同浓度醋酸体系中具有较好的预测性。     

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

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

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

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

分布器对鼓泡塔中气泡直径分布的影响

采用电导探针测定了冷态鼓泡塔中不同气速下的气泡直径及气含率的轴向分布,考察了分布板对鼓泡塔操作性能的影响.结果表明:随着开孔率的减小,从均匀鼓泡区到过渡区的转变提前;在均匀鼓泡区,开孔率对气泡直径影响较小;在过渡区,开孔率大的分布器形成的稳定气泡直径较小、气含率较大;分布板开孔直径越大,形成的初始气泡直径越大,但对轴向气泡直径分布的影响仅限于分布器区.包含分布器影响的气泡直径经验关联式为d/D=140.2Bo-0.5Ga-0.12Fr0.099(h/D)-0.15T-0.34(0.5 cm/s＜ug＜7 cm/s).     

静电效应对有无埋管气固鼓泡床内气泡特性的影响分析

采用基于双流体模型(TFM)耦合静电模型的方法,研究颗粒的静电对有无埋管气固鼓泡床内气固流动特性和气泡特性的影响。首先在无静电场存在的条件下,利用双流体模型对自由鼓泡床和埋管鼓泡床内的流动情况进行模拟并与实验结果进行对比;进一步耦合静电模型,考察静电对自由鼓泡床和埋管鼓泡床内床层的整体性质和气泡特性的影响。研究结果表明,在无静电场条件下采用双流体模型能较好地预测自由鼓泡床和埋管鼓泡床内的气固流动状况以及气泡的平均直径和气泡的上升速度。埋管的存在使鼓泡床内气固流动发生强烈扰动,并使气泡的平均直径和气泡的上升速度均呈振荡分布。静电的存在对自由鼓泡床和埋管鼓泡床内床层的平均固含率影响不大,但对气泡分布规律影响较大,使得自由鼓泡床内气泡数目减少,而埋管鼓泡床下部区域的气泡分布比较集中,上部有大气泡出现。     

散射管气体鼓泡塔流场分布特性及优化

鼓泡塔因其较好的气-液传质性能具有高污染物脱除效率,被广泛应用于生物化工和烟气处理等领域。鼓泡塔散射管气体分布器的几何尺寸和结构是影响相间传质效率的关键因素,优化塔内流场对于提高鼓泡塔内气-液两相间的传质效率至关重要。采用Fluent软件对有内构件散射管横向进气口式的鼓泡塔进行模拟研究,基于双流体方法和群体平衡模型(PBM)模型对鼓泡塔三维建模,采用一阶迎风差分格式离散,使用Phase Coupled Simple算法进行压力速度耦合。研究了散射管所在圆环直径d分别为0. 375D、0. 5D、0. 625D、0. 75D时(D为鼓泡塔直径),散射管进气口的布置对整体和局部气含率、液速和气泡尺寸等的影响。研究结果表明,随着散射管分布环直径的增大,整体气含率先增大后减少,平均气泡直径先减小后增大;当散射管所在圆环直径d=0. 5D时,鼓泡塔整体气含率和液相循环速度最大,平均气泡直径最小,鼓泡塔流场综合性能最好。     

鼓泡塔气液两相流不同曳力模型的数值模拟

采用欧拉-欧拉双流体模型对圆柱形鼓泡塔内气液两相流动进行了三维数值模拟. 通过UDF自定义程序对气相出口边界进行了速度修正,解决了模拟中自由区域内有漩涡的问题;分别使用单一气泡尺寸模型和群体平衡模型(PBM)计算气泡尺寸,并比较其对气含率分布的预测结果,分别采用Schiller-Naumann, Grace和Tomiyama曳力系数模型进行模拟. 结果表明,在全塔径均匀进气的简化条件下,单一气泡尺寸模型不适用,在合适的Hamaker数下,PBM模型中原用于颗粒计算的Abrahamson模型可计算气泡聚并速率;Tomiyama曳力模型耦合PBM模型可更好地描述塔内流动情况,并与文献值吻合良好. Schiller-Naumann模型所得平均气含率与实验值相差约40%,而Grace模型所得湍动耗散比Tomiyama的结果高14.5%,气含率分布与文献值相差16.3%.     

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

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

装配水平筛板的气升式反应器中气液传质特性的数值模拟

利用Turbulent–Lehr组合模型对装配水平筛板的气升式反应器进行了计算流体力学(CFD)模拟,研究水平筛板对气含率、气泡直径、体积传质系数(kLa)和气液流速的影响。结果表明,筛板对气相的囤积作用和对液相的阻碍作用增加了反应器的整体气含率。筛板对气相的二次均布作用减弱了筛板和液面之间区域的气泡聚并过程,筛板筛孔对气泡的破碎作用产生了大量小于初始直径的气泡,增加了气泡比表面积(a);筛板对液相的阻碍作用提高了筛板附近的气–液相流动速度差,从而提高了该区域的液膜传质系数(kL),强化了反应器内的气液传质效果。     

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.     

基于EMMS方法的鼓泡塔反应器CFD及群平衡模拟

能量最小多尺度（energy-minimization multi-scale,EMMS）方法已经被应用于气液体系中群平衡（population balance model,PBM）模型的改进。EMMS模型可计算气泡破碎聚并过程的能量,进而获得聚并速率的修正因子。应用这一模型对高气速鼓泡塔进行了模拟计算,并进一步对比了均一尺径模型、CFD-PBM模型以及CFD-PBM-EMMS模型的模拟结果与实验数据。结果表明,在高表观气速条件下,基于EMMS方法的群平衡模型可以更加准确地预测鼓泡塔中不同高度的气泡尺径分布和轴向液速,同时提高了对整体气含率和局部气含率的模拟准确性。在表观气速为0.16 m·s^-1和0.25 m·s^-1时,CFD-PBM-EMMS模型对气泡尺径分布的预测精度更高,同时整体气含率模拟的相对误差下降为5%和15%,局部气含率模拟平均相对误差下降为8%和17%。     

Numerical simulation of bubble columns flows: effect of different breakup and coalescence closures

Two-dimensional axisymmetric Eulerian/Eulerian simulations of two-phase (gas/liquid) transient flow were performed using a multiphase flow algorithm based on the finite-volume method. These numerical simulations cover laboratory scale bubble columns of different diameters, operated over a range of superficial gas velocities ranging from the bubbly to the churn turbulent regime. The bubble population balance equation (BPBE) is implemented in the two-fluid model that accounts for the drag force and employs the modified k-ε turbulence model in the liquid phase. Several available bubble breakup and coalescence closures are tested. Quantitative agreements between the experimental data and simulations are obtained for the time-averaged axial liquid velocity profiles, as well as for the kinetic energy profiles, only when model predicted breakup rate is increased by a factor of ten to match the coalescence rate. The calculated time-averaged gas holdup profiles deviate in shape from the measured ones and suggest that full three-dimensional simulation is needed. Implementation of BPBE leads to better agreement with data, especially in the churn-turbulent flow regime, compared to the simulation based on an estimated constant mean bubble diameter. Differences in the predicted interfacial area density, with and without BPBE implementation, are significant. The choice of bubble breakup and coalescence closure does not have a significant impact on the simulated results as long as the magnitude of breakup is increased tenfold.     

Numerical simulations of gas-liquid mass transfer in bubble columns with a CFD-PBM coupled model

Gas-liquid mass transfer in a bubble column in both the homogeneous and heterogeneous flow regimes was studied by numerical simulations with a CFD-PBM (computation fluid dynamics-population balance model) coupled model and a gas-liquid mass transfer model. In the CFD-PBM coupled model, the gas-liquid interfacial area a is calculated from the gas holdup and bubble size distribution. In this work, multiple mechanisms for bubble coalescence, including coalescence due to turbulent eddies, different bubble rise velocities and bubble wake entrainment, and for bubble breakup due to eddy collision and instability of large bubbles were considered. Previous studies show that these considerations are crucial for proper predictions of both the homogenous and the heterogeneous flow regimes. Many parameters may affect the mass transfer coefficient, including the bubble size distribution, bubble slip velocity, turbulent energy dissipation rate and bubble coalescence and breakup. These complex factors were quantitatively counted in the CFD-PBM coupled model. For the mass transfer coefficient k_l, two typical models were compared, namely the eddy cell model in which k_l depends on the turbulent energy dissipation rate, and the slip penetration model in which k_l depends on the bubble size and bubble slip velocity. Reasonable predictions of k_la were obtained with both models in a wide range of superficial gas velocity, with only a slight modification of the model constants. The simulation results show that CFD-PBM coupled model is an efficient method for predicting the hydrodynamics, bubble size distribution, interfacial area and gas-liquid mass transfer rate in a bubble column.     

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.     

基于气泡群相间作用力模型的加压鼓泡塔流体力学模拟

目前,多数文献报道了冷态加压湍动鼓泡塔内流动特征,并且通过实验数据回归相关经验关联式。然而,此类关联式适用范围有限,难以直接外推到工业鼓泡塔反应器条件。因此,在FLUENT平台上建立了基于气泡群相间作用力的、动态二维加压鼓泡塔计算流体力学模型。通过数值模拟考察了操作压力为0.5~2.0 MPa,表观气速为0.20~0.31 m·s^-1,内径0.3 m鼓泡塔内流场特性参数分布,并且与冷态实验数据进行比较。结果表明,采用修正后的气泡群曳力模型、径向力平衡模型以及壁面润滑力模型描述气泡群相间作用力,能够较为准确地反映平均气含率和气含率径向分布随操作压力和表观气速变化的规律。     

单孔及微孔曝气低气速鼓泡床内气泡行为比较

引言鼓泡床反应器被广泛应用于吸收、液相氧化、好氧生化等气液反应过程,气体在液相中的分散情况对鼓泡床的反应和传质特性都有很大影响.为了提高气液传质效率,增加生产强度,工业反应器很多都是在高气速下操作(Ug>0·05m·s-1),很多研究都集中在高气速湍动鼓泡区[1~3].但对有机