LIGGGHTS: Dense packing cannot be achieved in some conditions

Submitted by qfhou on Wed, 09/19/2012 - 05:54
Plain text icon cubic.txt3.12 KB
Image icon loose_packing.jpg1.75 MB
Plain text icon log.liggghts.txt16.07 KB
Image icon cubic.avi_.jpg16.99 MB

Problem: dense packing under gravity will not be formed even the particle velocities are nearly zero.
Conditions: insert particles more than a critical number.

tested versions 2.0.4, 2.1.1

I think this should be a bug where the gravity field is lost when the total particle number is larger than a critical value (I am doing test and will tell the critical value in the following given case).

Files attached: geometry (cubic.stl, extension has to be changed), input (in.packing, extension as well), and a result using liggghts 2.0.4.