wall/gran/hertz/history frictionflag

Submitted by SergeiD on Tue, 05/03/2011 - 14:50

Hello.
What is mean of the frictionflag in wall/gran style?
I thought if the frictionflag is zero then the wall without friction. But the test showed that it is not.
In test script a sphere moves along yplane with given initial velocity. If yplane is frictionless then sphere's velocity should not be changed. But it is changing, see attached figs. And in different ways depending on the value of the flag.

And if I need frictional and frictionless walls should I to define 2 atomtype?
Thanks.

Test script:

#Particle packing by insertion and successive growing of particles

atom_style granular
atom_modify map array
boundary p m p
newton off
#echo both

communicate single vel yes

units si

#Simulation domain
region reg block -0.1 0.1 0. 0.1 -0.1 0.1 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 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/hertz/history 1 0 #Hertzian without cohesion
pair_coeff * *

timestep 0.00001

#Wall
fix ywalls all wall/gran/hertz/history 1 0 yplane 0. NULL 1

#Sphere
create_atoms 1 single 0 0.005 0 units box
set atom 1 diameter 0.01 density 2500 vx 0.1

#Gravity
fix gravi all gravity 9.81 vector 0.0 -1.0 0.0

#region
group nve_group 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 100 post/dump.wallfric id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
run 50000

AttachmentSize
Image icon frictionflag=0.png16.5 KB
Image icon frictionflag=1.png9.32 KB
ckloss's picture

ckloss | Tue, 05/03/2011 - 14:59

>>If yplane is frictionless then sphere's velocity should not be changed

Probably we should rename "frictionflag" into "dampflag" - b/c that's what it does - remove tangential damping

The sphere will have a tangential overlap with the wall, which results in a tangential spring force, and energy stored in both the tangential spring and the rotary motion.
The oscillation you see in the first picture (frictionflag=0) is the energy being transferred from the tangential spring to the kinetic energy back and forth without any damping. If frictionflag=1, this oscillation is damped and a steady state is reached.

Christoph

tmhealy | Mon, 04/09/2012 - 05:37

Related to the original question, I'm running a simulation where I'd like to have some walls with friction and others without. Since I've now learned that setting the friction flag to 0 will not accomplish the goal of getting a frictionless wall, I'm wondering how I could get a will without friction. Do I need to define an third atom type and set its coefficients of friction with the other atoms to zero? Is there an easier way than creating this "dummy" atom type?

I would appreicate any and all help people could offer.

Thanks and Regards,
Tim