气-固流化床中颗粒的内循环流动

从多尺度范围考察了气－固流化床内尾涡颗粒流和乳化相颗粒流的运动规律,将分散的尾涡颗粒流和乳化相颗粒流连续介质化,从流变学角度定义了颗粒流粘度,用流体力学方法建立了内循环流动结构的多尺度、连续介质流模型,较好地揭示了颗粒循环流动的规律。实验观测支持模型预测结果。     

从质量流向漏斗流转变过程中的动力学分析

球床模块式高温气冷堆的堆芯是全陶瓷型包覆铀燃料制成的球形颗粒,与石墨颗粒混合堆积而成,堆芯颗粒流的流态取决于颗粒尺度的平移、旋转等动力学量,以及力链、涡旋等介尺度物理量。为了分析颗粒的平移、旋转等动力学量对颗粒流流态的影响。基于筒仓颗粒流的物理模型,首先开展了筒仓颗粒流流变过程的实验测量,并使用基于 Hertz-Mindlin和 RVD （relative velocity dependent）滚动摩擦接触模型的离散单元法 （distinct element method, DEM）,研究了锥形筒仓颗粒流流变过程中球形颗粒的动力学量。进一步,基于DEM计算结果进行分析,发现筒仓自上而下呈现出质量流向漏斗流过渡的混合流状态。在筒仓混合流的不同流型区域中,平移速度和旋转速度之间的相关性是相反的;颗粒间的相对切向运动较大的区域集中在漏斗流区域与边壁区域。了解筒仓流变过程中颗粒的动力学特征,有助于优化筒仓颗粒流动,并减少颗粒表面的磨损。     

基于EMMS范式的离散模拟及其化工应用

化工过程通常涉及化学、化工、过程系统工程3个层次,而每个层次又包含微尺度、介尺度和宏尺度,如化工层次的颗粒、颗粒团和反应器尺度。每个层次中的微尺度单元都自然适用离散模型,即通过跟踪每个单元的运动获得整个体系演化的宏观规律。但由于单元数量巨大,工程模拟往往依赖经过统计平均的连续介质模型。由此带来的精度问题,特别是忽略了介尺度结构的问题,随着对化工过程效率和绿色度等要求的提高而日渐突出。介绍了通过问题、模型、软件和硬件结构的一致性提升离散模拟的精度、能力和效率的方法、进展及其在复杂分子体系、颗粒流、气固流态化等方面的应用,展示了通过离散模拟实现虚拟过程工程的可能性。     

偏心楔形喂料斗卸料过程中颗粒流动特性

采用三维离散单元法,研究了偏心楔形喂料斗中不同粒径颗粒卸料过程的流场分布,建立了适用于偏心楔形喂料斗的整体流系数模型,分析了料斗卸料流型以及卸料流型与卸料质量流率的相关性,并通过实验验证了离散元模型的可靠性。结果表明：偏心楔形喂料斗内颗粒流场分布以靠近喂料管口区域为高速区,并呈辐射状朝远离垂直壁面端的料斗上部低速区过渡,高速区颗粒流场的整体一致性好,低速区颗粒流场呈局部涡流状;颗粒粒径增大,颗粒整体流动性变差,高速区域范围减小,低速区域范围增大,过渡区域变模糊;当颗粒粒径不大于10 mm时,整体流系数与颗粒粒径呈线性负相关,卸料质量流率与整体流系数呈线性正相关。     

气—固流化床颗粒的内循环流动

从多尺度范围考察了气－固流化床内尾涡颗粒流和乳化相颗粒流的运动规律，将分散的尾涡颗粒流和乳化相颗粒流连续介质化，从流变学角度定义了颗粒洗粘度，用流体力学方法建立了内循环流动结构的多尺度，连续介质流模型，较好地揭示了颗粒循环循环流动的规律。实验观测支持模型预测结果。     

基于颗粒尺度的离散颗粒传热模型

颗粒间传热在诸多工业过程中有着十分重要的作用。详细考虑颗粒间传热机理,对颗粒间各传热途径建模,包括颗粒内部导热、颗粒粗糙表面传热、颗粒表面气膜及接触颗粒间隙气膜传热,并与离散颗粒模型(DEM)耦合,建立颗粒尺度下离散颗粒传热模型。以固定床为对象,考察颗粒粒径、颗粒比热容、颗粒热导率及压缩负载对固定床有效传热系数的影响,并将本文计算值和文献的实验值及模型预测值对比,结果表明,该模型可定量预测固定床有效传热系数。本文建立的离散颗粒传热模型为合理预测颗粒体系内的传热提供了一种有效方法。     

提升管与下行床颗粒团聚行为的离散颗粒模拟

从微观机理出发,采用计算流体力学和离散单元方法（CFD-DEM）结合的模型对二维提升管和下行床气固流动体系进行了数值模拟。模拟选用了粒径为520 μm、密度为2620 kg·m-3的球形颗粒和周期性边界条件,展示了气固并流逆重力场和顺重力场运动的颗粒聚团瞬态图像,定性或半定量地揭示了两个不同体系的颗粒微观聚集行为。提升管中颗粒聚团较为严重,且表现明显的颗粒返混现象;下行床中的颗粒聚团比较松散,且具有与宏观流动相同的流速方向,几乎无颗粒返混。通过统计分析获得宏观时均流体力学行为,包括两相的相分布和速度分布,并与文献报道的实验现象进行定性的比较。     

单颗粒褐煤高温烟气干燥过程数值模拟

褐煤干燥对于提高其品质具有重要意义。为了模拟高温烟气干燥这一高温差、变温差非稳态传热传质过程中褐煤内水分蒸发过程,采用有限体积法建立了一维球坐标系下蒸发界面向内迁移的单颗粒褐煤干燥数学模型,并利用该模型分析了初始烟气温度和颗粒粒径对单个褐煤颗粒干燥特性的影响。模型模拟结果与实验结果对比表明二者变化趋势一致,所建模型能较好地反映出高温烟气干燥过程中褐煤内水分蒸发过程。结果表明,初始烟气温度越高,颗粒粒径越小,蒸发界面向内迁移速度越快,水分脱除越快,干燥时间越短;蒸发界面平均迁移速度均与初始烟气温度和颗粒粒径呈线性关系;在初始烟气温度700℃下,较短的停留时间使得颗粒表面温度未达到挥发分析出温度,本研究中不同粒径褐煤颗粒在干燥过程中基本没有挥发分的析出。     

颗粒形状对包衣设备内药片颗粒运动特性影响的数值模拟

为研究颗粒形状对包衣设备内药片颗粒运动特性的影响,基于离散单元法及自行编写的喷雾区颗粒检测算法,采用数值模拟的方法对五种不同形状(棒状、长椭球、扁椭球、双凸形和球形)的药片颗粒在包衣设备内的运动行为规律进行了研究。分析了颗粒形状对颗粒系统的能量、床面颗粒平动速度、颗粒温度及颗粒流在喷雾区域停留时间的分布及其相对标准差的影响。结果表明,颗粒形状对颗粒的平均动能、颗粒床面速度、颗粒温度、喷雾区域停留时间分布及颗粒间包衣均匀性有重要影响。除双凸形颗粒系统外,对于其他四种形状的颗粒系统,随着颗粒球形度增大,颗粒系统具有的动能、床面速度和颗粒温度均呈减小趋势。除棒状颗粒系统外,对于其余四种形状的颗粒系统,随着颗粒球形度增大,颗粒系统在包衣喷雾区域内平均停留时间减小,平均停留时间的相对标准差增大,包衣均匀性变差。与球形颗粒系统相比,非球形颗粒系统的包衣均匀性更好;药片颗粒形状对包衣设备内颗粒运动特性及颗粒之间的包衣均匀性有重要影响。     

气液交叉流自增湿促进PM2.5颗粒长大模型

研究气液交叉流阵列内传热传质自增湿形成过饱和条件下水蒸汽在PM_(2.5)颗粒表面经历异质核化凝结长大的机理,并建立颗粒长大计算模型,逐排计算,确定颗粒粒径变化量。模型分析与计算结果解释了自增湿由强到弱的阵列条件下,实验检测到的PM_(2.5)颗粒粒径增长速度先增大后减小的规律。在气体温度T_(in)=41.2℃,饱和度S_R=0.85,横掠温度T_w=5℃的一个单元(100排)降膜交叉阵列实验条件下,最大可达5.5μm/s的增长速度,且模型预测颗粒长大后粒径分布与实验结果吻合较好。     

Silo music and silo quake: granular flow-induced vibration

Acceleration and sound measurements during granular discharge from silos are used to show that silo music is a sound resonance produced by silo quake. In tall and narrow silos, the latter is produced by stick-slip friction between the wall and the granular material. For the discharge rates studied, the occurrence of flow pulsations is determined primarily by the surface properties of the granular material and the silo wall. The measurements show that the pulsating motion of the granular material drives the oscillatory motion of the silo.     

Prediction of particle motion in a two-dimensional bubbling fluidized bed using discrete hard-sphere model

The solids motion in a gas-solid fluidized bed was investigated using a discrete hard-sphere model. Detailed collision between particles and a nearest list method are presented. The turbulent viscosity of gas phase was predicted by subgrid scale (SGS) model. The interaction between gas and particles phases was governed by Newton's third law. The distributions of concentration, velocity and granular temperature of particles are obtained. The radial distribution function is calculated from the simulated spatio-temporal particle distribution. The normal and shear stresses of particles are predicted from the simulated instantaneous particle velocity. The pressure and viscosity of particles are obtained from both the kinetic theory of granular flow and the calculated stresses of particles. For elastic particles the individual lateral and vertical particle velocity distribution functions are isotropic and Maxwellian. The observed anisotropy becomes more pronounced with increasing degree of inelasticity of the particles.     

Numerical analysis of silo behavior using non-coaxial models

Granular solids in silos experience considerable principal stress rotations, which result in the non-coaxiality between principal stresses and plastic strain rates. This paper discusses the influences of the use of elastoplastic non-coaxial models for granular solids on predictions of wall pressure distributions in silos by using the finite element method. A well established non-coaxial model in geomechanics, the yield vertex model, is employed. Simulations are performed on a steep hopper characterized with a mass flow and a flat-bottomed silo with a semi-mass flow. The simulations indicate that the non-coaxiality does not influence predictions of wall pressures after filling. On the other hand, the predicted discharge wall pressures with non-coaxial considerations are larger than those without it. Its mechanism is discussed in this paper. The suppressed shear-dilatancy of granular solids in silos leads to a larger increase of normal stress with non-coaxial models.     

Application of particle image velocimetry (PIV) for deformation measurement during granular silo flow

The paper presents results of deformation measurements in dry cohesionless sand during free flow in small rectangular model silos using a non-invasive, indirect method called particle image velocimetry (PIV). It is an optical technique for measuring surface deformations from successive digital images. Laboratory model tests were performed with a mass and funnel flow silo to investigate the kinematics of the flowing sand. The measurements were carried out for granular flow in model silos without inserts and in a funnel flow silo equipped with three different types of inserts: cone-in-cone, inverted cone and double cone. The effect of the initial sand density and roughness of silo walls on the volumetric and deviatoric strain in sand was investigated. The accuracy of measurements was discussed. Advantages and disadvantages of PIV were outlined. The results were qualitatively compared with sand displacements obtained with colored sand layers and with the aid of X-rays (initially dense sand).     

A second-order moment method of dense gas-solid flow for bubbling fluidization

A gas-solid two-fluid model with the second-order moment method is presented to close the set of equations applied to fluidization. With the kinetic theory of granular flow, transport equations for the velocity moments are derived for the particle phase. Closure equations for the third-order moments of velocity and for the fluid-particle velocity correlation are presented. The former is based on a modified model with the contribution of the increase of the binary collision probability, and the latter uses an algebraic model proposed by Koch and Sangani [1999. Particle pressure and marginal stability limits for a homogeneous monodisperse gas-fluidized bed: kinetic theory and numerical simulations. Journal of Fluid Mechanics 400, 229-263]. Boundary conditions for the set of equations describing flow of particles proposed by Strumendo and Canu [2002. Method of moments for the dilute granular flow of inelastic spheres. Physical Review E 66, 041304/1-041304/20] are modified with the consideration of the momentum exchange by collisions between the wall and particles. Flow behavior of gas and particles is performed by means of gas-solid two-fluid model with the second-order moment model of particles in the bubbling fluidized bed. The distributions of velocity and moments of particles are predicted in the bubbling fluidized bed. Predictions are compared with experimental data measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168]. in the bubbling fluidized beds. The simulated second-order moment in the vertical direction is 1.1-2.5 [Muller, C.R., Holland, D.J., Sedeman, A.J., Scott, S.A., Dennis, J.S., Gladden, L.F., 2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and 1.1-4.0 [Yuu, S., Umekage, T., Johno, Y., 2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168] times larger than that in the lateral direction because of higher velocity fluctuations for particles in the bubble fluidized bed. The bubblelike Reynolds normal stresses per unit bulk density used by Gidaspow et al. [2004. Hydrodynamics of fluidization using kinetic theory: an emerging paradigm 2002 Flour-Daniel lecture. Powder Technology 148, 123-141.] are computed from the simulated hydrodynamic velocities. The predictions are in agreement with experimental second-order moments measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and fluctuating velocity of particles measured by Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168].     

Flow Channel Boundaries in Silos

The critical slip planes at the silo filling state are compared with the flow channel boundary during silo discharge for semi‐mass flows. The static critical slip planes are determined by using the dynamic programming method based on the stress field of granular solids stored in silos at the filling state. The flow channel boundary is estimated through the finite element analysis of the silo discharge. The results indicate that the critical slip line lies above the flow channel boundary. This characteristic can be attributed to the changeover of major principal stress directions of granular solids from the silo filling to the silo discharge. The analysis demonstrates that the silo wall friction tends to lift up the critical slip plane and flow boundary. A simple correlation is developed between the positions of critical slip planes and flow boundaries and is experimentally verified.     

Silo music — Mechanism of dynamic flow and structure interaction

This papers deals with the strong dynamic effects (called silo music) appearing during confined granular flow in the cylindrical silos. Silo experiments with dry cohesionless sand during gravitational outflow were performed in a cylindrical perspex model silo. During tests the wall accelerations and acoustic signals were recorded and the mode shapes of the silo structure were determined. In addition, experiments were performed with additional modifications of the silo structure and silo flow. A novel simple mechanism of the origin of silo music was proposed.     

Shape dependence of resistance force exerted on an obstacle placed in a gravity‐driven granular silo flow

Resistance force exerted on an obstacle in a gravity‐driven slow granular silo flow is studied by experiments and numerical simulations. In a two‐dimensional granular silo, an obstacle is placed just above the exit. Then, steady discharge flow is made and its flow rate can be controlled by the width of exit and the position of obstacle. During the discharge of particles, flow rate and resistance force exerting on the obstacle are measured. Using the obtained data, a dimensionless number characterizing the force balance in granular flow is defined by the relation between the discharge flow rate and resistance‐force decreasing rate. The dimensionless number is independent of flow rate. Rather, we find the weak shape dependence of the dimensionless number. This tendency is a unique feature for the resistance force in granular silo flow. It characterizes the effective flow width interacting with the obstacle in granular silo flow. © 2018 American Institute of Chemical Engineers AIChE J, 64: 3849–3856, 2018