yade-users team mailing list archive

[Jibril Coulibaly <question685707@xxxxxxxxxxxxxxxxxxxxx>]

New question #685707 on Yade:
https://answers.launchpad.net/yade/+question/685707

Dear Yade community,

I am trying to perform the seemingly simple computation of the permeability of a periodic assembly using the PFV module. From that previous question (https://answers.launchpad.net/yade/+question/662105), I gather that I need to define a macroscopic pressure gradient, as well as to fix the pressure in a given cell.

Trying to do so, Yade complains "setCellPImposed" and "setCellPressure" being out of range with max value 0 and returns wrong values for my boundary flows. Following is my input script, I would very much appreciate your help on figuring this out:

###############################################
from yade import pack,Vector3

O.periodic=True
O.cell.hSize=Matrix3(0.1,0,0, 0,0.1,0, 0,0,0.1)
sp=pack.SpherePack()
radius=5e-3
num=sp.makeCloud((0,0,0),(.1,.1,.1),radius,.2,200,periodic=True)
O.bodies.append([sphere(s[0],s[1]) for s in sp])

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb()],verletDist=.05*radius),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
	NewtonIntegrator(damping=0.),
	PeriodicFlowEngine(dead=0,label="flow")
]

O.dt=0.1e-8
O.dynDt=False

flow.defTolerance=0.3
flow.meshUpdateInterval=300
flow.useSolver=0 # Solver >0 not supported for PBC
flow.permeabilityFactor=1
flow.viscosity=1
flow.gradP=Vector3(1,0,0)
flow.setCellPImposed(0,True)
flow.setCellPressure(0,0)

O.run(1,1)
Qin = flow.getBoundaryFlux(0)
Qout = flow.getBoundaryFlux(1)
print("Qin=",Qin," Qout=",Qout)
###############################################

Thank you,
Jibril B. Coulibaly

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