yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #23166
Re: [Question #690803]: bndCondIsPressure
Question #690803 on Yade changed:
https://answers.launchpad.net/yade/+question/690803
Guilherme das Neves Seguro posted a new comment:
Hello mr Scholtès;
it is the same code I've shared in [1]. I'm having trouble in measuring
the permeability and I was concerned that it might be because of a wrong
setting the Pressure Condition.
The main code is shared in [2] (by now I corrected my mistake of making
it open to everyone!) and this current question is on this part of the
code:
PRESS=0.0
[...]
bottom=0
[then the code calls a secondary script and comes to]
### Flow Engine ###
flow=DFNFlowEngine(
isActivated=False
,useSolver=3 # (0: Gauss Seidel, 1: Taucs, 2: Pardiso, 3: CHOLMOD)
#,boundaryUseMaxMin=[0,0,0,0,0,0] # [left, right, bottom, top, back, front]: False means boundary made with walls
,bndCondIsPressure = [0,0,bottom,0,0,0] # bndCondIsPressure(=vector<bool>(6, false))
# bndCondIsPressure=[left, right, bottom, top, back, front]
# ,bndCondValue=[0,0,0,0,PRESS,0]
,bndCondValue=[0,0,PRESS,0,0,0] # bndCondValue(=vector<double>(6,0))
,permeabilityFactor=pFactor
,viscosity=visc
,fluidBulkModulus=KFluid
### DFN related
,clampKValues=False
,jointsResidualAperture=slotAperture
)
[fluid injection occurs and after that I need to impose the pressure on
y axis as I asked. So I change the parameters to]
bottom = 1
PRESS = 3e6
O.run(1,True)
# getBoundaryFlux get the total discharge [m3/s]
Qin = flow.getBoundaryFlux(2)
Qout = flow.getBoundaryFlux(3)
# if Qout is the total discharge, we can compute k=Q*nu*Length/(Area*(Pout-Pin))
# if Qout is the flux, we can compute k=Q*nu*Length/(Pout-Pin) -> getFlux gives total discharge -> Qout (m3/s)!
permeability = abs(Qout)*flow.viscosity*Y/(X*Z) # !!! if Pout=1, Pin=0
permeability2 = flow.averageVelocity()*flow.viscosity*Y # !!! if Pout=1, Pin=0
conductivity = permeability*DENS_FLUID*9.82/flow.viscosity # K=rho*g*k/nu
print "\n\nQin=",Qin,"\nQout=",Qout,"\nOBS: ARE THEY EQUAL? IF NOT => NO FLOW!\n"
print "\n\nPermeability [m2]=",permeability,"\nHydraulic conductivity [m/s]=",conductivity, '\n\nTHE END!\n'
[1] https://answers.launchpad.net/yade/+question/690573
[2] https://drive.google.com/drive/folders/1G-
rrLStPBemmjIBDlGO_ZVMupDyzDO5b?usp=sharing
Thank you very much for all your support, guys!
--
You received this question notification because your team yade-users is
an answer contact for Yade.