复合移动床电石反应器中颗粒运动的离散单元法模拟

针对复合移动床反应器内固体颗粒运动,采用离散单元法模型（DEM）考察布料器分别为扇形开口和矩形开口时,布料器转速和开口对颗粒运动的影响,并基于文献结果论证了本文模型的准确性。模拟结果表明：①对于不同布料器,颗粒在移动床中呈现平推流和汇聚流两种流动形态。②随布料器转速及开口的增加,颗粒质量通量非线性增加。③随布料器转速的增加,下落床径向上颗粒分布更均匀;随布料器开口的增大,下落床径向上颗粒分布范围变大,颗粒分布更均匀;对下落床径向上颗粒分布,布料器扇形开口时分布呈U形、矩形开口时分布呈M形。④沿反应器轴向向下,颗粒分布有均匀化趋势;扇形开口布料器对颗粒分布的离散系数大于1,矩形开口布料器对颗粒分布的离散系数约为0.5。     

梯度磁场下气固流化床中磁颗粒运动数值模拟

在外加竖直方向梯度磁场的气固鼓泡流化床中,考虑铁磁颗粒受到的梯度磁场力和颗粒间磁感应力,对气相采用流体力学方法(CFD),颗粒相采用离散元法(DEM),建立二维磁鼓泡流化床数学模型,模拟不同磁场强度下全磁颗粒圆形床料的气固流动过程,分析了不同磁场强度对磁流化床中气泡生长、颗粒运动、床层压降和磁颗粒受力的影响。研究结果表明在沿高度磁感应强度递减的梯度磁场中,磁颗粒在颗粒间磁感应力的作用下凝聚成链,破坏了大气泡的形成。随着磁场强度增加,颗粒扩散系数减小,颗粒间磁感应力和梯度磁场力增大;气体与颗粒相间作用力先减小、后增加;而颗粒接触力先增加、后减小。     

二维导流管喷动床离散颗粒动力学模拟

采用离散单元法(DEM)-计算流体力学(CFD)双向耦合数值方法对二维导流管喷动床进行了模拟,颗粒的运动通过DEM模型描述,而气体的运动用Navier-Stokes方程进行求解,气体和固体颗粒之间的相互作用通过曳力形式传递。文中将DEM和边界元方法(BEM)结合起来解决颗粒在具有复杂边界设备内的运动。通过采用BEM+DEM-CFD相结合的方法进行模拟计算,得到了喷动床的最小喷动速度,研究了不同表观气速下床内的流型,得到了二维导流管喷动床的床层压降与表观气速的关系,统计分析了喷射区、环隙区内颗粒的运动速度和相应的空隙率,全面地描述了二维导流管喷动床内的气固流动特征。     

颗粒移动床内不稳定运动的计算颗粒动力学模拟

采用考虑滚动摩擦的三方程离散单元法（DEM）模型对侧开孔的移动床中的颗粒流动进行了数值研究。结果表明,计算颗粒动力学(CGD)方法可对复杂颗粒系统内颗粒的运动行为进行准确的预测,包括时均速度场和脉动速度场。讨论了模型中颗粒摩擦参数的重要影响,并对颗粒流动表现出的间歇现象进行了分析。颗粒流动与流体流动有相似之处,都存在随机的脉动,但颗粒流的随机脉动机理与流体中的湍流机理有很大不同,颗粒流动会表现出很强的不连续性。     

水平气力输送管中颗粒浓度分布的数值模拟

以商业软件FLUENT为平台,通过UDF将DEM程序编译进FLUENT,利用DEM与CFD耦合的方法研究了水平管气力输送过程中颗粒粒径和密度对颗粒浓度分布的影响.结果表明,水平管道内颗粒浓度沿轴向方向不断降低,在管中心浓度逐渐扩散,使径向分布逐渐变得均匀;与双欧拉模型的计算结果不同,在管下部的高颗粒浓度区呈现出中心浓度高于边壁浓度的分布,而在管上部则表现为中心浓度低于边壁浓度的分布;随着颗粒密度及颗粒体积的增加,颗粒的跃动距离减少.     

移动床内颗粒共振行为的DEM模拟

移动床在各种过程工业中普遍存在。从颗粒尺度出发刻画移动床内颗粒物料的复杂流动行为，对于大型移动床反应器的设计、放大和优化具有重要的意义。本工作基于三维离散单元法(Discrete Element Method, DEM)数值模拟，考察了漏斗流和半整体流卸料流型下移动床内颗粒物料的运动特性，重点探讨了两种卸料流型下颗粒运动脉动特性的异同。首先通过对比流动区轮廓及其特征宽度随时间演化的模拟预测结果和实验测量结果，检验了DEM模拟结果的可靠性。DEM模拟结果表明，两种卸料流型下，流动区上部区域不同位置处颗粒平均轴向速度随时间的变化都呈现出显著的非随机波动，表现为颗粒平均轴向速度时间序列的离散傅里叶变换频谱图都出现了明显的特征峰。空间相关性分析结果表明，在流动区上部区域，不同位置处颗粒轴向速度随时间的波动呈现强烈的空间相关性，两种卸料流型下体系内都形成了局部共振。不同轴向位置处颗粒轴向速度的延迟相关性分析结果表明，这种共振行为源自于床层底部出口上方。对颗粒轴向速度和颗粒间接触力之间延迟相关性的分析结果显示，两者之间存在显著的相关性，且后者的波动先于前者的波动，说明模拟得到的共振现象可能源于自由...     

逆向气体射流对下行床颗粒混合的影响

下行床入口结构的研究一直被人们所重视。今在内径为0．192m的下行床中颗粒达到均匀分布的部位，沿床四周均布了三个45度方向逆流场气体射流入口，在此处设立气体入口可以使颗粒分散与气固快速接触不再同时进行。采用磷光颗粒示踪技术对下行床有逆流场射流气体存在时颗粒的轴径向混合行为进行了研究。这种逆流场射流气体对下行床颗粒的轴向混合行为无明显影响，在各操作条件下下行床内颗粒均能以接近平推流的方式运动；但该射流气体可以大大加强颗粒的径向混合，有利于气固接触，在下行床颗粒径向混合越差的操作条件下，射流气体对颗粒径向混合的影响效果越明显，下行床的这种入口结构具有良好的应用前景。     

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

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

下行床气粒流动行为的Eulerian-Lagrangian模拟

采用计算流体力学和离散单元方法耦合模型（CFD-DEM）对二维下行床内的气粒流动行为进行了全床数值模拟。模拟结果展示了下行床典型操作条件下特有的气固动态流动结构：沿流动方向存在明显的入口控制区、过渡区和（完全）发展区;颗粒聚团并不是出现在浓度相对较高的入口区,而是在过渡区之后的发展过程中逐渐形成较多的、松散的动态聚团结构。下行床发展段呈现典型的近壁浓环结构,这与实验结果基本一致。考察了颗粒之间以及颗粒与壁面之间的碰撞参数对下行床内气固流动结构的影响,发现完全弹性碰撞颗粒体系在入口区呈现最快速的颗粒分散;而对本文研究的操作条件,颗粒碰撞参数对发展段时均流体力学行为只产生轻微的影响。     

模拟颗粒流动的离散元方法及其应用

介绍了离散元 (DEM)方法的基本原理、颗粒运动控制方程和颗粒相互作用力的数学模型。综述了DEM在流化床和固定床反应器 ,以及一些单元操作如料仓卸料过程、混合过程等中的最新应用和研究结果 ,表明DEM能够反映过程的本质机理 ,可以利用基本的数据模拟复杂的颗粒流动系统。最后指出了DEM发展中亟待解决的问题     

Numerical simulation of particle motion in vibrated fluidized beds

Flow behavior of particles is simulated in a two-dimensional vibrated bubbling fluidized bed. The motion of particles is simulated by discrete element method (DEM). The distributions of velocity and concentration of particles are analyzed at the different amplitudes and frequencies of vibration. In the case with vertical vibration, the bed consists of three different regimes along bed height: a low particle concentration regime with a vibration gap near the bottom, a high concentration regime in the middle of the bed and a transition regime at the bed surface. The gas pressure losses are reduced with the vibration. The phase delay time between the time at the maximum pressure loss and the time at the maximum displacement of vibrated distributor in one cycle decreases with the increase of amplitudes. The contacting interaction between the vibrated distributor and particles increases as the distributor moves up, and trends to zero as the distributor moves down with the increase of frequency of vibration. The vibrating energy is transferred by collisions between the vibrated distributor and particles at the period of acceleration of the distributor.     

Numerical simulation of particle motion in vibrated fluidized bed

The particle motion in vibrated fluidized bed was examined by using numerical simulation. A two-dimensional fluidized bed was used, and discrete element method (DEM) was used as for calculating the particle motion in vibrated fluidized bed. Geldart group B particle, which is regarded to have no effect of agglomerate, was used as the fluidizing particle. The vibration directions (vertical and horizontal) and vibration parameters (frequency and amplitude of vibration) were changed.In the case with vertical vibration, large bubbles caused by vibration gap (defined as the gap between particle bed and wall caused by vibration) appear in bed. When the vibration strength is high and the vibration frequency is low, it is difficult to fludize the particle bed in the case with horizontal vibration because the vibration gap acts like the channel. When the vibration frequency increases at the same vibration strength, the effect of vibration on the particle motion becomes small. The effect of vibration on the pressure loss in particle bed is large in the case with horizontal vibration.     

Spouting of fine powder from vertically vibrated bed

Experimental and numerical studies of spouting of powder from vertically vibrated bed are performed. The powder flows out vigorously through a side orifice of the vertically vibrated vessel in which the powder is contained, provided that the vibration acceleration is greater than the gravitational one. The mechanism of the efflux and the relationship between the efflux rate and the vibration condition are examined. The interstitial air pressure in the bed is measured and is compared to the numerical analysis considering the relative motion of the powder bed to the vibrated vessel and percolation air flow. The outflow behavior of powder is observed by using a high-speed video camera. The experimental and the numerical results show that the powder spouts out intermittently with periodic air outflow, which is generated in response to the change of the gap between the bed and the vessel base. Furthermore, it is found that the efflux rate of powder is in proportion to the generated air pressure and in inverse proportion to the vibration frequency.     

Circulation of particles in a vibrated bed with an inner tube

The particle motion in a vibrated bed with an inner tube was simulated by the discrete element method (DEM). A height difference is observed in the vibrated particle bed between the interior and annulus of the inner tube. The bed height difference is strongly affected by the ratio of the cross section at the interior to that at the annulus of the inner tube. When the inner tube is immersed in the particle bed, the bed height difference causes the circulation of particles in the bed. The direction and velocity of particle circulation can be controlled by changing the inner tube diameter and the circulation velocity is also controlled with vibration conditions.     

Mixing in a vibrated granular bed: Diffusive and convective effects

In this study, Discrete Elementary Method (DEM) is employed to simulate the motion and mixing behavior of granular materials in a three-dimensional vibrated bed, which is energized by vertical sinusoidal oscillations under different vibrating conditions. With frictional sidewalls, the convection flow is a very important phenomenon in the vibrated granular bed. The influence of vibrating conditions, including vibration acceleration and frequency, on the formation of symmetric convection flow is investigated in this study. In order to characterize the convective flow and the diffusive motion of granular materials, the dimensionless convection flow rate, J_conv, and the vertical self-diffusion coefficient, D_yy, are defined, respectively. Péclet number, Pe, is employed to characterize the ratio of the convective flow to the diffusive motion in vertical direction. The role of Pe in the formation of symmetric convection flow is discussed in detail. Moreover, the top-bottom initial loading pattern of two groups of glass beads with different colors is employed to investigate mixing behavior of granular materials. The well-known Lacey index is employed as the mixing degree, M, to quantify the mixing quality. The mixing rate is calculated from a least-square fit using the time evolution of M. The simulation results demonstrate that the mixing rates increase with increasing J_conv and D_yy in exponential relations.     

Mixing in vibrated granular beds with the effect of electrostatic force

In this study, the DEM (Discrete Elementary Method) is used to simulate the behavior of granular mixing in vibrated beds. First, the velocity fields are simulated by the DEM model to examine the convective currents in a three-dimensional vibrated granular bed. Then, in order to characterize the effect of electrostatic force on the granular flow, the electrostatic number Es is defined as the ratio of the electrostatic force to the particle weight. Also, to quantify the quality of mixing, the mixing degree M is used by the well-known Lacey index. The top–bottom and the side–side initial loading patterns of two groups of glass bead with different colors are employed to investigate the influence of the convective currents on the granular mixing. To simplify the electrostatic effect, these two groups of glass beads are given opposite charges and the charged strength is assumed to be constant. The simulation results demonstrate that the granular temperatures increase linearly with the increasing Es number. Meanwhile, the mixing rate constants, calculated from the time evolutions of mixing degree, increase with the increasing Es number in power law relations. The role of granular temperature in the granular mixing is also discussed in this paper.     

振动流化床中浅床层时床层与分布板之间的弹性作用

采用A类颗粒作为实验粉体,通过测量床层对分布板的作用力,分析了振动流化床中分布板与床层的作用机理,提出了振动流化床中在浅床层共振时分布板与床层的弹性作用模型. 实验表明,分布板与床层之间是连续的弹性作用,周期性的机械振动通过分布板与床层之间的弹性作用以波动形式在床层中传播. 该模型计算结果与实验测定结果符合的较好.     

CFD study of hydrodynamics behavior of a vibrating fluidized bed using kinetic-frictional stress model of granular flow

The hydrodynamics of a vertically vibrating fluidized bed was studied using an Eulerian-Eulerian two-fluid model (TFM) incorporating the kinetic theory of granular flow and including the frictional stress effects. Influences of frictional stresses, vibration amplitudes and frequency on behavior of the particles were studied. In the case with vertical vibration, the numerical results showed three regions of solid concentration along the bed height: a low particle concentration region near the bottom of the bed, a high concentration region in the middle of the bed, and a transition region at top of the bed. The accuracy of results was found to be closely related to the inclusion of the frictional stresses. Ability of the two-fluid model for predicting the hydrodynamics of vibrating fluidized beds was discussed and confirmed.     

Testing the validity of the spherical DEM model in simulating real granular screening processes

We investigate the flow of a granular material over a vibrated horizontal screen. We perform a direct quantitative comparison, across a range of operating conditions, between laboratory scale experiments and simulations using the discrete element method (DEM). We test the extent to which the commonly employed DEM approximation of particles being spherical affects the ability of the model to realistically reproduce the behaviour of industrial screening systems where the particles are generally non-spherical in shape.The simulation geometry and input particle size distribution are set up to exactly match the experimental system, which consists of a horizontal screen with a wire mesh cloth onto which quarry rock is fed at a series of input flow rates. The screen is vibrated, causing the granular bed to flow over the deck and vertically stratify with finer material passing through the screen, where it is collected in a series of bins located along the length of the screen. The size distribution of the material flowing through each section of the screen is found by analyzing the contents of each collection bin.The best agreement is found for very low flow rates, where the vast majority of the below aperture size material is rapidly captured just after it enters the screen in both the simulation and experiment. At higher flow rates, significant quantitative errors are found with the over-prediction of the flow rate through the screen for near grate sized particles. This is attributed to the higher rate of percolation through the bed and the easier capture by the screen surface of the spherical shaped material. The near aperture sized spherical particles also show a very strong tendency to peg the screen, becoming trapped in the screen openings and limiting further flow through those parts the screen. The use of spherical particles in the DEM simulation of vibrating screens is therefore found to be inadequate for modelling realistic flow and separation of particles that are not actually spherical.     

Movement mechanisms of solid-like and liquid-like motion states in a vibrating granular bed

In this study we analyze the granular motion and the variation in energy on the surface of a vibrating bed as well as the flow regime inside the bed. The variation in the granular kinetic energy at the free surface of the bed under different vibration conditions produces the following motion states: solid-like state, transition region and liquid-like state. The granular kinetic energy is measured and found to increase slowly until the vibration acceleration reaches the first specific dimensionless vibration acceleration Γ_I. The granular kinetic energy increases sharply as the vibration acceleration continues to increase after exceeding Γ_I until reaching the second specific dimensionless vibration acceleration Γ_II. The granular kinetic energy increases with the vibration acceleration afterwards. The movement mechanism inside the granular bed is transformed from particle reorganization to surface convection as the bed state transfers from a solid-like state to a liquid-like state.