forked from lijiext/lammps
141 lines
5.4 KiB
Groff
141 lines
5.4 KiB
Groff
LAMMPS (15 Jan 2010)
|
|
# Pour granular particles into chute container, then induce flow
|
|
|
|
atom_style granular
|
|
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 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.14886 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.604955 6600
|
|
9000 1206 1243.5289 52.989305 6600
|
|
10000 1608 1403.0711 56.9465 6600
|
|
11000 1608 1225.8111 52.938639 6600
|
|
12000 1608 1180.1935 40.544973 6600
|
|
13000 2010 1328.8958 49.216222 6600
|
|
14000 2010 1178.4104 44.431482 6600
|
|
15000 2010 998.18775 45.388975 6600
|
|
16000 2412 1050.7228 44.704814 6600
|
|
17000 2412 1025.5915 40.030825 6600
|
|
18000 2412 919.51551 36.650545 6600
|
|
19000 2814 1006.5829 41.592546 6600
|
|
20000 2814 964.34863 40.248622 6600
|
|
21000 2814 755.00069 47.423444 6600
|
|
22000 2814 516.03089 37.374261 6600
|
|
23000 3000 412.07349 38.529474 6600
|
|
24000 3000 352.98971 25.413823 6600
|
|
25000 3000 224.4892 26.293648 6600
|
|
Loop time of 16.4411 on 1 procs for 25000 steps with 3000 atoms
|
|
|
|
Pair time (%) = 10.0191 (60.9391)
|
|
Neigh time (%) = 1.58554 (9.64374)
|
|
Comm time (%) = 0.336891 (2.04908)
|
|
Outpt time (%) = 0.072582 (0.441467)
|
|
Other time (%) = 4.42704 (26.9266)
|
|
|
|
Nlocal: 3000 ave 3000 max 3000 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Nghost: 675 ave 675 max 675 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Neighs: 13593 ave 13593 max 13593 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
|
|
Total # of neighbors = 13593
|
|
Ave neighs/atom = 4.531
|
|
Neighbor list builds = 1154
|
|
Dangerous builds = 0
|
|
|
|
unfix ins
|
|
fix 2 all gravity 1.0 chute 26.0
|
|
run 25000
|
|
Memory usage per processor = 9.91094 Mbytes
|
|
Step Atoms KinEng 1 Volume
|
|
25000 3000 224.4892 26.293648 6600
|
|
26000 3000 109.26365 18.701584 6600
|
|
27000 3000 105.11897 12.063416 6600
|
|
28000 3000 184.91344 10.98296 6600
|
|
29000 3000 297.3293 13.247594 6600
|
|
30000 3000 470.52936 16.158596 6600
|
|
31000 3000 679.00541 23.059875 6600
|
|
32000 3000 930.87174 27.785348 6600
|
|
33000 3000 1335.2762 36.678555 6600
|
|
34000 3000 1757.0303 41.554334 6600
|
|
35000 3000 2279.8561 52.203628 6600
|
|
36000 3000 2916.5804 59.514974 6600
|
|
37000 3000 3520.1004 65.363759 6600
|
|
38000 3000 4240.5237 74.969118 6600
|
|
39000 3000 5064.0821 78.067916 6600
|
|
40000 3000 5827.5787 100.77169 6600
|
|
41000 3000 6790.2326 112.89539 6600
|
|
42000 3000 7867.015 119.80944 6600
|
|
43000 3000 8898.1328 123.18993 6600
|
|
44000 3000 9963.0265 140.481 6600
|
|
45000 3000 11240.984 167.23621 6600
|
|
46000 3000 12461.112 171.90559 6600
|
|
47000 3000 13796.683 185.03283 6600
|
|
48000 3000 15197.085 209.71877 6600
|
|
49000 3000 16699.196 221.86029 6600
|
|
50000 3000 18215.728 216.18732 6600
|
|
Loop time of 32.2512 on 1 procs for 25000 steps with 3000 atoms
|
|
|
|
Pair time (%) = 21.9871 (68.1745)
|
|
Neigh time (%) = 2.34306 (7.26503)
|
|
Comm time (%) = 0.715272 (2.21782)
|
|
Outpt time (%) = 0.118213 (0.366539)
|
|
Other time (%) = 7.08756 (21.9761)
|
|
|
|
Nlocal: 3000 ave 3000 max 3000 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Nghost: 751 ave 751 max 751 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Neighs: 14038 ave 14038 max 14038 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
|
|
Total # of neighbors = 14038
|
|
Ave neighs/atom = 4.67933
|
|
Neighbor list builds = 802
|
|
Dangerous builds = 0
|