Dear all,
I cannot find the bug in my script. Please help me, because I don't understand, why there doesn't appear any repulsion of my particels after falling on the ground.
Here is my little script:
dimension 2
atom_style granular
atom_modify map array
boundary m m p
newton off
echo log
communicate single vel yes
units si
region reg block -0.03 0.03 -0.03 0.03 -0.0002 0.0002 units box
create_box 1 reg
neighbor 0.002 bin
neigh_modify delay 0
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 30e9
fix m2 all property/global poissonsRatio peratomtype 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.05
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
#New pair style
pair_style gran/hertz/history 1 0
pair_coeff * *
timestep 0.00001
fix xwalls all wall/gran/hertz/history 1 0 xplane -0.3 0.3 1
fix ywalls all wall/gran/hertz/history 1 0 yplane -0.3 0.3 1
fix gravi all gravity 9.81 vector 0.0 -1.0 0.0
#region and insertion
group nve_group region reg
fix ins nve_group pour/legacy 500 1 1 vol 0.7 1000 diam 0.03 0.03 dens 2500 2500 vel 0. 0. -0.3 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 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
run 1
dump dmp all custom 350 post/dump.oW.* id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
run 50000 upto
unfix ins
raguelmoon | Wed, 04/20/2011 - 19:43
Hi, Value of young's modulus
Hi,
Value of young's modulus is very high according to your time step.
Best,
Ram
Silias | Wed, 04/20/2011 - 20:04
Problem solved
Hi,
thanks, I checked timestep and young's modulus, but not the coefficient of restitution...
With 0.05 it was too low, so that there was no restitution - why has logic to be always so hard to understand... :/
kind regards :)
Silias