LAMMPS (7 Jul 2009) # Pour granular particles into chute container, then induce flow atom_style granular boundary p p fm newton off region reg block -10 10 -10 10 -0.5 16 units box create_box 1 reg Created orthogonal box = (-10 -10 -0.5) to (10 10 16) 2 by 2 by 1 processor grid neighbor 0.2 bin neigh_modify delay 0 # IMPORTANT NOTE: these values of stiffness (2000) and timestep (0.001) # are used in this example file to produce a quick simulation and movie. # More appropriate values for realistic simulations would be # k = 2.0e5 and dt = 0.0001, as in bench/in.chute. pair_style gran/hooke/history 2000.0 NULL 50.0 NULL 0.5 0 pair_coeff * * timestep 0.001 fix 1 all nve/sphere fix 2 all gravity 1.0 spherical 0.0 -180.0 fix zlower all wall/gran 2000.0 NULL 50.0 NULL 0.5 0 zplane 0.0 2000.0 region slab block -9.0 9.0 -9.0 9.0 10.0 15.0 units box fix ins all pour 3000 1 300719 vol 0.13 50 region slab Particle insertion: 402 every 3162 steps, 3000 by step 22135 compute 1 all erotate/sphere thermo_style custom step atoms ke c_1 vol thermo 1000 thermo_modify lost ignore norm no compute_modify thermo_temp dynamic yes dump id all atom 1000 dump.pour run 25000 Memory usage per processor = 3.14935 Mbytes Step Atoms KinEng 1 Volume 0 0 -0 0 6600 1000 402 779.98578 0 6600 2000 402 1424.4252 0 6600 3000 402 1425.0863 12.9927 6600 4000 804 1724.8627 47.36161 6600 5000 804 1599.9516 63.808734 6600 6000 804 1359.1774 63.495338 6600 7000 1206 1532.7896 68.602762 6600 8000 1206 1403.019 62.696462 6600 9000 1206 1246.5766 52.235655 6600 10000 1608 1384.2625 50.982242 6600 11000 1608 1243.1052 51.141884 6600 12000 1608 1195.3028 40.935179 6600 13000 2010 1328.4985 50.722781 6600 14000 2010 1198.4098 45.628483 6600 15000 2010 997.57313 41.390864 6600 16000 2412 1037.0841 43.223404 6600 17000 2412 1018.8984 38.835381 6600 18000 2412 930.69902 35.880622 6600 19000 2814 1007.138 37.703169 6600 20000 2814 967.57176 41.267685 6600 21000 2814 772.7741 44.763252 6600 22000 2814 525.78939 38.98055 6600 23000 3000 429.75161 34.235783 6600 24000 3000 364.4473 21.966186 6600 25000 3000 218.90078 24.596278 6600 Loop time of 12.4018 on 4 procs for 25000 steps with 3000 atoms Pair time (%) = 2.93148 (23.6375) Neigh time (%) = 0.440789 (3.55422) Comm time (%) = 5.24973 (42.3302) Outpt time (%) = 0.0958874 (0.773171) Other time (%) = 3.68395 (29.7049) Nlocal: 750 ave 757 max 742 min Histogram: 1 0 0 1 0 0 0 0 1 1 Nghost: 389.75 ave 396 max 382 min Histogram: 1 0 0 1 0 0 0 0 1 1 Neighs: 3579 ave 3625 max 3541 min Histogram: 1 0 0 0 2 0 0 0 0 1 Total # of neighbors = 14316 Ave neighs/atom = 4.772 Neighbor list builds = 1152 Dangerous builds = 0 unfix ins fix 2 all gravity 1.0 chute 26.0 run 25000 Memory usage per processor = 9.68014 Mbytes Step Atoms KinEng 1 Volume 25000 3000 218.90078 24.596278 6600 26000 3000 108.47989 19.098297 6600 27000 3000 95.772626 12.451488 6600 28000 3000 168.56408 11.2661 6600 29000 3000 266.55027 12.736385 6600 30000 3000 429.91963 16.099577 6600 31000 3000 653.50177 20.833152 6600 32000 3000 941.4538 27.486019 6600 33000 3000 1287.1593 35.913238 6600 34000 3000 1705.1323 44.1634 6600 35000 3000 2274.278 50.268183 6600 36000 3000 2864.078 55.78932 6600 37000 3000 3565.6924 69.196132 6600 38000 3000 4360.8456 82.984183 6600 39000 3000 5227.4355 101.09226 6600 40000 3000 6138.5626 102.21419 6600 41000 3000 7143.2929 117.41759 6600 42000 3000 8246.8072 123.70515 6600 43000 3000 9409.9225 134.14761 6600 44000 3000 10593.288 157.7653 6600 45000 3000 12018.948 178.91396 6600 46000 3000 13410.044 196.8436 6600 47000 3000 14750.988 191.71558 6600 48000 3000 16088.89 203.0547 6600 49000 3000 17512.26 218.17825 6600 50000 3000 18836.338 213.33912 6600 Loop time of 17.9768 on 4 procs for 25000 steps with 3000 atoms Pair time (%) = 6.41246 (35.6707) Neigh time (%) = 0.669206 (3.72261) Comm time (%) = 6.30764 (35.0876) Outpt time (%) = 0.143959 (0.800805) Other time (%) = 4.44356 (24.7183) Nlocal: 750 ave 757 max 747 min Histogram: 1 2 0 0 0 0 0 0 0 1 Nghost: 408.5 ave 414 max 402 min Histogram: 1 0 0 0 0 2 0 0 0 1 Neighs: 3665.75 ave 3717 max 3597 min Histogram: 1 0 0 0 1 0 0 0 1 1 Total # of neighbors = 14663 Ave neighs/atom = 4.88767 Neighbor list builds = 821 Dangerous builds = 0