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

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

On the prediction of gas hold-up in two-phase flow systems using an Euler–Euler model

We quantify the ability of the two-fluid Euler–Euler model to predict the overall gas hold-up during two-phase flow in vertical columns using a combination of experiments and simulations. Gas hold-up in a bubble column and gas hold-up in the less-frequently studied co-current flow are investigated. For homogeneous flow characterized by nearly uniform bubble size, Euler–Euler model predictions are within 10% of the experimental values for both modes of operation, if the bubble diameter supplied as input to the model is the average bubble diameter in the physical system. This also holds true for heterogeneous flow in bubble columns despite the presence of a broad distribution of bubble sizes, if turbulence and bubble swarm effects on momentum exchange between phases are properly accounted for. Swarm corrections adequate for bubble columns, are less successful for co-current heterogeneous flow, for which gas hold-up predictions are least accurate (average error of 22%).     

二维鼓泡床内气液流动特性实验与数值模拟

采用高速摄像法测量了0.20 m×0.02 m×2.00 m拟二维床内气泡尺寸分布和流型等变化规律,结果表明,随着表观气速的增大,鼓泡床内依次呈现均匀鼓泡区、过渡区和湍动区3种形式,以气泡个数概率表示的气泡尺寸分布呈对数正态分布。以计算流体力学软件ANSYS CFX 10.0为平台,采用k-ε湍流模型和GRACE曳力模型对气液鼓泡床内流体动力学行为展开了数值模拟,其结果与实验值比较吻合。研究表明,从多相流理论出发的计算流体力学模拟方法可以用来预报鼓泡床内流型过渡等流体动力学特性。     

Bubble drift velocity from the bed collapse technique in three-phase fluidized beds

Transient behavior of a bed collapsing after cut-off of gas supply into a three-phase fluidized bed was determined in a 0.21 m-diameter half-tube acrylic column having a test section 1.8 m high. The transient behavior of the bed collapse after cut-off of the gas supply to the beds was monitored by a video camera (30 frames/s). A theory was developed to account for the dynamic behavior of the bed collapse after the gas supply shut-off to three-phase fluidized beds. The bubble drift velocity was theoretically calculated by gas and liquid phase holdups at steady state condition. At a liquid velocity of 0.103 m/s and gas velocity of 0–0.023 m/s, bubble size was uniform in the dispersed bubble flow regime. However, as the gas velocity increased above 0.023 m/s, the discrete or coalesced bubble flow regime could be observed. The agreement between the predicted and experimental values is acceptable in the dispersed bubble flow regime, but the agreement becomes poorer with increasing gas velocity.     

群体平衡模型对复杂气液泡状流数值模拟

引言 气液泡状流广泛出现于航空航天、石油化工、核工程技术等领域,且多为复杂湍流泡状流.在工程领域中,准确预测这种复杂泡状流的含气率、气一液速度、气泡直径等参数对工业设备安全与优化分析十分重要.     

Bubble column research in Japan

Many experimental studies on the bubble column have been reported by Japanese researchers since around 1960. They include studies of bubble behaviour, bubble size distribution, transition from the homogeneous bubbly flow regime to the heterogeneous liquid circulation regime, liquid velocity distribution, longitudinal liquid mixing, hydrodynamic modelling, the gas holdup, and the volumetric coefficient of gas-liquid mass tranfer k_La. Studies covered various modified bubble columns, such as the airlift reactor with an external or internal loop, the packed bubble column, and others. Performance of three-phase bubble columns, which deal with suspensions or emulsions, and their use as bioreactors or chemical reactors were also studied.     

Bubble Size Distribution in Oil‐Based Bubble Columns

A practical population balance model was used to evaluate the bubble size distribution in a bubble column. In addition, the bubble size distribution in the bubble column was measured at different gas velocities by photography and analysis of the pictures. Four types of liquid, i.e., water and three petroleum‐based liquids, were used in the experiments. The gas phase was air. It was found that the existing models in the literature are not able to satisfactorily predict the experimentally measured bubble size distribution. The model can be corrected by applying a correction factor to the energy dissipation rate. The corrected model fits the experimental bubble size distribution considerably better than the existing models. The variation of this correction factor is reported for different systems at different gas velocities.     

On the development of flow pattern in a bubble column reactor: Experiments and CFD

A comprehensive analysis of the development of flow pattern in a bubble column reactor is presented here through extensive LDA measurements and CFD predictions. In the LDA measurements, the simultaneous measurements of 2D velocity-time data were carried out at several radial locations and many axial cross-sections of the column for two different spargers. The profiles of mean axial liquid velocity, fractional gas hold-up and bubble slip velocity showed excellent agreement between the predictions and the experimentally measured values. The experimental results showed that the mean tangential velocity varies systematically in the radial as well as along the axial co-ordinates. The turbulence parameters viz. turbulent kinetic energy, energy dissipation rate and eddy diffusivity were also analysed. The estimated values of local energy dissipation rate obtained using eddy isolation model were used for establishing the energy balance in the column. The experimental data were used for the estimation of normal and shear stress profiles. For the case of single point sparger, just above the sparger region, the bubble plume was seen to have a strong tangential component of motion thereby yielding higher gas hold-up slightly away from the centre. This visual observation was well captured in profiles of all the hydrodynamic parameters obtained from the experimental data. CFD simulations of the mean velocities, gas hold-up and turbulent kinetic energy compared well with the experimental results.     

Stability analysis of bubble columns: Predictions for regime transition

A criterion for the transition from the homogeneous to the heterogeneous regime in a bubble column is developed based on the theory of linear stability. Hydrodynamics of bubble column is described by two-fluid model incorporating the interphase forces like drag force and added mass force. Added mass force affects the hydrodynamics of gas-liquid flows significantly and is formulated by taking into account the bubble deformation. A proper understanding of the nature of gas-liquid interface (clean or contaminated) is desired for the reliable predictions of the added mass coefficient. Data from the literature on the transition in bubble columns is critically analyzed. A good agreement has been obtained between the experimental transition gas hold-up and the predictions of the same obtained by the theory developed in this work.     

鼓泡塔列管传热系数的测量与估算

在φ500 mm×5000 mm的冷模实验装置中,使用自制传热探头,对无内构件的空塔和安装31根竖直换热列管的鼓泡塔内列管传热系数进行了测量.实验表明,列管传热系数随表观气速的增加而增大,传热系数沿径向呈抛物型分布,垂直列管内构件的加入使得传热系数的径向分布变得更为陡峭.基于表面更新理论,结合鼓泡塔内气含率和液速分布的测量及计算结果,提出了计算传热系数的数学模型.该模型既可以用于空塔的局部传热系数与平均传热系数计算,也可以用于安装列管束的局部传热系数与平均传热系数计算.模型计算值与实验数据符合良好,最大相对误差为5.62％.     

Radial Profiles of Gas Bubble Behavior in a Gas‐Liquid Bubble Column Reactor under Elevated Pressures

The effects of operating conditions on radial variation of gas holdups, bubble swarm rising velocity, and Sauter mean diameter were investigated in a bubble column reactor under elevated pressures using a conductivity probe method. Air served as gas phase and tap water as liquid phase with varying gas velocity and pressure. All three parameters increased constantly with higher superficial gas velocity. Maximum holdup, velocity, and Sauter mean diameter were found at the center of the cross section. Two different cases for Sauter mean diameter distribution were observed. The gas holdups increase continuously with higher system pressure, but decrease for bubble swarm rising velocity and Sauter mean diameter. According to experimental results, an empirical correlation of the gas holdup profiles is presented.     

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

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

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.     

筛板上液体横向流动时气泡的形成特性

An apparatus,desinged to simulate bubbling of a sieve tray operated in froth regime,was employed. Bubble contact angles in and above the incipient weeping regimer for an air-water-plexiglas system were investigated. The influence of both liquid cross-flow and gas up-flow upon bubble contact angles was examined. A model considering the influence of liquid cross-flow was developed to predict bubble size from a sieve hole in froth operation regime.The comparison shows that the bubble sizes predicted by the present model are consistent with our experimental values and the available published experimental data.     

Solid Particle Effects in Centi-scale Slurry Bubble Columns

Both solid particles and column diameter affect the gas holdup and flow regimes in slurry bubble columns, but investigations of the combined effects are not to be found. This study shows the simultaneous impacts on the overall gas holdup and flow regime transitions and determines the dominant effects in slurry bubble columns on the centi-scale containing solid particle concentrations up to 20 vol %. Additional tomography measurements are presented to visualize the gas phase flow and the spatial gas phase distribution in the column.     

Bubble flow in trickle beds: investigations using resistive sensors

The local hydrodynamics of co-current gas-liquid down-flow through porous media are investigated in a quasi two-dimensional regular arrangement by means of a network of resistive sensors. The investigations are focused on a liquid continuous flow regime, the dispersed bubble flow, where the gas is divided into elongated bubbles. Due to the variation of the local flow channel orientation and the local void fraction, the average bubble velocity strongly depends on the local geometry. The flow is more coherent in vertical constrictions, compared to all other types of sites, this is probably due to bubble stagnation in the flow channel enlargements. At a given liquid superficial velocities and for sufficiently high superficial gas velocities, the average bubble size is independent of the gas flow rate; it is of the order of magnitude of the volume of the enlargements of the porous medium. The maximum bubble size is about three times its average size, corresponding thus to the coalescence of three average sized bubbles.     

湍动浆态床流体力学研究（Ⅲ）垂直列管内构件的影响

浆态床反应器内部往往安装有密集的竖直列管换热内构件,有关列管束对流场分布的影响还少有研究。本文在Φ500 mm×5000 mm的大型冷模实验装置中测定了安装不同密度列管束时的速度分布和气含率分布。实验表明,列管的存在一方面会显著提高浆料轴向速度,促进大尺度流体循环,另一方面也抑制了液体与气泡的径向湍动,使速度和气含率的径向分布更不均匀,造成液体与气体轴向返混加剧,增大了有列管束的浆态床反应器的放大风险。在低气速湍动鼓泡条件下,列管加入造成的“烟囱效应”将更为显著。     

Studies on hydrodynamics of turbulent slurry bubble column（Ⅲ）Effect of vertical pipe bundles

浆态床反应器内部往往安装有密集的竖直列管换热内构件,有关列管束对流场分布的影响还少有研究。本文在Φ500 mm×5000 mm的大型冷模实验装置中测定了安装不同密度列管束时的速度分布和气含率分布。实验表明,列管的存在一方面会显著提高浆料轴向速度,促进大尺度流体循环,另一方面也抑制了液体与气泡的径向湍动,使速度和气含率的径向分布更不均匀,造成液体与气体轴向返混加剧,增大了有列管束的浆态床反应器的放大风险。在低气速湍动鼓泡条件下,列管加入造成的“烟囱效应”将更为显著。     

Flow Regime Identification in a Slurry Bubble Column from Gas Holdup and Pressure Fluctuations Analysis

Analysis of gas holdups and pressure fluctuations are conducted in a slurry bubble column to study changes in flow regimes and comparisons are made with solid‐free bubble column. The pressure fluctuations are measured by fast response pressure transducers mounted on the column wall in the distributor and bulk regions. Air, tap water and 35 micron glass beads are used as the gas, liquid and solid phases respectively. Statistical analysis of pressure fluctuation data combined with gas holdup analysis provided information about flow regime transition and interesting insights into bubble size distributions and changes in flow structure.     

Population balance modelling for bubbly flows with heat and mass transfer

Population balance equations combined with a three-dimensional two-fluid model are employed to predict bubbly flows with the presence of heat and mass transfer processes. Subcooled boiling flow belongs to this specific category of bubbly flows is considered. The MUSIG (MUltiple-SIze-Group) model implemented in CFX4.4 is further developed to account for the wall nucleation and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements near the test channel exit is made for the radial distribution of the bubble Sauter diameter, void fraction, interfacial area concentration and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Additional comparison was also performed against existing boiling model in CFX4.4 and the modified model developed in our previous work (Int. J. Heat Mass Transfer 45 (2002) 1197). Good agreement is better achieved with the local radial bubble Sauter diameter, void fraction, interfacial area concentration and liquid velocity profiles against measurements using the newly formulated MUSIG boiling model over the simpler boiling models. However, significant weakness of the model is still evidenced in the prediction of the vapour velocity. Work is in progress to circumvent the deficiency of the model by the consideration of additional momentum equations or an algebraic slip model to account for bubble separation.