yade-users team mailing list archive

[Mithushan Soundaranathan <question701228@xxxxxxxxxxxxxxxxxxxxx>]

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

Description changed to:
Hi,

I am measuring pore radius and pore volume using the TwoPhaseFlowEngine() for my packing of sphere (using ymport),  and initialization the twophaseflow. I am getting following error message when I run my code:
"negative volume for an ordinary pore (temp warning, should still be safe)"
*** stack smashing detected ***: terminated,

than Aborted (core dumped), and the whole thing crashes.

I was also trying with include TwoPhaseEngine in O.engine, and used
savePoreNetwork to extract pore data. With this method I am getting an
empty dataset for poredata.

Best regards,
Mithu

Here is my code:

import os
import numpy as np
import pandas as pd
from yade import pack, ymport
from yade import utils, plot, timing

O=Omega()
path_data='/home/mithushan/Swelling/PH101_20_1mm_removed_die/Particle_pos'
filename='PH101_particle_pos_20_1_2.0000009999999997.txt'

a=ymport.text(os.path.join(path_data,filename), shift=Vector3(0, 0, 0), scale=1.0)
O.bodies.append(a)

#O.engines
O.engines = [
  TwoPhaseFlowEngine(dead=0,label="flow"),
]

flow.getCellVoidVolume(100)
flow.getCellInSphereRadius(100)
flow.savePoreNetwork()


#Initialization
unsat = TwoPhaseFlowEngine()
unsat.initialization()
unsat.getCellVoidVolume(100)
unsat.getCellInSphereRadius(100)
unsat.savePoreNetwork()

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