Simulation Study of Granular Compaction Dynamics under Vertical Tapping

We study by numerical simulation the compaction dynamics of frictional hard disks in two dimensions, subjected to vertical shaking. Shaking is modeled by a series of vertical expansions of the disk packing, followed by dynamical recompression of the assembly under the action of gravity. The second phase of the shake cycle is based on an efficient event−driven molecular−dynamics algorithm. We analyze the compaction dynamics for various values of friction coefficient and coefficient of normal restitution. We find that the time evolution of the density is described by ρ(t)=ρ∞ − ρEα[−(t/τ)α], where Eα denotes the Mittag−Leffler function of order 0