Question #688685 on Yade changed:
https://answers.launchpad.net/yade/+question/688685
Status: Answered => Open
Othman Sh is still having a problem:
Hi Robert,
I have added a line to my code to save vtk (flow.saveVtk()) and run Paraview
but then I get the error (Triangulation does not exist. Sorry.). Then I noticed
I didn't have flow.emulateAction() in my code so probably that's why there is
no triangulation? When I add emulateAction() I get this:
negative volume for an ordinary pore (temp warning, should still be safe)
*** stack smashing detected ***: <unknown> terminated
Aborted (core dumped)
Another question, can I use periodic flow engine for non-periodic packing?
Thanks,
Othman
--------------------------------------------
# -*- coding: utf-8 -*-
import numpy as np
from yade import pack, ymport, plot
radiuscyl=.08
heightcyl=.25
dP=4e3 #Pa
visc=1e-3 #Pa.sec
density=1000 #kg/m3
g=9.81 #m/s2
O.materials.append(FrictMat(young = 5e10, poisson = 0.15,frictionAngle =
atan(.2), density=1920))
############################ spheres #############################
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(.4,.4,2),rMean=0.0083,rRelFuzz=0.1)
#### cylinder extraction
pred=pack.inCylinder((.3,.3,0),(.3,.3,heightcyl),radiuscyl)
spFilter=filterSpherePack(pred,sp,Material=Material, returnSpherePack=True)
spFilter.toSimulation()
yade.qt.View()
Height=max(utils.aabbDim())
#Fix all particles in their positions. No deformation
for i in O.bodies:
i.state.blockedDOFs='xyzXYZ'
########## Activate flow engine and set boundary conditions in order to get
permeability #####
flow=FlowEngine()
#flow.dead=0
flow.useSolver=3
flow.viscosity=visc
flow.boundaryUseMaxMin=[1,1,1,1,1,1]
flow.bndCondIsPressure=[0,0,0,0,1,1] #[xmin,xmax,ymin,ymax,zmin,zmax]
flow.bndCondValue=[0,0,0,0,0,dP]
flow.emulateAction()
flow.saveVtk()
#Permeability calculations
Q=max(flow.averageVelocity())
Permeability=abs((density*g*Q*Height)/dP)
print ("average velocity: ",Q, "Permeability: ",Permeability)
[1] https://answers.launchpad.net/yade/+question/685707
--
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
