yade-users team mailing list archive

[Andres Escobar <question697783@xxxxxxxxxxxxxxxxxxxxx>]

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

Hello, 

I recently added periodic boundaries to my simulation, in which I have a granular flow on an inclined plane. For this reason I have a reservoir that is  filled with spheres so my incoming flow is steady. Previously I used the following boxfactory command on the engines to generate spheres as observed here:

BoxFactory(maxParticles=-1,center=(lr/2,0,(hr+hc)/2),extents=(lr/2,wc/2,(hc-hr)/2),rMin=d-disp,rMax=d+disp,vMin=0.5,vMax=0.5,vAngle=0,massFlowRate=1000,normal=(0,0,-1),materialId=sphereMat,stopIfFailed=False,silent=True,maxAttempt=3,label='BoxFact')

Yet I obtain the following error:
In [1]: <FATAL ERROR> ThreadRunner:35 void yade::ThreadRunner::run(): Exception occured: 
InsertionSortCollider::probeBoundingVolume: handling periodic boundary not implemented.

I skipped this problem by creating a function that checks if the reservoir is getting empty every few steps, if it's empty I would manually add more. However, since I need to ensure a constant flow I need to add a huge quantity of spheres(Otherwise this would influence the results), making the script very inefficient. 

Based on this, does an alternative to this command exists that could work in periodic conditions? 
Thank you very much in advance for any help or suggestions that could exist! 





If it helps, I assume my current version of yadedaily wouldn't be a problem. This is my version as for the startup message:
Welcome to Yade 20210306-5036~03fcf2d~focal1 
Using python version: 3.8.5 (default, May 27 2021, 13:30:53) 

Since the error message involves "insertionSortcollider" here's my engines setup:

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()],allowBiggerThanPeriod=True),
   InteractionLoop(
#  handle sphere+sphere and box+sphere collisions
   [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
   [Ip2_ViscElMat_ViscElMat_ViscElPhys()],
   [Law2_ScGeom_ViscElPhys_Basic()]),
   NewtonIntegrator(gravity=(g*sin(slope),0,-g*cos(slope)), damping=0), 
]

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