 Free cooling
N= 30, V= 1 x 1 = 1
n = N/V = 30, diameter d= 0.1, phi = n Pi d^2 /4 = 0.235619
T_0= 0.89744, g(d)= 1.53509, tau_b= 0.0992593, tau_e= 0.0646604, t_e= 1.36127
Restitution: 0.9
Elastic cutoff: 1e-09
Random positions at beginning
1 seconds of elastic thermalization Free cooling
N= 1000, V= 1 x 1 = 1
n = N/V = 1000, diameter d= 0.01, phi = n Pi d^2 /4 = 0.0785398
T_0= 0.951452, g(d)= 1.13726, tau_b= 0.0289202, tau_e= 0.0254296, t_e= 0.535361
Restitution: 0.9
Elastic cutoff: 1e-09
Random positions at beginning
1 seconds of elastic thermalization Gravity (no vibration)
N= 30, V= 1 x 1 = 1
n = N/V = 30, diameter d= 0.1, phi = n Pi d^2 /4 = 0.235619
T_0= 0.89744, g(d)= 1.53509, tau_b= 0.0992593, tau_e= 0.0646604, t_e= 0.718449
Restitution: 0.8
Elastic cutoff: 1e-09
Gravity: -10
Solid top with restitution 1
Solid bottom with restitution 1
Random positions at beginning
1 seconds of elastic thermalization
(collapse near t=1.5) Gravity (with too small vibration)
N= 30, V= 1 x 1 = 1
n = N/V = 30, diameter d= 0.1, phi = n Pi d^2 /4 = 0.235619
T_0= 0.89744, g(d)= 1.53509, tau_b= 0.0992593, tau_e= 0.0646604, t_e= 0.718449
Restitution: 0.8
Elastic cutoff: 1e-09
Gravity: -10
Solid top with restitution 1
Solid bottom with restitution 1
Vibrating bottom with amplitude 0.00625 and omega 50
Random positions at beginning
1 seconds of elastic thermalization
(collapse near t=1.5) Gravity (with "good" vibration)
N= 30, V= 1 x 1 = 1
n = N/V = 30, diameter d= 0.1, phi = n Pi d^2 /4 = 0.235619
T_0= 0.89744, g(d)= 1.53509, tau_b= 0.0992593, tau_e= 0.0646604, t_e= 0.718449
Restitution: 0.8
Elastic cutoff: 1e-09
Gravity: -10
Solid top with restitution 1
Solid bottom with restitution 1
Vibrating bottom with amplitude 0.05 and omega 50
Random positions at beginning
1 seconds of elastic thermalization Gravity (with "good" vibration)
N= 100, V= 1 x 1 = 1
n = N/V = 100, diameter d= 0.02, phi = n Pi d^2 /4 = 0.0314159
T_0= 1.09365, g(d)= 1.05127, tau_b= 0.134874, tau_e= 0.128296, t_e= 1.42551
Restitution: 0.8
Elastic cutoff: 1e-09
Gravity: -10
Solid top with restitution 1
Solid bottom with restitution 1
Vibrating bottom with amplitude 0.01 and omega 50
Random positions at beginning
1 seconds of elastic thermalization Vibration without gravity
N= 200, V= 50 x 50 = 2500
n = N/V = 0.08, diameter d= 2, phi = n Pi d^2 /4 = 0.251327
T_0= 1, g(d)= 1.58792, tau_b= 1.76309, tau_e= 1.11032, t_e= 23.3751
Restitution: 0.9
Elastic cutoff: 1e-09
Solid top with restitution 1
Solid bottom with restitution 1
Vibrating bottom with amplitude 0.477465 and omega 6.28319
Random positions at beginning
0 seconds of elastic thermalization

BIG animation, go here (nb: 55 MB, you are warned).
Free cooling (from homogeneous to shear/clusters)
N= 1024, V= 32 x 32 = 1024
n = N/V = 1, diameter d= 0.3, phi = n Pi d^2 /4 = 0.0706858
T_0= 0.959087, g(d)= 1.1221, tau_b= 0.960163, tau_e= 0.855682, t_e= 4.56364
Restitution: 0.5
Elastic cutoff: 1e-09
Random positions at beginning
100N collisions of elastic thermalization