DSMC-LES方法研究循环流化床内颗粒团聚物的流动特性

采用直接模拟Monte Carlo方法法DSMC）模拟颗粒间的碰撞,采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟（LES）研究气相湍流.单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定.数值模拟垂直管内气固两相上升流动,对管内气相速度和颗粒相速度、浓度以及聚团流动进行分析.研究平均单个颗粒团聚物的存在时间、颗粒团聚物的时间份额和颗粒团聚物的生成频率分布特性,模拟结果与文献的实验结果基本吻合.     

DSMC-LES方法数值模拟鼓泡流化床内气泡和颗粒流动行为

基于稠密气体分子动力学和气固两相流体动力学,建立流化床稠密气固两相离散颗粒运动-碰撞解耦模型,采用直接模拟蒙特卡罗方法(DSMC)模拟颗粒间的碰撞,采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟(LES)研究气相湍流,单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定.数值模拟流化床中颗粒流动以及气泡的生成、长大和破碎过程,获得颗粒轴向和径向速度的概率密度分布,及颗粒浓度分布.计算结果表明床内气泡的形成造成床内颗粒的循环,使得流化床内颗粒具有不同的轴向和径向脉动速度,颗粒分速度分布近似服从高斯分布.颗粒温度随颗粒浓度增加,达到最大值后,随颗粒浓度增大而下降.流化床颗粒浓度脉动主要是低频部分,高频分量较低,表明在流化床内颗粒浓度脉动频率低,能量高,颗粒浓度脉动主频率为0.4～1.0Hz,其值与Pain et al.获得的颗粒浓度脉动主频率基本吻合.     

Mixing and combustion in a shallow coal-limestone fluidized bed combustor

A model based on the Monte Carlo approach was developed to simulate the mixing and combustion behavior of a shallow coal-limestone fluidized bed combustor. The model involved the coupling of two sub-models: a combustion sub-model based on the two-phase concept of fluidization and a mixing sub-model based on our previously developed dynamic mixing model. The combustion sub-model considered both the volatile and char combustion. It assumed that the combustor consisted of three distinct phases, i.e., jet, bubble and emulsion, with combustion occurring only in the emulsion phase. The mixing sub-model considered the upward or downward movement of a coal particle in the bed as being governed by certain probability laws; these laws were, in turn, affected by the bubbling hydrodynamics. In all, the combustor simulation model took into consideration the effects of coal feed rate, coal size distribution, limestone size, air flow rate and combustor temperature on the combustor behavior. The simulation results included the dynamic response of coal concentration profile, coal size distribution, coal particle elutriation rate as well as the mixing status between the coal and limestone particles.     

Analysis of cyclic combustion of solid fuels

The cyclic nature of coal particles combustion results from the movement of loose material in the flow contour of the circulating fluidized bed (CFB): the combustion chamber, the cyclone, the downcomer.The experimental results proved that the cyclic change of the oxygen concentration around coal particles, led to the vital change of both mechanism and combustion kinetics. The mathematical model of the process of coal combustion has been scientifically described whose original concept is based on the allowance for cyclic changes of concentrations of oxygen around the char particle. It enables the prognosis for change of the surface and the centre temperatures and a mass loss of the char particles during the cyclic combustion. It allows to appoint mass-rate of combustion of a char particle in the above conditions.     

Modeling and experimental studies on single particle coal devolatilization and residual char combustion in fluidized bed

Single particle devolatilization followed by combustion of the residual coal char particle has been analyzed in a batch-fluidized bed. The kinetic scheme with distributed activation energy is used for coal devolatilization while multiple chemical reactions with volume reaction mechanism are considered for residual char combustion. Both the models couple kinetics with heat transfer. Finite Volume Method (FVM) is employed to solve fully transient partial differential equations coupled with reaction kinetics. The devolatilization model is used to predict the devolatilization time along with residual mass and particle temperature, while the combined devolatilization and char combustion model is used to predict the overall mass loss and temperature profile of coal. The computed results are compared with the experimental results of the present authors for combustion of Indian sub-bituminous coal (15% ash) in a fluidized bed combustor as well as with published experimental results for coal with low ash high volatile matter. The effects of various operating parameters like bed temperature, oxygen mole fraction in bulk phase on devolatilization time and burn-out time of coal particle in bubbling fluidized bed have been examined through simulation.     

Comparison of fibre optical measurements and discrete element simulations for the study of granulation in a spout fluidized bed

Spout fluidized beds are frequently used for the production of granules or particles through granulation. The products find application in a large variety of applications, for example detergents, fertilizers, pharmaceuticals and food. Spout fluidized beds have a number of advantageous properties, such as a high mobility of the particles, which prevents undesired agglomeration and yields excellent heat transfer properties.A discrete element model is used describing the dynamics of the continuous gas phase and the discrete droplets and particles. For each element momentum balances are solved. The momentum transfer among each of the three phases is described in detail at the level of individual elements.The results from the discrete element model simulations are compared with local measurements of time time-averaged particle volume fractions as well as particle velocities by using a novel fibre optical probe in a fluidized bed of 400 mm I.D. Simulations and experiments were carried out for three different cases using Geldart B type aluminium oxide particles: a freely bubbling fluidized bed; a spout fluidized bed without the presence of droplets and a spout fluidized bed with the presence of droplets. It is found that the experimental and numerical results agree in a qualitative manner.It is demonstrated how the discrete element model can be used to obtain information about the interaction of the discrete phases, i.e. the growth zone in a spout fluidized bed. Additional analysis of the numerical results indicates that liquid breakthrough does not take place for the studied test case.     

典型尺寸燃煤颗粒富氧燃烧特性及燃烧本征动力学研究

利用微型流化床加热速度快、温度分布均匀以及气体近平推流等优势,在直径20 mm自动控温的微型流化床反应分析仪中研究了粒度分布为1.7～3.35 mm和0.12～0.23 mm两种典型尺寸燃煤颗粒在790～900℃温度范围内的富氧燃烧行为。通过快速响应过程质谱对燃烧产生的烟气进行实时监测,成功地识别和记录了粗颗粒燃烧过程中经历的挥发分燃烧和原位新生半焦燃烧两个主要阶段。挥发分析出速度最快,然后快速燃烧,而半焦燃烧速度较慢。相比之下,细颗粒燃烧的这两个阶段具有几乎相同的速率,因而相互耦合而难以区分。根据实验结果,挥发分析出和燃烧为快速反应,煤颗粒燃烧过程速率受原位新生半焦燃烧过程控制。进一步研究了挥发分和原位新生半焦燃烧动力学行为,获得其本征动力学的活化能分别为107.2和143.9 kJ/mol。     

Determination of PVC powder drying kinetics at particle scale: Experimental study and modeling

An original experimental method is used to determine drying kinetic at particle scale. The particle scale kinetics was obtained by immersion of a small mass of wet polyvinyl chloride (PVC) particles (cake) in a batch dense fluidized bed containing inert hot particles (glass bead). The results are summarized here and prove clearly that the PVC drying is controlled by a competition between internal and external transfers. The drying kinetic was described by a particle scale model taking into account the convective–diffusive (mass transfer) and the convective–evaporative (heat transfer) phenomena. To validate this model with the experimental data, the experimental fluidized bed dryer is modeled following two different approaches: a perfect stirred reactor model and a 3D numerical simulation using the multiphase flow code NEPTUNE_CFD. The aim of this 3D simulation is to simulate the phenomena occurring, at local scale, in a dense fluidized bed dryer and to show the limitations of the perfect stirred reactor model.     

Depth of combustion inside char from evidence of oxidized pyrite

In coal combustion processes the rate controlling mechanism varies from external diffusion and pore diffusion to chemical kinetics. The particle zone where combustion takes place therefore changes from the outer layer, for large particles at high temperatures, to the total internal pore volume, for small particles at low temperatures. Partial penetration of oxygen occurs in intermediate cases. Recent publications report on fluidized bed experiments where, according to the model used, combustion takes place at or near the outer surface of the particle (‘shrinking particle model’) in a narrow boundary layer. Knowledge of the depth of this layer could contribute to combustion modelling. This paper shows that during fluidized bed combustion at 900 °C oxygen penetrates into the char to a depth of 50–100 μm. This is concluded from the width of the zone where pyrite particles in the char are oxidized. The presence of open pores may increase the depth of the internal combustion layer up to several hundreds of micrometres.     

鼓泡流化床内高灰煤气化过程的颗粒尺度特性研究

基于多相流体质点网格方法(MP-PIC)对高灰煤在三维鼓泡流化床气化过程进行了数值模拟研究。在欧拉-拉格朗日框架下将气相和固相分别视作连续介质和离散相处理。首先,将模拟得到的出口处气体组分结果与实验数据进行对比,实验数据与模拟结果具有良好的一致性。其次,研究了煤颗粒在气化炉内的温度、传热系数、速度和停留时间,从颗粒尺度揭示了鼓泡流化床气化炉内的颗粒分布特性和气固流动特征。结果表明:在气化炉入口附近煤颗粒与床层温差最大,传热系数最大;由于流化床内强非线性的气固流动,床中煤温度和传热系数的空间分布不均匀;煤颗粒和床料的瞬时速度具有稳定的波动幅度,其中垂直方向速度波动最明显,且煤颗粒的瞬时速度比床料的瞬时速度略大;由于颗粒间的剧烈碰撞,延长了煤颗粒停留时间。此外,对鼓泡流化床中煤气化过程颗粒尺度的研究,有助于深入了解固体颗粒的流动行为以及气固相相互作用特性,对鼓泡流化床反应器的设计优化具有重要意义。     

Modeling the thermochemical degradation of biomass inside a fast pyrolysis fluidized bed reactor

A fast pyrolysis process in a bubbling fluidized bed has been modeled, thoroughly reproduced and scrutinized with the help of a combined Eulerian/Lagrangian simulation method. The 3‐D model is compared to experimental results from a 100 g/h bubbling fluidized bed pyrolyzer including such variables as particle composition at the outlet and gas/vapor/water yields as a function of fluidization conditions, biomass moisture concentrations, and bed temperatures. Multiprocessor simulations on a high‐end computer have been carried out to enable the tracking of each of the 0.8 million individual discrete sand and biomass particles, making it possible to look at accurate and detailed multiscale information (i.e., any desired particle property, trajectory, particle interaction) over the entire particle life time. The overall thermochemical degradation process of biomass is influenced by local flow and particle properties and, therefore, accurate and detailed modeling reveals unprecedented insight into such complex processes. It has been found, that the superficial fluidization velocity is important while the particle moisture content is less significant for the final bio‐oil yield. © 2011 American Institute of Chemical Engineers AIChE J, 58: 3030–3042, 2012     

Model for biomass char combustion in the riser of a dual fluidized bed gasification unit: Part II — Model validation and parameter variation

The two-phase combustion model for biomass char combustion in a riser of a dual fluidized bed gasification unit that has been presented in part I is validated using the data obtained from the 8 MWth dual fluidized bed reactor at Guessing/Austria. The model is capable of calculating the average temperatures in all zones, the gas phase composition, solid hold up, char feed rates and air ratio. The model predictions for the temperature profile along the riser and for the exiting gas composition are in good agreement with the measured values. The simulation results show that the residual char from the gasifier is only partly converted in the riser for char particles larger than 0.6 mm. Un-combusted char is circulated back into the gasification reactor. Parameter variations show that the exact location where additional liquid fuels are introduced in the middle zone of the riser does not affect the global behaviour of the combustion reactor. Based on the simulation results it is proposed that external supply of char (additional) may be a very effective method for reducing producer gas recycling to the riser, which is currently necessary to obtain the desired gasification temperatures.     

Validation of the similar particle assembly (SPA) model for the fluidization of Geldart's group A and D particles

The similar particle assembly (SPA) model proposed for large‐scale discrete element method simulation was validated. The SPA model describes the scaling law for the motion of particles in an assembly, which is derived from the equation of motion of a particle. In the SPA model, particles with similar physical or chemical properties are represented by a single particle, which is called the representative particle. The local assembly of particles in the original bed is supposed to be maintained in the bed with the representative particles. The SPA model was validated by numerical simulations of a 2D fluidized bed of noncohesive and cohesive particles. Using the SPA model, the flow behavior of the representative particles of a certain magnification can be made similar to that of the original particles with significantly less computational load. A good approximation of bubble size distribution between the SPA bed and the original bed was obtained. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

气固流化床内鼓泡行为的离散颗粒数值模拟

采用离散颗粒数值模型探讨了气固流化床内多个气泡的运动规律. 数值计算考虑了颗粒与气体的相互作用以及颗粒间的相互作用,颗粒与气体的相互作用采用Gidaspow阻力修正公式,颗粒间的相互作用则采用硬球碰撞模式描述. 给出了一种在二维非结构网格下精确计算颗粒体积分数及两相全隐耦合求解颗粒速度的方法,同时,在硬球模型中对颗粒碰撞事件采用散列表的处理方法,大大提高了计算效率. 数值计算结果表明,由于多个气泡的相互作用而导致气泡的横向运动显著,气泡运动轨迹与气泡形状依赖于流化床内静压力的分布,气泡总是选择压力小、压差梯度大的路径上升,这一现象与气液分离流中气泡的运动规律类似.     

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 gas-solid 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.     

Combustion of single char particles in a turbulent fluidized bed

The combustion of single bituminous char particles (4-12 mm diameter) was studied in a turbulent fluidized bed operated at 1098 K using air as the fluidising medium. Results indicated that particles burn with constant density following a shrinking sphere model. Burning rates are much higher than those observed in a bubbling fluidized bed. The rate of transfer of oxygen to the particle surface is also higher than that observed in bubbling beds. A model is proposed to calculate the Sherwood numbers of the burning carbon particles. Experimental values of the Sherwood numbers agree well with those predicted from the model.     

有机硅单体合成的气固流化床的三维数值模拟

为了探索有机硅单体合成气固流化床内硅粉颗粒的流化特性,作者利用计算流体力学CFD软件,采用双欧拉气固两相流模型及SIMPLE算法,模拟了三维的气固流化床内硅粉颗粒的流化特性;分析了气泡生成、长大和破裂的过程,及不同床层高度的固体颗粒运动速度矢量图,不同床层高度处横截面颗粒体积分数变化。结果表明:三维模拟能直观的表现颗粒在流化床中的流化状态,为工业生产及应用提供了有效的依据。     

Modelling of dense gas-particle flow in a circulating fluidized bed by Distinct Cluster Method (DCM)

Computational Fluid Dynamics (CFD) is a powerful tool to study the dense gas-solid flow in a circulating fluidized bed. Most of the existing methods focus on the microscopic properties of individual particle. Therefore, the simulation scale is significantly limited by the huge number of individual particles, and so far the numbers of particles in most of the reported simulations are less than 10⁵. The hydrodynamics behaviour of particle clustering in a dense gas-solid two-phase flow has been verified by several experimental results. The Distinct Cluster Method (DCM) was proposed in this paper by studying the macroscopic particle clustering behaviour, and comprehensive models for cluster motion, collision, break-up, and coalescence have been well developed. We model the dense two-phase flow field as gas-rich lean phase and solid-rich cluster phase. The particle cluster is directly treated as one discrete phase. The gas turbulent flow is calculated by Eulerian approach, and the particle behaviour is studied by Lagrangian approach. Using the proposed method, a three-dimension dense gas-particle two-phase flow field in a circulating fluidized bed with square-cross-section, with particle number up to 7.162 × 10⁷ are able to be numerically studied, on which few results have been reported. Details on instantaneous and time-averaged distributions are obtained. Developing process of non-uniform particle distribution is visualized. These results are in agreements with experimental observations, which justified the feasibility of using the DCM method to model and simulate dense gas-solid flow in a circulating fluidized bed with large number of particle numbers.     

A descriptive model of carbon attrition in the fluidized combustion of a coal char

A model is proposed which takes into account the interaction between combustion in the pores, weakening of the mechanical structure of the particle and attrition during the fluidized combustion of a char. A constitutive relationship between the degree of enfeeblement of the outer region of particles and the rate of regression of the surface is introduced to characterize the behaviour of the char in respect to combustion assisted attrition. An orthogonal collocation technique is used to solve model equations.Results obtained in the batch fluidized bed combustion of the char of a bituminous coal, by means of a 40 mm ID laboratory unit especially equipped for time resolved collection of attrited carbon, are used to test the model. Attrited carbon rates and carbon combustion efficiencies measured in experiments carried out with oxygen concentration in fluidizing gas changing from 21 to 0.75% show a fair agreement with values given by model calculations.     

Modelling of circulating fluidized bed combustion of a char

The combustion of a char in the 41 mm ID riser of a laboratory circulating fluidized bed combustor has been investigated at different air excesses and rates of solids (char and sand) circulating in the loop. Riser performance was characterized by an axial oxygen concentration profile as well as by the overall carbon content and particle size distribution. The proposed model accounts for carbon surface reaction, intraparticle and external diffusion, and attrition. External diffusion effects were relevant in the riser dense region where char was potentially entrapped in large clusters of inert solids. Experimental data and results of the model calculations are in satisfactory agreement.