快速颗粒流碰撞应力的动理模型

颗粒碰撞是快速颗粒流的动量传递,能量传递及耗散的主要机制。因此,碰撞应力是快速颗粒流研究的主要内容。本文基于大小均匀、光滑圆球颗粒的非弹性碰撞模型,类比气体分子动理论中的处理方法得出快速颗粒流积分形式的碰撞应力,能量传递通量及能量耗散率等关系。在颗粒速度分布函数的一级近似式符合Maxwell分布的假设下,对碰撞应力积分得出相应于Boltzmann方程二级近似解的碰撞应力表达式。在简单快速颗粒流动条件下,本文的理论结果与实验资料符合较好。     

Gravity-driven dry granular slow flows down an inclined moving plane: a comparative study between two concepts of the evolution of porosity

Two concepts in modeling the effects of the evolution of porosity in dry granular flows are investigated to illuminate their performance and limitations. To this end, the thermodynamic analysis, based on the Müller-Liu entropy principle, and the quasi-linear theory, are employed to deduce the ultimate constitutive models and the restrictions on their thermodynamic consistencies. The models are employed to study an isothermal dry granular slow flow down an inclined moving plane, of which the results are compared with the experimental outcomes. Results show that, while the two models deliver appropriate equilibrium expressions of the Cauchy stress tensor for compressible grains, the model in which the evolution of porosity is treated kinematically yields a spherical stress tensor for incompressible grains. Only the model with a dynamic evolution of porosity can give rise to a non-spherical stress tensor at equilibrium. Moreover, whilst the former model can better capture the characteristics of flows with slow to moderate speeds, the latter model is more able to describe the features of very rapid flows like avalanches. The present study illustrates the essential difference between the two concepts in modeling the effects of the evolution of porosity, and can be extended for further studies on other microstructural effects in granular flows.     

Granular materials: Bridging damaged solids and turbulent fluids

Granular materials exhibit abundant dissipations due to fluctuations in both granular motions and configurations (i.e., granular skeleton) evolutions. Twin granular temperatures T_k and T_e are introduced accounting for two types of fluctuations, and the so-called twin granular temperatures theory is established as an extension of granular solid hydrodynamics. By using simulations, the nonaffine deformations in a 2D assembly are simulated by using discrete element methods. By analogizing with microdamages in deformed solids, double scalar damage variables, D_P and D_q, are proposed to describe the deformed granular solid under triaxial compressions. Granular flows are found intrinsically turbulent due to the presence of T_k and the Naiver–Stokes equation is obtained for granular flows.     

颗粒流动力学及其离散模型评述

颗粒流是由众多颗粒组成的具有内在相互作用的非经典介质流动. 自然界常见颗粒流都是密集流, 颗粒间接触形成力链, 诸多力链相互交接构成支撑整个颗粒流重量和外载荷的网络, 其局部构型及强度在外载荷下演化, 是颗粒流摩擦特性和接触应力的来源.本文介绍球形颗粒间无粘连作用时的Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论. Campbell依据是否生成较为稳定的力链把颗粒流分为弹性流和惯性流两大类, 其中弹性-准静态流和惯性-碰撞流分别对应准静态流和快速流, 作为两种极端流动情况通常处理成连续体, 分别采用摩擦塑性模型和动理论予以描述, 但是表征接触力链的颗粒弹性参数并不出现这两个模型和理论框架中, 如何进一步考虑颗粒弹性参数将非常困难. 目前离散动力学方法逐渐成为复现其复杂颗粒流动现象、提取实验不可能获得的内部流动信息进而综合起来探索颗粒流问题的一种有效工具, 其真实性强于连续介质理论的描述. 软球模型对颗粒间接触力简化处理, 忽略了切向接触力对法向接触力及其加载历史的依赖, 带来了法向和切向刚度系数如何标度等更艰难的物理问题, 但由于计算强度小而广泛应用于工程问题中. 硬球模型不考虑颗粒接触变形, 因而不能描述颗粒流内在接触应变等物理机理, 仅适用于快速颗粒流, 这不仅仅是由于两体碰撞的限制. 因此基于颗粒接触力学的离散颗粒动力学模型是崭新的模型,适用于准静态流到快速流整个颗粒流态的模拟, 可以细致考虑接触形变及接触力的细节,建立更为合理的颗粒流本构关系, 进而有力的促进颗粒流这一非经典介质流动的研究.      

Nonlinear elasto-plastic model for dense granular flow

This work proposes a model for granular deformation that predicts the stress and velocity profiles in well-developed dense granular flows. Recent models for granular elasticity [Jiang, Y., Liu, M., 2003. Granular elasticity without the Coulomb condition. Phys. Rev. Lett. 91, 144301] and rate-sensitive fluid-like flow [Jop, P., Forterre, Y., Pouliquen, O., 2006. A constitutive law for dense granular flows. Nature 441, 727] are reformulated and combined into one universal elasto-plastic law, capable of predicting flowing regions and stagnant zones simultaneously in any arbitrary 3D flow geometry. The unification is performed by justifying and implementing a Kröner–Lee decomposition, with care taken to ensure certain continuum physical principles are necessarily upheld. The model is then numerically implemented in multiple geometries and results are compared to experiments and discrete simulations.     

Simulations of granular flow along an inclined plane using the Savage–Hutter model

In this work a velocity-dependent friction is introduced into a depth-averaged Savage–Hutter dynamical model for shallow granular flows. The process of granular material flowing along an inclined plane and then depositing on a horizontal plane is simulated. The surface profiles and evolution of various types of energy are investigated and compared when using the standard Coulomb-type friction versus velocity-dependent friction. Interestingly, there is a small difference between the two different types of friction.     

Simulation of mass balance behavior in a large-scale circulating fluidized bed reactor

We determine using a compound model the influence of the mass of granular matter on the behavior of a supercritical circulating fluidized bed (CFB) reactor. Population balance enables a stationary-regime modeling of the mass flow of granular matter inside a CFB unit in a large-scale. The simulation includes some important dynamic processes of gas-particle flows in fluidized bed such as attrition, fragmentation, elutriation, and fuel combustion. Numerical calculations with full boiler loading were performed of operational parameters such as furnace temperature, furnace pressure, feeding materials mass flows, and excess air ratio. Furthermore, three bed inventory masses were adopted as experimental variables in the simulation model of mass balance. This approach enables a sensitivity study of mass flows of granular matter inside a CFB facility. Some computational results from this population balance model obtained for a supercritical CFB reactor are presented that show consistency with the operational data for large-scale CFB units.     

On the equilibrium response of different modeling approaches for the porosity variation in dry granular matter

A selection of models for the variation in porosity in dry granular flows is investigated and compared on the basis of thermodynamic consistency to illustrate their performance and limitations in equilibrium situations. To this end, the thermodynamic analysis, based on the Müller–Liu entropy principle, is employed to deduce the ultimate constitutive equations at equilibrium. Results show that while all the models deliver appropriate equilibrium expressions of the Cauchy stress tensor for compressible grains, the model in which the variation in porosity is treated kinematically yields a spherical stress tensor for incompressible grains. Only the model in which the variation in porosity is modeled by a dynamic equation can give rise to a non-spherical stress tensor at equilibrium. The present study illuminates the validity and thermodynamic justification of the two modeling approaches for the porosity variation in dry granular matter.     

基于应力波动的修正非局部流变模型

基于Pouliquen 提出的非局部流变模型,考虑颗粒流中某个位置重新排列引起的自激发过程,详细分析颗粒介质中应力波动幅值的概率密度分布形式以及剪切速率与体积分数的耦合作用,提出一种修正的非局部流变模型. 采用此修正非局部流变模型对斜面剪切颗粒流的流动特性进行了预测,颗粒流动的临界厚度、平均流动速度及剪切速率廓线的预测结果与实验结果定量吻合. 此修正模型的提出为复杂的密集颗粒流的描述和表征提供了一种新的研究思路.      

A mathematical model for the waste region of a long-wall mineworking

The mechanical behaviour of the broken rock material occupying the waste region of a long-wall mineworking is investigated. A mathematical model is proposed, developed from physical considerations, which is based upon the continuum, plane strain, rigid-plastic theory for double-shearing flows of compressible granular materials. The model comprises the stress equilibrium equations, the Coulomb yield criterion, kinematic equations due to Mehrabadi and Cowin (1978) and the continuity equation, together with the waste region geometry and the boundary conditions imposed upon the field variables. Statical indeterminancy is a property of the formulation and this results in a connection between the stress and velocity fields through the boundary conditions. A strategy is presented for the iterative construction of the stress and velocity fields and approximations to solutions of the equations of the model are obtained numerically.     

Towards a constitutive law for the unsteady contact stress in granular media

In order to build a unified modelling for granular media by means of Eulerian averaged equations, it is necessary to study two contributions in the effective Cauchy stress tensor: the first one concerns solid and fluid matter, including contact and collisions between grains; the second one focuses on the random movements of grains and fluid, similar to Reynolds stress for turbulent flows. It is shown that the point of view of piecewise continuous media already used for two phase flows allows one to derive a constitutive equation for the first contribution, under the form of a partial differential equation. Similarly as for the Reynolds stress in turbulent flows, this equation cannot be written only in terms of averaged quantities without adequate approximations. The structure of the closed equation is discussed with respect to irreversible thermodynamics, and in connection with some already existing models. It is emphasised that numerical simulations by the discrete elements method can be used in order to validate these approximations, through numerical experiments statistically considered. Finally an extension of this approach to the second contribution of the effective Cauchy stress tensor is briefly discussed, showing how the modelling of both contributions have to be coupled.      

The influence of the drag force due to the interstitial gas on granular flows down a chute

Fully-developed steady flow of granular material down an inclined chute has been a subject of much research interest, but the effect of the interstitial gas has usually been ignored. In this paper, new expressions for the drag force and energy dissipation caused by the interstitial gas (ignoring the turbulent fluctuations of the gas phase) are derived and used to modify the governing equations derived from the kinetic theory approach for granular–gas mixture flows, where particles are relatively massive so that velocity fluctuations are caused by collisions rather than the gas flow. This new model is applied to fully-developed, steady mixture flows down an inclined chute and the results are compared with other simulations. Our results show that the effect of the interstitial gas plays a significant role in modifying the characteristics of fully developed flow. Although the effect of the interstitial gas is less pronounced for large particles than small ones, the flowfields with large particles are still very different from granular flows which do not incorporate any interactions with the interstitial gas.     

Concurrent use of finite element and finite volume methods for shallow water flows in locally 1‐D channel networks

A numerical model is developed for shallow water equation in locally 1‐D channel networks. The model concurrently uses the standard Galerkin finite element method for the continuity equation and the finite volume method with an upwind scheme for the momentum equation. The surface gradient method is consistently employed. A minimum treatment is given for channel junctions so that application to multiply connected channels do not require any special consideration The model is capable of computing different types of transcritical flows, wet and dry flows, and flows with complex source terms. Standardized test problems and laboratory experimental data are used for verifying the model. Applicability of the models is validated in a multiply connected channel network draining hydromorphic farmlands located in a West African savanna, and Manning's roughness coefficient is identified, so that the steady solution is consistent with field observations. Unsteady simulation demonstrates that the model is capable of stably reproducing shifts of hydraulic jumps in flows of sub‐millimeter water depths. Copyright © 2011 John Wiley & Sons, Ltd.     

颗粒柱塌落中的尺寸效应和瞬态流变性研究

颗粒材料是自然界和工程中广泛存在和普遍应用的材料,泥石流、滑坡和混凝土等均可视为颗粒材料。颗粒材料研究有助于更好地控制相关自然灾害或利用其某些特性。对颗粒材料柱体的塌落动力学研究不仅可以方便理解颗粒材料在瞬态流动时的流变性,还可以引申到泥石流等岩土材料的运动与堆积形态。本文利用扩展多面体离散单元法对颗粒堆积柱的塌落进行了细致的研究,探索高宽比、摩擦系数以及颗粒柱相对尺寸等对柱体塌落的影响。对颗粒集合进行分网并分析每个网格内的应力与应变率之间的关系,讨论其瞬态本构关系。相关研究对于深入理解颗粒材料重力流的动力学性质以及颗粒集合体的堆积形态具有重要意义。     

颗粒流的动力学模型和实验研究进展

首先叙述了颗粒流研究的背景和基本概念,接着概述了颗粒流动力学研究的三种基本数学力学模型:颗粒动理论(kinetic theory)模型、摩擦塑性模型和离散元模拟模型.介绍了它们的基本原理,并分析了各个模型的优缺点和适用范围.随之,对颗粒流的实验研究情况作了简要的介绍,列举了研究中几种常用的试验方法和观测手段,以及观测到的一些典型现象.最后,我们简述了颗粒流研究的主要困难,并提出了某些有待解决的研究课题.      

Flow through isotropic granular porous media

A generalization of the Navier-Stokes equation is developed to include laminar flow through a rigid isotropic granular porous medium of spatially varying permeability. The model is based on a theory of interspersed continua and the mean geometrical properties of an idealized granular porous microstructure. The derived momentum transport equations are applicable to granular porous media over the entire porosity range from zero through unity. No restriction with respect to flow velocity is imposed, except for the assumption of laminar flow within the pores. The results provide useful and versatile equations and substantiate many of the empirical equations currently in use. One of the major advantages of the generalized momentum equation is its adaptability to numerical simulation.     

Simulated configurational temperature of particles and a model of constitutive relations of rapid-intermediate-dense granular flow based on generalized granular temperature

In gas–solid flat-base spout bed with a jet, the flow of particles must go through an intermediate regime where both kinetic/collisional and frictional contributions play a role. In this paper, the statistical framework is proposed to define the generalized granular temperature which sums up the configurational temperature and translational granular temperature. The configurational temperature, translational and rotational granular temperatures of particles are simulated by means of CFD-DEM (discrete element method) in a 3D flat-base spout bed with a jet. The configurational temperatures of particles are calculated from instantaneous overlaps of particles. The translational and rotational granular temperatures of particles are calculated from instantaneous translational and angular velocities of particles. Roughly, the simulated translational and rotational granular temperatures increase, reach maximum, and then decrease with the increase of solids volume fractions. However, the configurational temperature increases with the increase of solids volume fractions. At high solid volume fraction, the predicted configurational temperatures are larger than the translational and rotational granular temperatures, indicating that the rate of energy dissipation do contributes by contact deformation of elastic particles. The generalized granular temperature is proposed to show the relation between the variance of the fluctuation velocity of deformation and the variance of the translational fluctuation velocity of particles. The constitutive relations of particle pressure, viscosity, granular conductivity of fluctuating energy and energy dissipation in rapid-intermediate-dense granular flows are correlated to the generalized granular temperature. The variations of particle pressure, shear viscosity, energy dissipation and granular conductivity are analyzed on the basis of generalized granular temperature in a flat-base spout bed with a jet. The axial velocities of particles predicted by a gas–solid two-fluid model of rapid-intermediate-dense granular flows agree with experimental results in a spout bed.     

Solving the Savage–Hutter equations for granular avalanche flows with a second‐order Godunov type method on GPU

In this article, we apply Davis's second‐order predictor‐corrector Godunov type method to numerical solution of the Savage–Hutter equations for modeling granular avalanche flows. The method uses monotone upstream‐centered schemes for conservation laws (MUSCL) reconstruction for conservative variables and Harten–Lax–van Leer contact (HLLC) scheme for numerical fluxes. Static resistance conditions and stopping criteria are incorporated into the algorithm. The computation is implemented on graphics processing unit (GPU) by using compute unified device architecture programming model. A practice of allocating memory for two‐dimensional array in GPU is given and computational efficiency of two‐dimensional memory allocation is compared with one‐dimensional memory allocation. The effectiveness of the present simulation model is verified through several typical numerical examples. Numerical tests show that significant speedups of the GPU program over the CPU serial version can be obtained, and Davis's method in conjunction with MUSCL and HLLC schemes is accurate and robust for simulating granular avalanche flows with shock waves. As an application example, a case with a teardrop‐shaped hydraulic jump in Johnson and Gray's granular jet experiment is reproduced by using specific friction coefficients given in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

An integrated smoothed particle hydrodynamics model for complex interfacial flows with large density ratios

In this paper, an integrated smoothed particle hydrodynamics (SPH) model for complex interfacial flows with large density ratios is developed. The discrete continuity equation and acceleration equation are obtained by considering the time derivative of the volume of particle and Eckart's continuum Lagrangian equation. A continuum surface force model is used to meet the fact that surface force may not be distributed uniformly on each side of the interface. An improved boundary condition is imposed to model wall free-slip and no-slip condition for interfacial flows with large density ratios. Particle shifting algorithm (PSA) is added for interfacial flows by imposing the normal correction near the interface, called as Interface-PSA. Then four representative numerical examples, including droplet deformation, Rayleigh-Taylor instability, dam breaking, and bubble rising, are presented and compared well with reference data. It is demonstrated that inherent interfacial flow physics can be well captured, including surface tension and the dynamic evolution of the complex interfaces.     

THE KEY SCIENTIFIC PROBLEMS IN THE EULERIAN MODELING OF LARGE-SCALE MULTI-PHASE FLOWS——DRAG MODEL

描述相间相互作用的曳力模型是决定流态化模拟成败的关键科学问题,核心是如何表征稠密非均匀流动中曳力的一般规律. 首先综述了现有曳力模型的发展现状和优缺点,然后讨论了流态化系统能量分析方法（energy minimization multiscale method, EMMS）,以及在此基础上建立的曳力模型及其发展历程,重点分析了EMMS 理论在代表非均匀流动特征的颗粒团尺寸、内部固含率等关键参数上的缺陷、文中的修正方法和实验检验结果. 最后,讨论了曳力模型的发展方向即普适性问题.