yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #10045
[Bug 1233727] [NEW] Shop::unbalancedForce is not valid for capillary law
Public bug reported:
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
** Affects: yade
Importance: Undecided
Status: New
--
You received this bug notification because you are a member of Yade
developers, which is subscribed to Yade.
https://bugs.launchpad.net/bugs/1233727
Title:
Shop::unbalancedForce is not valid for capillary law
Status in Yet Another Dynamic Engine:
New
Bug description:
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
To manage notifications about this bug go to:
https://bugs.launchpad.net/yade/+bug/1233727/+subscriptions
Follow ups
References