Particles penetrating walls

Submitted by ecmo on Thu, 03/22/2012 - 16:28

I could not resolve particle penetration to any walls. What mistake did I make? Below is a code that I worked on. This code will be coupled with the CFDEM.
Thanks in advance for your help.

===
echo both

atom_style granular
atom_modify map array
communicate single vel yes

#boundary f f p
boundary f f f
newton off

units si

region reg block -0.0002 0.0002 -0.0002 0.0002 0 0.0004 units box
create_box 1 reg

neighbor 0.00008 bin
neigh_modify delay 0 binsize 0.00008

#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 5.e7
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5

#pair style
pair_style gran/hertz/history 1 0 #Hertzian without cohesion
pair_coeff * *

#timestep, gravity
timestep 0.000001
#timestep 0.00001
#timestep 0.0000001

fix gravi all gravity 98.1 vector 0.0 0.0 1.0

#walls
fix xwalls all wall/gran/hertz/history 1 0 xplane -0.0002 0.0002 1
fix ywalls all wall/gran/hertz/history 1 0 yplane -0.0002 0.0002 1
fix zwalls all wall/gran/hertz/history 1 0 zplane 0 0.0004 1

#create single partciles
create_atoms 1 single 0.00015 0.00015 0.00001 units box
create_atoms 1 single 0.00015 0.00005 0.00001 units box
set atom 1 diameter 0.000008 density 1.5 vx 0 vy 0 vz 0
set atom 2 diameter 0.000008 density 1.5 vx 0 vy 0 vz 0

#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere #wenn das ausgeblendet, dann kein vel update

dump dmp all custom 100 post/dump.liggghts_init id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

run 100000 upto #80000 upto
===

ckloss's picture

ckloss | Thu, 03/22/2012 - 23:37

Hi,

have you tried the following:

+ use another unit system (e.g. cgs, your particle diameter is really small)
+ use fix check/timestep/gran to check your time-step - I guess it is too large for the particle size

Cheers, Christoph