Collision statistics in sheared inelastic hard spheres.
Collision statistics in sheared inelastic hard spheres.
复制标题
剪切非弹性硬球中的碰撞统计。
DOI：
--
复制
发表时间：
2009
期刊：
影响因子：
--
通讯作者：
L. Lue
中科院分区：
文献类型：
--
作者：
M. N. Bannerman;T. E. Green;P. Grassia;L. Lue
采用非平衡态分子动力学模拟和直接模拟Monte Carlo方法研究了剪切非弹性硬球系统的动力学行为。在分子动力学模拟中，Lees-Edwards边界条件用于施加剪切。选择模拟箱的尺寸，以确保系统是均匀的，并且均匀地施加剪切。各种系统属性进行监测，包括单粒子速度分布，颗粒温度，应力张量，碰撞率，碰撞之间的时间。单粒子速度分布被认为是同意合理的各向异性高斯分布，只有轻微的人口过剩的高速尾巴。速度分布是强烈的各向异性，特别是在较低的密度和较低的恢复系数值，在剪切方向上的最大方差。剪切非弹性硬球系统的压缩因子的密度依赖性与弹性硬球流体的压缩因子的密度依赖性非常相似。当系统变得更无弹性时，掠射碰撞开始超过更直接的正面碰撞。对碰撞之间时间分布的研究表明，在高度非弹性系统中，粒子所经历的碰撞是强相关的。模拟数据的比较与直接的Monte Carlo模拟的Enskog方程。还包括基于Enskog方程的Montanero [J. Fluid Mech.389，391（1999）]动力学模型的结果。在一般情况下，良好的协议被发现为高密度，弱非弹性系统。
The dynamics of sheared inelastic-hard-sphere systems is studied using nonequilibrium molecular-dynamics simulations and direct simulation Monte Carlo. In the molecular-dynamics simulations Lees-Edwards boundary conditions are used to impose the shear. The dimensions of the simulation box are chosen to ensure that the systems are homogeneous and that the shear is applied uniformly. Various system properties are monitored, including the one-particle velocity distribution, granular temperature, stress tensor, collision rates, and time between collisions. The one-particle velocity distribution is found to agree reasonably well with an anisotropic Gaussian distribution, with only a slight overpopulation of the high-velocity tails. The velocity distribution is strongly anisotropic, especially at lower densities and lower values of the coefficient of restitution, with the largest variance in the direction of shear. The density dependence of the compressibility factor of the sheared inelastic-hard-sphere system is quite similar to that of elastic-hard-sphere fluids. As the systems become more inelastic, the glancing collisions begin to dominate over more direct, head-on collisions. Examination of the distribution of the times between collisions indicates that the collisions experienced by the particles are strongly correlated in the highly inelastic systems. A comparison of the simulation data is made with direct Monte Carlo simulation of the Enskog equation. Results of the kinetic model of Montanero [J. Fluid Mech. 389, 391 (1999)] based on the Enskog equation are also included. In general, good agreement is found for high-density, weakly inelastic systems.