Hey everyone,
I have a simulation of a packing which works flowless on an older version of Liggghts but i only get 2 atoms inserted in the Version 3.0.x.
Here is the code:
#Simple chute wear test
atom_style granular
atom_modify map array
boundary p f f
newton off
processors 2 2 1
communicate single vel yes
units si
region reg block -0.0015 0.0015 -0.02 0.02 0.0 0.055 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.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
fix m5 all property/global characteristicVelocity scalar 2.
#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#---------% WALLS %-----------------------------------#
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane -0
fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane -0.02
fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 0.02
#fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane -0.0015
#fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 0.0015
fix zwallInsert all wall/gran model hertz tangential history primitive type 1 zplane 0.1
#---------% MESHES %-----------------------------------#
fix PlungerCAD all mesh/surface file meshes/Plunger_Face.stl type 1
#-------------distributions for insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.0002
fix pts2 all particletemplate/sphere 15485867 atom_type 1 density constant 2500 radius constant 0.004
fix pdd1 all particledistribution/discrete 32452843. 2 pts1 0.3 pts2 0.7
#----------------% OTHER REGIONS - INSERTION REGION %----------------------------#
group all region reg
region bc cylinder z 0.0 0.0 0.015 0.03 0.09 units box
region bc2 block -0.0015 0.0015 -0.02 0.02 0.03 0.05 units box
#----------------% INSERTING PARTICLES %------------------#
fix ins all insert/pack seed 32452867 distributiontemplate pdd1 insert_every once overlapcheck yes vel gaussian 0 0 -.4 0 0 .2 volumefraction_region 0.6 region bc2
#fix ins all insert/rate/region seed 5330 distributiontemplate pdd1 nparticles 17000 particlerate 70000 insert_every 100 overlapcheck yes vel gaussian 0 0 -.4 0 0 .2 region bc2
#-------------------apply nve integration to all particles that are inserted as single particles
#fix integr all nve/sphere
fix 1 all nve/sphere
#output settings, include total thermal energy
compute FC all pair/gran/local id pos force
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 100 post/dump*.Packing_3 id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dumpstl all mesh/stl 100 post/dump*.stl
dump fc all local/gran/vtk 100 post/forcechain*.vtk FC
cfdemLiggghtsPar in.Packing_3 4
#insert particles
run 18000 upto
unfix ins
#------------ Starting movement of the mesh ----------------#
fix Plunger_Wall all wall/gran model hooke tangential history mesh n_meshes 1 meshes PlungerCAD
fix MovePlunger all move/mesh mesh PlungerCAD linear 0.0 0.0 -0.3
run 24000 upto
run 27000 upto
Here is the command that i use to run the simulation:
cfdemLiggghtsPar in.Packing_3 4
Can someone detect the problem please?
Best Regards