竖直上升气液泡状流数学模型封闭的研究

竖直上升管气液两相流广泛应用于相变传热、核反应堆等工业过程。本文以竖直上升气液两相流为研究对象,运用欧拉双流体模型,针对表观液速为0.45m/s、表观气速分别为0.015m/s和0.1m/s的泡状流数值模拟过程中的升力、壁面润滑力、湍流扩散力、气泡诱导湍流（BIT）等封闭模型,开展数值模拟比较研究。模拟发现：①低气速泡状流中,升力和壁面润滑力的同时加入能够改善壁面附近的气含率,气泡在这两个力作用下在径向上达到一个相对平衡,得到与实验气含率类似的壁面峰,模拟的液相速度较合理;低气速时,BIT的影响可以忽略。②高气速泡状流中,BIT对气-液两相流的模拟结果影响比较明显,湍动耗散源项的加入能使液速分布的模拟结果得到改善,Troshko模型相对Sato模型更能反映气泡诱导湍流对液相湍流的作用。③高气速时升力的引入使气含率产生壁面峰,加入湍流扩散力能使峰值略微降低,但仍没有解决高气速时引入升力出现的气含率壁面峰问题,说明在径向上湍流扩散力还不足以抵抗升力。     

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

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

大管径垂直管道内高压汽液混合流动的数值研究

采用流体体积模型(VOF)对高压环境下190 mm大管径垂直管内水蒸汽-水混合流动进行数值研究。数值计算得到了5.07,10.13与17.23 MPa高压下大管径垂直管内汽液流型分布图及搅混流态的相分布图和速度场分布,并与常压下的计算结果进行对比,以研究压力环境带来的影响。数值结果表明,高压环境下大管径垂直管内的流型图与Hewitt和Roberts流型图的吻合度较差。高压环境下大管径垂直管内没有出现雾状流;泡状流和搅混流的发生区域扩大;弹状流的发生区域被压缩得很小;环状流的变化最小。随着压力的增大,大管径垂直管内汽液搅混流中界面波的高度有所降低;液膜铺展在壁面的面积扩大;系统的稳定性提高。速度场分布是管道中心处速度较大;近壁面处速度场发生振荡;壁面处速度迅速减小至零。随着压力的增大,近壁面处速度场振荡的紊乱程度减轻。根据汽液两相流动特性,分析了发生以上现象的原因。     

水平管泡状流相分布特性

<正>气液两相泡状流的主要特征是连续液相中携带散布其中的细小气泡,气泡的存在不仅对气液两相流的传热、传质及阻力特性有很大的影响,而且对两相流动的稳定性也有很大的影响.前人有关泡状流的研究大多偏重于两相流的平均参数,对于两相流局部统计参数如局部空隙率等参数的变化规律是近年来两相流研究的新趋势.对于垂直管内的流动已经积累了相当数量的数据,而同样有广泛应用的水平管内相分布规律还知之甚少.气泡对气液两相流的传热、传质及流动结构影响机理的研究必须以了解相分布及气泡的局部统计参数为前提,同时对相分布特性的深人研究也为气液两相流的数学模型化提供实验依据.本文以空气、水为工质,研究水平管内气液两相流的相分布特性,给出了典型泡状流的时域信号图,研究了相分布随气液两相流量的规律变化,并与前人的有关结果进行了比较.     

湍流泡状流混合层流动的PIV测量

在单相实验的基础上,使用PIV对不同位置注入气泡的泡状流湍流混合层流动进行实验研究。混合层高低侧流速比为4∶1,基于两股流体速度差和管道水力直径的Reynolds数范围为4400～158400。与单相相比,气液两相混合层流动脉动速度的测量结果表明,在低Reynolds数时气泡的加入会增强混合层的速度脉动,但随着Reynolds数的增大气泡的加入反而会减弱速度脉动。雷诺应力集中在隔板下游一个较窄的区域内,并随Reynolds数的增大而增大。从流动横截面上雷诺应力的分布可以看出,气泡的加入减小了混合层内的平均雷诺应力,并且与单相相比泡状混合层流动在同一横截面上的雷诺应力分布曲线会产生波动。     

基于丝网探针的螺旋管内气液两相流气泡行为研究

核反应堆蒸汽发生器的传热面由螺旋管束组成。螺旋管的三维螺旋结构使得泡状流和塞状流等气液两相流中的气泡在重力、离心力和浮力等作用下在管道内部呈现不对称的相分布状态,两相滑移速度增大,显著影响换热性能并导致DNB型传热恶化难以预测。实验介质为空气-水,结合自主开发的电导式丝网探针技术并发展先进的数据后处理算法,实现了复杂流场的三维时空重构和离散气泡粒径的精细测量,获得了螺旋管内泡状流和塞状流的截面空泡分布规律。基于研究结果,可根据气泡分布规律对螺旋管道的几何结构进行调整以避免传热恶化,为螺旋管式蒸发器的安全设计提供了基础数据和优化思路。     

基于旋涡强度法的湍流泡状流混合层流动

用旋涡强度法对PIV测得的湍流泡状流混合层流动进行研究,混合层高低侧流速之比为4∶1,基于速度差和管道水力直径的Reynolds数范围为4400～132000。气泡分别从隔板尾部混合层的中心起始位置与低速侧注入流体,直径为0.5～2 mm。利用旋涡强度法分析流场清楚地显示出排除了剪切作用的涡结构主要集中在混合层中心的锥形区域,随着Reynolds数的增大,旋涡强度的值不断增大而锥形区域不断变窄。在同一Reynolds数下,沿主流方向旋涡强度先增大、后减小。从混合层中心注气时气泡主要分布在隔板下游旋涡强度值较大的区域,而低速侧注入的气泡主要分布在低速侧,但由于涡结构的卷吸作用,低速侧注入的气泡也会进入隔板下游的中心区域。与单相流动相比,低Reynolds数情况下气泡的加入可以加剧流场本身的变化,而高Reynolds数时则对旋涡强度产生的影响较弱。     

微通道内气-液两相传质过程行为及其应用

微通道内气-液两相体系中Taylor流和泡状流具有气泡尺寸均一、停留时间分布窄、可调控性强和比表面积高等优点,具有广泛的应用前景。从Taylor气泡和泡状气泡的传质过程出发,系统综述了微尺度下气泡的溶解规律、传质过程机理和传质/溶解模型等方面的研究进展,并介绍上述流型在反应或过程强化、基础物性及动力学数据测量和微纳材料合成方面的应用。最后总结并展望了技术领域的研究难点与研究方向。     

微通道内气-液两相传质过程行为及其应用

微通道内气-液两相体系中Taylor流和泡状流具有气泡尺寸均一、停留时间分布窄、可调控性强和比表面积高等优点,具有广泛的应用前景。从Taylor气泡和泡状气泡的传质过程出发,系统综述了微尺度下气泡的溶解规律、传质过程机理和传质/溶解模型等方面的研究进展,并介绍上述流型在反应或过程强化、基础物性及动力学数据测量和微纳材料合成方面的应用。最后总结并展望了技术领域的研究难点与研究方向。     

基于丝网探针的螺旋管内气液两相流气泡行为研究

核反应堆蒸汽发生器的传热面由螺旋管束组成。螺旋管的三维螺旋结构使得泡状流和塞状流等气液两相流中的气泡在重力、离心力和浮力等作用下在管道内部呈现不对称的相分布状态,两相滑移速度增大,显著影响换热性能并导致DNB型传热恶化难以预测。实验介质为空气-水,结合自主开发的电导式丝网探针技术并发展先进的数据后处理算法,实现了复杂流场的三维时空重构和离散气泡粒径的精细测量,获得了螺旋管内泡状流和塞状流的截面空泡分布规律。基于研究结果,可根据气泡分布规律对螺旋管道的几何结构进行调整以避免传热恶化,为螺旋管式蒸发器的安全设计提供了基础数据和优化思路。     

加热上升管内过冷流动沸腾数值模拟

采用计算流体动力学（CFD）程序CFX4.4对加热上升管内过冷流动沸腾工况下气水两相流动局部两相流参数（空泡份额和汽泡尺寸）进行了数值模拟。对数值差分方法、相关模型（界面力和气泡诱导的紊流）和汽泡尺寸进行了敏感性分析。空泡份额分布计算结果与实验结果比较表明,在低空泡份额工况下,两者符合较好,在高空泡份额工况下两者存在一定偏差,并且气相速度和汽泡尺寸的计算结果不理想。计算结果与实验结果之间的差异说明程序模型对于加热上升管内过冷流动沸腾模拟并不完善,建立更为合理的汽泡尺寸模型,考虑汽泡的合并和撕裂是必要的。     

Simulation of oscillatory baffled column: CFD and population balance

In the present work, an attempt has been made to combine population balance and a CFD approach for simulating the flow in oscillatory baffled column (OBC). Three-dimensional Euler-Euler two-fluid simulations are carried out for the experimental data of Oliveira and Ni [2001. Gas hold-up and bubble diameter in a gassed oscillatory baffled column. Chemical Engineering Science 56, 6143-6148]. The experimental data include the average hold-up profile and bubble size distribution in the OBC. All the non-drag forces (turbulent dispersion force, lift force) and the drag force are incorporated in the model. The coalescence and breakage effects of the gas bubbles are modeled according to the coalescence by the random collision driven by turbulence and wake entrainment while for bubble breakage by the impact of turbulent eddies. Predicted liquid velocity and averaged gas hold-up are compared with the experimental data. The profile of the mean bubble diameter in the column and its variation with the superficial gas velocity is studied. Bubble size distribution obtained by the model is compared with the experimental data.     

注气孔位置对文丘里管式微气泡发生器成泡特性的影响分析

尽管文丘里管式微气泡发生器的注气口位置会对气泡在文丘里流道内的碎化特征产生直接影响,但迄今缺乏针对性的深入研究。通过可视化实验方法,对比分析了注气口分别位于喉管处（结构1型）和进水管处（结构2型）时的气液流型、气泡破碎特征以及成泡特性。实验表明,气、液相流量对结构1型微气泡发生器内的气液流型影响显著,初始成泡区域随液相流量增加,环状流或泡状流向弹状流转变,而随气相流量增加则由泡状流或弹状流向环状流转变;结构2型微气泡发生器则在此过程中始终为泡状流,其对操作工况的适应范围大于结构1型。在相同工况下,结构1型微气泡发生器的成泡Sauter平均粒径小于结构2型,但随着液相Reynolds数的增大,二者间的成泡平均粒径差值随之减小。分析原因是由于弹状流流型下,延伸至扩张段区域的弹型泡的表面积更大,能量转化率更高,气泡界面失稳碎化的程度更显著。随着液相Reynolds数的增大,初始成泡体积减小,湍流破碎机理作用占据主导,掩盖了由于界面失稳引起的气泡破碎。结构1型微气泡发生器的成泡能耗高于结构2型,并且随液相Reynolds数的增大,两者之间的差值随之增大。综合来看,结构2型微气泡发生器能够在低能耗下实现高效成泡,面向工程应用将更具优势。     

鼓泡塔中气泡尺寸分布和局部气含率研究

在内径为0.38 m的鼓泡塔中采用双电导探针法对不同通气速率下的气泡尺寸分布和局部气含率进行了实验研究,分析了气泡尺寸的概率密度分布。结果表明:气泡尺寸随轴向高度的增加而增大,随径向距离增加而减小;鼓泡塔中气液流动可分为过渡流域和充分发展流域,在过渡流域气含率随轴向高度增加而增大,在充分发展流域气含率趋于均值,径向局部气含率分布呈抛物线型下降。高气速下气泡尺寸概率密度分布比低气速下宽,且随轴向高度的增加分布变宽。     

Regime transition in viscous and pseudo viscous systems: A comparative study

A comprehensive quantitative study on the effect of liquid viscosity (1 ≤ µ_L ≤ 1149 mPa‐s) on the local flow phenomena of the gas phase in a small diameter bubble column is performed using ultrafast electron beam X‐ray tomography. The internal dynamic flow structure and the bubble size distribution shows a dual role of the liquid viscosity on the hydrodynamics. Further, the effect of solid concentration (C_s = 0.05, 0.20) on the local flow behavior of the gas phase is studied for the pseudo slurry viscosities similar to the liquid viscosities of the gas–liquid systems. The effects of liquid and pseudo slurry viscosities on flow structure, bubble size distribution, and gas phase distribution are compared. The bubble coalescence is significantly enhanced with the addition of particles as compared to the system without particles for apparently same viscosity. The superficial gas velocity at which transition from homogeneous bubbly to slug flow regime occurs is initiated by the addition of particles as compared to the particle free system for apparently same viscosity. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3079–3090, 2014     

Phase distribution in dispersed liquid–liquid flow in vertical pipe: Mean and turbulent contributions of interfacial force

We applied an Eulerian–Eulerian two‐fluid model on an upward dispersed oil–water flow in vertical pipe with 80 mm diameter and 2.5 m length. The numerical profiles of the radial distribution of the oil drops at 1.5 m from the inflow are compared to the experimental data of Lucas and Panagiotopoulos (Flow Meas Instrum. 2009;20:127–135) This article analyzes the roles of turbulence and interfacial forces on the phase distribution phenomenon. In liquid–liquid flow the relative velocity is low and the distribution of the dispersed phase is mainly governed by the turbulence. This work highlights the important role of the turbulent contribution obtained by averaging the added mass force on the radial distribution profiles of the oil drops. The numerical results present improved profiles of the dispersed phase comparing to the experimental data when this turbulent contribution is taken into account in the momentum balance. © 2017 American Institute of Chemical Engineers AIChE J, 63: 4214–4223, 2017     

Mean and turbulent fluctuating velocities in oil-water vertical dispersed flows

Experiments of oil-water upward and downward flows have been carried out in a 38 mm ID pipe to investigate the modifications of turbulent flow characteristics by the presence of dispersed phase, i.e., mean and turbulent velocity profile of the continuous phase and mean velocity profile of the dispersed phase. Results for both oil-in-water (o/w) and water-in-oil (w/o) dispersions are presented. In o/w upward flow, the axial mean velocity profiles are found to be flatter than in single-phase flow and then change to centre peaked as the input oil fraction increases; a flatter profile is seen in w/o upward flow. In downward flow, the presence of oil drops always tends to flatten the continuous phase velocity profile in o/w dispersions, while a slightly centre peaked profile is observed in all cases of w/o systems. For both upward and downward flows, the presence of the dispersed phase tends to flatten the turbulence intensity profile and to result in a more uniform distribution of the turbulent energy over the pipe cross-section. It is also found that turbulence is more likely to be enhanced in the pipe centre area, where the volume fraction and the size of the dispersed phase are larger, while suppressed in the area close to the wall. Turbulence intensity is increased with mixture velocity and is slightly higher in upward than in downward flows. The current study suggests that local dispersed phase fraction and size as well as dispersed phase velocity seem to affect turbulence characteristics in oil-water flows. Previous models based on particle-laden flows for the prediction of turbulence enhancement or suppression were examined and agreement was found to depend on the type of dispersion (i.e., whether oil or water constitute the continuous phase).     

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.     

A dual-scale turbulence model for gas-liquid bubbly flows☆

A dual-scale turbulence model is applied to simulate cocurrent upward gas–liquid bubbly flows and validated with available experimental data. In the model, liquid phase turbulence is split into shear-induced and bubble-induced turbulence. Single-phase standard k-εmodel is used to compute shear-induced turbulence and another transport equation is added to model bubble-induced turbulence. In the latter transport equation, energy loss due to interface drag is the production term, and the characteristic length of bubble-induced turbulence, simply the bubble diameter in this work, is introduced to model the dissipation term. The simulated results agree well with experimental data of the test cases and it is demonstrated that the proposed dual-scale turbulence model outperforms other models. Analysis of the predicted turbulence shows that the main part of turbulent kinetic en-ergy is the bubble-induced one while the shear-induced turbulent viscosity predominates within turbulent vis-cosity, especially at the pipe center. The underlying reason is the apparently different scales for the two kinds of turbulence production mechanisms:the shear-induced turbulence is on the scale of the whole pipe while the bubble-induced turbulence is on the scale of bubble diameter. Therefore, the model reflects the multi-scale phe-nomenon involved in gas–liquid bubbly flows.