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.

_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp

Reply via email to