hydraulic gradient with PFV-DEM

Asked by azim

Dear all,
I have a non-periodic simulation with a cube filled with particles(V=c^3). I want to know the way we can exert hydraulic gradient to our model!! (DEM-PFV)
any help will be appreciated.
Azim

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Robert Caulk
Solved:
Last query:
Last reply:
Revision history for this message
Jérôme Duriez (jduriez) said :
#1

Hi,

Did you try using attributes such as FlowEngine.bndCondValue [*] ?
You may also check the example script [**] (mentioned in the above mentioned FlowEngine doc)

[*] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngine.bndCondValue
[**] https://github.com/yade/trunk/blob/master/examples/FluidCouplingPFV/oedometer.py

Revision history for this message
azim (mirzavand) said :
#2

Hi Jérôme,
thank for your reply.
I imposed pressure in [**] as:
flow.bndCondIsPressure=[0,0,1,1,0,0]
flow.bndCondValue=[0,0,a,0,0,0]
for having gradient equal to: i=a/(c*waterUnitWeight) where c is the length in the direction of flow.
and ran simulation for thousands of iterations but i can't see the movement of particles!
Am I doing something wrong??
Azim

Revision history for this message
Jérôme Duriez (jduriez) said :
#3

Your commands seem correct to me but I'm not really familiar with FlowEngine..

Anyway, on a general note, it is rare there is much too see in YADE qt.View() and I would rely on other indicators to check for a correct simulation definition.

Why not look instead at the plots defined in [**] ?

Revision history for this message
Bruno Chareyre (bruno-chareyre) said :
#4

Hi, the question in unclear.
I would say it is already answered by the example script pointed out by Jerôme [**], where hydraulic gradient is imposed.
Bruno

Revision history for this message
azim (mirzavand) said :
#5

Hi Jérôme,
I already looked at [**]. they prove the good performance of DEM+PFV.(Thanks to Bruno and Emanuele for this powerful engine).
but my question is about the movement of particles. what is the easiest way to check it??

Revision history for this message
azim (mirzavand) said :
#6

Hi Bruno,
thanks for your answering my question.
I will be more grateful if you answer my following questions.
1) is my expression (in #2) right for imposed gradient??
2) should we impose macroscopic hydraulic gradient? or Do we need to do something_ like calibration_ as we do for soil particles properties?? because my particle do not move even with many large gradients.
3) how can we watch flow in PARAVIEW??( if my question is irrelevant please tell me to open a new thread.)

Revision history for this message
Best Robert Caulk (rcaulk) said :
#7

Hello Azim,

1) Yes, you are imposing a gradient in the Y direction with your commands.

2) Flow engine boundary condition assignments are designed to impose "macroscopic hydraulic gradients" on cubical specimens. If you want to use this hydraulic gradient to calibrate the macroscopic permeability of your sample, you will want to adjust the permeability factor [1] and estimate the corresponding macroscopic permeability.

3) In an attempt to keep this relevant to this thread - if you want to see the hydraulic gradient in paraview, you will want to export pressures with flow.saveVTK [2].

Cheers,

Robert

[1]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngineT.permeabilityFactor
[2]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngineT.saveVtk

Revision history for this message
Chareyre (bruno-chareyre-9) said :
#8

To compliment 2) by Robert: if permeability factor is set negative it will
be used as (approximate) Darcy permeability. Otherwise permeability is
predicted based on micrstructure and fluid viscosity, then multiplied by
permeability factor in case it is not 1.
B

Le jeu. 9 août. 2018 17:17, Robert Caulk <
<email address hidden>> a écrit :

> Question #671365 on Yade changed:
> https://answers.launchpad.net/yade/+question/671365
>
> Status: Open => Answered
>
> Robert Caulk proposed the following answer:
> Hello Azim,
>
> 1) Yes, you are imposing a gradient in the Y direction with your
> commands.
>
> 2) Flow engine boundary condition assignments are designed to impose
> "macroscopic hydraulic gradients" on cubical specimens. If you want to
> use this hydraulic gradient to calibrate the macroscopic permeability of
> your sample, you will want to adjust the permeability factor [1] and
> estimate the corresponding macroscopic permeability.
>
> 3) In an attempt to keep this relevant to this thread - if you want to
> see the hydraulic gradient in paraview, you will want to export
> pressures with flow.saveVTK [2].
>
> Cheers,
>
> Robert
>
> [1]
> https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngineT.permeabilityFactor
> [2]
> https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FlowEngineT.saveVtk
>
> --
> You received this question notification because your team yade-users is
> an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>
>

Revision history for this message
azim (mirzavand) said :
#9

Hi Robert,
I really appreciate your help.
>> you will want to export pressures with flow.saveVTK [2].
I want to watch flow with grains together. VTK shows a solid with pressures, it is a very good engine, but not my need.
Thanks again
Azim

Revision history for this message
azim (mirzavand) said :
#10

Hi Bruno,
I think permeability set to 1 (Default) is my choice.
thanks for your helpful comments.
Azim

Revision history for this message
azim (mirzavand) said :
#11

Thanks Robert Caulk, that solved my question.