ball mill with walls made out of particles

#Contact model example

atom_style	granular
atom_modify	map array
boundary	f f p	#periodic in z
newton		off

communicate	single vel yes

units		si

region		reg block -0.51 0.51 -0.51 0.51 -0.01 0.51 units box
#create_box	2 reg

neighbor	0.02 bin
neigh_modify	delay 0

read_data wall.pour

#Material properties required for new pair styles

fix 		m1 all property/global youngsModulus peratomtype 1.e7 1e7
fix 		m2 all property/global poissonsRatio peratomtype 0.3 0.3
fix 		m3 all property/global coefficientRestitution peratomtypepair 2 0.2 0.2 0.2 0.2
fix 		m4 all property/global coefficientFriction peratomtypepair 2 0.3 0.3 0.3 0.3
fix 		m5 all property/global coefficientRollingFriction peratomtypepair 2 0.3 0.3 0.3 0.3


#New pair style
pair_style gran model hertz tangential history rolling_friction epsd2
pair_coeff	* *

timestep	1e-5

#output settings, include total thermal energy
thermo_style	custom step atoms spcpu
thermo			1000
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes
shell mkdir post
dump			dmp all custom 1000 post/dump*.liggghts id type x y z vx vy vz fx fy fz omegax omegay omegaz radius mass


group walls type 1
fix wall walls freeze
#rotate the drum
fix rot walls move rotate 0.0 0.0 0.0 0.0 0.0 1.0 5.0 #Px Py Pz Rx Ry Rz period

neigh_modify exclude type 1 1 #walls can not collide with walls

create_atoms 2 single 0 0 0.25 #ony free atom in center of system
group free type 2
set group free diameter 0.04 density 2800



#apply gravity and nve integration to free particles
fix			gravi free gravity 9.81 vector 0.0 -1.0 0.0
fix			integr free nve/sphere

run			100000 upto