纳米流体强化气液传质研究进展

简要阐述了气相沉积法、共混法和分散法3种纳米流体的制备方法以及物理法、化学法两种纳米流体的分散技术。重点综述了纳米流体强化气液传质过程以及强化机理方面的最新研究成果。分析给出了纳米流体强化气液传质的几点原因：掠过效应、抑制气泡聚并机理、边界层混合机理、渗透机理以及多个影响因素相互关联作用。并预测出可能成为研究热点并有助于统一的强化理论表述提出的4个研究方向：影响纳米颗粒强化气液传质的各种因素相互耦合作用;纳米流体对氨、CO2、CO、O2和水蒸气以外的其它气体的物理化学吸收,进一步提出纳米颗粒强化气液传质普适性模型;搜集纳米颗粒影响气液界面附近的浓度分布;速度分布的微观信息以及纳米颗粒与气液界面相互作用的研究。     

加入第二液相强化气液传质

以水吸收二氧化碳为研究体系,考察了第二液相(有机相)的加入对气液传质的增强作用.运用经验关联式,对气泡平均直径进行了预测,应用渗透模型计算并讨论了第二液相的加入对气液接触界面积、体积传质系数的影响.实验结果表明:不同搅拌速度对气泡大小影响不同,适宜的搅拌速度有利于传质;少量的第二液相的加入减少了气泡聚并现象,增大了气液接触界面积.气液传质增强效果用增强因子表示.     

微通道内气-液弹状流动及传质特性研究进展

气-液弹状流,又称Taylor流,是一种以长气泡和液弹交替形式流动的流动形态。微通道内气-液弹状流因其气泡与液弹尺寸分布均一、停留时间分布窄、径向混合强等优点,是一种适于强化气-液反应的理想流型。本文首先介绍了微通道内气泡的生成机理、气泡和液弹长度,以及气泡生成阶段的传质特征。其次系统综述了主通道中弹状流动及传质过程的研究进展,包括气泡形状与液膜厚度、液弹内循环和泄漏流特征、气-液传质系数的测量与预测,以及物理与化学吸收过程中的传质特性等方面内容。最后阐述了当前研究的不足并展望了气-液弹状流的研究方向。     

吸收过程的界面传质机理

以分子热力学为基础 ,对气体吸收过程进行了理论分析 ,导出了传质通量的数学表达式。根据该文分析 ,气液界面传质的源动力来自界面处气液两相的不平衡 ,即只要有传质发生 ,液相界面处的浓度就不会达到与气相呈平衡的浓度。对于气相阻力可以忽略的吸收过程 ,两相传质速率的大小主要取决于液相的溶质界面浓度和液膜厚度 ,影响溶质界面浓度和液膜厚度的主要因素是近界面液相侧的流场分布。利用近界面浓度与液膜厚度的激光测定结果 ,计算了甲醇、乙醇及正丙醇吸收CO2 气泡的传质通量 ,计算结果与实验值吻合良好     

鼓泡破泡一体化高效精馏塔盘流动特性与CFD模拟

就痕量精馏中塔板传质效率低、需强化气液传质的问题,研究者提出了新型鼓泡破泡一体化高效精馏塔盘,通过在筛板上泡沫层高度范围内设置一层破泡装置,打破大气泡,减小气泡体积,强制界面进行更新,从而提高传质效率。采用双欧拉模型分别对鼓泡破泡一体化塔盘和筛板进行了气液流场的数值模拟,并对模型进行了验证。对比两种塔板的计算结果可以看出：在相同操作条件下,破泡装置将大气泡破裂成无数小气泡,使高气含率区域面积较普通筛板进一步增大,且气含率梯度变化更均匀;增加破泡装置后,在相同气速条件下气泡上升速度下降,气体在液层中的滞留时间延长,使鼓泡层高度增加,可显著提高传质效率,且降低了气体雾沫夹带量;破泡装置还明显改善了气相的纵向分布,气含率由塔板底部向上逐渐增大且存在明显分界;破泡装置附近湍动较剧烈,气泡破碎喷出的气体会进一步撕裂液膜,气体破碎作用会抑制气泡聚并,促进界面的快速更新更有利于传质过程的进行。研究结果可对工业塔板设计和优化提供指导。     

Lightnin静态混合器内气液两相混合与传质强化特性

采用计算流体力学(CFD)耦合群体平衡模型(PBM)研究Lightnin静态混合器(LSM)内气泡分散特性，使用不同的破碎核和聚并核函数系统研究雷诺数(Re)、气相体积分数(α_d)和元件数量对气泡分散行为以及混合效率的影响，采用气液界面积和体积传质系数(k_La)量化LSM和Kenics静态混合器(KSM)内气泡破碎性能以及传质速率，基于变异系数(CoV)和流体微元拉伸率分析LSM和KSM的分布混合性能和分散混合性能。结果表明：Luo聚并模型和Prince模型高估了LSM内气泡的聚并效率，通过Luo-Turbulent模型计算的气泡尺寸与实验数值具有很好的一致性；随着Re和α_d的增大，LSM内气泡的分散行为被进一步强化，在高Re的条件下增大元件数量可以显著提高气液传质效率；CoV曲线表明，LSM具有比KSM更好的分布混合性能，3个LSM元件可以保证混合程度大于95%。LSM的分散混合效率是KSM的1.06～1.16倍；基于压力波动信号时间序列，确定了LSM中流型从过渡区向非均匀区转变的临界表观速度。     

微通道内气液两相流行为研究进展

综述了微通道内的气液两相流行为及传质特性。在微通道内流型一般分为泡状流、弹状流、环状流和弹状-环状流,没有分层流。气液传质效率比常规尺度中的提高了2～3个数量级。讨论了气泡对气液两相流的影响及其生成、生长和聚并规律。介绍了微通道内气液两相流的计算机模拟结果。从实验、理论和数值模拟3个方面对微通道内气液两相流的研究和应用前景进行了展望。     

基于气泡数密度模型的气体穿越液池过程气泡特性数值模拟

采用考虑了气泡破碎和聚并的平均气泡数密度(ABND)输运方程,并与计算流体力学(CFD)中的湍流双流体模型相结合,建立CFD-ABND耦合计算模型。运用该模型对含固体颗粒气体通过冷却管穿越液池过程中液池内的气泡尺寸和气液界面浓度进行数值模拟。定量获得了液池内气泡尺寸和气液界面浓度分布,并分析了气速变化对其影响规律。结果表明:所建立的CFD-ABND模型能够对液池中的气泡尺寸及气液界面浓度分布等气泡特性进行较好预测;在靠近冷却管外壁面附近区域形成较大尺寸气泡和较高气液界面浓度;液池内隔板的存在有助于气液扰动,使液池内总体的气泡尺寸得到有效降低及气液界面浓度得到提升。     

ILAR中醇及电解质对非牛顿流体中气-液传质的影响

在-气升式内环流反应器中试验考察了非牛顿流体羧甲基纤维素钠(CMC)中的气泡聚并现象以及表面活性物质对液相体积传质系数的影响。结果表明,非牛顿流体中气-液传质效率随黏度的增加而降低,其原因是黏度增加使Taylor泡的尾流趋于稳定,降低了液相扰动,气泡间易聚并,从而气-液传质效率低。向非牛顿流体中添加醇类物质会影响气-液传质行为,对于聚合物含量低的流体,添加微量醇可以促进气-液传质,聚合物含量高的非牛顿流体,微量醇的加入反而不利于气-液间传质过程。非牛顿流体在ILAR上升管中的气含率随着黏度的增加变化不大,而下降管中的气含率有所提高。     

鼓泡破泡一体化高效精馏塔盘流动特性与CFD模拟

就痕量精馏中塔板传质效率低、需强化气液传质的问题,研究者提出了新型鼓泡破泡一体化高效精馏塔盘,通过在筛板上泡沫层高度范围内设置一层破泡装置,打破大气泡,减小气泡体积,强制界面进行更新,从而提高传质效率。采用双欧拉模型分别对鼓泡破泡一体化塔盘和筛板进行了气液流场的数值模拟,并对模型进行了验证。对比两种塔板的计算结果可以看出:在相同操作条件下,破泡装置将大气泡破裂成无数小气泡,使高气含率区域面积较普通筛板进一步增大,且气含率梯度变化更均匀;增加破泡装置后,在相同气速条件下气泡上升速度下降,气体在液层中的滞留时间延长,使鼓泡层高度增加,可显著提高传质效率,且降低了气体雾沫夹带量;破泡装置还明显改善了气相的纵向分布,气含率由塔板底部向上逐渐增大且存在明显分界;破泡装置附近湍动较剧烈,气泡破碎喷出的气体会进一步撕裂液膜,气体破碎作用会抑制气泡聚并,促进界面的快速更新更有利于传质过程的进行。研究结果可对工业塔板设计和优化提供指导。     

幂律流体中平行上升双气泡间相互作用研究

采用流体体积(VOF)法对非牛顿幂律流体中平行上升双气泡间相互作用进行了模拟研究.结果表明:剪切变稀效应及气泡周围流场结构对气泡间相互作用的贡献分别随着气泡间距增大而减小,气泡间相互作用主要受气泡间流体涡旋结构控制.当初始间距很小时,气泡平行聚并发生;随着初始间距增大,双气泡行为转变为相互影响的平行上升;当初始间距增大...     

Interfacial area density in bubbly flow

The interfacial area per unit volume is one of the key parameters in bubbly flow. Momentum, mass and energy transfer occur through the interface between the phases. The functionality of two phase reactors with bubbly flow depends mainly on these three transfer processes. Thus, the design process of a reactor requires the prediction of interfacial area density. In the present work a simple equation for the interfacial area density is derived from the population balance, taking into account the events of coalescence and bubble break-up for each bubble fraction. The system of partial integro-differential equations is simplified. Since the integrals in these equations complicate a numerical treatment. This reduces the balance to one single partial differential equation. An approximate analytical solution is given. If the resulting equation is applied to large gas fluxes, the instability of the coalescence process causes large bubbles and gas plugs to develop. From the instability the volume fraction of the large bubbles and gas plugs may be predicted. Additives may hinder the coalescence process. Experiments show that coalescence hindrance changes the coalescence kernel only by a factor. Calculations are done for bubble columns and vertical pipe flow.     

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

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

Physical explanation of the empirical coefficients of gas-liquid mass transfer equations

Equipment design and scale up is one of the biggest problems that chemical engineers face. A lot of research has been done at laboratory and pilot plant scale. However, experimental data are, many times, useless in equipment scale up because they overlook the processes involving bubbles. Therefore, a new approach considering the hydrodynamics of the bubbles is proposed to explain and understand the mass transfer rates. Semi-theoretical equations have been developed for bubble columns and stirred tanks based on hydrodynamic considerations of the flow, the bubbles and dispersions in the tank. These equations are able to explain the physical meaning, identify the effect of the scale on mass transfer rates and even predict the coefficients of the empirical correlations for k_La based on the hydrodynamic processes that bubbles experience. It can be concluded that the proportional constant of the correlations for k_La depends on the bubble size, on the physical and transport properties and on the size of the tank because it affects the mixing. Meanwhile, the exponents related to the power input and the superficial gas velocity depend on bubble break up and coalescence, and the dispersions generated. In the case of stirred tanks, the physical effect of the impeller on the bubbles also plays an important role on the exponents.     

Three‐dimensional numerical simulation of coalescence and interactions of multiple horizontal bubbles rising in shear‐thinning fluids

The dynamics of multiple horizontal bubbles rising from different orifice arrangements in shear‐thinning fluids was simulated numerically by three‐dimensional Volume of Fluid method. The effects of bubble size, rheological properties of shear‐thinning fluids, and orifice structure arrangements on multiple bubbles interaction and coalescence were analyzed, and the mechanisms of bubble coalescence and breakup were fully discussed and elucidated. The variation of bubble rising velocity during coalescence process and freely rising processes for different orifice arrangements was also deeply investigated. The critical initial horizontal intervals for coalescence of multiple horizontal bubbles with various orifice arrangements were attained by simulation, which could serve as the critical criterion of bubble coalescence or noncoalescence. Furthermore, the critical bubble interval was predicted based on the film drainage model, the prediction accords well with the simulation result and is quite conducive for the design and optimization of perforated gas–liquid contact equipment. © 2015 American Institute of Chemical Engineers AIChE J, 61: 3528–3546, 2015     

An improved subgrid scale model for front-tracking based simulations of mass transfer from bubbles

Gas–liquid bubble column reactors are often used in industry because of their favorable mass transfer characteristics. The bubble mass boundary layer in these systems is generally one order of magnitude thinner than the momentum boundary. To resolve it in simulations, a subgrid scale model will account for the sharp concentration variation in the vicinity of the interface. In this work, the subgrid scale model of Aboulhasanzadeh et al., Chem Eng Sci, 2012, 75:456–467 embedded in our in-house front tracking framework, has been improved to prevent numerical mass transfer due to remeshing operations. Furthermore, two different approximations of the mass distribution in the boundary layer have been tested. The local and global predicted Sherwood number has been verified for mass transfer from bubbles in the creeping and potential flow regimes. In addition, the correct Sherwood number has been predicted for free rising bubbles at several Eötvös and Morton numbers with industrial relevant Schmidt numbers (10³–10⁵).     

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.     

高气液密度比的传热相变复合模型

为研究两相流动中热量传递机制,基于格子Boltzmann热模型及大密度比模型,将相变源项引入到控制两相密度分布函数中,来描述温度场对气液相变的影响,提出了一个可以描述气液密度比达到2778的传热相变复合模型。通过对压力速度分布函数的回归修正克服了气液密度比过大造成的数值不稳定问题。模拟了溴化锂水溶液中双气泡的上升运动过程及周围的温度场分布,研究发现：双气泡上升时,碰撞前上方气泡温度高于下方气泡,碰撞时,两气泡间液桥打开,发生热量传递,气泡内部温度变得均匀;双气泡体积先减小再增大,碰撞时体积达到最大值,在融合成一个气泡后体积逐渐缩小,最终趋于稳定;初始气泡的体积越大,气泡上升过程中的速度越大。