Discussion Forum

Please help the particle-particle interaction is not working.
Pasted below is the LAMMPS code.


# granular simulation through silo

units cgs
atom_style sphere
dimension 2
boundary f fm p
newton off
comm_modify vel yes


lattice sq 0.05
region p1 block 0 30 0 50 -0.5 0.5
region p2 block 1 28 16 50 -0.5 0.5
region p3 block 14 16 14 15 -0.5 0.5
region p4 block 1 28 1 14 -0.5 0.5
region p34 union 2 p3 p4

create_box 2 p1

create_atoms 1 box
delete_atoms region p34

group crystal region p1
group delcrystal region p34
group crystal1 region p2
group net subtract crystal delcrystal crystal1
#group net include 1

neighbor 0.2 bin
neigh_modify delay 0

#set type 1 crystal1

pair_style gran/hooke/history 1000000 NULL 1000.0 NULL 0.8 1
pair_coeff * *

fix 1 crystal1 nve/sphere
fix 2 crystal1 gravity 981 spherical 0.0 -180.0
#fix 4 net setforce 0 0 0
fix xwalls crystal1 wall/gran 2000000 NULL 500.0 NULL 0.3 1 &
xplane 0 30
fix ywalls crystal1 wall/gran 2000000 NULL 500.0 NULL 0.3 1 &
yplane 16 NULL

#region slab block 5 25 18 47 -0.5 0.5
#fix ins 1 pour 1000 1 4767548 vol 0.4 10 &
# diam one 0.05 region slab

fix 3 crystal1 enforce2d


timestep 0.0001

compute 1 crystal1 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 crystal atom 250 dump.lammpstrj

#dump 2 all image 500 image.*.jpg type type &
# zoom 1.6 adiam 1.5
#dump_modify 2 pad 5

#dump 3 all movie 500 movie.mpg type type &
# zoom 1.6 adiam 1.5
#dump_modify 3 pad 5

run 25000