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

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

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

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

滑动速度对铝合金-石墨复合材料磨损特性的影响

作者研究了滑动速度对铸造铝合金-石墨颗粒复合材料在干摩擦下滑动摩擦磨损特性的影响,指出这种复合材料在与钢干摩擦滑动时,摩擦系数和磨损速率有一个共同的滑动速度临界值。在临界值以下的摩擦系数和磨损速率都比较低且变化不大,当超过临界滑动速度时则二者都明显升高。文章还就滑动速度影响摩擦磨损性能的机理和对摩方式可能产生的影响进行了分析。     

颗粒流润滑的现状和展望

颗粒流润滑的思路是将固体材料以颗粒（粉末）状态直接导入摩擦副,使摩擦间隙充满固体颗粒,利用微小颗粒的摩擦、变形、碰撞、挤压和滑滚等微观运动,减少作相对运动两表面微凸体的直接接触,从而达到减少摩擦和保护表面的目的.总结了颗粒流润滑在航天磁悬浮备用轴承、径向轴承、制动器、磨削加工、压铸加工等工程领域中的研究和应用现状.分析和对比了已有的颗粒流润滑理论如类流体润滑理论、碰撞动理学润滑理论、离散单元和元胞自动机理论,指出摩擦间隙中颗粒流行为和颗粒接触模型是颗粒流润滑理论的关键问题.最后,从多尺度特性、非连续特性、环保等角度对颗粒流润滑在未来的进一步发展进行了展望.     

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

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

颗粒材料柱体崩塌物质点法数值模拟

为了研究颗粒材料崩塌的运动规律和堆积特性,采用物质点法对颗粒材料柱体崩塌试验进行数值模拟,并将模拟结果与试验结果进行对比验证。对颗粒材料柱体崩塌过程中颗粒的流动特性(滑动距离、堆积高度、速度、能量和动能通量的演化)进行了分析。进一步探究了颗粒材料柱体高宽比对颗粒流动能通量的影响,从而反应颗粒材料柱体崩塌过程中颗粒流的破坏能力。颗粒材料柱体高宽比越大,颗粒材料柱体外侧边缘颗粒速度越大,其溃散的程度更加强烈,并且滑动距离和动能均在增大。对于动能通量的分布,水平方向越靠近初始颗粒材料柱体,动能通量越大。     

锐缘硬质轮在金属表面的滚动与滑动

在工程应用的很多摩擦对偶中，往往有一方是锐形的，因此，研究锐形滑块的滚动与滑动具有重要的实用意义。以锐缘硬质轮在挤压铝和软钢两种平面试样表面的滚动和滑动为研究对象，考察了两种摩擦过程中材料的受力状况和变形状况，结果表明，要在同种金属试样表面造成相同宽度的槽痕，滚动法所需要的垂直负荷比滑动法所需要的大，然而滚动法所需要的切向牵引力却比滑动法所需要的小．对两种摩擦过程中界面上的应力与材料的强度性能之关系的分析讨论指出，滚动试验中接触界面上的正应力相应于变形材料的屈服应力，因而滚动试验可以作为材料硬度的一种测量方法；滑动试验中接触界面上的剪应力相应于变形材料的流动剪应力；但滚动试验中界面剪应力与材料的流动剪应力之间，以及滑动试验中界面正应力与材料的屈服应力之间却都没有直线关系．     

粗糙颗粒动理学及稠密气固两相流动的数值模拟

考虑颗粒碰撞过程中摩擦作用,给出了粗糙颗粒碰撞动力学.引入颗粒相拟总温来表征颗粒平动和转动脉动能量的特征.基于气体分子运动论,建立颗粒碰撞中平动和旋转共同作用的粗糙颗粒动理学,给出了颗粒相压力和黏度等输运参数计算模型.运用基于颗粒动理学的欧拉-欧拉气固两相流模型,数值模拟了流化床内气体颗粒两相流动特性,分析了颗粒旋转流动对颗粒碰撞能量交换和耗散的影响.模拟得到的流化床内径向颗粒浓度和提升管内颗粒轴向速度与他人实验结果相吻合.模拟结果表明随着颗粒浓度的增加,颗粒相压力和能量耗散逐渐增加,而颗粒拟总温先增加后下降.随着颗粒粗糙度系数的增加,床内平均颗粒相拟总温和能量耗散增加,表明颗粒旋转产生的摩擦将导致颗粒旋转脉动能量的改变,影响床内气体-颗粒两相宏观流动特性.      

明渠层流近床面颗粒拖曳力研究

泥沙颗粒受到的拖曳力是泥沙运动的主要驱动力,而当前应用于计算流体力学-离散颗粒法(CFD-DPM)耦合模型进行水沙运动模拟的泥沙颗粒拖曳力公式均没有考虑明渠流底床边壁作用的影响。求解不可压缩Navier-Stokes方程,对明渠层流不同雷诺数条件下床面附近不同高度处颗粒所受拖曳力进行了模拟,根据模拟结果变化规律,提出了综合考虑床面和水流惯性对标准拖曳力影响的修正拖曳力计算公式。与常用的单颗粒标准拖曳力公式和考虑遮蔽效应的多颗粒拖曳力公式相比,采用本文修正公式得到的水沙作用力更接近高精度数值解,应用于CFD-DPM输沙模拟获得的输沙结果与输沙率公式结果一致,应用分析表明输沙模拟应当采用粗糙底床边界。     

低浓度固液两相流的颗粒相动理学模型

用广义Fokker-Planck扩散模型描述液相湍动对颗粒的挟带作用，用修正的BGK模型描述粒间碰撞效应，建立了封闭的颗粒相PDF输运方程．运用Chapman-Enskog迭代法求得方程的二阶近似解，获得颗粒相脉动速度二阶矩和三阶矩闭合关系．模型与颗粒流模型相容，与液相湍流闭合模型是否相容依赖于扩散模型的具体形式，并据此比较了不同的涡一颗粒作用模型．模型与二维明渠流轻质沙和天然沙试验资料符合很好．表明细小粒径颗粒能够充分跟随水流运动；大粒径颗粒的相间平均速度差和壁面滑移速度明显，近壁区内的颗粒沿流向和垂向脉动强度都可能大于水流，并存在一定程度的颗粒碰撞效应．     

硬球加强模型在CFD-DEM耦合计算中的验证与分析

针对CFD-DEM耦合计算中,颗粒计算时间步的选取影响颗粒碰撞计算精度和效率的问题。本文引入插值算法,将动量定理求解颗粒碰撞前后速度进行加权平均;根据弹性理论计算得到颗粒碰撞力,进行动力学方程求解;通过速度收敛准则修正初值速度并自动调整迭代求解次数,提出一种计算精度不受计算时间步长影响,无需对碰撞过程进行精细描述的高效率和高精度的加强硬球模型。对两个颗粒匀和变速碰撞算例进行数值模拟,碰撞后速度、碰撞力和碰撞时间与理论计算误差小于4%,与采用软球碰撞模型的DEM方法相比,颗粒碰撞计算精度不受计算时间步长影响,计算效率提高36.3%和36.8%。对单个颗粒在静水中沉降进行数值模拟,计算步长取10 s～5 s,颗粒与壁面即可得到精确解,计算效率提高33.5%。通过压力损失实验验证了该模型能够准确计算颗粒体积分数小于12%条件下两相流的压力损失。     

聊聊动态塑性和黏塑性

固体力学研究者致力于具有应力-应变本构关系（以下简称为形变型本构关系）的变形体的力学响应研究,而流体力学研究者致力于具有应力-应变率本构关系（以下简称为流动型本构关系）的流动体的力学响应研究。当涉及结构和材料的动态塑性时,到底应该用“塑性变形”还是“塑性流动”来表示？本文从宏观塑性本构理论和微观位错动力学机理两个角度,分别讨论并指出塑性本构关系属于流动型黏塑性率相关本构关系,且同时适用于加载和卸载;因而不应该用应力-应变图来描述塑性加-卸载过程。弹塑性本构关系则是一种形变型和流动型本构关系的耦合。

     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

土的工程力学性质的颗粒流模拟

基于颗粒流理论,引入不同的颗粒接触连接本构模型,分别建立了砂土和粘性土的颗粒流模型.通过颗粒流数值模型试验,对砂土和粘性土的室内平面应变试验及其剪切带形成和发展进行了数值模拟,分别对比了不同围压下颗粒流试样与室内试验的应力应变关系曲线,基本再现了砂土和粘性土试样应力-应变关系.通过砂土和粘性土PFC试样剪切带模拟表明,当围压较小时试样内部颗粒位移量小而且分布范围较广,当围压增大时,试样内部颗粒位移量也增大,而且发生较大位移颗粒的分布范围趋于集中,同时随着围压的增大试样内部形成明显的剪切带.无论砂土还是粘性土的PFC试样,随着围压的增加剪切带的形状趋于集中,而且剪切带宽度在减小.在围压很小时,试样内形成大的破坏区域,在围压较大时出现明显的线破坏区.这些规律基本与室内试验结果相似.     

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θmodel of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θmodel in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θmodel has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θmodel was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.     

Acoustomechanical constitutive theory for soft materials

Acoustic wave propagation from surrounding medium into a soft material can generate acoustic radiation stress due to acoustic momentum transfer inside the medium and material, as well as at the interface between the two. To analyze acoustic-induced deformation of soft materials, we establish an acoustomechanical constitutive theory by com-bining the acoustic radiation stress theory and the nonlinear elasticity theory for soft materials. The acoustic radiation stress tensor is formulated by time averaging the momen-tum equation of particle motion, which is then introduced into the nonlinear elasticity constitutive relation to construct the acoustomechanical constitutive theory for soft materials. Considering a specified case of soft material sheet subjected to two counter-propagating acoustic waves, we demonstrate the nonlinear large deformation of the soft material and ana-lyze the interaction between acoustic waves and material deformation under the conditions of total reflection, acoustic transparency, and acoustic mismatch.     

Rheological aspects of dense lignite–water suspensions; structure development on consecutive flow loops

Aspects of dense lignite–water slurries (LWS) rheology were investigated using controlled stress and controlled strain rheometers with parallel disks and Couette geometries. During the preparation of the slurries, the achieved solids volume fractions were up to 0.425 and the particle size distributions were polydispersed with sizes up to 300 μm. In the ascending parts of consecutive flow loops, a slope transition of the flow curve was observed and studied in relation to the solids volume fraction. The obtained results with the different geometries and rheometers were qualitatively the same. By following the model proposed by Cheng (Rheol Acta 42:372–382, 2003) for thixotropic fluids, and taking into account the yield stress appearance, a suitable correlation for LWS is proposed, which is consistent with the experimental flow curves.     

A dynamic flow rule for viscoplasticity in polycrystalline solids under high strain rates

For visco-plasticity in polycrystalline solids under high strain rates, we introduce a dynamic flow rule (also called the micro-force balance) that has a second order time derivative term in the form of micro-inertia. It is revealed that this term, whose physical origin is traced to dynamically evolving dislocations, has a profound effect on the macro-continuum plastic response. Based on energy equivalence between the micro-part of the kinetic energy and that associated with the fictive dislocation mass in the continuous dislocation distribution (CDD) theory, an explicit expression for the micro-inertial length scale is derived. The micro-force balance together with the classical momentum balance equations thus describes the viscoplastic response of the isotropic polycrystalline material. Using rational thermodynamics, we arrive at constitutive equations relating the thermodynamic forces (stresses) and fluxes. A consistent derivation of temperature evolution is also provided, thus replacing the empirical route. The micro-force balance, supplemented with the constitutive relations for the stresses, yields a locally hyperbolic flow rule owing to the micro-inertia term. The implication of micro-inertia on the continuum response is explicitly demonstrated by reproducing experimentally observed stress–strain responses under high strain-rate loadings and varying temperatures. An interesting finding is the identification of micro-inertia as the source of oscillations in the stress–strain response under high strain rates.     

Probability density functions of contact forces for cohesionless frictional granular materials

A theory is developed for the probability density functions of contact forces for cohesionless, frictional granular materials in quasi-static equilibrium. This theory is based on a maximum information entropy principle, with an expression for information entropy that is appropriate for granular materials. Entropy is maximized under the constraints of a prescribed stress and that the normal component of the contact force is compressive and that the tangential component of the contact force is limited by Coulomb friction. The theory results in a dependence of the probability density function for the tangential contact forces on the friction coefficient. The theoretical predictions are compared with results from discrete element simulations on isotropic, two-dimensional assemblies under hydrostatic stress. Good qualitative agreement is found for means and standard deviations of contact forces and the shape of the probability density functions, while the quantitative agreement is fairly good. Discrepancies between theory and simulations, such as the difference in shape of the probability density function for the normal force and the observed dependence on elastic properties of the exponential decay rate of tangential forces, are attributed to the fact that the method does not take into account any kinematics, which are essential in relation to elastic effects.     

A CLOSED SYSTEM OF EQUATIONS FOR DENSE TWO-PHASE FLOW AND EXPRESSIONS OF SHEARING STRESS OF DISPERSED PHA’E AT A WALL

In this paper, each of the two phases in dense two-phase flow is considered as continuous medium and the fundamental equations for two-phase flow are described in Eulerian form. The generalized constitutive relation of the Bingham fluid is applied to the dispersed phase with the analysis of physical mechanism of dense two-phase flow. The shearing stress of dispersed phase at a wall is used to give a boundary condition. Then a mathematical model for dense two-phase flow is obtained. In addition, the expressions of shearing stress of dispersed phase at a wall is derived according to the fundamental model of the factional collision between dispersed-phase particles and the wall.