What is wrong with this script file?

Submitted by pnkumar on Mon, 12/21/2015 - 08:35

I am trying to run this script file but it is showing that velocity projection on face normal is less than 1e-3.

Script File

#Simple chute wear test

atom_style granular
atom_modify map array
boundary f f f
newton off

communicate single vel yes

units si

region reg block 0. 0.3 -0.05 0.05 0.0 0.2 units box # 30cm x 10cm x 20cm
create_box 1 reg

neighbor 0.003 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5

#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *

timestep 0.00001

fix gravi all gravity 9.81 chute 25

#the chute
fix cad all mesh/surface file meshes/trial.stl type 1
fix inface all mesh/surface file meshes/trial2.stl type 1
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad

#distributions for insertion
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.002
fix pdd1 all particledistribution/discrete 1 1 pts1 1.0

#region and insertion
group nve_group region reg
region bc block 0 0.02 -0.02 0.02 0.05 0.15 units box

#particle insertion
fix ins nve_group insert/stream seed 5330 distributiontemplate pdd1 &
nparticles 6000 massrate 0.1 insert_every 1000 overlapcheck yes all_in no vel constant 0.0 0.0 -1.0 &
insertion_face inface

#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 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

#insert the first particles so that dump is not empty
run 1
dump dmp all custom 1000 post/dump*.txt x y z vx vy vz
#dump_modify dmp sort 1
#dump_modify dmp sort 1 element S1 every 1000

dump pic all image 1000 dump.*.ppm type type size 1024 1024 zoom 1.5 adiam 0.0015 adiam 0.0025
#dump stress all mesh/gran/VTK 1000 post/dump*.vtk stress wear cad
#dump mov all movie 1000 movie.m4v type type size 1280 720

#insert particles
run 10000 upto
unfix ins

and I attached the mesh file which is in txt extension just change the extension to stl.

AttachmentSize
Plain text icon inclined plane stl file282.02 KB
Plain text icon insertion face stl file16.44 KB