Specify region- or group- dependant coefficients of friction

Submitted by alexander.polson on Fri, 01/14/2011 - 23:34

Hi Everybody

I wanted to set the friction to different values, depending on where the particle is in the model, e.g. in a hopper, and in the process below the hopper.

My approach was to:
- define a region, REG_hopper
- associate a group of particles with this region
- set the friction first for "all", and then for the group using "fix property/global"

The steps below shows only the applicable commands:
region REG_hopper block x_lo x_hi ... z_lo z_hi units box
...
group GRP_hopper region REG_hopper
...
fix m1 all property/global coefficientFriction peratomtypepair 1 ${mu_1}
fix m2 GRP_hopper property/global coefficientFriction peratomtypepair 1 ${mu_2}

When I follow the above, I get:
ERROR: Error in fix property/global. There is already a fix that registers a variable of the same name

Is there a better way to do this?

Thank you

Alexander

ckloss's picture

ckloss | Sat, 01/15/2011 - 03:24

Sorry, I don't think there is any general way to do it out-of the box right now. If you need such a feature it would have to be coded. Currently, the group is ignored for fix property/global.

The only thing you could do is split the simulation up into two parts (if that is feasible), i.e. let the hopper discharge be the boundary condition for the process below. Then you could choose the properties differently

Cheers,
Christoph