Alternative permeability values during the compaction
Hello everyone,
Hope you are doing well.
Regarding my previous question about the wrong permeability values during the compaction [1] which is not solved yet, I tried to import the packing in each step and use flowengine function to obtain the permeability.
But then I have received alternating values as well which is weird because the porosity decreases by passing the time.
I tried to follow the oedometer.py [2] for the modeling.
Should I change any boundary condition?
I am grateful for any hints or insights into this issue.
Here is my code:
###
from __future__ import print_function
from builtins import range
from yade import pack, export, ymport
young=1e6
compFricDegree = 3 # initial contact friction during the confining phase
finalFricDegree = 30 # contact friction during the deviatoric loading
mn,mx=Vector3(
O.materials.
O.materials.
walls=aabbWalls
wallIds=
spheres = ymport.
O.bodies.
O.engines = [
]
#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead = 0
flow.defTolerance = 0.3
flow.meshUpdate
flow.useSolver = 3
flow.permeabili
flow.viscosity = 10
flow.bndCondIsP
flow.bndCondValue = [0, 0, 1, 0, 0, 0]
flow.boundaryUs
O.dt = 0.1e-3
O.dynDt = False
O.run(1, 1)
Qin = flow.getBoundar
Qout = flow.getBoundar
permeability = abs(Qin) / 1.e-4 #size is one, we compute K=V/∇H
print("Qin=", Qin, " Qout=", Qout, " permeability=", permeability)
###
Many thanks in advance.
Best regards,
Soheil
[1] https:/
[2] https:/
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: