LAMMPS (5 Oct 2016) # Pour 2d granular particles into container dimension 2 atom_style sphere boundary f fm p newton off comm_modify vel yes region reg block 0 100 0 50 -0.5 0.5 units box create_box 1 reg Created orthogonal box = (0 0 -0.5) to (100 50 0.5) 1 by 1 by 1 MPI processor grid neighbor 0.2 bin neigh_modify delay 0 # IMPORTANT NOTE: these values of stiffness (4000) 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 = 4.0e5 and dt = 0.0001, as in bench/in.chute (if it were Hertzian). pair_style gran/hertz/history 4000.0 NULL 100.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 xwalls all wall/gran hertz/history 4000.0 NULL 100.0 NULL 0 1 xplane 0 100 fix ywalls all wall/gran hertz/history 4000.0 NULL 100.0 NULL 0 1 yplane 0 NULL region slab block 1.0 99.0 30 34.5 -0.5 0.5 units box fix ins all pour 1000 1 4767548 vol 0.4 10 diam range 0.5 1.0 region slab Particle insertion: 224 every 3000 steps, 1000 by step 12001 fix 3 all enforce2d 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 250 dump.pour #dump 2 all image 500 image.*.jpg type type # zoom 1.6 adiam 1.5 #dump_modify 2 pad 5 #dump 3 all movie 500 movie.mpg type type # zoom 1.6 adiam 1.5 #dump_modify 3 pad 5 run 25000 Neighbor list info ... 2 neighbor list requests update every 1 steps, delay 0 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 1.2 ghost atom cutoff = 1.2 binsize = 0.6 -> bins = 167 84 2 Memory usage per processor = 2.75186 Mbytes Step Atoms KinEng c_1 Volume 0 0 -0 0 5000 1000 224 201.77464 0 5000 2000 224 372.42868 0 5000 3000 224 599.60221 0 5000 4000 448 1084.3752 0 5000 5000 448 1592.4543 0 5000 6000 448 1763.3153 37.274939 5000 7000 672 1805.206 89.331853 5000 8000 672 1778.0015 111.58381 5000 9000 672 1592.6805 97.550311 5000 10000 896 1631.962 76.905078 5000 11000 896 1590.1527 68.644552 5000 12000 896 1564.4707 60.119562 5000 13000 1000 1490.072 56.982972 5000 14000 1000 1266.4669 57.579371 5000 15000 1000 1034.1055 46.079374 5000 16000 1000 748.70283 43.933154 5000 17000 1000 462.65166 37.522016 5000 18000 1000 327.40481 29.450285 5000 19000 1000 223.24319 31.981144 5000 20000 1000 54.792042 23.449101 5000 21000 1000 25.805484 16.352528 5000 22000 1000 20.617765 11.756965 5000 23000 1000 16.115382 7.9435503 5000 24000 1000 11.971003 5.9411126 5000 25000 1000 10.404117 4.3811155 5000 Loop time of 3.00132 on 1 procs for 25000 steps with 1000 atoms Performance: 719682.852 tau/day, 8329.663 timesteps/s 99.9% CPU use with 1 MPI tasks x no OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 1.2554 | 1.2554 | 1.2554 | 0.0 | 41.83 Neigh | 0.4288 | 0.4288 | 0.4288 | 0.0 | 14.29 Comm | 0.0078034 | 0.0078034 | 0.0078034 | 0.0 | 0.26 Output | 0.0003016 | 0.0003016 | 0.0003016 | 0.0 | 0.01 Modify | 1.177 | 1.177 | 1.177 | 0.0 | 39.21 Other | | 0.132 | | | 4.40 Nlocal: 1000 ave 1000 max 1000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 2374 ave 2374 max 2374 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 2374 Ave neighs/atom = 2.374 Neighbor list builds = 2089 Dangerous builds = 0 Total wall time: 0:00:03