yade-users team mailing list archive

[rhaven <question670610@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,
I am trying to tessellate particle geometry generated with Yade. For this purpose I am using the TwoPhaseFlowEngine. 
When I run the TwoPhaseFlowEngine on geometry generated directly from makeCloud everything runs smoothly. However, when I try to run the TwoPhaseFlowEngine on geometry which I have generated to simulate an aggregate. I receive a segmentation fault. In order to track down this issue, I have written the code below. 
This code reads a packing of 113 spheres. If I run the TwoPhaseFlowEngine at this point, there is no error, however looking at the exported pore network, only 8 pores are generated by the tesselation, and they are all located at the same location (see here: http://pasteall.org/1012264 )

In order to try and dig into the issue. I wrote the following code. First I load the sphere geometry. I try to run the tesselation. I then save the individual sphere locations and radii and add them back in one at a time to the simulation. At every step I try to run the tesselation. In this case, when I add the 5th particle I receive a Segmentation fault (core dumped) and the contents of PoreBodyLocation.txt is a list of 0,0,0 's . 

from yade import export,ymport

sp = ymport.textClumps('/tmp/clump0.txt')
print('Number of particles = ' , len(O.bodies))

TwoPhase=TwoPhaseFlowEngine()
TwoPhase.initialization()
TwoPhase.savePoreNetwork('networkTest')
#raw_input('Pause')

spheres = []
for b in O.bodies:
	if b.isClump == False:
		spheres.append([b.state.pos, b.shape.radius])

O.resetThisScene()

for pos,r in spheres:
	O.bodies.append(sphere(pos,r))
	print('Number of particles = ' , len(O.bodies))
	#raw_input('Pause')
	TwoPhase.initialization()
	TwoPhase.savePoreNetwork('networkTest')
	
try:
	from yade import qt
	qt.View()
except:
	pass

and clump0.txt at http://pasteall.org/1012270

Does anyone know the reason for this segmentation fault, and or the reason for the pore locations = 0,0,0 ?

many thanks in advance
Jesse 


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