89 lines
3.2 KiB
Groff
89 lines
3.2 KiB
Groff
LAMMPS (1 Oct 2006)
|
|
# Pour 2d granular particles into container
|
|
|
|
dimension 2
|
|
atom_style granular
|
|
boundary p fm p
|
|
newton off
|
|
|
|
region reg block 0 100 0 50 -0.5 0.5 units box
|
|
create_box 1 reg
|
|
Created box = (0 0 -0.5) to (100 50 0.5)
|
|
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/hertzian 2000.0 50.0 0.0 0
|
|
timestep 0.001
|
|
|
|
fix 1 all nve/gran
|
|
fix 2 all gravity spherical 90.0 -90.0
|
|
fix xwalls all wall/gran xplane 0 100 50 0
|
|
fix ywalls all wall/gran yplane 0 NULL 50 0
|
|
|
|
region slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
|
|
fix ins all insert 1000 1 4767548 vol 0.4 10 diam 0.5 1.0 region slab
|
|
Particle insertion: 224 every 3000 steps, 1000 by step 12001
|
|
|
|
fix 3 all enforce2d
|
|
|
|
thermo_style granular
|
|
thermo 1000
|
|
thermo_modify lost ignore
|
|
|
|
dump id all atom 100 dump.pour.2d
|
|
run 25000
|
|
Memory usage per processor = 0.328064 Mbytes
|
|
Step Atoms Trans-KE Rot-KE Volume
|
|
0 0 0 0 5000
|
|
1000 224 348.20549 0 5000
|
|
2000 224 650.98655 0 5000
|
|
3000 224 1056.0039 0 5000
|
|
4000 448 1941.6817 0 5000
|
|
5000 448 2869.092 0 5000
|
|
6000 448 3310.8798 0 5000
|
|
7000 672 3505.2381 0 5000
|
|
8000 672 3325.8942 0 5000
|
|
9000 672 3165.2505 0 5000
|
|
10000 896 3170.2667 0 5000
|
|
11000 896 2944.0565 0 5000
|
|
12000 896 2862.9537 0 5000
|
|
13000 1000 2862.2101 0 5000
|
|
14000 1000 2638.5235 0 5000
|
|
15000 1000 2218.2739 0 5000
|
|
16000 1000 1583.5051 0 5000
|
|
17000 1000 963.53755 0 5000
|
|
18000 1000 666.8307 0 5000
|
|
19000 1000 339.73655 0 5000
|
|
20000 1000 100.52829 0 5000
|
|
21000 1000 56.112393 0 5000
|
|
22000 1000 47.013733 0 5000
|
|
23000 1000 32.615049 0 5000
|
|
24000 1000 25.853047 0 5000
|
|
25000 1000 17.216223 0 5000
|
|
Loop time of 22.8255 on 1 procs for 25000 steps with 1000 atoms
|
|
|
|
Pair time (%) = 9.2328 (40.4494)
|
|
Neigh time (%) = 2.19741 (9.627)
|
|
Comm time (%) = 0.301715 (1.32183)
|
|
Outpt time (%) = 0.97282 (4.26198)
|
|
Other time (%) = 10.1208 (44.3398)
|
|
|
|
Nlocal: 1000 ave 1000 max 1000 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Nghost: 27 ave 27 max 27 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
Neighs: 2441 ave 2441 max 2441 min
|
|
Histogram: 1 0 0 0 0 0 0 0 0 0
|
|
|
|
Total # of neighbors = 2441
|
|
Ave neighs/atom = 2.441
|
|
Neighbor list builds = 1758
|
|
Dangerous builds = 0
|