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

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

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)

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)

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

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

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).

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