Simulation study of the effect of wall roughness on the dynamics of granular flows in rotating semicylindrical chutes

A discrete element model (DEM) is used to investigate the behavior of spherical particles flowing down a semicylindrical rotating chute. The DEM simulations are validated by comparing with particle tracking velocimetry results of spherical glass particles flowing through a smooth semicylindrical chute at different rotation rates of the chute. The DEM model predictions agree well with experimental results of surface velocity and particle bed height evolution. The validated DEM model is used to investigate the influence of chute roughness on the flow behavior of monodisperse granular particles in rotating chutes. To emulate different base roughnesses, a rough base is constructed out of a square close packing of fixed spherical particles with a diameter equal to, smaller, or larger than the flowing particles. Finally, the DEM model is used to study segregation in a binary density mixture for different degrees of roughness of the chute. © 2015 American Institute of Chemical Engineers AIChE J, 61: 2117–2135, 2015     

Bulk coarse particle arching phenomena in a moving bed with fine particle presence

Bridging or arching of flowing solids particles is a serious hazard in the operation of moving bed systems. The mechanics of the arching has been extensively analyzed in the context of particle discharge from a hopper with conical geometry by considering the particulate layer stress distribution. However, bridging can also occur in a moving bed system with cylindrical geometry during the continuous mass flow of solids particles. Experimental work conducted in this study reveals that the appearance of solids bridging is normally accompanied by the presence of fine particles in the coarse moving particles as well as by the countercurrent interstitial gas flow. In this study, a stress analysis of the layered particles distributed in a cylindrical, vertical moving bed that flows downward opposing to upward flow of the interstitial gas is developed to quantify the bridging phenomenon. The analysis takes into account of the effects of presence of fine powder in the coarse particle flows and properties, such as particle‐size distribution, bed voidage, and interstitial gas flow rate. The experimental validation of the present stress analysis for moving bed systems with varied fine and coarse particle concentration distributions, and interstitial gas velocities is also conducted. The stress distributions of the particles under flowing and arching conditions are obtained. An arching criterion is formulated, which indicates that the critical radius of the standpipe to avoid arching phenomenon is only related to the property of the bulk solids in the present geometric configuration of the flow system. © 2014 American Institute of Chemical Engineers AIChE J, 60: 881–892, 2014     

Bubble dynamics in a 3‐D gas–solid fluidized bed using ultrafast electron beam X‐ray tomography and two‐fluid model

Bubble characteristics in a three‐dimension gas‐fluidized bed (FB) have been measured using noninvasive ultrafast electron beam X‐ray tomography. The measurements are compared with predictions by a two‐fluid model (TFM) based on kinetic theory of granular flow. The effect of bed material (glass, alumina, and low linear density polyethylene (LLDPE), d_p ～1 mm), inlet gas velocity, and initial particle bed height on the bubble behavior is investigated in a cylindrical column of 0.1‐m diameter. The bubble rise velocity is determined by cross correlation of images from dual horizontal planes. The bubble characteristics depend highly upon the particle collisional properties. The bubble sizes obtained from experiments and simulations show good agreement. The LLDPE particles show high gas hold‐up and higher bubble rise velocity than predicted on basis of literature correlations. The bed expansion is relatively high for LLDPE particles. The X‐ray tomography and TFM results provide in‐depth understanding of bubble behavior in FBs containing different granular material types. © 2014 American Institute of Chemical Engineers AIChE J, 60: 1632–1644, 2014     

颗粒斜槽流的实验研究与数值模拟

从实验和理论两个方面对颗粒慢速斜槽流进行了研究。建立颗粒慢速斜槽流的实验装置 ,采用示踪颗粒法测定表面速度 ;通过测量表面速度和流层厚度 ,初步分析了流率及壁面状况对流动的影响。用有限元法对粗糙表面上的慢速斜槽流进行了数值模拟     

Effect of particle size on flow and mixing in a bladed granular mixer

A number of studies have modeled flow and mixing of granular materials using the discrete element method (DEM). In an attempt to reduce computational costs, many of these DEM studies model particles larger than the actual particle size without investigating the implications of this assumption. Using DEM, the influence of the modeled particle size on flow and mixing in a bladed granular mixer is studied. The predicted flow microdynamics, including mixing rates, are strongly dependent on the particle diameter. The effect of particle size on macroscopic advective flow also is significant, particularly for dilute flow regions. These results suggest that the influence of particle size needs to be taken into consideration when using larger particles in DEM mixing simulations. To guide scale‐up efforts, particle‐size‐based scaling relationships for several key flow measurements are presented. © 2014 American Institute of Chemical Engineers AIChE J, 61: 46–57, 2015     

颗粒流拟流体的本构关系

针对慢速密集颗粒流的特点,采用拟流体的方法,类比具有屈服应力的非牛顿流体,给出了适用于慢速密集颗粒流的本构方程表达形式,建立了用粗糙底面斜槽测量颗粒流当量粘性的数学模型,并以小麦颗粒为例,通过粗糙底面斜槽实验得出了小麦颗粒流的屈服应力及其本构方程的关键系数. 将实验得出的小麦颗粒流本构关系系数作为其当量粘性的特征值对斜面流进行解析计算,体积流量的模型预测值与实验值的相对误差在15%以内.     

Direct numerical simulations of dynamic gas‐solid suspensions

Direct numerical simulation results for gas flow through dynamic suspensions of spherical particles is reported. The simulations are performed using an immersed boundary method, with careful correction for the grid resolution effect. The flow systems we have studied vary with mean flow Reynolds number, solids volume fraction, as well as particle/gas density ratio. On the basis of the simulation results, the effect of particle mobility on the gas‐solid drag force is analyzed and introduced into the existing drag correlation that was derived from simulations of stationary particles. This mobility effect is characterized by the granular temperature, which is a result of the particle velocity fluctuation. The modified drag correlation is considered so‐far the most accurate expression for the interphase momentum exchange in computational fluid dynamics models, in which the gas‐solid interactions are not directly resolved. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1958–1969, 2016     

大涡模拟方法模拟鼓泡床内气固两相流动

采用大涡模拟（LES）方法模拟气相湍流,颗粒动理学方法考虑颗粒相碰撞产生的动量和能量传递和耗散,采用颗粒相大涡模拟方法（LESp）模拟颗粒脉动导致的能量耗散,同时考虑介观尺度对颗粒相压力的影响,建立了气体-颗粒LES-θ-LESp双流体模型,研究鼓泡流化床内气固两相流动的特性。数值模拟与文献实测颗粒速度和实测颗粒浓度结果具有相同的变化趋势。     

Simulation of particles and gas flow behavior in a riser using a filtered two-fluid model

Flow behavior of gas and particles is predicted by a filtered two-fluid model by taking into the effect of particle clustering on the interphase momentum-transfer account. The filtered gas–solid two-fluid model is proposed on the basis of the kinetic theory of granular flow. The subgrid closures for the solid pressure and drag coefficient (Andrews et al., 2005) and the solid viscosity (Riber et al., 2009) are used in the filtered two-fluid model. The model predicts the heterogeneous particle flow structure, and the distributions of gas and particle velocities and turbulent intensities. Simulated solids concentration and mass fluxes are in agreement with experimental data. Predicted effective solid phase viscosity and pressure increase with the increase of model constant c_g and c_s. At the low concentration of particles, simulations indicate that the anisotropy is obvious in the riser. Simulations show the subgrid closures for viscosity of gas phase and solid phase led to a qualitative change in the simulation results.     

Investigating the dynamics of cylindrical particles in a rotating drum using multiple radioactive particle tracking

The behavior of granular flows inside rotating drums is an ongoing area of research. Only a few studies have investigated non‐spherical particles despite the fact that particle shape is known to have a significant impact on flow behavior. In addition, the experimental techniques limit the interpretation of the results of these studies. In this work, we compared the flow behavior of cylindrical and spherical particles using the multiple radioactive particle tracking technique to capture the positions and orientations of cylindrical particles simultaneously. We analyzed two important components of the transverse flow dynamics, that is, the boundary between the active and passive layers, and the velocity profile on the free surface. For the cylindrical particles, two general models are proposed to calculate the velocity profiles on the free surface and the effective particle sizes in the active and passive layers. © 2016 American Institute of Chemical Engineers AIChE J, 62: 2622–2634, 2016     

Simulation of the granular flow of cylindrical particles in the rotating drum

This study aims at unveiling the effect of particle shape on granular flow behavior. Discrete element method is used to simulate cylindrical particles with different aspect ratios in the rotating drum operating in the rolling regime. The results demonstrate that the cylindrical particles exhibit similar general flow patterns as the spherical particles. As the aspect ratio of the cylindrical particles increases, the active‐passive interfaces become steeper, and the number fraction, solid residence time, and collision force in the active region decreases. The mechanism underlying the difference is the preferential orientation, with particles of greater aspect ratios increasingly orientating their longitudinal axes perpendicular to the drum length. Also, particle alignment in the active region is more uniform than that in the passive region. The results obtained in this work provide new insights regarding the impact of particle shape on granular flow in the rotating drum. © 2018 American Institute of Chemical Engineers AIChE J, 64: 3835–3848, 2018     

Rapid Gravity Flow of a Granular Medium

An experimental study is performed to analyze the equation of state of a granular medium, which relates its dilatancy, shear rate, and pressure by an analogy between a granular medium under rapid shear and a dense gas. A method is developed for the contactless determination of the solid-phase concentration distribution in the two-dimensional flow of a granular medium using x-ray analysis. The method is used to explore the concentration profiles in rapid gravity flows of ceramic particles down a rough incline. The adequacy of the equation of state of a granular medium under rapid shear and for unordinary properties of the gravity flow down a rough incline is directly confirmed.     

粗糙颗粒动理学及流化床内气固流动的数值模拟

基于气体分子动理学和颗粒动理学理论,考虑颗粒旋转流动对颗粒碰撞能量交换和耗散的影响,建立粗糙颗粒动理学。采用Chapman-Enskog颗粒速度分布函数,提出了颗粒相应力、热通量和颗粒碰撞能量耗散计算模型。采用欧拉-欧拉气固双相流模型,数值模拟鼓泡流化床内气体-颗粒两相流动特性。模拟结果得到了床内颗粒相速度和脉动速度分布,与Yuu等实验结果相吻合。分析不同的切向弹性恢复系数对颗粒相拟总温的变化规律,结果表明在低颗粒浓度时颗粒拟总温随切向弹性恢复系数而增加。     

The effect of particle shape on simple shear flows

Simple shear flows, (without gravity force and implemented using periodic boundary conditions or in Couette flow configurations with gravity) have been the subject of study using DEM simulation for more than two decades. Earlier studies explored the effect of attributes such as shear rate, particle size and domain scale on the distribution of the particles in the flow, velocity profiles and the stress distributions. These studies were conducted using simple shapes for the particles such as spheres. In recent years, the importance of particle shape on flow has been recognized in a range of industrial application including mixing, comminution, hopper discharge and chute flows. In this paper, we return to the simple shear flows and quantitatively explore the effect of particle shape on velocity, volume fraction, granular temperature and stress distributions across the channel. Particle shape is found to sharply increase the strength of the material making it stronger and harder to shear. The generation of particle spin throughout the flow of non-circular particles leads to high granular temperatures, dilative pressures and lower solid fractions in the core of the flow. For aspect ratios between 0.6 and 0.5, a transition in the effective behaviour of the wall boundary conditions is identified. The connections of shape to spin, to granular temperature, to bulk flow changes are elaborated.     

Regime transition in viscous and pseudo viscous systems: A comparative study

A comprehensive quantitative study on the effect of liquid viscosity (1 ≤ µ_L ≤ 1149 mPa‐s) on the local flow phenomena of the gas phase in a small diameter bubble column is performed using ultrafast electron beam X‐ray tomography. The internal dynamic flow structure and the bubble size distribution shows a dual role of the liquid viscosity on the hydrodynamics. Further, the effect of solid concentration (C_s = 0.05, 0.20) on the local flow behavior of the gas phase is studied for the pseudo slurry viscosities similar to the liquid viscosities of the gas–liquid systems. The effects of liquid and pseudo slurry viscosities on flow structure, bubble size distribution, and gas phase distribution are compared. The bubble coalescence is significantly enhanced with the addition of particles as compared to the system without particles for apparently same viscosity. The superficial gas velocity at which transition from homogeneous bubbly to slug flow regime occurs is initiated by the addition of particles as compared to the particle free system for apparently same viscosity. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3079–3090, 2014     

Numerical simulation of the transient temperature distribution inside moving particles

This work is devoted to the two‐dimensional (2D) numerical simulation of heat and fluid flow by granular mixing in a horizontally rotating kiln. The heat and fluid flow in the gas phase are solved directly using a fixed Eulerian grid. At the same time, the particle dynamics and their collisions are solved on a Lagrangian grid. The no‐slip boundary condition on the particle surface is implemented using the fictitious boundary method. The heat transfer inside the particles is calculated utilising two models: the first is the direct solution of the energy conservation equation in Lagrangian and Eulerian space and the second is the so‐called linear model, which assumes a homogeneous distribution of the temperature inside each particle. Numerical simulations showed that if the thermal diffusivity of the gas phase significantly exceeds the same parameter of the particles, the linear model over‐predicts the heating rate of the particles. The analysis of the time‐averaged flow field inside the kiln showed that in the gas phase a double vortex structure is formed which increases the convective heat transfer in the upper part of the particulate bed. The influence of the particle size, the angular velocity of the drum and the fluid on the heating rates of particles is studied and discussed.     

Three-dimensional, rapid shear flow of particles with continuous size distributions

Rapid granular flows occur in nature and industry and often contain particles of many sizes. Over the last two decades, significant theoretical and experimental effort has been directed at rapid granular flows with monodisperse or binary particle-size distributions. In contrast, the behavior of rapid granular flows with more than two particle sizes has received only limited attention. The particle-size distributions in many natural and industrial granular flows may be represented as continuous distributions (e.g., Gaussian or lognormal distributions), providing incentive for the investigation of rapid granular flows with these particle-size distributions. As an extension of previous work for two-dimensional simulations of rapid shear flows with Gaussian and lognormal particle-size distributions, this work is directed at three-dimensional flows with continuous size distributions. Event-driven, discrete-particle (“molecular dynamic”) simulations are employed for the three-dimensional simple shear flow of smooth, spherical particles with Gaussian and lognormal particle-size distributions. The results parallel those found previously in two dimensions and demonstrate the effect of distribution width on the stress tensor and granular energy.     

Flow of gas and particles in a bubbling fluidized bed with a filtered two-fluid model

Numerical simulations of gas-particles flow in a bubble fluidized bed with two large eddy simulations of gas and solid phases are presented. For gas phase and solid phase, the sub-grid scale model for the viscosity is based on the Smagorinsky form. The sub-grid model for the particle pressure proposed by Igci et al. (2008) is modified by replacing the minimum fluidization velocity. The collisional interaction of particles is considered by the kinetic theory of granular flow. Flow behavior of gas and particles is performed by means of these two sub-grid scale models. The subgrid closure for the particle phase viscosity and pressure led to a qualitative change in the simulation results. Predictions are compared with experimental data measured by Yuu et al. (2000) and Taghipour et al. (2005) in the bubbling fluidized beds. The distributions of concentration and velocity of particles are predicted in the bubbling fluidized bed. The predicted filtered particle phase pressure increases and the filtered particle phase viscosity decreases with the increase of particle concentration. The qualitative importance of the model constant c_s of particles is demonstrated.