Homec4science
Diffusion lammps (lammm-master-deprecated) Browse

lammps/examples/pour/log.15May15.pour.2d.molecule.g++.4
42ff3a7abafflammm-master-deprecated

Last Change
Raw File

log.15May15.pour.2d.molecule.g++.4
View Options

LAMMPS (30 Apr 2015)
# Pour 2d granular particles into container
dimension 2
atom_style sphere
atom_modify map array
boundary f fm p
newton off
comm_modify vel yes cutoff 2.5
fix prop all property/atom mol ghost 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)
4 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 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
molecule object molecule.vshape
Read molecule object:
5 atoms with 1 types
0 bonds with 0 types
0 angles with 0 types
0 dihedrals with 0 types
0 impropers with 0 types
fix 3 all rigid/small molecule mol object
0 rigid bodies with 0 atoms
2.23607 = max distance from body owner to body atom
# insure region size + molecule size does not overlap wall
region slab block 3.0 97.0 30 34.5 -0.5 0.5 units box
fix ins all pour 500 0 4767548 vol 0.8 10 region slab mol object rigid 3
Particle insertion: 26 every 3000 steps, 500 by step 57001
fix 4 all enforce2d
compute 1 all erotate/sphere
compute Tsphere all temp/sphere
thermo_style custom step atoms ke c_1 vol
thermo_modify lost ignore norm no temp Tsphere
compute_modify Tsphere dynamic yes
thermo 1000
#dump id all atom 100 tmp.dump
#variable colors string # "red green blue yellow white # purple pink orange lime gray"
#variable mol2 atom mol%10
#dump 2 all image 250 image.*.jpg v_mol2 type # zoom 1.6 adiam 1.5
#dump_modify 2 pad 5 amap 0 10 sa 1 10 ${colors}
#variable colors string # "red green blue yellow white # purple pink orange lime gray"
#variable mol3 atom mol%10
#dump 3 all movie 250 movie.mpg v_mol3 type # zoom 1.6 adiam 1.5
#dump_modify 3 pad 5 amap 0 10 sa 1 10 ${colors}
run 25000
Neighbor list info ...
2 neighbor list requests
update every 1 steps, delay 0 steps, check yes
master list distance cutoff = 1.2
Memory usage per processor = 2.72041 Mbytes
Step Atoms KinEng 1 Volume
0 0 -0 0 5000
1000 130 259.24123 2.3773387e-30 5000
2000 130 465.39775 9.5010816e-30 5000
3000 130 736.55426 2.1365846e-29 5000
4000 260 1274.295 6.045774e-05 5000
5000 260 1855.4991 0.0021687846 5000
6000 260 1936.8934 6.9705094 5000
7000 390 1909.1519 11.402993 5000
8000 390 1958.2549 10.910611 5000
9000 390 1968.6293 14.080709 5000
10000 520 2278.4559 8.2512574 5000
11000 520 1952.4436 6.8525602 5000
12000 520 1391.9826 3.6742003 5000
13000 650 1608.8924 2.2123232 5000
14000 650 1562.872 2.9306477 5000
15000 650 1548.2428 2.4495128 5000
16000 780 1712.8658 2.2773941 5000
17000 780 1649.1423 4.9833774 5000
18000 780 1421.6645 4.0064802 5000
19000 910 1574.37 2.4833854 5000
20000 910 1354.8838 5.1248395 5000
21000 910 1156.7385 2.5239056 5000
22000 1040 1233.7916 2.9364978 5000
23000 1040 1150.2238 2.5059515 5000
24000 1040 1224.372 1.6145537 5000
25000 1170 1202.871 2.6775319 5000
Loop time of 1.09127 on 4 procs for 25000 steps with 1170 atoms
Pair time (%) = 0.230632 (21.1343)
Neigh time (%) = 0.070335 (6.44524)
Comm time (%) = 0.0610123 (5.59095)
Outpt time (%) = 0.00059557 (0.0545758)
Other time (%) = 0.728695 (66.7749)
Nlocal: 292.5 ave 304 max 278 min
Histogram: 1 0 1 0 0 0 0 0 0 2
Nghost: 42.5 ave 58 max 27 min
Histogram: 1 1 0 0 0 0 0 0 1 1
Neighs: 427 ave 456 max 388 min
Histogram: 1 0 0 1 0 0 0 0 0 2
Total # of neighbors = 1708
Ave neighs/atom = 1.45983
Neighbor list builds = 1732
Dangerous builds = 0
c4science · Help