沉浸管式流化床的颗粒尺度模拟

为模拟具有复杂几何结构的气固流动系统,文中将计算流体力学和离散单元法与边界元方法结合起来,对沉浸管式流化床内颗粒及气泡的运动行为进行了数值模拟。模拟计算得到的瞬态流型图揭示了气泡绕流沉浸管束时出现的合并和破碎状态及颗粒群的详细运动行为,发现床内气固二相的流动受到沉浸管束存在的显著影响。当颗粒及气相的流动受到沉浸管的阻碍而绕管流动过程中气泡会发生变形,变得扭曲狭长且易被撕碎。同时颗粒与管道壁面碰撞会造成气固二相复杂的动态运动形式,床内的管道大部分时间会被气穴包围,将严重阻碍管道与颗粒之间的传热。     

Particle-scale simulation of fluidized bed with immersed tubes

In order to simulate gas-solids flows with complex geometry, the boundary element method was incorporated into the implementation of a combined model of computational fluid dynamics and discrete element method. The resulting method was employed to simulate hydrodynamics in a fluidized bed with immersed tubes. The transient simulation results showed particle and bubble dynamics. The bubble coalescence and break-up behavior when passing the immersed tubes was successfully predicted. The gassolid flow pattern in the fluidized bed is changed greatly because of the immersed tubes. As particles and gas are come in contact with the immersed tubes, the gas bubbles will be deformed. The collisions between particles and tubes will make the tubes surrounded by air pockets most of the time and this is unfavorable for the heat transfer between particles and tubes. __________ Translated from Chemical Engineering, 2007, 35(11): 21–24 [译自:化学工程]     

Computational fluid dynamics simulations of heat transfer between the dense‐phase of a high‐temperature fluidized bed and an immersed surface

The transient process of heat transfer between a high‐temperature emulsion packet and the wall of an immersed surface is simulated using computational fluid dynamics (CFD). From these simulations, the total heat transfer coefficient and its radiant contribution due to the emulsion (dense) phase are evaluated. The results are compared with experimental data (Ozkaynak et al., “An experimental investigation of radiant heat transfer in high temperature fluidized beds,” in Fluidization IV, 1983:371–378) and with predicted values from the generalized heterogeneous model (GHM), (Mazza et al., “Evaluation of overall heat transfer rates between bubbling fluidized beds and immersed surfaces,” Chem Eng Commun., 1997;162:125–149). The CFD simulations are in good agreement with both, experimental data and theoretical GHM predictions and provide a reliable way to quantify the studied heat transfer process. Also, the GHM is validated as a practical tool to this end. © 2011 American Institute of Chemical Engineers AIChE J, 58: 412–426, 2012     

CFD–DEM simulation of tube erosion in a fluidized bed

The erosion of the immersed tubes in a bubbling‐fluidized bed is studied numerically using an Eulerian–Lagrangian approach coupling with a particle‐scale erosion model. In this approach, the motion of gas and particles is simulated by the CFD–DEM method, and an erosion model SIEM (shear impact energy model) is proposed to predict the erosion of the tubes. The model is validated by the good agreement of the simulation results and previous experimental data. By analyzing the simulation results, some characteristics of the tube erosion in the fluidized bed are obtained, such as the distribution of the erosion rate around the tube, the variation of the erosion rate with the position of the tube, the effect of the friction coefficient of particles on the erosion, the relationship between the maximum and the average erosion rate, etc. The microscale behavior of particles around the tubes is also revealed and the linear relationship between the erosion and the shear impact energy is confirmed by the simulation results and experiment. The agreement between simulation and experiment proves that the microscale approach proposed in this article has high accuracy for predicting erosion of the tubes in the fluidized bed, and has potential to be applied to modeling the process in other chemical equipment facing solid particle erosion. © 2016 American Institute of Chemical Engineers AIChE J, 63: 418–437, 2017     

竖直管气固鼓泡流化床传热机理的CPFD模拟

针对细颗粒气固鼓泡流化床中床料与竖直传热管壁面间的传热行为,在前期实验的基础上,采用计算颗粒流体力学(CPFD)方法从颗粒在传热壁面更新的角度,深入分析了传热特性与壁面气固流动行为之间的关联性。结果表明,模拟得到的传热管壁面颗粒更新通量和基于颗粒团更新模型的颗粒团平均停留时间均能很好解释实验测得的传热系数变化规律,这证实颗粒团更新是影响传热过程的控制性因素。模拟还发现随加热管从床层中心向边壁的移动,加热管周向方向上颗粒更新通量和传热系数的不均匀性都呈增大趋势。随着表观气速的增大,气泡行为导致床层颗粒内循环流率增大,这是导致颗粒团在加热管壁面上的更新频率增大以及床层与壁面间传热系数增大的根源。     

埋管流化床内不同粒径颗粒传热行为的欧拉-拉格朗日模拟研究

将有限元方法、基于非结构化网格的计算流体力学方法与离散单元法结合,建立了CFD-DEM-FEM耦合方法,并在此基础上采用k-ε湍流模型及考虑颗粒间及气固间作用的多向耦合传热模型,对埋管流化床内的流动和传热行为进行了模拟和分析.通过计算结果从微观尺度探讨了埋管流化床内的传热机制,分析了影响床内传热的关键因素,得到了换热管道周围固含率和传热系数的分布规律,考察了颗粒直径对埋管周围传热系数的影响.数值模拟结果与实验数据基本一致,证实了CFD-DEM-FEM耦合方法模拟复杂气固流动和传热的可行性和准确性.     

Mixing behaviors of wet granular materials in gas fluidized bed systems

The discrete element method combined with computational fluid dynamics was coupled with a capillary liquid bridge force model for computational studies of mixing behaviors in gas fluidized bed systems containing wet granular materials. Due to the presence of strong capillary liquid bridge forces between wet particles, relative motions between adjacent particles were hindered. There was a high tendency for wet particles to form large aggregates within which independent motions of individual particles were limited. This resulted in much lower mixing efficiencies in comparison with fluidization of dry particles. Capillary liquid bridge forces were on average stronger than both fluid drag forces and particle–particle collision forces and this accounted for the difficulty with which individual particles could be removed and transferred between aggregates. Such exchange of particles between aggregates was necessary for mixing to occur during fluidization of wet granular materials but required strong capillary liquid bridge forces to be overcome. © 2013 American Institute of Chemical Engineers AIChE J, 59: 4058–4067, 2013     

不同入口结构下行床内气固流动及混合行为的CFD-DEM模拟

采用计算流体力学和离散单元模型（CFD-DEM）耦合一种简单的气固催化反应模型对具有不同入口结构的二维下行床内的气粒流动和混合行为进行全床数值模拟。模拟得到了不同入口结构下行床内的多尺度气固运动状态、全床的固含、速度及反应生成物浓度分布,以及气体和颗粒在下行床内的停留时间分布,发现入口结构对反应器内的流动、混合和气固接触效率起着关键性的作用,入口气体和颗粒的不均匀分布将导致下行床内气体停留时间的宽分布以及气固接触效果的恶劣。     

Particle scale study of heat transfer in packed and bubbling fluidized beds

The approach of combined discrete particle simulation (DPS) and computational fluid dynamics (CFD), which has been increasingly applied to the modeling of particle‐fluid flow, is extended to study particle‐particle and particle‐fluid heat transfer in packed and bubbling fluidized beds at an individual particle scale. The development of this model is described first, involving three heat transfer mechanisms: fluid‐particle convection, particle‐particle conduction and particle radiation. The model is then validated by comparing the predicted results with those measured in the literature in terms of bed effective thermal conductivity and individual particle heat transfer characteristics. The contribution of each of the three heat transfer mechanisms is quantified and analyzed. The results confirm that under certain conditions, individual particle heat transfer coefficient (HTC) can be constant in a fluidized bed, independent of gas superficial velocities. However, the relationship between HTC and gas superficial velocity varies with flow conditions and material properties such as thermal conductivities. The effectiveness and possible limitation of the hot sphere approach recently used in the experimental studies of heat transfer in fluidized beds are discussed. The results show that the proposed model offers an effective method to elucidate the mechanisms governing the heat transfer in packed and bubbling fluidized beds at a particle scale. The need for further development in this area is also discussed. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

垂直管内三相流化床沸腾传热特性

研究了三相流化床沸腾传热的特性和影响传热系数的诸因素。在传热过程中，由于固体粒子的存在，强化了传热。以玻璃球粒子为固相的三相流化床沸腾传热系数，是相同条件下汽液两相流沸腾传热的二倍。以铜粒子为固相的三相流化床沸腾传热系数，是相同条件下汽液两相流沸腾传热系数的３倍。     

振动流化床中大颗粒与水平管局部传热特性

在二维振动流化床中,以平均粒径1.83 mm的玻璃珠为物料,研究了大颗粒与水平管间局部传热规律;考察了气速、振动频率等因素对局部传热系数的影响,同时与小米和小玻璃珠实验结果进行对比。结果表明:大颗粒与小颗粒局部传热系数有很大差异;对于大颗粒,低速下局部传热系数随振动频率的增大先增加后减小,高速下局部传热系数随着振动频率的增加而降低;一定振动频率下,气速小时局部传热系数在60°左右达到最大,气速逐渐增加后,其最大值向90°转移。通过实验数据得到了计算大颗粒与水平管局部传热系数的关联式,计算值与实验值吻合较好,误差在±20%范围内。结果可为带浸没水平管的振动流化床设计和研究提供参考。     

带浸没管的惰性粒子振动流化床传热特性

采用带浸没加热管的惰性粒子振动流化床对膏状物料干燥进行了实验研究。考察了加料速率、进气温度、进气速度、加热管功率、振动强度等参数对床温和体积传热系数的影响,得出了计算体积传热系数的关联式。结果表明,在流化床中增设振动和浸没加热管装置,能大大强化传热传质,体积传热系数随加料量、振动强度、加热管功率、进风速度的增加而增大,随进气温度的增加而减小。其结果对惰性粒子流化床干燥器的设计和改进具有重要的指导意义。     

An efficient and reliable predictive method for fluidized bed simulation

In past decades, the continuum approach was the only practical technique to simulate large‐scale fluidized bed reactors because discrete approaches suffer from the cost of tracking huge numbers of particles and their collisions. This study significantly improved the computation speed of discrete particle methods in two steps: First, the time‐driven hard‐sphere (TDHS) algorithm with a larger time‐step is proposed allowing a speedup of 20–60 times; second, the number of tracked particles is reduced by adopting the coarse‐graining technique gaining an additional 2–3 orders of magnitude speedup of the simulations. A new velocity correction term was introduced and validated in TDHS to solve the over‐packing issue in dense granular flow. The TDHS was then coupled with the coarse‐graining technique to simulate a pilot‐scale riser. The simulation results compared well with experiment data and proved that this new approach can be used for efficient and reliable simulations of large‐scale fluidized bed systems. © 2017 American Institute of Chemical Engineers AIChE J, 63: 5320–5334, 2017     

网格尺寸对循环流化床内气固两相流动的影响

将基于能量最小多尺度方法(EMMS)的曳力模型耦合到双流体模型中,并针对循环流化床内的气固两流动进行了模拟研究。采用全滑移壁面边界条件处理颗粒相,考察了3种网格尺度对轴向空隙率和出口颗粒循环量等气固流动特性的影响。计算结果表明,应用EMMS曳力模型处理相间作用力,同时在采用全滑移壁面边界条件处理颗粒相时,双流体模型能够正确预测轴向空隙率分布。采用网格尺寸为2.325 mm×20 mm时,模拟结果和实测数据吻合较好,表明在循环流化床的数值模拟中选择恰当的网格尺度是极为重要的。     

A CFD model for predicting the heat transfer in the industrial scale packed bed

Compared to the traditional lumped-parameter model,computational fluid dynamics (CFD) attracted more attentions due to facilitating more accurate reactor design and optimization methods when analyzing the heat transfer in the industrial packed bed.Here,a model was developed based on the CFD theory,in which the heterogeneous fluid flow was resolved by considering the oscillatory behavior of voidage and the effective fluid viscosity.The energy transports in packed bed were calculated by the convection and diffusion incorporated with gaseous dispersion in fluid and the contacting thermal conductivity of packed particles in solids.The heat transfer coefficient between fluid and wall was evaluated by considering the turbulence due to the packed particles adjacent to the wall.Thus,the heat transfer in packed bed can be predicted without using any adjustable semi-empirical effective thermal conductivity coefficient.The experimental results from the literature were employed to validate this model.     

气固流化床外取热器内流动和换热特性分析

外取热器是维持催化裂化反应-再生系统热平衡和保持装置平稳运行的关键设备之一。外取热器的优化设计和合理调控,要求深入理解外取热器内的流动特性、换热特性及两者之间关系。在一套大型冷模热态实验装置上,分别考察了表观气速、颗粒质量流率对换热管附近的局部固含率和气泡频率、床层与换热管间传热系数的影响。结果表明:增加表观气速可以降低局部固含率、增加局部气泡频率、强化床层与换热管间换热;随着颗粒质量流率增加,局部固含率和局部气泡频率均增加;在较低表观气速下,增加颗粒质量流率不利于换热,而在较高表观气速下,传热系数随颗粒质量流率逐渐增加。不同流型下,气固流动特性对换热特性的影响不同。在鼓泡床流型下,过高的局部固含率不利于颗粒在换热表面的更新,增加换热管附近的局部气泡频率可以明显强化换热;而在湍流床流型下,换热管附近的局部固含率和气泡频率的增加,均使传热系数逐渐增大。建立了针对不同流型的换热经验关联式,预测值与实验值的平均相对偏差分别为6.9%和1.3%。     

Computational study of heat transfer in a bubbling fluidized bed with a horizontal tube

A combined approach of discrete particle simulation and computational fluid dynamics is used to study the heat transfer in a fluidized bed with a horizontal tube. The approach is first validated through the good agreement between the predicted distribution and magnitude of local heat transfer coefficient with those measured. Then, the effects of inlet fluid superficial velocity, tube temperature and main particle properties such as particle thermal conductivity and Young's modulus are investigated and explained mechanistically. The relative importance of various heat transfer mechanisms is analyzed. The convection is found to be an important heat transfer mode for all the studied conditions. A large convective heat flux corresponds to a large local porosity around the tube, and a large conductive heat flux corresponds to a large number of particle contacts with the tube. The heat transfer is enhanced by the increase of particle thermal conductivity while it is little affected by Young's modulus. Radiative heat transfer becomes increasingly important as the tube temperature is increased. The results are useful for temperature control and structural design of fluidized beds. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Effect of annular fins on heat transfer of a horizontal immersed tube in bubbling fluidized beds

Experiments were conducted in a bubbling air-fluidized bed to investigate the effect of annular fins of constant thickness on heat transfer. Steady state time averaged local heat transfer coefficient measurements were made by the local thermal simulation technique in a cold bubbling fluidized bed (90 mm ID, 260 mm tall) with horizontally immersed tube initially with no fin and then with three fixed annular fins of constant thickness. Silica sand of mean particle diameter 307 μm and 200 μm were used as the bed materials. The superficial velocity of air was from minimum fluidization conditions, u_mf, to approximately 3 × u_mf. The results indicate that, although the heat transfer coefficient falls with the use of fins, the total heat transfer rises as a result of the greater surface area. Increasing the particle diameter reduces the heat transfer coefficient not only for unfinned horizontal tube but also for annular finned horizontal tube at the same conditions of fluidized bed. Based on the experimental data, correlations are proposed for predicting heat transfer coefficient from fluidized bed to horizontally immersed tubes with and without fins.