Alvarez-Hamelin J Ignacio, Puglisi Andrea
CONICET and Facultad de Ingeniería, Universidad de Buenos Aires, Paseo Colón 850, Buenos Aires, Argentina.
Phys Rev E Stat Nonlin Soft Matter Phys. 2007 May;75(5 Pt 1):051302. doi: 10.1103/PhysRevE.75.051302. Epub 2007 May 2.
We address the problem of recollisions in cooling granular gases. To this aim, we dynamically construct the interaction network in a granular gas, using the sequence of collisions collected in an event driven simulation of inelastic hard disks from time 0 until time t . The network is decomposed into its k -core structure: particles in a core of index k have collided at least k times with other particles in the same core. The difference between cores k+1 and k is the so-called k -shell, and the set of all shells is a complete and nonoverlapping decomposition of the system. Because of energy dissipation, the gas cools down: its initial spatially homogeneous dynamics, characterized by the Haff law, i.e., a t{-2} energy decay, is unstable toward a strongly inhomogeneous phase with clusters and vortices, where energy decays as t{-1} . The clear transition between those two phases appears in the evolution of the k -shells structure in the collision network. In the homogeneous state the k -shell structure evolves as in a growing network with a fixed number of vertices and randomly added links: the shell distribution is strongly peaked around the most populated shell, which has an index k{max} approximately 0.9(d) with (d) the average number of collisions experienced by a particle. During the final nonhomogeneous state a growing fraction of collisions is concentrated in small, almost closed, communities of particles: k{max} is no more linear in (d) and the distribution of shells becomes extremely large developing a power-law tail approximately k{-3} for high shell indexes. We conclude proposing a simple algorithm to build a correlated random network that reproduces, with few essential ingredients, the whole observed phenomenology, including the t{-1} energy decay. It consists of two kinds of collisions (links): single random collisions with any other particle and long chains of recollisions with only previously encountered particles. The algorithm disregards the exact spatial arrangement of clusters, suggesting that the observed stringlike structures are not essential to determine the statistics of recollisions and the energy decay.
我们研究了冷却颗粒气体中的再碰撞问题。为此,我们在颗粒气体中动态构建相互作用网络,利用从时间0到时间t的非弹性硬球事件驱动模拟中收集的碰撞序列。该网络被分解为其k -核结构:索引为k的核中的粒子与同一核中的其他粒子至少碰撞了k次。k + 1核与k核之间的差异就是所谓的k -壳层,所有壳层的集合是系统的一个完整且不重叠的分解。由于能量耗散,气体冷却:其初始的空间均匀动力学,以哈夫定律为特征,即能量按t^(-2)衰减,对于具有团簇和涡旋的强非均匀相是不稳定的,在该相中能量按t^(-1)衰减。这两个相之间的明显转变出现在碰撞网络中k -壳层结构的演化中。在均匀状态下,如果顶点数量固定且随机添加链接,k -壳层结构的演化类似于一个增长的网络:壳层分布在人口最多的壳层周围强烈峰值,其索引k_max约为0.9(d),其中(d)是粒子经历的平均碰撞次数。在最终的非均匀状态下,越来越多的碰撞集中在小的、几乎封闭的粒子群落中:k_max不再与(d)呈线性关系,壳层分布变得极其分散,对于高壳层索引发展出近似k^(-3)的幂律尾部。我们通过提出一种简单算法来构建相关随机网络,该算法用很少的基本要素重现了整个观察到的现象学,包括t^(-1)能量衰减,从而得出结论。它由两种碰撞(链接)组成:与任何其他粒子的单次随机碰撞以及仅与先前遇到的粒子的长链再碰撞。该算法忽略了团簇的确切空间排列,这表明观察到的线状结构对于确定再碰撞统计和能量衰减并非必不可少。