yade-users team mailing list archive

[Robert Caulk <question698778@xxxxxxxxxxxxxxxxxxxxx>]

Question #698778 on Yade changed:
https://answers.launchpad.net/yade/+question/698778

Robert Caulk posted a new comment:
Hey Zoheir,

I don't think you are missing anything.  It seems that the pseudocode I
wrote here is incomplete. Looking a bit closer at the source, I am
unsure we can get the facet fluxes in python with the current wrapping.
Do you *really* need this in python? The only reason I suggested this
pseudocode in python was so that you could double-check the function
that I wrote in C++ [1]. Is [1] not working for you? Apparently you said
"Now, the influx and outflux are the same." which tells me there is no
need to be double-checking this in python.

Please clarify your needs. If you for some reason need the facet fluxes
to be wrapped in python, then I can do that - but it needs justification
first.

For some reason, I started chasing geese again on this problem, before
realizing that this was in fact solved, and not a problem, when you said
that you got influx and outflux to be the same using
flow.getCavityFlux()....oh well, that is my fault for not rereading the
entire thread. Anyways, during that goose chase, I visualized your
cavity with flow.saveVTK(). The cavity appears to be scattered and
incorrectly assigned, based on the VTK. Although, perhaps you have an
updated script where you made sure the cavity was properly defined. If
you did not change the cavity assignment, I would not expect any cavity
fluxes to be properly estimated.

-rc

[1]https://gitlab.com/yade-
dev/trunk/-/blob/master/lib/triangulation/FlowBoundingSphere.ipp#L1310

Cheers,

Robert

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.