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.     

Impulsive motion of a suspension: Effect of antisymmetric stresses and particle rotation

A two-phase continuum theory (two-fluid model) for a suspension of rigid spherical particles in a Newtonian fluid is applied to investigate theoretically the flow induced by impulsive motion of an infinite flat plate. Consideration of rotational intertia of the particles gives rise to an antisymmetric part of the volume averaged stress tensor of the continuous phase. The influence of particle rotation and of antisymmetric stresses of the continuous phase, which depend on the relative rotational motion between the particles and the ambient fluid, on the motion of each phase and on the skin friction is examined.Approximate solutions to the equations, corresponding to the physical situation of large and small particle slip, are obtained by power series expansions for small and large times.     

Normal stress effects in the gravity driven flow of granular materials

In this paper, we study the fully developed gravity-driven flow of granular materials between two inclined plates. We assume that the granular materials can be represented by a modified form of the second grade fluid where the viscosity depends on the shear rate and volume fraction and the normal stress coefficients depend on the volume fraction. We also propose a new isotropic (spherical) part of the stress tensor which can be related to the compactness of the (rigid) particles. This new term ensures that the rigid solid particles cannot be compacted beyond a point, namely when the volume fraction has reached the critical/maximum packing value. The numerical results indicate that the newly proposed stress tensor has obvious and physically meaningful effects on both the velocity and the volume fraction fields.     

Volume-weighted mixture theory for granular materials

In the present work we treat granular materials as mixtures composed of a solid and a surrounding void continuum, proposing then a continuum thermodynamic theory for it. In contrast to the common mass-weighted balance equations of mass, momentum, energy and entropy for mixtures, the volume-weighted balance equations and the associated jump conditions of the corresponding physical quantities are derived in terms of volume-weighted field quantities here. The evolution equations of volume fractions, volume-weighted velocity, energy, and entropy are presented and explained in detail. By virtue of the second law of thermodynamics, three dissipative mechanisms are considered which are specialized for a simple set of linear constitutive equations. The derived theory is applied to the analysis of reversible and irreversible compaction of cohesionless granular particles when a vertical oscillation is exerted on the system. In this analysis, a hypothesis for the existence of a characteristic depth within the granular material in its closely compacted state is proposed to model the reversible compaction.     

Microcontinuum derivation of Goodman–Cowin theory for granular materials

This study presents a derivation of the Goodman–Cowin (GC) equation using the microcontinuum field theory. Through the decomposition of various microcontinuum field quantities into the straining, dilatant, and rotational parts, a microcontinuum can be classified into seven subclasses. One of the subclasses, called a microdilatation continuum, is introduced when only the dilatant motion in a macroelement is taken into account. The balance equation of equilibrated force in the GC theory can be derived while introducing the equilibrated intrinsic body force in the energy balance equation of the microdilatation continuum. The internal length of granular materials, appearing in the modified GC equation, is interpreted as the gyration radius of a macroelement. This study also obtains the evolution equation of the internal length from the microcontinuum point of view.      

A tensor method for the derivation of the equations of rigid body dynamics

A tensor method for the derivation of the equations of rigid body dynamics,based onthe concepts of continuum mechanics,is presented.The formula of time derivative of theinertia tensor with zero corotational rate is used to prove the equivalences of five methods,namely,Lagrange’s equations,Nielsen’s equations,Gibbs-Appell’s equations,Kane’sequations and the generalized momentum type of Kane’s equations.Some differentialidentities on angular velocity and angular acceleration are given.     

Micro-macro scale instability in 2D regular granular assemblies

Instability and stress–strain behavior were investigated for 2D regular assemblies of cylindrical particles. Biaxial shear experiments were performed on three sets of assemblies with regular, albeit increasingly defective structures. These experiments revealed unique instability behavior of these assemblies. Continuum models for the assemblies were then constructed using the granular micromechanics approach. In this approach, the constitutive equations governing the behavior of inter-particle contacts are written in local or microscopic level. The behavior of the RVE is then retrieved by using either kinematic constraint or least squares (static constraint) along with the principle of virtual work to equate the work done by microscopic force–displacement conjugates to that of the macroscopic stress and strain tensor conjugates. The ability of the two continuum approaches to describe the measured stress–strain behavior was evaluated. The continuum models and the local constitutive laws were used to perform instability analyses. The onset of instability and orientation of shear band was found to be well predicted by the instability analyses with the continuum models. Further, macro-scale instability was found to correlate with the instability of inter-particle contacts, although with some variations for the two modeling approaches.     

一种改进的颗粒移动床μ(I)拟流体模型及应用

颗粒移动床在工业领域应用广泛, 发展实用可靠的颗粒移动床模型具有理论和应用价值. 本文基于颗粒流μ(I)模型, 补充局部颗粒体积分数与颗粒局部压力和局部颗粒流密度的关系式, 将移动床内密集颗粒处理成可压缩拟流体, 建立了颗粒流单相可压缩流μ(I)模型, 并建立了颗粒流?壁面摩擦条件, 在计算中对颗粒流拟黏度和拟压力项进行正则化处理. 采用上述模型与方法对3种典型散料在移动床缩口料仓内的流动进行模拟, 与实验对比, 得到了玻璃珠、刚玉球和粗沙的μ(I)模型参数, 分析了3种不同散料在料仓内的颗粒速度、体积分数等分布特性, 模拟结果较好地揭示了料仓内不同物料的整体流和漏斗流特性; 进而以玻璃珠为例, 对移动床颗粒单管绕流流动进行了模拟, 所得结果合理揭示了管流附近的流动特性. 计算结果表明, 对于本文的计算工况, 颗粒体积分数变化最大范围为0.510 ~ 0.461, 绝大部分区域流动惯性数小于0.1, 改进的单相μ(I)模型能合理预测出密集颗粒流移动床内的流动特性, 方法可行且较多相流算法能明显减小计算量.      

Some Properties of a New Model for Slow Flow of Granular Materials

Some properties of a new continuum model for the bulk flow of a dense granular material in which neighbouring grains are in contact for a finite duration of time and in which the contact force is non-impulsive – the so called slow flow regime – are presented. The model generalises both the plastic potential and double-shearing models and contains an additional kinematic quantity – the intrinsic spin. The stress tensor is, in general, non-symmetric and separate yield conditions govern translational and rotational yield. We consider homogeneous, quasi-static loadings for the symmetric part of the stress and dynamic loading for the anti-symmetric part of the stress. A solution for the stress state in terms of a single parameter, namely the major principal direction of the symmetric part of the stress, is presented. This direction itself is determined by a consideration of the flow equations in the context both dilatant and isochoric simple shear flows. These simple flows are used to complete the characterisation of the relationship between the anti-symmetric part of the stress and the intrinsic spin.     

Couette flow of granular materials

The flow of granular materials between rotating cylinders is studied using a continuum model proposed by Rajagopal and Massoudi (A method for measuring material moduli for granular materials: flow in an orthogonal rheometer, DOE/PETC/TR90/3, 1990). For a steady, fully developed condition, the governing equations are reduced to a system of coupled non-linear ordinary differential equations. The resulting boundary value problem is non-dimensionalized and is then solved numerically. The effect of material parameters, i.e., dimensionless numbers on the volume fraction and the velocity fields are studied.     

Assessing continuum postulates in simulations of granular flow

Continuum mechanics relies on the fundamental notion of a mesoscopic volume “element” in which properties averaged over discrete particles obey deterministic relationships. Recent work on granular materials suggests that a continuum law may be inapplicable, revealing inhomogeneities at the particle level, such as force chains and slow cage breaking. Here, we analyze large-scale three-dimensional discrete-element method (DEM) simulations of different granular flows and show that an approximate “granular element” defined at the scale of observed dynamical correlations (roughly three to five particle diameters) has a reasonable continuum interpretation. By viewing all the simulations as an ensemble of granular elements which deform and move with the flow, we can track material evolution at a local level. Our results confirm some of the hypotheses of classical plasticity theory while contradicting others and suggest a subtle physical picture of granular failure, combining liquid-like dependence on deformation rate and solid-like dependence on strain. Our computational methods and results can be used to guide the development of more realistic continuum models, based on observed local relationships between average variables.     

纯滚动圆轮静摩擦力方向的确定

给出了匀质圆轮受力偶或水平力作用在水平地面上作纯滚动时, 其静摩擦力方向的确定方 法. 只受力偶作用作纯滚动的圆轮, 其摩擦力的方向仅取决于外力偶矩的转向, 判断依据为 摩擦力在该圆轮的平面运动中起阻碍圆轮的转动和为圆轮的平移提供主动力两相对的作用. 当圆轮只受水平力作用作纯滚动时, 用刚体平面运动微分方程和动能定理两种方法对摩擦力 方向进行了推导判定, 摩擦力的方向不仅与水平力的作用位置相关, 还与圆轮的大小和相对 质心的惯性半径等因素相关.     

Virtual experimentation through 3D full-loop simulation of a circulating fluidized bed

Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-scale （EMMS） model was used to simulate a semi-industry scale circulating fluidized bed （CFB）. Three-dimensional （3D）, time-dependent simulation of a full-loop CFB revealed that the axial profiles of cross-sectionally averaged solid volume fraction, and the radial profiles of solid axial velocity and solid volume fraction were in reasonable agreement with experimental data. Based on this agreement, database derived from experiments not yet accomplished was replenished with such simulations, and fluid regime diagrams and pressure balance around the CFB loop were derived accordingly. This work presents an integrated viewpoint on CFB and unfolds a fresh paradigm for CFB modeling, which can be expected to help resolve certain issues long in dispute but bard for experiments.     

Virtual experimentation through 3D full-loop simulation of a circulating fluidized bed

Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-male (EMMS) model was used to simulate a semi-industry scale circulating fiuidized bed (CFB).Three-dimensional(3D), time-dependent simulation of a full-loop CFB revealed that the axial profiles of cross-sectionally averaged solid volume fraction,and the radial profiles of solid axial velocity and solid volume fraction were in reasonable agreement with experimental data.Based on this agreement,database derived from experiments not yet accomplished was replenished with such simulations, and fluid regime diagrams and pressure balance around the CFB loop were derived accordingly. This work presents an integrated viewpoint on CFB and unfolds a fresh paradigm fur CFB modeling, which can be expected to help resolve certain issues long in dispute but hard for experiments.     

Averaged topological equations for dispersed two-phase flows

A general relationship between the volume fraction and the specific interfacial area for averaged dispersed two-phase flows is proposed. This relationship, expressed as a basic set of two scalar evolution equations and two vectorial non-uniformity state equations, is an analytical result obtained by a systematic approach using the derivatives of some generalized functions and a local volume-averaging technique. The proposed set of equations was expressed for measurable macroscopic parameters of the system and has the same generality as the averaged transport equations of two-phase flows. By combination of the basic set of equations, called the averaged topological equations (ATEs), second-order ATEs for the volume fraction were found. The second-order ATEs were expressed both by a Lagrangian formulation and by a Eulerian formulation. The importance and physical meaning of the ATEs developed in this study were clarified within the framework of the theory of kinematic waves.     

Micromechanical study of plasticity of granular materials

Plastic deformation of granular materials is investigated from the micromechanical viewpoint, in which the assembly of particles and interparticle contacts is considered as a mechanical structure. This is done in three ways. Firstly, by investigating the degree of redundancy of the system by comparing the number of force degrees of freedom at contacts with the number of governing equilibrium equations; Secondly, by determining the spectrum of eigenvalues of the stiffness matrix for the structure that is represented by the particles and their contacts; Thirdly, by investigating the evolution with imposed strain of the continuum elastic stiffness tensor of the system. It is found that, with increasing imposed strain, the degree of redundancy rapidly evolves towards a state with small redundancy, i.e. the system becomes nearly statically determinate. The spectrum of the system shows many singular and near-singular modes at peak shear strength and at large strain. The continuum elastic stiffness tensor becomes strongly anisotropic with increasing imposed strain and shows strong non-affinity of deformation. The assumption of a constant and isotropic elastic stiffness tensor in elasto-plastic continuum constitutive relations for granular materials is generally incorrect. Overall, the plastic continuum behaviour of granular materials originates from the plastic frictional behaviour at contacts and from damage in the form of changes in the contact network.     

A higher order subdomain method for finding local stress fields in composites

The problem of finding local and volume averaged stresses in a two-dimensional heterogeneous solid is formulated in terms of fundamental point load solutions (Green's function) leading to singular integral equations. The resulting equations are solved approximately using a subdomain method in which closed form solutions for a rectangular subdomain are obtained and utilized to find the full field solution. Previously, closed form solutions for a rectangular subvolume had been found, but only for the case of an assumed constant strain. In the present paper the solution is obtained for a quadratic form which includes not only the usual constant term but also linear and quadratic terms. The advantages of using the higher order solutions is illustrated by finding the local field in a periodic composite with square fibers. The numerical solution takes less than 90 CPU s on a workstation. The solution yields average properties independent of the reference modulus as would be expected for an accurate solution of the singular integral equation and the effective transverse modulus vs volume fraction is close to that from Christensen's model developed for round fibers.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅰ) —MICROPOLAR CONTINUA

IntroductionThemicropolarcontinuumtheoryisatypicalandwideapplicatedtheoryinthegeneralizedcontinuumtheoriesandwasgenerallyrecognizedlongago .Manymonographsandalargenumberofscientificpapersconcerningthisfieldarepublished .Wehaverestudiedtheexistingpolarcontinuumtheoriesandfoundthatsomebasicbalancelawsandprinciplesofthemareincompleteandthereexistsometheoreticaldefectswhichshouldberemoved .Forcontrastandclaritywenowlistthetraditionalbasicbalancelawsandequationsformicropolarcontinuumtheoryasfollows…     

A rigorous derivation of the bioheat equation for local tissue heat transfer based on a volume averaging theory

A general three-dimensional bioheat equation for local tissue heat transfer has been derived with less assumptions, exploiting a volume averaging theory commonly used in fluid-saturated porous media. The volume averaged energy equations obtained for the arterial blood, venous blood and tissue were combined together to form a single energy equation in terms of the tissue temperature alone. The resulting energy equation turns out to be remarkably simple as we define the effective thermal conductivity tensor, which accounts not only for the countercurrent heat exchange mechanism but also for the thermal dispersion mechanism. The present equation for local tissue heat transfer naturally reduces to the Weinbaum-Jiji equation for the unidirectional case.