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

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

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

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

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

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

颗粒体流-固特性的实验研究方法评述

颗粒材料是一种常见的复杂多体物质,其形态各不相同,颗粒之间通过相互作用使其表现出十分复杂的流体或固体特性。因此结合理论、实验及模拟的方法对颗粒材料进行研究,可以进一步发现颗粒的物理规律并对其进行表征,同时也有助于学者对颗粒材料的特性形成更加系统的认识。本文对颗粒材料在不同物质状态下表现出的不同流态特性与多尺度特点进行简要介绍,同时回顾了颗粒材料流-固特性的研究现状与2种状态下的实验研究方法,概括了相关测量方法的优缺点。     

椭圆颗粒在剪切流中旋转特性的数值研究

研究颗粒在流体剪切作用下的运动特性是理解和预测颗粒悬浮流流动行为的关键.当流体的惯性不能忽略时,颗粒的运动往往变得非常复杂.本文采用格子Boltzmann方法对中等雷诺数下椭圆颗粒在剪切流中的旋转运动进行了模拟.首先,研究了雷诺数(0Re 170)的影响,结果表明当雷诺数低于临界值时,颗粒以周期性的方式旋转,角速度最小时对应的长轴方向随着雷诺数的增大而逐渐远离水平方向,而且这一倾角与雷诺数呈分段线性关系;当雷诺数大于临界值时,椭圆形颗粒最终保持静止状态,且静止时的转角与雷诺数呈幂函数关系,雷诺数越大,转角越小,椭圆的长轴越远离水平位置.其次,研究了椭圆颗粒的长短轴之比α(1α10)的影响,结果表明颗粒旋转的周期与α呈幂函数关系,α越大,颗粒旋转周期越小.此外,当α超过临界值时,颗粒也在水平位置附近保持静止状态,此时的转角与α也呈幂函数关系,α越大,转角越小.研究还发现,当雷诺数较大时椭圆颗粒在旋转过程中会产生过冲现象.     

考虑颗粒碰撞的多重Monte Carlo算法

从减少计算代价和改进碰撞算法出发, 提出了考虑颗粒碰撞的多重Monte Carlo算法, 它采用直接模拟Monte Carlo算法来考虑颗粒碰撞, 并与求解颗粒拉氏Langevin方程的Monte Carlo算法耦合起来, 跟踪比实际颗粒数目小得多的虚拟颗粒. 提出了时间步长选定标准、虚拟碰撞伙伴所在控制容积的判断准则、颗粒碰撞发生的判 断准则、虚拟碰撞伙伴的选择、基于随机碰撞角度的碰撞动力学, 构成了考虑颗粒碰撞的完整多重Monte Carlo算法. 对理想工况的细微颗粒流和粗重颗粒流进行了数值模拟, 颗粒碰撞率的模拟结果与理论分析解和DNS结果均符合很好, 颗粒场演变的细节信息, 如时间平均和特定时刻的颗粒数密度, 速度和颗粒湍动能等, 均与DNS结果符合很好. 数值模拟结果证明该算法不仅具有较低的计算代价, 而且能够达到足够的计算精度.     

可燃颗粒悬浮流中的正激波

本文讨论了激波后可燃颗粒点火、燃烧的机制和波后流场的松弛结构,结果表明,点火阶段颗粒外层浓度梯度极大,燃烧速率由非均相化学反应速率决定,而在燃烧阶段,外层的浓度梯度极大,燃烧速率由气体的扩散速率所决定。在考虑管壁效应时,对守恒方程中的能量耗散项作了说明。     

基于颗粒流的混凝土材料数值实验研究

混凝土材料是一种典型的非均质材料,其力学行为和破坏过程很复杂.本文利用PFC颗粒流软件的内嵌fish语言编写了骨料生成程序、边界条件控制程序,程序模拟了MTS伺服功能,实现了虚拟的混凝土单轴压缩试验.相对于一般有限元程序投放的骨料,文中骨料可以破裂.在生成的数值试件基础上,设计了单级配和全级配混凝土单轴压缩试验,并对其破裂形态、裂纹扩展过程、应力应变曲线和破裂过程能量变化规律进行了研究.利用该虚拟试验平台,可对混凝土的动、静力学特性进行系统研究.     

颗粒流润滑的现状和展望

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

水沙流中离散颗粒流动应力关系

在已有的基于膨胀体模型、动理论模型和连续介质理论得到的关系式基础上,分析了水沙流中颗粒流动的应力关系应由摩擦应力、滑动应力和碰撞应力3部分组成,并将颗粒弹性恢复系数引入,对颗粒碰撞项进行了修正。根据颗粒碰撞时间间隔T_c和颗粒弛豫时间T_e的对比,提出了颗粒流动进入碰撞惯性区的临界值,即当T_c/T_e<0.02时,水沙流中颗粒流动碰撞作用占优,可忽略摩擦应力和滑动应力,并将提出的临界值与实验资料进行了比较。利用修正后的水沙流中颗粒流动应力关系,对缓槽水沙流动进行了模拟计算,获得了与实测资料一致的结果。     

关于二相流、多相流、多流体模型和非牛顿流等概念的探讨

本文分析了单相流、二相流和多相流等概念上的差异,也分析了单流体模型、双流体模型和多流体模型等概念上的差异,指出前面三种概念是按流动介质的客观物理构成划分的,而后者是按主观采用的研究方法划分的．目前这些概念在使用中存在一些混乱,如二相流与多相流,多相流与多流体模型等．本文还研究了扩散模型、非牛顿流模型和颗粒流模型等,指出前两种模型在分类上属于单流体模型,分析了非牛顿流模型、扩散模型和双（多）流体模型的特点和应用范围,最后,以泥石流为例讨论了以上概念的应用．      

颗粒流润滑系统力传输行为的试验研究

为了研究颗粒流润滑状态下的力传输行为,设计了一个剪切平行板间颗粒流润滑的试验装置.利用该试验装置,研究了颗粒材料、颗粒尺寸、剪切速度和颗粒湿度对颗粒流润滑状态下力传输行为的影响.此外,试验结果主要考虑了测量得到的法向力和切向力,这些力是从下摩擦副通过颗粒润滑介质传递到上摩擦副的.研究结果表明,测量得到的法向力和切向力随着剪切速度和颗粒湿度的增大而减小,随着颗粒尺寸的增大而增大.颗粒流润滑比干摩擦具有较好的润滑效果,比油润滑具有更好的承载性能.颗粒流润滑状态下的力传输路径主要决定于颗粒润滑介质之间的力链.试验结果揭示了颗粒流润滑时的力传递机制并可以为颗粒流润滑的工业应用提供指导.     

壁面展向震荡诱导颗粒湍槽流减阻的直接数值模拟研究

对湍槽流的减阻研究具有科学意义和工程应用价值,已有大量研究表明向单相湍流中添加离散物质是一种有效的被动减阻方法.相比于被动减阻技术,主动减阻技术如壁面震荡减阻的可控性更高,近年来也得到广泛的关注,但对于壁面展向震荡诱导减阻的研究主要针对单相湍槽流,还未见有相关研究将这一手段用于含颗粒湍槽流的减阻.因此,文章采用直接数值模拟方法开展了壁面展向震荡诱导颗粒湍槽流减阻的机理研究.一方面关注壁面震荡对颗粒湍槽流的调制效果及机理.另一方面关注颗粒和震荡对单相湍槽流的耦合减阻效应.结果表明:壁面震荡可以达到有效减阻,存在最优震荡周期使减阻率达到最大,且最优震荡周期与单相流结果相近.在相同体积分数下,施加壁面震荡的小颗粒湍槽流减阻效果更好.相比于单相湍槽流,当震荡周期小于最优周期时,震荡和颗粒的耦合效应对减阻率的额外贡献较小且可能为负,当大于最优周期时额外贡献逐渐增大,对整体减阻率的占比最高可达10%左右.     

颗粒在黏弹性流体扩张和收缩流中沉降的格子Boltzmann模拟分析

对于Oldroyd-B型黏弹性流体,本文应用格子Boltzmann方法(LBM),实现了流体在二维1:3扩张流道及3:1收缩流道中流动的数值模拟,获得了黏弹性流体在扩张和收缩流道中的流场分布．结合颗粒的受力和运动规则,基于点源颗粒模型,数值分析了颗粒在扩张流和收缩流中的沉降过程和特征,讨论了颗粒相对质量和起始位置以及雷诺数Re和威森伯格数Wi对颗粒沉降特征的影响．结果表明,颗粒相对质量和起始位置以及Re对颗粒沉降轨迹和落点影响较大,而Wi的影响则较小．     

两种颗粒湍流扩散模型数值模拟气液两相流泄漏扩散的比较

采用相间耦合的欧拉-拉格朗日方法,模拟了装有液化气（丙烷）的容器出现小孔或裂缝时,发生泄漏后的气液两相扩散过程。分别应用随机轨道模型及颗粒群模型来考察湍流对液滴扩散的影响,并与确定性轨道进行了比较。给出了数学物理模型,计算结果与实验数据进行了对比。结果表明：采用随机轨道模型能较好地描述液滴的湍流扩散,适用于有液相蒸发的两相流扩散问题。     

测量两相流中颗粒特性的相位多普勒方法理论模型

本文在几何光学近似Ｍｉｅ散射理论的基础上建立了一种计算激光多普勒信号中相左与粒径关系的数学模型。对几种常用的相对折射率情况下的相差一粒径特性进行了计算,分析并给出了光学系统的最佳配置。计算结果与Ｍｉｅ散射理论的计算结果进行了对比,有限好的一致性。     

床面固体颗粒圆柱绕流的试验现象观察

本文介绍了床面固体颗粒随水流绕过圆柱体时，将在圆柱周围的床面上形成一个无粒子运动区的试验现象。水槽试验结果表明，当固体颗粒的粒径减小时，无粒子区的范围将增大；无粒子区的范围随圆柱直径的增大而增大；水流条件的变化直接影响着床面固体颗粒的运动情况，同无粒子区的形成、消失及范围大小有密切的关系。根据试验资料，结合量纲分析，建立了无粒子区的无量纲经验关系式     

纳米颗粒多相流体动力学研究及应用

纳米颗粒多相流研究是目前多相流研究中新的研究方向及重点发展领域. 为探索纳米尺度多相流相间作用机理及内部存在机制,采用理论分析及数值计算手段,对一般动力学方程的封闭处理、颗粒碰撞率宏观模型的有效构建、颗粒凝并系统动力学演变特性的机理分析、非稀相问题碰撞率的求取、双变量问题求解方法的建立以及一些实际应用进行了系统研究,提出了新的针对纳米尺度颗粒动力学演变的一般动力学方程求解方法,并将其应用于实际工业过程问题的研究. 该文对上述研究工作进行了综述.     

中等Reynolds数平板绕流的动理论分析

郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边界层理论中的经典工作.许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估,评估的依据是Janour与Schaaf和Sherman的实验数据.本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边界层理论中的经典工作.许多研究者对平板绕流阻力系数的郭水怀公式以及后续工作进行了评估,评估的依据是Janour与Schaaf和Sherman的实验数据.本文的动理论分析和计算表明:Schaaf和Sherman在低亚声速条件下(0.16相似文献   

周期荷载作用下煤岩声发射特征的颗粒流模拟

为研究周期荷载应力水平对煤样声发射特征的影响,采用PFC数值软件开展了 3种不同应力水平的等幅周期荷载数值模拟试验,分析了周期荷载应力水平对煤样破坏循环次数、声发射计数及损伤演化特征的影响.研究结果表明:周期荷载作用下,煤样破坏过程中的声发射活动呈现初始、相对平静和活跃三阶段演化规律,且在煤样破坏前的周期荷载卸载阶段及...