Alternative permeability values during the compaction

Asked by Soheil Safari

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(0,0,0),Vector3(1,1,2) # corners of the initial packing

O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)

spheres = ymport.text('filename.txt', shift=Vector3(0, 0, 0), scale=1.0)
O.bodies.append(spheres)

O.engines = [
        FlowEngine(dead=1, label="flow"), #introduced as a dead engine for the moment, see 2nd section
]

#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead = 0
flow.defTolerance = 0.3
flow.meshUpdateInterval = 200
flow.useSolver = 3
flow.permeabilityFactor = 1
flow.viscosity = 10
flow.bndCondIsPressure = [0, 0, 1, 1, 0, 0]
flow.bndCondValue = [0, 0, 1, 0, 0, 0]
flow.boundaryUseMaxMin = [0, 0, 0, 0, 0, 0]
O.dt = 0.1e-3
O.dynDt = False

O.run(1, 1)
Qin = flow.getBoundaryFlux(5)
Qout = flow.getBoundaryFlux(4)
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://answers.launchpad.net/yade/+question/701670
[2] https://gitlab.com/yade-dev/trunk/-/blob/master/examples/FluidCouplingPFV/oedometer.py

Question information

Language:
English Edit question
Status:
Open
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:

Can you help with this problem?

Provide an answer of your own, or ask Soheil Safari for more information if necessary.

To post a message you must log in.