yade-users team mailing list archive

[Roxana Saghafian Larijani <question704035@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi all,
I am trying to continue my simulation in a rotating drum after I stop it. To do this I save my first simulation using O.save at the end of my first script. I also export vtu files for visualisation in Paraview. This works correctly. Then I O.load the saved file (test.bz2) at the beginning of the second script. The simulation goes on and gives me vtu files. However, in the vtu files some points are Nan and I do not get a proper visualisation in Paraview (it does not show the particles in paraview). I was wondering if you could help me with this issue. 

I have attached my two scripts here as well :
Script 1:

from yade import ymport
from yade import pack
from yade import utils, plot,wrapper
from yade import export



nCG=1
fr = 0.38
rho = 3700
En = 0.22
Et = 0.22
poi=0.3
yoM=200e6
r = 0.0005        


Gamma = 0.072    
Theta = 0
vB = 6.03  * 1e-9 * (nCG**3)
Eta= 0.001         
CapType="Rabinovich"
Tc=0.00045

##defining material

mat=O.materials.append(
        ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM, Vb=vB, gamma=Gamma, eta=Eta, theta=Theta, Capillar=True, CapillarType=CapType,en=En, et=Et,tc=Tc)
)
mat2=O.materials.append(
        ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM, Vb=vB, gamma=Gamma, eta=0, theta=90, Capillar=True, CapillarType=CapType, en=En, et=Et,tc=Tc)
)


#defining the spheres

sp=pack.SpherePack()
sp.makeCloud((-0.035,-0.035,-0.01),(0.035,0.035,0.01),psdSizes=[0.001,0.002,0.003,0.004], psdCumm=[0.0,0.8,0.95,1.0])
sp.toSimulation(material=mat)
Nprtcl=len(O.bodies)
print(Nprtcl)
Tt= utils.PWaveTimeStep()
O.dt = 0.5*Tt
print(O.dt)

#liquidMigration
VV=0.03
Vmin=0.0
for s in O.bodies:
        if not type(s.shape)==wrapper.Sphere:
                continue
        s.state.Vf=VV * (4/3) * 3.14*(s.shape.radius)**3
        s.state.Vmin=Vmin
##

Drum=geom.facetCylinder(material=mat2,center=(0.0,0.0,0.0),radius=0.05,height=0.03,orientation=Quaternion(Vector3(0,0,1),(pi/2.0)))
walls = O.bodies.append(Drum)


##engine
O.engines = [
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Facet_Aabb()]),
        InteractionLoop(
                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
                [Ip2_ViscElCapMat_ViscElCapMat_ViscElCapPhys()],
                [Law2_ScGeom_ViscElCapPhys_Basic()],
                
        ),
       
        NewtonIntegrator(gravity=[0, -9.8, 0]),
        RotationEngine(ids=walls,rotationAxis=[0,0,1],rotateAroundZero=True, zeroPoint=[0.0,0.0,0.0], angularVelocity=2.0),
        LiqControl(label='lqc'),
        VTKRecorder(iterPeriod=1000,recorders=['spheres','facets','colors'],fileName='test-')]




O.run( 500000,True)

O.save('test.bz2')



Script 2:

from yade import ymport
from yade import pack
from yade import utils, plot,wrapper
from yade import export


O.load('test.bz2')
Nprtcl=len(O.bodies)
print(Nprtcl)


O.run(15000,True)

###
BTW, I am using this version yade-2022-04-13.git-5bd3ade.

Thanks in advance for your respopnse!

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