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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp