Homec4science
Diffusion lammps (master) Browse

lammps/examples/pour/log.pour.2d.9Jan12.linux.1
5026f4b7886bmaster

Last Change
Raw File

log.pour.2d.9Jan12.linux.1
View Options

LAMMPS (10 Jan 2012)
# Pour 2d granular particles into container
dimension 2
atom_style sphere
boundary f fm p
newton off
communicate single 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 4000.0 NULL 100.0 NULL 0 1 xplane 0 100
fix ywalls all wall/gran 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 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
run 25000
Memory usage per processor = 9.58 Mbytes
Step Atoms KinEng 1 Volume
0 0 -0 0 5000
1000 224 185.87695 0 5000
2000 224 348.08417 0 5000
3000 224 565.2162 0 5000
4000 448 1042.9716 0 5000
5000 448 1542.0614 0 5000
6000 448 1766.5136 41.192729 5000
7000 672 1846.4596 67.352592 5000
8000 672 1814.3429 92.898167 5000
9000 672 1662.3648 111.88805 5000
10000 896 1654.7992 91.210682 5000
11000 896 1539.2643 72.174442 5000
12000 896 1479.3613 60.304466 5000
13000 1000 1464.6099 54.360683 5000
14000 1000 1357.776 50.981395 5000
15000 1000 1129.2816 50.644045 5000
16000 1000 795.55925 47.94934 5000
17000 1000 448.80612 42.4465 5000
18000 1000 320.15236 34.231318 5000
19000 1000 150.74762 35.789901 5000
20000 1000 49.636313 19.317565 5000
21000 1000 24.980341 13.517257 5000
22000 1000 17.164368 9.7652983 5000
23000 1000 15.131329 7.6229469 5000
24000 1000 13.291283 5.8822567 5000
25000 1000 12.481182 3.8348041 5000
Loop time of 2.97817 on 1 procs for 25000 steps with 1000 atoms
Pair time (%) = 1.23279 (41.3941)
Neigh time (%) = 0.302755 (10.1658)
Comm time (%) = 0.0113213 (0.380143)
Outpt time (%) = 0.0628569 (2.11059)
Other time (%) = 1.36845 (45.9494)
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: 2377 ave 2377 max 2377 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 2377
Ave neighs/atom = 2.377
Neighbor list builds = 2001
Dangerous builds = 0
c4science · Help