流化床密相区颗粒扩散系数的CFD数值预测

应用离散颗粒模型直观获得颗粒运动情况,并从单个颗粒和气泡作用的角度分析颗粒运动和混合,证实气泡在床层中上升、在床层表面爆破以及气泡上升引起的乳化相下沉运动对颗粒混合起关键作用。应用基于颗粒动理学的双流体模型系统地对床宽分别为0.2、0.4、0.8 m的二维流化床在鼓泡区和湍动区的气固两相流动行为进行数值模拟。受离散颗粒模型启发,在双流体模型计算结果基础上,引入理想示踪粒子技术计算床内平均颗粒扩散系数。计算结果表明,颗粒横向扩散系数(D_x)总体上随流化风速增大而增大,但受床体尺寸影响较大;颗粒轴向扩散系数随流化风速增大而增大,受床体尺寸影响较弱。文献报道的密相区颗粒横向扩散系数分布在10^-4～10^-1 m²·s^-1数量级。本文提出的计算方法在数量级上与文献实验结果吻合,表明在大尺寸流化床且高流化风速下,颗粒横向扩散系数远大于小尺寸鼓泡流化床,为不同研究者实验结果的分歧提供了理论依据,也为预测大型流化床内颗粒扩散速率提供了放大策略。     

在具有锥形料腿的循环流化床中流化CaCO3超细颗粒

根据粘附性颗粒在流化过程中形成的聚团具有较宽粒径分布并因此导致大聚团在流化床中沉积和死床的问题,提出了循环流化床的锥形回料系统设计. 该回料系统包括两部分：锥形料腿和带辅助进气的V型阀. 实验证明,锥形料腿通过提供变化的表观流化气速,克服了流化聚团沉积死床等现象;而V型阀的辅助进气,对于保证V型阀顺利输送粘附性颗粒具有关键性作用. 借助这种回料系统,实现了高粘附性超细CaCO3颗粒在循环流化床的稳定快速流化. 从提升管内部拍摄的照片显示,尽管提升管采用较高的流化气体速度,但超细CaCO3颗粒仍然是以聚团的形式被流化. 对在提升管不同高度采集的聚团分析表明,处于快速流化状态的CaCO3聚团的直径远小于传统流化床中聚团的直径,并且在提升管高度方向聚团直径没有较大的变化. 同时实验还显示,提升管轴向空隙率呈S型分布,而径向则体现环-核结构,具有典型的快速床特征.     

气固并流下行床气体扩散行为的研究

采用氢气稳态示踪方法在内径140mm的气固并流下行循环流化床中对气体扩散行为进行了实验研究.实验结果表明:下行床中气体扩散行为可用二维拟均相模型进行描述,其气体的径向扩散系数与气速、固体循环量及颗粒密度的关系可用下列准数关联式表示Pe_r=4.35×10_(-3)Re~(0.95)ε~(-73.4) 1>ε>0.99而下行床中气体轴向扩散系数要比提升管中小1个数量级以上.     

气固流化床内宽筛分硅粉颗粒流化特性的数值模拟

为了探索气固流化床内宽筛分硅粉颗粒的流化特性,作者利用计算流体力学CFD软件,采用Eulerian气固多相流模型及SIMPLE算法,模拟了二维气固流化床内不同粒级硅粉颗粒在不同操作件下的气固流化特性;分析了气泡生成、长大和破裂的过程,研究了床内气固两相的流动特性.结果表明:模拟计算值与实验值吻合较好,最大相对误差为10...     

颗粒相壁面条件对非球形颗粒流动影响的数值模拟

采用双流体模型对设置竖直隔板的气固密相流化床中非球形颗粒的运动进行了模拟,颗粒形状的影响由相间曳力模型考虑,重点考察壁面处颗粒边界条件的影响。同时进行了实验室规模三维流化床的流化实验,以验证模型的有效性。通过压降轴向分布、颗粒浓度径向分布以及物料出口处颗粒质量流率功率谱估计等定量分析,结果表明：对不设置内构件的自由床,壁面反射系数对系统宏观流动特性影响较小,而对壁面处局部颗粒运动影响较大;对壁面面积大幅增加的内构件床,壁面反射系数可显著改变气体和颗粒的运动特征,取值需控制在适当范围内。     

喷动流化床流动特性数值模拟

为了捕捉喷动流化床中微观层次上的颗粒运动信息,建立了基于CFD的二维非稳态喷动流化床欧拉-欧拉两相流模型。分析了不同流化气速对喷动流化床气固流动特性的影响,即不同工况下的炉内压力降、颗粒浓度、床内空隙率分布、气体速度分布和固体颗粒速度分布。数值模拟研究结果表明:随流化气速的增大,压降和炉内平均空隙率逐渐增大,密相床层高度逐渐增加,沿着轴向方向的气体流量增大,喷动气的射流深度逐渐增加,同时射流半径也逐渐增加。     

颗粒相壁面条件对非球形颗粒流动影响的数值模拟

采用双流体模型对设置竖直隔板的气固密相流化床中非球形颗粒的运动进行了模拟,颗粒形状的影响由相间曳力模型考虑,重点考察壁面处颗粒边界条件的影响。同时进行了实验室规模三维流化床的流化实验,以验证模型的有效性。通过压降轴向分布、颗粒浓度径向分布以及物料出口处颗粒质量流率功率谱估计等定量分析,结果表明:对不设置内构件的自由床,壁面反射系数对系统宏观流动特性影响较小,而对壁面处局部颗粒运动影响较大;对壁面面积大幅增加的内构件床,壁面反射系数可显著改变气体和颗粒的运动特征,取值需控制在适当范围内。     

摇摆流化床的气固流动特性

采用二维床及D类玻璃珠颗粒,在表观气速U_g=0.267~0.978 m/s、摇摆幅值Θ=5°~15°、摇摆周期T=8~20 s的实验条件下,对摇摆流化床内气固流动过程及气体通过流化床的时均总压降进行了研究,并通过与常规直立床和倾斜床进行对比,分析了床体摇摆对气固流动的影响。结果表明,在平均角速度ω_ave>2(°)/s的条件下,当初始装料量和表观气速相同时,气体通过摇摆流化床的时均总压降低于直立床,高于相同最大倾角时的倾斜床;惯性力所产生的压降在0.15 kPa以下,其对床层压降的影响较小,床体倾斜导致气体向边壁区域聚集是影响摇摆流化床内气固流动特性的主要因素,由此导致床内存在固定床和下行移动床状态的非流化区域,使得处于流化区域的颗粒量减少,同时还降低了流化床层在竖直方向的静压。非流化区域的存在还会造成流化区域的气速高于直立床表观气速,两者表观气速之比为1.04~1.49。     

气固流化床多传感器信号时间序列分析

利用光纤、压力传感器同时获取气固流化床不同位置颗粒浓度、压差、压力信号时间序列：运用吸引子和涨落复杂性参数分析三种信号时间序列，研究细颗粒在低气速下的流化状况。实验结果证实了Takens的相空间嵌入原理，且在流化床由散式流化到聚式流化过渡的过程中，三种信号的复杂性分析都体现床内运动存在涨落过程，这与Wackerbauer研究的结论是吻合的。同时结果表明不同信号涨落复杂性参数随气速的变化趋势是一致的，有理由认为涨落复杂性参数对应的仅是系统本身的一定运动状况，且和轴向位置关系不大。     

液-固流化床中颗粒分级和混合的研究

本文从分析固体颗粒在液一固流化床中受力运动出发,提出了颗粒分级速度的解析表达式,从而建立起混合颗粒分级床中颗粒浓度分布的数学模型。在φ70mm的圆筒形流化床中以水作为流化介质,测定了双组分混合颗粒的轴向浓度分布,由此确定了团体颗粒的轴向扩散系数。在此基础上,依据因次分析法和实验数据推导出扩散系数与诸影响因素间的准数关联式。研究表明:床层在空隙度为0.75～0.85的范围内混合最好,混合颗粒中细颗粒的扩散能力比粗颗粒的扩散能力大,散式床的分级效果比聚式床好。     

微型流化床内气体返混

微型流化床基础和应用在近几年受到越来越多的关注。针对微型流化床对气固反应分析的应用要求,利用脉冲示踪法研究了内径10 mm和21 mm两种尺寸微型流化床中的气体返混特性,具体考察了管内径、颗粒静床高度、床料颗粒粒径和气体表观流速对气体返混程度的影响。结果表明:随着床内径、颗粒静床高度和表观气速的减小和床料颗粒粒径的增大,气体在床内的返混程度减小。使用粒径约270 μm粗颗粒时,两种床径的浅层微型流化床中的气体返混程度都较小,对应的Peclet数在27以上,证明了床内气体流动接近平推流,从而为利用微型流化床最小化气体返混对反应测试的影响,获得近本征反应动力学参数提供了流动特性的保障。     

RADIAL DISPERSION AND BUBBLE CHARACTERISTICS IN THREE-PHASE FLUIDIZED BEDS

The effects of gas and liquid velocities, liquid viscosity and particle size on the radial dispersion coefficient of liquid phase (D_r) and the bubble properties in three-phase fluidized beds have been determined. A new flow regime map based on the drift flux theory in three-phase fluidized beds has been proposed.

In three-phase fluidized beds, D, increases with increasing gas velocity in the bubble coalescing and in the slug flow regimes, but it decreases in the bubble disintegrating regime. The coefficient exhibits a maximum value in the bed of small particles with increasing liquid velocity at lower gas velocities. However, it increases with increasing liquid velocity at higher gas velocities. In two and three-phase fluidized beds of larger particles (6,8 mm), D_r exhibits a maximum value with an increase in liquid viscosity at lower gas velocities, but it increases at higher gas velocities. The mean bubble chord length and its rising velocity increase with increasing gas velocity and liquid viscosity. However, the bubble chord length decreases with an increase in liquid velocity and it exhibits a maximum value with increasing particle size in the bed. The radial dispersion coefficients in the bubble coalescing and disintegrating regimes of three-phase fluidized beds in terms of the Peclet number in the present and previous studies have been well represented by the correlations based on the concept of isotropic turbulence theory.     

快速流化床中气体轴向混合特性

在内径90mm、高8m的快速循环流化床试验装置上,用非理想脉冲示踪、双探头检测、计算机实时采样分析的方法,对湍流流态化、快速流态化和气力输送不同流动状态的气体混合程度分别进行了测量,结果表明均存在不同程度的轴向混合.轴向扩散系数主要与由气体速度和固体循环量决定的床层空隙度有关;用拟均相一维扩散模型对结果进行处理,时域拟合法求取模型参数,并将轴向扩散系数进行了关联,得到如下统一关联式D_e=0.1953ε~(-4.1197)     

Deterministic chaos: a new tool in fluidized bed design and operation

Deterministic chaos theory offers new and useful quantitative tools to characterize the non-linear dynamic behaviour of fluidized beds. The dimension and entropy of the fluidized bed's strange attractor can be used for various purposes, such as the classification of fluidization regimes or fluidized bed scale-up. This is illustrated by experimental and model simulation examples of deterministic chaotic behaviour in ambient gas-solids fluidized beds of Geldart B particles. It is shown that the Kolmogorov entropy is dependent on, amongst other parameters, the gas velocity and the bed aspect ratio. In dimensionless scaling of fluidized bed reactors this type of relationship can probably be of use in establishing full dynamic similarity.     

A two-dimensional gas fluidized bed for hydrodynamic and elutriation studies

A two-dimensional fluidized bed having the dimensions of 52.1 cm (20.5 in) by 2.00 cm (0.787 in) is designed and tested for its use in hydrodynamic and elutriation studies. The fluidization column is provided with a calming section and freeboard which are 45.7 cm (18.0 in) and 129.5 cm (51.0 in) high respectively. A porous distributor plate is provided whose pressure drop is found to vary linearly with air velocity in the range of current interest. Fluidization experiments with three sand particles (788, 488 and 167 μm), glass beads (427 μm), millet (2064 μm) and green peas (4578 μm) are reported. Bed expansion and bubble growth characteristics are examined in some detail. Variations of bed height and pressure drop with fluidization velocity are analyzed to establish bed voidage as a function of gas velocity, and minimum fluidization velocity. The latter is also measured for three particles in a 0.305 m square fluidized bed. These studies reveal that two-dimensional fluidized beds are great tools for making novel qualitative investigations for mechanistic details of processes taking place in three-dimensional fluidized beds. Currently, investigations are underway for elutriation phenomenon.     

柱形流化床传热特性的数值模拟

对Shedid等搭建的圆柱体流化床采用欧拉?欧拉法进行三维数值模拟,考察了颗粒球形度、表观进气速度和床料初始堆积高度对流化床内垂直加热壁面与流动床料之间对流传热特性的影响,采用有效导热系数分别计算气相和固相的对流传热系数。结果表明,随表观进气速度增大,流化床内颗粒物料湍流运动加剧,加热壁面平均温度和流体平均温度下降,壁面流体间传热平均温度差减小,壁面流体间对流传热系数增大;随初始床料高度增加,流化床内颗粒与加热壁面的接触面积增大,导致固相平均对流传热系数增大。     

气-固微型流化床压降特性及最小流化速度的实验研究

在内径3～20 mm的4个气-固微型流化床中,分别考察了A类和B类两种类型颗粒的流化特性,同时研究了床几何结构、操作条件、物相性质等各因素对其最小流化速度的影响.结果 表明,气-固微型流化床中的床层压降特性与颗粒类型密切相关,不同的流动状态下两种类型颗粒的流动特性存在显著地差异.在固定床阶段,与B类颗粒相比,A类颗粒与...     

Gas mixing in the cocurrent downflow circulating fluidised bed

An experimental study on the gas dispersion behaviour of a cocurrent downflow gas-solids suspension in a 140 mm i.d. circulating fluidised bed (CDCFB) using the steady-state tracer method is presented. The influence of gas velocity, solids circulating rate and particle density on radial gas dispersion has been examined. Gas dispersion can well be described by an eddy diffusion mechanism and a proposed two dimensional dispersed plug-flow model can fit the experimental data very well. Correlations of the radial diffusion coefficient were obtained. It is found that the axial diffusion coefficient obtained in the CDCFB is much lower than that in conventional circulating fluidised beds.     

Experimental investigation of particle contact time at the wall of gas fluidized beds

The contact time of particles at the walls of gas fluidized beds has been studied using a radioactive particle tracking technique to monitor the position of a radioactive tracer. The solids used were sand or FCC particles fluidized by air at room temperature and atmospheric pressure at various superficial velocities, covering both bubbling and turbulent regimes of fluidization. Based on the analysis of tracer positions, the motion of individual particles near the walls of the fluidized bed was studied. The contact time, contact distance and contact frequency of the particles at the wall were evaluated from these experimental data. It was found that in a bed of sand particles, the mean wall contact time of the fluidized bed of sand particles decreases by increasing the gas velocity in the bubbling and increases in the turbulent fluidization. In other words, the particle-wall contact time is minimum at the onset of turbulent fluidization in the bed of sand particles. However, the mean wall contact time is almost constant in both regimes of fluidization in the bed of FCC particles. All the existing models in the literature predict a decreasing contact time when the gas velocity in the bed is increased. It was also shown that the contact distance increases monotonously by increasing the gas velocity in the bed of sand particles, while it is almost constant for the bed of FCC particles. Contact frequency has a trend similar to that of the contact time for both sand and FCC particles.     

快速流态化气固两相间的动量交换

本文根据一维定常态流动模型,对快速流化床内气固两相间的相互作用及其机理进行了研究.结果表明,在快速流态化条件下,颗粒总是趋于聚集,以减小气固两相间的相互作用力,从而使曳力系数c_D小于单颗粒标准曳力系数c_(DS)(c_D/c_(DS)<1.0).c_D/c_(DS)不仅与(?)有关,而且受气固流动状况以及颗粒物性、床层直径等因素的影响.通过对大量数据的分析,得到预测曳力系数的经验关联式(平均相对偏差小于5%)c_D/c_(DS)=1.685(?)~(0.253)(Re_r/Re_t)~(-1.213)(d_p/D)~(0.105)