准静态颗粒流流动规律的热力学分析.pdfVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 3 (2014) 034702 准静态颗粒流流动规律的热力学分析∗ 刘中淼 孙其诚 宋世雄 史庆藩 1)(北京理工大学物理学院, 北京 100081) 2)(清华大学水沙科学与水利水电工程国家重点实验室, 北京 100084) ( 2013 年8 月20 日收到; 2013 年11 月5 日收到修改稿) 本文分析了颗粒流的介观结构及其特征, 提出了颗粒流的双颗粒温度概念 kin 和 conf , 表征颗 粒无序运动和构型无序演化的程度; 进而作为非平衡变量, 与经典非平衡热力学(classical irreversible thermodynamics, CIT) 变量共同构成颗粒流的热力学状态变量集, 确定了颗粒流的能量转换规律和熵产生 率等, 发展了颗粒流双颗粒温度(two granular temperate, TGT) 模型. 以体积恒定的简单剪切准静态颗粒流 为例, 结合离散元模拟(discrete element method, DEM), 确定了双颗粒温度模型所需的材料参数, 分析了颗 粒流发展段的规律和稳恒段的有效摩擦系数. 关键词: 颗粒流, 颗粒温度, 非平衡热力学, 能量耗散 PACS: 47.57.Gc, 05.70.Ln, 81.70.Pg, 61.85.+p DOI: 10.7498/aps.63.034702 面, 颗粒动理学(granular kinetics)较为成熟, 它仅 1 引 言 适用于稀疏颗粒流动, 对于准静态流和慢速流还缺 乏理论研究. 另外, 中国是世界上地质灾害最为严 颗粒流是大量离散固体颗粒流动形成的宏观 重的国家之一, 特别是“5.12” 汶川地震之后, 震区 现象, 其核心研究内容是颗粒流动的本构关系. 它 数千平方公里范围内山体松动、岩石破碎, 已发生 一般只考虑颗粒群在空气中运动, 这时可认为粒间 多起特大规模(亦即体积大于 m ) 散体灾害, 极 没有流体; 或颗粒在某一种流体中运动, 但颗粒之 具破坏力. 地质散体的蠕变、失稳起动和高速远程 间的作用占优势, 而流体作用可以忽略不计. 人们 的流动规律, 尚缺乏深入认识; 初步分析表明, 地质 通常把颗粒流划分为准静态颗粒流、慢速流和快速 灾害的这些力学过程可分别对应颗粒流的三种流 流等三种流态, 并分别采用弹塑性模型、黏塑性模 态, 因此对颗粒流流态的连贯研究能够为地质灾害 型和非牛顿流体模型来分别建立这三种流态的本 防治提供理论基础. 构关系, Forterre 等人提出了惯性系数的概念, 建 颗粒流的结构可以划分为三种尺度: 微观尺度 1 立了统一描述颗粒流本构关系的唯像模型 . 颗粒 (颗粒内部的分子和原子), 介观尺度(颗粒及颗粒 流流变关系分别涉及固体力学和流体力学的内容, 团簇) 和宏观尺度(亦即颗粒体系), 每种尺度上人 由于学科的精细划分, 使得颗粒流力学性质的连贯 们关心的问题不同. 比如在微观尺度上, 关心能量 研究格外困难, 比如颗粒恢复系数是影响快速流流 的耗散过程, 弹性的产生机理;