颗粒相流动边界条件对其浓度分布的影响机制

运用描述颗粒相内部作用力的离散介质动力理论和“双流体”欧拉方法对循环流化床中的稠密气固两相流动进行了数值模拟,并专门研究了颗粒边界条件对其在近壁区浓度分布的影响,探讨了产生畸变的机制,比较和检验了相关的理论模型。研究表明,近壁区颗粒轴向速度梯度大小,和在那里局部加密的计算网格,对颗粒浓度的径向分布有强烈的影响,导致浓度分布畸变的原因是不合理的边界条件过高估计了当地颗粒的脉动强度,提出选择边界条件的原则是允许颗粒在壁面上有一定程度的滑移。     

离散元法铁粉末压制中粒径分布对力链演化机制的影响

为阐明粒径分布对铁粉压制中体系内部细观力学行为的影响, 基于离散元理论, 通过改变铁粉颗粒粒径分布建立压制模型, 结合力链提取方法, 通过对力链空间分布、力链数目、力链长度和力链方向性的分析, 探究粒径分布对力链演化的影响机理. 研究结果表明: 不同粒径分布的粉体压制时形成的力链空间分布具有差异, 粒径分布范围越小, 形成的力链分布越集中, 反之, 粒径分布范围越大, 形成的力链分布越松散且均匀; 在粉末压制时, 粒径分布对力链数目也有影响, 具体表现为随着粉体的粒径分布范围变大, 力链总数逐渐减少; 粉体的粒径分布对颗粒形成短力链的数目起着显著影响, 而对力链长度的影响较为有限; 随着粒径分布范围的增大, 力链的方向由均匀分布逐渐集中在特定角度方向, 表现出一定各向异性, 形成的交叉力链网络结构有利于提高粉体致密化程度. 本文为从粉体粒径分布影响层面拓展粉末压制细观力学理论提供基础, 亦为进一步结合粉体粒径分布及体系内力链演变过程改善粉末致密化行为提供指导.      

离散颗粒流动堆积行为离散元模拟及实验研究

工程应用中存在许多颗粒流动堆积问题.首先设计了一系列测量方法,通过大量的实验和统计分析,得到了颗粒的多种物理参数.并以工程中高炉炉顶称量料罐为背景,采用离散元方法模拟不同物理条件下离散颗粒的流动堆积行为,得出料罐内颗粒系统中颗粒之间力的分布不均匀,而且强力链分布主要与料罐的左下壁方向平行.同时,设计并制作了具有多参数调节的离散颗粒料罐实验模型,进行了相应的物理实验,实测结果与数值模拟吻合良好.     

钟摆状态下湿颗粒坍塌流动行为及其动力相似性

基于建立的湿颗粒离散动力学模型,本文系统研究了钟摆状态下湿颗粒柱在重力驱动下的坍塌流动过程,主要考虑了颗粒粒径、液体表面张力系数和液体含量等参数对系统坍塌流动模式和动力学行为的影响。研究发现,在湿颗粒系统中,颗粒粒径和液体表面张力系数会改变颗粒间的毛细力大小,引起系统发生不同的坍塌流动模式,而液体含量仅定量影响颗粒坍塌后的堆积形态。在此基础上,进一步探讨了不同模式下系统坍塌流动行为与模型参数的相关性,发现无量纲Bond数是决定钟摆状态下湿颗粒物质坍塌流动动力学行为的本质因素。     

颗粒物质在竖直管中下行蠕动流动的实验测量

颗粒流蠕动行为是颗粒物质在竖直管中流动时常见的一种流动现象,其产生机理较复杂。为此,本文在在内径为150mm、高为5000mm的竖直管实验装置上,以FCC催化剂为固体颗粒物料,采用PV6型颗粒速度测量仪,测量不同颗粒流率下竖直管中的颗粒下行蠕动流动速度以及颗粒固含率,系统地考察了颗粒物质在竖直管中下行流动时的蠕动流动特性及产生机理。实验结果表明,颗粒物质在竖直管中下行流动时的流动行为可划分为两种形式。在颗粒流率较小时,颗粒物质下行速度呈现脉冲式变化,有速度停滞,可称之为蠕动I型流动。随着颗粒流率的增加,颗粒下行速度停顿消失,但仍是起伏变化,为蠕动II型流动。当颗粒流率增大到一定值后,颗粒物质下行蠕动行为消失,转变为流化流动。颗粒物质下行的蠕动行为是出口区颗粒成拱与崩塌、颗粒与器壁滑动摩擦和颗粒力链作用的综合反映。     

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

卡鲁塞尔氧化沟反应器三维流场体视PIV测量

采用体视PIV(Stereoscopic PIV)技术对卡鲁塞尔氧化沟模型直道、弯道及曝气叶轮段等处三维全场流速进行测量,克服了传统接触式点测量方法无法获取全场流动同步信息的缺陷,分析了不同位置处的流动结构,并对纵、横、垂三向流速沿程分布特点进行了研究和比较.结果表明,纵、垂两向的流动分布是决定沟内水力特性的主要因素;横、垂两向的流动是决定污泥沉积位置的主要因素;外沟靠近曝气叶轮直道段的流速分布上大下小,使低速区底部易发生污泥沉积;外沟远离曝气叶轮直道段流速分布上小下大,利于防止泥水分离;弯道段受横比降和横向环流的影响使内侧容易形成低速区或停滞区而发生污泥沉积;倒伞型叶轮使其附近流场在垂向呈螺旋形分布,利于气、液、固三相的均匀混合.     

明渠湍流边界层中颗粒的运动与分布

利用直接数值模拟、点球浸入边界法和颗粒离散元法相结合的方法, 模拟了颗粒在明渠湍流边界层中的运动, 并对颗粒的瞬时位置进行了Voronoi 分析, 定量研究了颗粒在湍流边界层中的运动和分布规律. 研究发现:颗粒的输运对湍流的统计特征有影响, 其运动与近壁区湍流拟序结构密切相关, 在"喷发"结构作用下被带离壁面, 在"扫掠" 结构和自身重力作用下回到壁面; 在湍流边界层中, 颗粒倾向于聚集在低流速带, 呈条带状分布;颗粒在大部分时间处于"簇"状态, 偶尔跳跃到"空" 状态, 但能够很快回到邻近低速区域.      

纳米W粉冲击烧结的分子动力学模拟

粉末冲击烧结是制备高品质W的一种有效方法，而分子动力学方法在尺度极小、过程迅速的数值模拟上有着独特的优势。因此运用分子动力学方法，结合W的嵌入原子势，对常温下的纳米W粉末的冲击烧结过程进行模拟，得到颗粒微观压实过程图、体系速度分布云图、p-U_p、T-U_p、T-p曲线以及径向分布函数。研究了不同颗粒速度及产生的射流对纳米W粉末冲击烧结影响，分析了微观冲击烧结机理。结果表明，低速冲击条件下（500 m/s以下），纳米颗粒无法压实。高速条件下（1 000 m/s及以上），颗粒能获得致密化很高的压实。颗粒间的相互挤压造成的高应力使颗粒表面的原子发生流动变形，原子向颗粒间空隙流动，形成压实。颗粒间产生的射流以及高速冲击导致的颗粒熔化，均促进烧结获得致密度更高的烧结体。     

矸石山剖面颗粒分布规律的实验研究

为研究矸石山内部剖面上不同粒径矸石的分布规律,首先利用因次分析方法,将影响矸石山堆积过程的10个物理力学参数简化为3个无因次量,并依此结果确定实验的相似判据与相似比,进而采用相似模拟方法对矸石山的整个堆积过程进行还原与再现。实验结果表明:矸石山顶部多分布为小粒径矸石,底部多分布大粒径矸石,中部区域由小粒径与中等粒径、中等粒径与大粒径交替发生尖灭而形成倾斜的锯齿状条纹结构;每条条纹存在2个临界角:一个是矸石颗粒能保持稳定的角度,约为24°～26°;另一个是矸石颗粒停止下滑的角度,约为34°～36°。     

Mixing and segregation of solid particles in a conical spouted bed: Effect of particle size and density

In this work, the mixing and segregation of binary mixtures of particles with different sizes and densities in a pseudo-2D spouted bed were studied experimentally. A binary mixture of solid particles including sand, gypsum, and polyurethane was used. To determine the particles mass fraction, and their mixing and segregation in the bed, an image-processing technique was developed and used. Important hydrodynamic parameters, such as the axial and radial segregation profiles of the solid particles, were measured. The effects of air velocity, particle size, and particle mass fraction were also evaluated. The flow regime in the spouted bed and the time required for reaching the equilibrium state of the solid particles were discussed. The results showed that the segregation of solid particles and the time to equilibrium both decreased when the air velocity increased to much larger than the minimum spouting velocity. The axial segregation increased with the diameter ratio of the particles. Upon completion of the test, coarse particles were concentrated mainly in the spout region, while fine particles were aggregated in the annulus region. Examination of the flow pattern in the spouted bed showed that the particles near the wall had longer flow paths, while those near the spout region had shorter flow paths.     

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.     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况，本文以颗粒长径比为2，颗粒之间的初始距离ΔS_Py=4D为例，基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布，在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下，研究了颗粒的受力和运动的变化，探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明，双柱状颗粒在线性剪切流场中易向速度大的流体区域运动；前后两颗粒运动状态和轨迹不同，颗粒之间距离较近时，曳力会产生较大的波动；只有当颗粒在壁面附近时，滞后颗粒才能追上领先颗粒，两颗粒发生牵引、翻滚和分离过程。     

Experimental and computational studies on flow behavior of gas-solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally.Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695--0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

回转筒体轴向窜动的力学性能分析及控制

基于长期生产实践，针对回转窑筒体在运转时发生持续下滑，引起系统频发故障，造成停机的问题，应用弹性滑动理论和摩擦传动理论分析论证了窑体下滑的原因；给出了窑体轴向滑动速度的计算方法；介绍了控制窑体轴向窜动的方法。     

Multi-fluid modeling of density segregation in a dense binary fluidized bed

This paper presents simulation results of the density segregation in a dense binary gas fluidized bed using a multi-fluid model from Chao et al. (2011). The segregation behavior of two types of particles with approximately same particle diameters and different particle densities was studied and validated using the experimental data from Formisani et al. (2008). Some detailed information regarding the gas, particle velocity profiles, the distributions of the particle volume fractions and the flotsam-to-total particle volume fraction ratios is presented. The simulation results show that the simulated axial average flotsam-to-total particle volume fraction ratio distribution agrees reasonably with the experimental data of Formisani et al. (2008). The binary particle velocities are closely coupled though the segregation exists. The segregation behavior and the particle velocity profiles are superficial gas velocity dependent. The number and distribution of particle velocity vortices change dramatically with superficial gas velocity: at a comparatively low superficial gas velocity, the particles mainly segregate axially, and at a comparatively high superficial gas velocity, the particles segregate both axially and radially.     

Local solids flow structure in a countercurrent liquid-upward and solids-downward fluidized bed (CCLSFB) system

The local solids holdup and local particle velocity in a Countercurrent Liquid-upward and Solids-downward Fluidized Bed (CCLSFB) were investigated in details using optical fiber probes with two different models in a Plexiglas column of 1.5 m in height and 7.0 cm in inner diameter. A new flow regime map including fluidized bed, transition, and flooding regimes was established. The axial solids holdup distribution is almost uniform at low liquid velocity and/or solids flowrate and becomes less uniform with higher solids holdup at the top of the column after the operating liquid velocity is reaching the flooding velocity. The radial solids holdup profile is also nearly flat with a slightly lower solids holdup in the near-wall region at low liquid velocity and solids flowrate but becomes nonuniform as the operating liquid velocity approaches the flooding velocity. Two equations were also proposed to correlate radial local solids holdups. The descending particle velocity in CCLSFB increases with the decrease of the liquid velocity and the increase of the solids flowrate. A generally uniform particle velocity distribution was found in the axial direction, as well as in the radial direction except for a small decrease near the wall. These results on the local solids flow structure would provide basic information and theoretical supports for the design and industrial application of CCLSFB.     

Investigation on gas-solids hydrodynamics and matching enhancement in feedstock injection zone of double-level nozzle riser

With the development of current energy economy, it is necessary to improve the product distribution of fluid catalytic cracking process, which is achieved by a riser reactor with double-level of nozzles. The new riser is constructed by adding a level of secondary nozzle 0.5 m below the main nozzle of traditional riser. This paper investigates the gas-solids flow and oil-catalyst matching feature based on the optical fiber and tracer technologies. According to the distribution of solids holdup, particle velocity and dimensionless jet concentration, the feedstock injection zone can be divided into the upstream flow control region, the main flow control region, and the secondary flow control region in the radial direction. The size of the regions is changed by the jet gas velocity and axial height. There is a poor match of secondary nozzle jet to particles below the main nozzle. The jet gas from secondary nozzles can improve the matching effect of oil-catalyst near the wall and reduce the probability of coking above the main nozzle.