Yade

How can I get local fluxes in flowEngine ?

Asked by Antoine Wautier

Hi,

I am interested in getting acces to the local fluxes (or equivalent radii of the constriction according to the PFV model) between two cells in FlowEngine. Is there a way to do that ?

Thanks in advance
Antoine

PS: I would like in the end to get an estimation of the fluid velocity between two cells of the triangulation defined with a FlowEngine. I found I can get the mean relative fluid velocity around a given grain with flow.avFlVelOnSph(id) however I am interested in characterizing the pore network, so I would prefer to get a pore centered quantity.

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Bruno Chareyre
Solved:
Last query:
Last reply:
Revision history for this message
Best Bruno Chareyre (bruno-chareyre) said : 		#1

There is a c++ function which assign the average velocity of each cell.
FlowBoundingSphere.ipp:96
void FlowBoundingSphere<Tesselation>::averageRelativeCellVelocity()

It has no python wrapping though, but maybe you can manage to do that? (I'll not be able to help on that at the moment)
Note that velocity is not a variable of the problem, only fluxes are really computed. This velocity is a post-processed quantity.
Also note that it is the velocity relative to the solid phase, so that it is usually weird (since one intuitively expect to see absolute velocities in a fixed reference frame).

Bruno

Revision history for this message
Antoine Wautier (antoine-wautier) said : 		#2

Thanks for your explanations.

I am not very used to C++ for now, so I will not be able to do the python wrapping soon, but I may try later. I'll keep this in mind.

Revision history for this message
Antoine Wautier (antoine-wautier) said : 		#3

Thanks Bruno Chareyre, that solved my question.