Shop::unbalancedForce is not valid for capillary law

[Christian Jakob]

Hi,

I create my granular material via gravitational settlement process. As abort criterion I usually use unbalancedForce [1]. If capillary law is deactivated, this works well. Unbalanced force is very high at the begining and slowly decreases, when particles settle/calm down.
If capillary law is activated, unbalanced force is around 0,015 all the time. So some calculations are aborted too early and some will never be aborted, because crterion (< 0,01) is never reached. So I think unbalanced force is not valid in this case.

In the code [2] I see, that the force f that acts on body b is used to determine F_unb. Capillary force is included in f, but I think it should not be considered in this case, right?

Has someone an idea how to fix this?

[1] https://github.com/yade/trunk/blob/master/pkg/dem/Shop.cpp#L159
[2] https://github.com/yade/trunk/blob/master/pkg/dem/Shop.cpp#L171

[Bruno Chareyre (bruno-chareyre)]

Closing the bug. Please re-open and provide example script if the problem persists.
UnbalancedForce is defined exactly the same way with and without capillary forces, since the input values are the total forces on each body independently of their origin.