forked from lijiext/lammps
141 lines
5.4 KiB
Groff
141 lines
5.4 KiB
Groff
LAMMPS (10 Jan 2012)
|
|
# Pour granular particles into chute container, then induce flow
|
|
|
|
atom_style sphere
|
|
boundary p p fm
|
|
newton off
|
|
communicate single vel yes
|
|
|
|
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)
|
|
1 by 1 by 1 MPI 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 = 9.46169 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.8532 47.325647 6600
|
|
5000 804 1598.6681 64.779741 6600
|
|
6000 804 1364.8231 61.75911 6600
|
|
7000 1206 1519.7287 69.028648 6600
|
|
8000 1206 1409.6162 62.604938 6600
|
|
9000 1206 1243.5344 53.025705 6600
|
|
10000 1608 1400.9581 57.452979 6600
|
|
11000 1608 1225.6036 52.631505 6600
|
|
12000 1608 1184.2735 45.623535 6600
|
|
13000 2010 1298.3639 50.266254 6600
|
|
14000 2010 1177.6607 41.899866 6600
|
|
15000 2010 986.18539 45.283438 6600
|
|
16000 2412 1032.5845 41.218715 6600
|
|
17000 2412 1024.6469 39.790283 6600
|
|
18000 2412 913.80876 40.351496 6600
|
|
19000 2814 1010.144 39.468447 6600
|
|
20000 2814 953.92496 40.689106 6600
|
|
21000 2814 762.62053 42.228087 6600
|
|
22000 2814 517.63012 38.535819 6600
|
|
23000 3000 424.8833 34.836058 6600
|
|
24000 3000 363.8636 21.972506 6600
|
|
25000 3000 210.08035 24.819644 6600
|
|
Loop time of 9.40419 on 1 procs for 25000 steps with 3000 atoms
|
|
|
|
Pair time (%) = 5.45799 (58.0379)
|
|
Neigh time (%) = 1.0322 (10.9759)
|
|
Comm time (%) = 0.250612 (2.66489)
|
|
Outpt time (%) = 0.0380349 (0.404446)
|
|
Other time (%) = 2.62535 (27.9168)
|
|
|
|
Nlocal: 3000 ave 3000 max 3000 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Nghost: 670 ave 670 max 670 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Neighs: 13595 ave 13595 max 13595 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
|
|
Total # of neighbors = 13595
|
|
Ave neighs/atom = 4.53167
|
|
Neighbor list builds = 1152
|
|
Dangerous builds = 0
|
|
|
|
unfix ins
|
|
fix 2 all gravity 1.0 chute 26.0
|
|
run 25000
|
|
Memory usage per processor = 10.2619 Mbytes
|
|
Step Atoms KinEng 1 Volume
|
|
25000 3000 210.08035 24.819644 6600
|
|
26000 3000 114.38566 17.150932 6600
|
|
27000 3000 104.67608 12.322911 6600
|
|
28000 3000 183.35846 11.491313 6600
|
|
29000 3000 304.77431 13.371079 6600
|
|
30000 3000 483.13283 15.775148 6600
|
|
31000 3000 698.65983 22.583167 6600
|
|
32000 3000 997.91493 29.132887 6600
|
|
33000 3000 1381.6582 36.071846 6600
|
|
34000 3000 1807.7933 39.766028 6600
|
|
35000 3000 2308.8013 49.301118 6600
|
|
36000 3000 2872.3529 55.638219 6600
|
|
37000 3000 3528.2521 68.672321 6600
|
|
38000 3000 4243.7668 78.210271 6600
|
|
39000 3000 5050.5727 86.453776 6600
|
|
40000 3000 5777.2063 90.664206 6600
|
|
41000 3000 6697.0059 105.74557 6600
|
|
42000 3000 7649.5485 121.09489 6600
|
|
43000 3000 8695.8277 129.0904 6600
|
|
44000 3000 9843.446 146.82081 6600
|
|
45000 3000 10998.776 143.2563 6600
|
|
46000 3000 12138.363 170.04079 6600
|
|
47000 3000 13412.445 174.29258 6600
|
|
48000 3000 14479.945 179.49372 6600
|
|
49000 3000 15616.206 182.18045 6600
|
|
50000 3000 16633.42 195.0737 6600
|
|
Loop time of 18.5002 on 1 procs for 25000 steps with 3000 atoms
|
|
|
|
Pair time (%) = 12.0726 (65.2567)
|
|
Neigh time (%) = 1.54973 (8.37683)
|
|
Comm time (%) = 0.540591 (2.92209)
|
|
Outpt time (%) = 0.0597088 (0.322747)
|
|
Other time (%) = 4.27755 (23.1217)
|
|
|
|
Nlocal: 3000 ave 3000 max 3000 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Nghost: 747 ave 747 max 747 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Neighs: 13861 ave 13861 max 13861 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
|
|
Total # of neighbors = 13861
|
|
Ave neighs/atom = 4.62033
|
|
Neighbor list builds = 807
|
|
Dangerous builds = 0
|