I am generating the particles in box and let them settle down. Once the particle settle down, I unfix the gravity earlier defined and using fix/freeze comand to stick the particle to the base so that I can have bumpy base. Now I tilt the gravity to 25 degree but freeze layer is also flowing Can any one please tell what's wrong with the my code?
#Simple chute wear test
atom_style granular
atom_modify map array
boundary p p f
newton off
communicate single vel yes
units si
region reg block 0.0 0.1 0.0 0.1 0.0 0.15 units box
create_box 1 reg
neighbor 0.001 bin
neigh_modify delay 0
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 5.e9
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.88
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
fix m5 all property/global characteristicVelocity scalar 0.0508
#New pair style
pair_style gran model hooke ktToKnUser on tangential_damping off
pair_coeff * *
timestep 0.00001
fix gravi all gravity 10 vector 0.0 0.0 -1.0
#the plane
fix zwall all wall/gran model hooke tangential history primitive type 1 zplane 0.0
#distributions for insertion
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 8000 radius constant 0.002
fix pdd1 all particledistribution/discrete 1 1 pts1 1.0
#region and insertion
region bottom block 0.0 0.1 0.0 0.1 0.0 0.004 units box
region flow block 0.0 0.1 0.0 0.1 0.004 0.15 units box
group flow1 region flow
group bottom1 region bottom
group nve_group region reg
#particle insertion
fix ins nve_group insert/rate/region seed 5101 distributiontemplate pdd1 &
nparticles 10000 massrate 10 insert_every 100 overlapcheck yes all_in no vel constant 0.0 0.0 0.0 &
region reg
#apply nve integration to all particles that are inserted as single particles
fix integr nve_group nve/sphere
#output settings, include total thermal energy
compute 1 all ke
thermo_style custom step atoms cpu ke
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
#insert the first particles so that dump is not empty
run 1
dump dmp1 all custom 100 post/dump*.xyz id type x y z vx vy vz
dump_modify dmp1 sort 1
dump pic all image 500 images/img*.ppm type type size 960 960 zoom 1.5 adiam 0.004
#insert particles
run 90000 upto
unfix ins
unfix gravi
fix bumpy bottom1 freeze
fix gravi2 flow1 gravity 10 chute 25
run 400000 upto
Regards
pnkumar
PG Student, Chemical Engineering
I.I.T. Kanpur