Particle–wall interaction in entrained-flow slagging coal gasifiers: Granular flow simulation and experiments in a cold flow model reactor

Entrained-flow slagging coal gasifiers display large conversion efficiencies and small levels of unconverted carbon at the exhaust. Both features are apparently at odds with the fairly small “space-time” of the particle-laden gas feeding, as compared with the time scale of heterogeneous gasification of carbon. This apparent inconsistency can be explained by considering that fuel residence times are longer than the “space-time” due to segregation of fuel particles in the near-wall region of the gasifier. Segregation is promoted by swirl flow, by particle–wall interaction as the wall is covered by a molten layer of slag and by the establishment of a dense-dispersed flow of granular solids in the proximity of the wall.This study presents results of granular flow simulations of the interaction of a dense-dispersed particle flow with the confining wall. Simulations consider that both the particles and the wall may be either “sticky” or “non sticky”, based on the prevailing elastic vs plastic behaviour upon collision. The effect of drag forces exerted on particles by swirled gas flow is simulated in a simplified manner. Particle–particle collisions are modelled with a Hertzian approach that includes torque and cohesion. The extent and time scale of segregation of a lump of particles loaded into a cylindrical vessel are assessed. Results clearly indicate the different structure of the layer of particles establishing at the wall surface in the different interaction regimes.Results of simulations are qualitatively compared with results of an experimental campaign performed in a reactor representing a cold flow model of the entrained-flow gasifier, where solid, molten or semi-molten particles have been simulated by atomized wax as surrogate material. Altogether, the results confirm the importance of the particle–particle and particle–wall micromechanical interactions for a correct prevision of the segregation of fuel particles in entrained-flow slagging gasifiers.     

CFD–DEM simulation of particle deposition characteristics of pleated air filter media based on porous media model

In this study, the three-dimensional physical model of pleated air filtration media was simplified to porous media model, and the calculation parameters of porous media were obtained based on experimental data. The model of V-shaped pleated air filter media is constructed, the height of the media pleat is 50 mm and the pleat thickness is 4 mm, the pleat angle is 3.7°. The Hertz-Mindlin contact model was modified by Johnson Kendall Roberts (JKR) adhesion contact model. The deposition process of particles in media was simulated based on computational fluid dynamics (CFD) theory and discrete element method (DEM). Results show that the CFD–DEM coupling method can be effectively applied to the macro research of pleated air filter media. The particles will form dust layer and dendrite structure on the fiber surface, and the dust layer will affect the subsequent air flow organization, and the dendrite structure will eventually form a “particle wall”. The formation of the “particle wall” will prevent the particles from moving further in the fluid domain, which makes area of pleated angle become the “low efficiency” part about the particle deposition. Compared with area of pleated angle, the particles are concentrated in the opening area and the middle area of the pleated to agglomerate and deposit.     

Onset velocity of circulating fluidization and particle residence time distribution: A CFD–DEM study

Until now, the onset velocity of circulating fluidization in liquid–solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity, and is reported to be only dependent on the liquid and particle properties. This study presents a new approach to calculate the onset velocity using CFD–DEM simulation of the particle residence time distribution (RTD). The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity. Our results are in reasonable agreement with experimental data. The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter. A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity. The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity, but influence the particle RTD, showing some humps and trailing. The particle RTD is found to be related to the particle trajectories, which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

Effects of supply air temperature and inlet location on particle dispersion in displacement ventilation rooms

The effects of supply temperature and vertical location of inlet air on particle dispersion in a displacement ventilated (DV) room were numerically modeled with validation by experimental data from the literature. The results indicate that the temperature and vertical location of inlet supply air did not greatly affect the air distribution in the upper parts of a DV room, but could significantly influence the airflow pattern in the lower parts of the room, thus affecting the indoor air quality with contaminant sources located at the lower level, such as particles from working activities in an office. The numerical results also show that the inlet location would slightly influence the relative ventilation efficiency for the same air supply volume, but particle concentration in the breathing zone would be slightly lower with a low horizontal wall slot than a rectangular diffuser. Comparison of the results for two different supply temperatures in a DV room shows that, although lower supply temperature means less incoming air volume, since the indoor flow is mainly driven by buoyancy, lower supply temperature air could more efficiently remove passive sources (such as particles released from work activities in an office). However, in the breathing zone it gives higher concentration as compared to higher supply air temperature. To obtain good indoor air quality, low supply air temperature should be avoided because concentration in the breathing zone has a stronger and more direct impact on human health.     

Joint PDF Closure of Turbulent Premixed Flames

In this paper, a novel model for turbulent premixed combustion in the corrugated flamelet regime is presented, which is based on transporting a joint probability density function (PDF) of velocity, turbulence frequency and a scalar vector. Due to the high dimensionality of the corresponding sample space, the PDF equation is solved with a Monte-Carlo method, where individual fluid elements are represented by computational particles. Unlike in most other PDF methods, the source term not only describes reaction rates, but accounts for “ignition” of reactive unburnt fluid elements due to propagating embedded quasi laminar flames within a turbulent flame brush. Unperturbed embedded flame structures and a constant laminar flame speed (as expected in the corrugated flamelet regime) are assumed. The probability for an individual particle to “ignite” during a time step is calculated based on an estimate of the mean flame surface density (FSD), latter gets transported by the PDF method. Whereas this model concept has recently been published [21], here, a new model to account for local production and dissipation of the FSD is proposed. The following particle properties are introduced: a flag indicating whether a particle represents the unburnt mixture; a flame residence time, which allows to resolve the embedded quasi laminar flame structure; and a flag indicating whether the flame residence time lies within a specified range. Latter is used to transport the FSD, but to account for flame stretching, curvature effects, collapse and cusp formation, a mixing model for the residence time is employed. The same mixing model also accounts for molecular mixing of the products with a co-flow. To validate the proposed PDF model, simulation results of three piloted methane-air Bunsen flames are compared with experimental data and very good agreement is observed.     

Effect of source location on particle dispersion in displacement ventilation rooms

A proposed computer model for predicting aerosol particle dispersion in indoor spaces was validated with experimental data found in the literature,and is then used to study the effect of the area and point source locations on particle dispersion in displacement ventilation (DV) rooms.The results show that aerosol source location has a strong impact on the spatial distribution and removal rate of indoor particles.Particle removal performance depends strongly on ventilation efficiency and particle deposition rate on indoor surfaces.Important consideration for both relative ventilation efficiency and deposition rate consists of the position of the aerosol source relative to the main airflow pattern and the occupied zone.     

Effect of source location on particle dispersion in displacement ventilation rooms

A proposed computer model for predicting aerosol particle dispersion in indoor spaces was validated with experimental data found in the literature, and is then used to study the effect of the area and point source locations on particle dispersion in displacement ventilation （DV） rooms. The results show that aerosol source location has a strong impact on the spatial distribution and removal rate of indoor particles. Particle removal performance depends strongly on ventilation efficiency and particle deposition rate on indoor surfaces. Important consideration for both relative ventilation efficiency and deposition rate consists of the position of the aerosol source relative to the main airflow pattern and the occupied zone.     

亚微米颗粒在汇作用下运动机理的实验研究

当前,城市空气质量的不断恶化,引起了公众的普遍性关注.空气中的悬浮颗粒物,是城市大气环境重要污染源之一,其分布、运动及扩散规律已成为科学领域的研究热点.与连续流体不同,大气中的悬浮颗粒物是离散的,确定颗粒运动的模型是研究大气细微颗粒污染问题的关键.本文拟研究小空间静稳空气中亚微米级颗粒在汇作用下的运动规律,并构建其运动模型.在密闭实验空间中通过燃烧生成亚微米颗粒,利用静电吸附装置模拟颗粒汇,并通过粒子图像测速(particle image velocimetry,PIV)实验和激光多普勒测速仪(lasser Doppler velocimeter,LDV)实验技术测量分析不同空间内亚微米颗粒在大气中的热运动速度和在汇作用下的运动规律,并推导出颗粒物的速度分布经验公式.结果显示:粒子在汇作用下的运动与连续流体汇运动规律类似,但在小空间内颗粒的运动不满足流体连续方程;说明在无气流夹带输运情况下,利用汇作用及颗粒的扩散而发展的颗粒净化技术是可行的.     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

A stochastic approach for the simulation of collisions between colloidal particles at large time steps

This paper presents a new approach for the detection and treatment of colloidal particle collisions. It has been developed in the framework of Lagrangian approaches where a large number of particles is explicitly tracked. The key idea is to account for the continuous trajectories of both colliding partners during a time step that is not restricted. Unlike classical approaches which consider only the distances between a pair of particles at the beginning and at the end of each time step (or assume straight-line motion in between), we model the whole relative, and possibly diffusive, trajectory. The collision event is dealt with using the probability that the relative distance reaches a minimum threshold (equal to the sum of the two particle radii). In that sense, the present paper builds on the idea of a previous work. However, in this first work, the collision event was simulated with a simplified scheme where one of the collision partners was removed and re-inserted randomly within the simulation domain. Though usually applied, this treatment is limited to homogeneous situations. Here, an extension of the stochastic model is proposed to treat more rigorously the collision event via a suitable evaluation of the time and spatial location of the collision and an adequate calculation of subsequent particle motion. The resulting collision kernels are successfully compared to theoretical predictions in the case of particle diffusive motion. With these promising results, the feasibility of simulating the collisional regime over a whole range of particle sizes (even nanoscopic) and time steps (from a ballistic to a purely diffusive regime) with a numerical method of reasonable computational cost has been confirmed. The present approach thus appears as a good candidate for the simulation of the agglomeration phenomenon between particles also in complex non-homogeneous flows.     

Onset velocity of circulating fluidization and particle residence time distribution:A CFD-DEM study

Until now,the onset velocity of circulating fluidization in liquid-solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity,and is reported to be only dependent on the liquid and particle properties.This study presents a new approach to calculate the onset velocity using CFD-DEM simulation of the particle residence time distribution(RTD).The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity.Our results are in reasonable agreement with experimental data.The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter.A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity.The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity,but influence the particle RTD,showing some humps and trailing.The particle RTD is found to be related to the particle trajectories,which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

粗糙颗粒动理学及稠密气固两相流动的数值模拟

考虑颗粒碰撞过程中摩擦作用,给出了粗糙颗粒碰撞动力学.引入颗粒相拟总温来表征颗粒平动和转动脉动能量的特征.基于气体分子运动论,建立颗粒碰撞中平动和旋转共同作用的粗糙颗粒动理学,给出了颗粒相压力和黏度等输运参数计算模型.运用基于颗粒动理学的欧拉-欧拉气固两相流模型,数值模拟了流化床内气体颗粒两相流动特性,分析了颗粒旋转流动对颗粒碰撞能量交换和耗散的影响.模拟得到的流化床内径向颗粒浓度和提升管内颗粒轴向速度与他人实验结果相吻合.模拟结果表明随着颗粒浓度的增加,颗粒相压力和能量耗散逐渐增加,而颗粒拟总温先增加后下降.随着颗粒粗糙度系数的增加,床内平均颗粒相拟总温和能量耗散增加,表明颗粒旋转产生的摩擦将导致颗粒旋转脉动能量的改变,影响床内气体-颗粒两相宏观流动特性.      

大速差射流燃烧室中烟煤与贫煤燃烧的数值模拟

本文基于颗粒相的轨道模型,对大速差射流燃烧室中烟煤粉与贫煤粉的二维流动,混合及燃烧进行了数值模拟,模拟结果从两相耦合的角度,阐明了煤粉颗粒在燃烧室中运动的规律,煤粉与大速差射流诱导的中心气体逆流之间的混合及其对煤粉火焰稳定的影响,指出此种燃烧室中煤粉火焰稳定的回流区燃烧机理,气相流场及回流区的预报结果与实验符合良好。     

Modelling the bed characteristics in fluidised-beds for top-spray coating processes

A particle sub-model describing the bed characteristics of a bubbling fluidised bed is presented. Atomisation air, applied at high pressures via a nozzle positioned above the bed for spray formation, is incorporated in the model since its presence has a profound influence on the bed characteristics, though the spray itself is not yet considered. A particle sub-model is developed using well-known empirical relations for particle drag force, bubble growth and velocity and particle distribution above the fluidised-bed surface. Simple but effective assumptions and abstractions were made concerning bubble distribution, particle ejection at the bed surface and the behaviour of atomisation air flow upon impacting the surface of a bubbling fluidised bed. The model was shown to be capable of predicting the fluidised bed characteristics in terms of bed heights, voidage distributions and solids volume fractions with good accuracy in less than 5 min of calculation time on a regular desktop PC. It is therefore suitable for incorporation into general process control models aimed at dynamic control for process efficiency and product quality in top-spray fluidised bed coating processes.     

Cold flow experiments in an entrained flow gasification reactor with a swirl-stabilized pulverized biofuel burner

Short particle residence time in entrained flow gasifiers demands the use of pulverized fuel particles to promote mass and heat transfer, resulting high fuel conversion rate. The pulverized biomass particles have a wide range of aspect ratios which can exhibit different dispersion behavior than that of spherical particles in hot product gas flows. This results in spatial and temporal variations in temperature distribution, the composition and the concentration of syngas and soot yield. One way to control the particle dispersion is to impart a swirling motion to the carrier gas phase. This paper investigates the dispersion behavior of biomass fuel particles in swirling flows. A two-phase particle image velocimetry technique was applied to simultaneously measure particle and gas phase velocities in turbulent isothermal flows. Post-processed PIV images showed that a poly-dispersed behavior of biomass particles with a range of particle size of 112–160 µm imposed a significant impact on the air flow pattern, causing air flow decelerated in a region of high particle concentration. Moreover, the velocity field, obtained from individually tracked biomass particles showed that the swirling motion of the carrier air flow gives arise a rapid spreading of the particles.     

DEM simulation of particle percolation in a packed bed

The phenomenon of spontaneous particle percolation under gravity is investigated by means of the discrete element method. Percolation behaviors such as percolation velocity,residence time distribution and radial dispersion are examined under various conditions. It is shown that the vertical velocity of a percolating particle moving down through a packing of larger particles decreases with increasing the restitution coefficient between particles and diameter ratio of the percolating to packing particles. With the increase of the restitution coefficient,the residence time and radial dispersion of the percolating particles increase. The packing height affects the residence time and radial dispersion. But,the effect can be eliminated in the analysis of the residence time and radial dispersion when they are normalized by the average residence time and the product of the packing height and packing particle diameter,respectively.In addition,the percolation velocity is shown to be related to the vertical acceleration of the percolating particle when an extra constant vertical force is applied. Increasing the feeding rate of percolating particles decreases the dispersion coefficient.     

DEM simulation of particle percolation in a packed bed

The phenomenon of spontaneous particle percolation under gravity is investigated by means of the discrete element method. Percolation behaviors such as percolation velocity, residence time distribution and radial dispersion are examined under various conditions. It is shown that the vertical velocity of a percolating particle moving down through a packing of larger particles decreases with increasing the restitution coefficient between particles and diameter ratio of the percolating to packing particles. With the increase of the restitution coefficient, the residence time and radial dispersion of the percolating particles increase. The packing height affects the residence time and radial dispersion. But, the effect can be eliminated in the analysis of the residence time and radial dispersion when they are normalized by the average residence time and the product of the packing height and packing particle diameter, respectively. In addition, the percolation velocity is shown to be related to the vertical acceleration of the percolating particle when an extra constant vertical force is applied. Increasing the feeding rate of percolating particles decreases the dispersion coefficient.     

Importance of unsteady contributions to force and heating for particles in compressible flows. Part 2: Application to particle dispersal by blast waves

Particle dispersal by blast waves is an interesting phenomenon. A model problem, i.e., a sudden release of a compressed gas–particle mixture contained in a spherical container, is employed to investigate the fundamental physics of particle dispersal. The problem is simulated by the multiphase flow models proposed in Part 1 of this article that include unsteady contributions in momentum and energy exchange between gas and particles. At early times, when particles are accelerated in the expansion fan, unsteady force and heating contributions are much larger than the corresponding quasi-steady contributions. Consequently, neglecting unsteady contributions leads to significant errors in particle front location (the boundary of the particle cloud). The complex wave interactions in the flow have strong influence on the particle motion. As a result, the particle motion is a non-monotonic function of particle density or diameter and the evolution of particle concentration is non-uniform and unsteady.     

Discrete element-embedded finite element model for simulation of soft particle motion and deformation

The motion and deformation of soft particles are commonly encountered and important in many applications. A discrete element-embedded finite element model (DEFEM) is proposed to solve soft particle motion and deformation, which combines discrete element and finite element methods. The collisional surface of soft particles is covered by several dynamical embedded discrete elements (EDEs) to model the collisional external forces of the particles. The particle deformation, motion, and rotation are independent of each other in the DEFEM. The deformation and internal forces are simulated using the finite element model, whereas the particle rotation and motion calculations are based on the discrete element model. By inheriting the advantages of existing coupling methods, the contact force and contact search between soft particles are improved with the aid of the EDE. Soft particle packing is simulated using the DEFEM for two cases: particle accumulation along a rectangular straight wall and a wall with an inclined angle. The large particle deformation in the lower layers can be simulated using current methods, where the deformed particle shape is either irregular in the marginal region or nearly hexagonal in the tightly packed central region. This method can also be used to simulate the deformation, motion, and heat transfer of non-spherical soft particles.     

Solids behavior in dilute zone of a CFB riser under turbulent conditions

The behavior of the solid phase in the upper zone of a circulating fluidized bed riser was studied using a phase Doppler anemometer. Glass particles of mean diameter 107 μm and superficial gas velocities U_g covering the turbulent and the beginning of the fast fluidization regime were investigated. Three static bed heights were tested. Ascending and descending particles were found co-existing under all operating conditions tested, and at all measurement locations. Superficial gas velocity proved/happened to have a larger effect on descending particles at the wall and on ascending particles in the central region. Transversal particle velocities in both directions (toward the center and toward the wall) behaved relatively equivalently, with only slight difference observed at the wall. However, observation of the number of particles moving in either transversal direction showed a change in bed structure when increasing U_g. Furthermore, a balance was constantly observed between the core zone and the annulus zone where the mutual mass transfer between these two zones occurred continuously. Transition from a slow to a fast particle motion was accompanied by a transition to high levels of velocity fluctuations, and was found corresponding to the appearance of significant solid particle flow rate.