← Back to team overview

yade-users team mailing list archive

Re: [Question #232941]: Capillary pressure between two particles

 

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

Seungcheol Yeom posted a new comment:
Hello all,

First of all, thanks for all your response. I have learned a lot from you.
By the way, I have modified the script as shown below:

 from yade import plot,qt

r = 1e-4 #particle radius
h = 1e-5 #praticle distance

#create two sphere paticles#
O.bodies.append([
   utils.sphere(center=(0,0,0),radius=r,fixed=False),
   utils.sphere((0,0,2*r+h),r)
])

#define engines#
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom()],         
      [Ip2_FrictMat_FrictMat_CapillaryPhys()], 
      [Law2_ScGeom_FrictPhys_CundallStrack(neverErase=True)]   
   ),
   Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=10000,label='Cap'),
   NewtonIntegrator(damping=0.4,gravity=(0,0,0)),
   PyRunner(iterPeriod=100,command='addPlotData()') 
]

Cap.createDistantMeniscii=True
O.run(1,True)
Cap.createDistantMeniscii=False

def addPlotData():
    plot.addData(i=O.iter,b0displ=O.bodies[0].state.displ(),b1displ=O.bodies[1].state.displ())

plot.plots={'i':('b0displ')}
plot.plot()
plot.live=True
plot.autozoom=True
O.dt=0.5*PWaveTimeStep()
qt.View()
O.saveTmp()

When I ran it, I am not able to see any displacements based on the plot. I hoped I can see the displacement as a function of the number of iteration but it seems something went wrong. 
Am I making something wrong?

Thanks for your help!

Seungcheol

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.