Submitted by RocheDEM on Wed, 07/24/2013 - 18:57

I typically use the fix_nve integration (only updates position and linear velocity) for my granular simulations as this gives a good approximation of blocky (angular) material behaviour. I am wondering if the sjkr cohesion model can be used with this integrator or do I need to use fix_nve_sphere to get cohesion?

In the mean time, I am making a little test model.

Thank you,


aaigner's picture

aaigner | Thu, 07/25/2013 - 17:03


If you use the pair_gran styles, you must use atom style sphere. Using the fix_nve integration the angular velocity, which is introduced with the atom style sphere, of all particles is ignored. All particles have zero angular velocity for all time. That's good for your blocky material. You could also use a rolling friction model with a high coefficient.

The sjkr model adds only a normal, attractive force proportional to the (simplified) contact area times the cohesive energy.

To put it in a nutshell, yes you can.