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

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

采用电阻层析成像技术测量气液加压鼓泡塔局部气泡参数

在内径0.3 m,高6.6 m的加压气液鼓泡塔反应器中,采用电阻层析成像技术(ERT)研究了空气-水体系中气泡群平均上升速度、局部气含率及其径向分布。在表观气速0.119~0.312 m s 1,压力0.5~2.0 MPa,考察了表观气速、压力对气泡群上升速度、局部气含率及其径向分布的影响。实验结果表明,鼓泡塔中局部气含率随着表观气速与压力的增大而增大,其径向分布呈现出中心高边壁低的分布特征,但整个截面的分布并非严格对称,在r/R=0~0.3,气含率变化较小,且极大值出现在该范围内;气泡群的局部上升速度随着表观气速的增大而增大,但是随着压力的增大而减小。     

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

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

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

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

鼓泡塔中离子液体-空气两相流的CFD-PBM耦合模拟

针对咪唑类离子液体介质,采用Euler-Euler双流体模型与群平衡模型（PBM）耦合的方法,引入由实验结果拟合获得的适用于该介质的气液相间曳力系数模型,对内径0.203 m、高2 m的鼓泡塔中离子液体-空气两相流进行计算流体力学模拟,研究了不同表观气速下塔内气液两相速度场分布、气含率和气泡尺寸分布等流体动力学性质。与现有的相间曳力系数模型Schiller-Naumann模型模拟结果对比,采用本文模型得到的气含率与实验值吻合更好,气泡尺寸分布与实验结果一致。     

逆流鼓泡塔气含率分布的试验测量与CFD模拟

利用热膜测速仪测得了气液逆流鼓泡塔内不同表观气速、表观液速和径向位置下的气液信号,采用改进的阈值法进行分析,得到塔内气含率的径向分布。结果表明气含率在各个截面上都是从塔中心到塔壁逐渐减小;同时利用计算流体力学方法对气液逆流的鼓泡塔内的气液两相流动进行了模拟,计算了不同气速和不同液速下的气含率,计算结果与试验数据吻合较好。     

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

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

基于径向力平衡的鼓泡塔二维流体力学模型

提出了一种鼓泡塔二维轴对称流体力学模型,模型中将气泡所受的升力以及湍动扩散力作为形成塔内气含率稳定分布的主要机制.采用Fluent 6.3流体力学软件求解模型,能得到稳定的二维流场,气含率与液速分布与实验值吻合良好,模型能准确反映表观气速(0.12～0.62 m·s^-1)以及塔径(ø200 mm、ø500 mm、ø800 mm)对流型的影响.利用该模型对更大直径鼓泡塔的流动参数进行了预测,结果与文献给出的经验关联式相符.     

超声场中鼓泡塔内气泡直径分布特征研究

利用自主研发的超声鼓泡塔和双电导探针气泡特征参数测量仪测定了不同频率组合下超声场辐射对鼓泡塔内气泡Sauter直径随表观气速和鼓泡塔径向及轴向位置的分布情况。实验结果表明,在距气体分布板轴向距离250 mm、塔中心处气泡Sauter直径由于超声作用而显著减小,最大减幅达55%;在超声场作用下随着表观气速增大,气泡Sauter直径呈增大趋势,当表观气速大于0.06 m/s时,气泡Sauter直径趋于不变;大气泡在塔中心区域聚集倾向明显,小气泡分布在塔壁面附近,气泡Sauter直径(在径向位置)呈现中心峰分布趋势;随着距气体分布板轴向距离增加,气泡Sauter直径增大,在多频混响声场中,气泡Sauter直径增速小于单频声场中增速,即在多频混响超声场中易形成小气泡。实验结果证实,超声场有利于减小气泡体积,增大气泡比表面积,这对于加强气液相间传质、提高传质效果会有显著作用。     

高气速下鼓泡塔中气含率分布的测定

在内径476 mm的鼓泡塔内用压差法测定全塔平均气含率与表观气速的关系,进一步利用响应特性良好的双电导探针,考察了不同气速下局部气含率的分布规律。实验结果表明,利用探针法计算得到的全塔平均气含率值与压差法测定值平均误差仅为4.5%,表明探针法测量局部气含率的可靠性良好。实验还表明在高气速下,除分布板影响区外,局部气含率均类似抛物线型分布;随着气速增加,气含率分布趋于陡峭。以实验为依据,拟合了不同气速下(0.05~1.0 m/s)鼓泡塔中局部气含率的关联式,认为塔内局部气含率与径向位置、表观气速和塔径等因素有关。     

基于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%。     

含垂直列管束内构件鼓泡塔的CFD模拟

针对含有密集垂直列管束内构件的气液鼓泡塔,在两相Euler二维轴对称k-ε模型中,分别考虑气相和液相受到的列管阻力。通过引入相应的动量源、湍动源以及耗散源建立带列管内构件的鼓泡塔二维CFD模型。模型能清晰、准确地描述带列管束鼓泡塔中气液流动的特征:"烟囱效应"以及分布器影响区延长。计算得到的气含率以及液速的二维分布在宽泛的表观气速(0.12～0.62 m·s^-1)范围内与实验值相符。     

Experimental and numerical investigations of scale-up effects on the hydrodynamics of slurry bubble columns

Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid concentrations(0–30 vol%). Radial profiles of time-averaged gas holdup, axial liquid velocity, and turbulent kinetic energy were measured by using in-house developed conductivity probes and Pavlov tubes. Effects of column diameter, superficial gas velocity, and solid concentration were investigated in a wide range of operating conditions. Experimental results indicated that the average gas holdup remarkably increases with superficial gas velocity, and the radial profiles of investigated flow properties become steeper at high superficial gas velocities. The axial liquid velocities significantly increase with the growth of the column size, whereas the gas holdup was slightly affected. The presence of solid in bubble columns would inhibit the breakage of bubbles, which results in an increase in bubble rise velocity and a decrease in gas holdup, but time-averaged axial liquid velocities remain almost the same as that of the hollow column. Furthermore, a 2-D axisymmetric k–ε model was used to simulate heterogeneous bubbly flow using commercial code FLUENT 6.2. The lateral lift force and the turbulent diffusion force were introduced for the determination of gas holdup profiles and the effects of solid concentration were considered as the variation of average bubble diameter in the model. Results predicted by the CFD simulation showed good agreement with experimental data.     

Hydrodynamics and mass transfer in a bubble column with an organic liquid

Hydrodynamic and mass transfer data in Soltrol-130 (a mixture of C₉+ iso-paraffins) were measured in a 0.305 m diameter bubble column. The gas holdup structure (i.e., the contributions to holdup from the small and large bubble fractions of the dispersion) for this hydrocarbon liquid in the churn turbulent flow regime was analyzed using the dynamic gas disengagement technique. The validity of the assumption of axially uniform gas holdup structure was checked. Literature correlations were found inadequate to explain the observed gas holdup and the volumetric mass transfer coefficients for Soltrol-130. The volumetric mass transfer coefficient per unit volume of large bubbles is shown to be independent of superficial gas velocity for the fully developed churn turbulent regime. The present hydrodynamic and mass transfer data in the churn turbulent regime should be useful in the design and scale-up of bubble columns used in organic process industries.     

湍动浆态床流体力学研究（Ⅰ）气含率及其分布规律

对高气速、高固含率、大塔径条件下的湍动浆态床平均气含率和气含率径向分布进行了实验测定,结合工业实验数据,归纳出可用于工业条件的气含率计算关联式,给出了简化的流体力学模型用于气含率分布的模拟。结果表明,浆态床气含率将随塔径增加而降低,固含率与塔径之间存在交互影响;同时,气含率的径向分布也随气速和塔径的增大而改变,存在明显的放大效应,简化模型能够较好地模拟实验结果。     

Flow regimes and mass transfer in counter-current bubble columns

Counter current bubble columns have the feature that specific gas-liquid interfacial area and gas holdup are larger than those for standard and cocurrent bubble columns. In this study, three different flow regimes, churn-turbulent flow, bubble flow and bubble down-flow, have been observed in a counter-current bubble column and correlations of gas holdup and volumetric liquid-phase mass transfer coefficient have been proposed as functions of operating variables such as the superficial velocities of gas and liquid, the gas-liquid slip velocity and the liquid properties.     

Hydrodynamics of air–kerosene bubble column under elevated pressure in homogeneous flow regime

A multiphase computational fluid dynamics(CFD) model coupled with the population balance equation(PBE) was developed in a homogeneous air–kerosene bubble column under elevated pressure(P). The specific pressure drop(DP/L), gas holdup(a_G), and Sauter mean diameter(d_(32)) were experimentally measured in the bubble column with 1.8 m height and 0.1 m inner diameter, which was operated at a superficial gas velocity of 12.3 mm·s~(-1), and P = 1–35 bar(1 bar = 10~5 Pa). A modified drag coefficient model was proposed to consider the effect of bubble swarm and pressure on hydrodynamics of the bubble column.The Luo breakage model was modified to account for liquid density, viscosity, surface tension and gas density. The DP/L, a_G, and d_(32) obtained from the CFD model were compared with experimental data,and the gas density-dependent parameters of the CFD model were identified. With increasing P from 1 to 35 bar, the aGvaried from 5.4% to 7.2% and the d_(32) decreased from 2.3 to 1.5 mm. The CFD-PBE model is applicable to predict hydrodynamics of pressurized bubble columns for gas–organic liquid in the homogeneous regime.     

Hydrodynamics of a Tapered Bubble Column

The hydrodynamic behavior of a single‐stage tapered bubble column using an air‐water two‐phase system has been studied. The experimental results indicate that the holdup increases with increasing superficial gas velocity and bubble slip velocity, while it remains constant with increasing superficial liquid velocity. The gas flow rate has a subtle effect on pressure drop owing to the dynamic pressure recovery stemming from the increase in flow area in the axial direction. The results further suggest that the tapered bubble column shows higher holdup with lower energy dissipation than the conventional bubble columns under similar hydrodynamic conditions. The experimental values of the holdup are in excellent agreement with the well‐known Akita and Yoshida correlation available in the existing literature. Also, the performance of the tapered system has been shown to be much better than that of conventional columns under similar conditions in water/alkaline scrubbing of fly ash and SO₂ either alone or in combination.     

SEBS高黏度溶液气液鼓泡塔的流体力学研究

针对SBS加氢反应器开发与设计,以SEBS-1650己烷溶液为液相,采用差压法和床层塌落法研究了气液鼓泡塔中高黏度溶液的流体力学行为,考察了黏度对低表面张力溶液的气含率、大小气泡气含率、大小气泡上升速度和比表面积等因素的影响。结果表明,随黏度增加,大气泡增多,气含率明显降低,塔内流型处于湍流区;由床层塌落曲线确定鼓泡塔内存在三种类型的气泡：大气泡、小气泡及细小气泡,随黏度增加,小气泡与细小气泡逐渐减少;黏度对大小气泡的上升速度略有影响,比表面积随黏度增加而明显降低。根据实验结果给出了大小气泡气含率与平均气含率的计算公式。