← Back to team overview

yade-dev team mailing list archive

[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