Need help in cohesion effect

Submitted by laudari on Thu, 10/21/2021 - 06:53

Hello everyone,

I have a question related to LIGGGHTS.

I am trying to visualize the effect of cohesion in my simulation. Since my particles are not interacting to each other and they are falling to inclined plane one by one. I was hoping to see the effect of cohesion between particles and inclined plane (especially on particles bouncing position) but I couldn't see. That means cohesion model doesn't make any effect on bouncing position.

Right now I am using these commands

```
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.4
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.6
fix m4 all property/global coefficientFriction peratomtypepair 1 0.6
fix m5 all property/global k_finnie peratomtypepair 1 1.0
fix m6 all property/global coefficientRollingFriction peratomtypepair 1 0.02
fix m7 all property/global characteristicVelocity scalar 2. # Required for hooke model
fix m8 all property/global cohesionEnergyDensity peratomtypepair 1 300000
```

And cohesion model is

```
pair_style gran model hertz tangential history cohesion sjkr #Hertzian with cohesion"
```

Is there any idea to bring some effects of cohesion between particles and plane (On particles bouncing position)?

Any idea or suggestion would be appreciated.

Thank you

Thank you

Daniel Queteschiner | Thu, 10/21/2021 - 12:07

The pair_style command only defines interaction between particles. Did you specify the cohesion model for the fix wall/gran as well?