Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Need help regarding 2d granular simulation through silo in LAMMPS

Please login with a confirmed email address before reporting spam

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



0 Replies Last Post 2017年3月26日 GMT-4 17:14
COMSOL Moderator

Hello Vivek Pankaj

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.