yade-users team mailing list archive

[Zhicheng Gao <question697642@xxxxxxxxxxxxxxxxxxxxx>]

Question #697642 on Yade changed:
https://answers.launchpad.net/yade/+question/697642

    Status: Solved => Open

Zhicheng Gao is still having a problem:
Dear Chareyre,
I modified the code in Section B in the way you suggested, but the same error still occurred. Thank you for your attention to this matter, I am honored and grateful for your support and eagerly look forward to your thoughts
This is the revised code in section B：
#B. Activate flow engine
flow.dead=0
flow.defTolerance=-1
flow.meshUpdateInterval=-1
flow.useSolver=0
flow.permeabilityFactor=1
flow.viscosity=10
flow.gradP=Vector3(0,1,0)
flow.duplicateThreshold=0.3
flow.updateTriangulation=True
O.run(2,1)
csdList=flow.getConstrictionsFull()
print(len(O.bodies),len(csdList),'finished')
flow.dead=1
print(O.bodies[10].shape.radius)
for i in range(10,90):
    O.bodies.erase(i)
print(len(O.bodies))
O.engines=O.engines[0:3]+[PeriodicFlowEngine(dead=1,label='flow2')]+O.engines[4:]
O.run(1000,True)
flow2.dead=0
flow2.defTolerance=-1
flow2.meshUpdateInterval=-1
flow2.useSolver=0
flow2.permeabilityFactor=1
flow2.viscosity=10
flow2.gradP=Vector3(0,1,0)
flow2.duplicateThreshold=0.3
flow2.updateTriangulation=True

O.run(2,1)
csdList1=flow2.getConstrictionsFull()
print(len(csdList1))

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