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

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

加压气液鼓泡塔的CFD数值模拟与ERT实验验证

在内径0.3 m,高6.6 m的加压鼓泡塔内,采用计算流体力学(CFD)数值模拟与气泡群平衡模型(PBM)耦合法进行塔内流体力学模拟,并将数值模拟结果与基于电阻层析成像技术(ERT)的实验结果对比分析,将通过ERT实时采集的横截面气含率分布和时间序列图与模拟结果进行比较。结果表明:ERT技术测量结果与CFD计算结果吻合良好,能很好地表示鼓泡塔内气液流动状态,进一步表明ERT技术对加压鼓泡塔内气液两相流进行可视化与实时测量是可行的。     

基于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-PBM耦合模型模拟气液鼓泡床的通用性研究

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

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

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

不同黏度下气液体系流体力学行为的PBM模拟

基于Eulef-Euler双流体模型,分别采用均一气泡尺寸和PBM方法对鼓泡床内气液两相流体系中的流体力学行为进行了数值模拟,经与实验结果比较表明,采用PBM方法获得的模拟结果更为准确,表明了PBM方法的优越性.采用PBM模型模拟了不同黏度气液体系中的流体力学行为,获得了包括平均及局部气含率分布、气泡尺寸分布,循环液速和大小气泡的空间分布等流体力学信息,并将模拟结果与相应操作条件下的实验结果进行了比较,吻合较好,说明PBM方法不仅可以用于鼓泡床反应器内气液两相流体系的模拟研究,还能对不同液相黏度下气液体系中的流体力学行为作出正确的预测和描述,为多相流反应器的设计和放大提供了有益的参考.     

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

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

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

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

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

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

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

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

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

Local hydrodynamics of gas–liquid-nanoparticles three-phase fluidization

The laser Doppler anemometer (LDA) and conductivity probes were used for measuring the local hydrodynamic performances such as gas holdup and liquid velocity in a lab-scale gas–liquid–TiO₂ nanoparticles three-phase bubble column. Effects of operating parameters on the local gas holdup and liquid velocity were investigated systematically. Experimental results showed that local averaged axial liquid velocity and local averaged gas holdup increased with increasing superficial gas velocity but decreased with increasing TiO₂ nanoparticles loading and the axial distance from the bottom of the bubble column. A three-dimensional computational fluid dynamic (CFD) model was developed in this paper to simulate the structure of gas–liquid–TiO₂ nanoparticles three-phase flow in the bubble column. The time-averaged and time-dependent predictions were compared with experimental data for model validation. A successful prediction of instantaneous local gas holdup, gas velocity, and liquid velocity were also presented.     

筛板结构对多级鼓泡塔中流体力学参数的影响

在外径0.3 m的多级鼓泡塔中采用双探针电导探头考察了筛板结构对鼓泡塔中筛板上下两处流体力学参数的影响,并与无筛板鼓泡塔中流体力学参数进行了对比。研究结果表明:带筛板的鼓泡塔中气含率径向分布规律和无筛板的鼓泡塔相类似,筛板上下二侧的平均气含率随表观气速的增加而增加;在鼓泡塔轴向位置上,筛板下侧的气含率普遍比筛板上侧气含率要高;筛板的孔分布方式对于气含率径向分布有着一定影响,但孔径大小对此影响较小。     

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

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

CFD‐Facilitated Flow Field Analysis of Bubble Columns with Concentric Plates for Biological Applications

As the hydrodynamic behavior in bubble columns is difficult to characterize, computational fluid dynamics (CFD) is a useful alternative tool for research and design. An experimental and computational analysis of the macromixing and gas holdup of a bubble column indicates that 3D simulations with CFD using an Eulerian‐Eulerian approximation yield results of the overall velocity field and gas holdup distribution that are suitable for engineering design purposes. Particularly, CFD simulations uncover that the inclusion of concentric solid plates into a bubbling column increase the gas holdup by 79 % and the mixing time by 48 % when compared with a column without plates operating at similar superficial gas velocities.     

Gas holdup in pulp fibre suspensions: Gas voidage profiles in a batch-operated sparged tower

Gas holdup in a semi-batch operated slurry (pulp fibre suspension) bubble column was investigated for two pulp types (softwood and hardwood kraft pulps) over a range of suspension mass concentrations (C_m=0–9% by mass) and superficial gas velocities (U_g=0.0027–0.027 m/s). Three techniques were used: height difference between gassed and ungassed operation; pressure difference as a function column height; and electrical resistance tomography (ERT). Depending on the technique used the average, axial and radial holdup profiles could be determined. In the pulp suspensions, the ERT determined gas holdups correlated well with those determined using the differential height method. In water, the ERT determined gas holdups were significantly lower, but the agreement was significantly improved by increasing the background conductivity by adding 1 g/L salt to the water. This, however, reduced the overall gas holdup due to the effect of the electrolyte on bubble coalescence. Other differences between the three measurement techniques were attributed to limitations in the detection methods and the averaging procedures used to compare results.The presence of pulp fibres reduced gas-holdup at all gas flow rates and suspension concentrations studied and is attributed to increased bubble coalescence which increases bubble size and consequently bubble rise velocity through the suspension. Gas holdup (as determined by ERT) increased with column height. The radial gas profiles were non-uniform and more peaked than the corresponding water profiles. At low suspension concentrations this was attributed to asymmetric suspension recirculation within the column. As suspension concentration increased, channels formed in the suspension with the average void fraction leveling off to a plateau.     

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.