yade-users team mailing list archive

[Soheil Safari <question701670@xxxxxxxxxxxxxxxxxxxxx>]

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

Soheil Safari gave more information on the question:
Dear Bruno,

I tried to fix the DOFs and set all velocities to zero in the
permeability definition but still I give alternative values !!!

###
def get_permeability():
    for b in O.bodies:
        b.state.blockedDOFs = "xyzXYZ"
        b.state.vel = (0,0,0)
        Qin = flow.getBoundaryFlux(2)
        Qout = flow.getBoundaryFlux(3)
        permeability = abs(Qin) / 1.e-4  #size is one, we compute K=V/∇H
        print("Qin=", Qin, " Qout=", Qout, " permeability=", permeability)
        return permeability
###

I am grateful for any hints or insights into this issue.

Many thanks in advance.

Best regards,
Soheil

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