Capillary forces between a network of particles

Asked by Vinícius Godim

The example capillaryLaplaceYoung has two particles interacting via capillary force. But how about simulating a network of spheres? I tried

pred = pack.inSphere((0,0,0),1)
sp=pack.randomDensePack(pred,spheresInCell=100,radius=.1,rRelFuzz=0)
O.bodies.append(sp)

O.engines=[ForceResetter()
           ,InsertionSortCollider([Bo1_Sphere_Aabb()])
           ,InteractionLoop(
               [Ig2_Sphere_Sphere_ScGeom()],
               [Ip2_FrictMat_FrictMat_CapillaryPhys()],
               [Law2_ScGeom_FrictPhys_CundallStrack(neverErase=1)]
                           )
           ,Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=1)
           ,NewtonIntegrator()
           ,GlobalStiffnessTimeStepper()
          ]

Nonetheless, any value however low I insert as an argument in capillaryPressure leads to a repulsive system while the console displays the following error:

In [1]: ERROR /build/yade-fDuCoe/yade-2018.02b/pkg/dem/Law2_ScGeom_CapillaryPhys_Capillarity.cpp:147 action: No meniscus found at a contact. capillaryPressure may be too large wrt. the loaded data files.

Is this expected? Do you believe there is a simple way to perform this sort of simulation?

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Vinícius Godim
Solved:
Last query:
Last reply:
Revision history for this message
Chareyre (bruno-chareyre-9) said :
#1

I don't understand your question. Could you clarify?

Le mer. 19 sept. 2018 18:08, Vinícius Godim <
<email address hidden>> a écrit :

> New question #674144 on Yade:
> https://answers.launchpad.net/yade/+question/674144
>
> The example capillaryLaplaceYoung has two particles interacting via
> capillary force. But how about simulating a network of spheres? I tried
>
> pred = pack.inSphere((0,0,0),1)
> sp=pack.randomDensePack(pred,spheresInCell=100,radius=.1,rRelFuzz=0)
> O.bodies.append(sp)
>
> O.engines=[ForceResetter()
> ,InsertionSortCollider([Bo1_Sphere_Aabb()])
> ,InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom()],
> [Ip2_FrictMat_FrictMat_CapillaryPhys()],
> [Law2_ScGeom_FrictPhys_CundallStrack(neverErase=1)]
> )
> ,Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=1)
> ,NewtonIntegrator()
> ,GlobalStiffnessTimeStepper()
> ]
>
> Nonetheless, any value however low I insert as an argument in
> capillaryPressure leads to a repulsive system while the console displays
> the following error:
>
> In [1]: ERROR
> /build/yade-fDuCoe/yade-2018.02b/pkg/dem/Law2_ScGeom_CapillaryPhys_Capillarity.cpp:147
> action: No meniscus found at a contact. capillaryPressure may be too large
> wrt. the loaded data files.
>
> Is this expected? Do you believe there is a simple way to perform this
> sort of simulation?
>
>
>
>
> --
> 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 address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>
>

Revision history for this message
Jérôme Duriez (jduriez) said :
#2

Hi,

Did you achieve getting a capillary force between 2 particles (i.e. did you correctly "install" the capillary files ?) ? In the spirit of examples/capillaryLaplaceYoung/capillaryBridge.py

Note that there also exists examples/capillaryLaplaceYoung/CapillaryPhys-example.py which deals with a packing ("network of spheres")

Revision history for this message
Vinícius Godim (viniciusgodim) said :
#3

I modified the capillaryBridge.py file in its directory and tried to run the simulation, from which I got that error.
I tried CapillaryPhys-example.py. It seems to be similar to what I am aiming at, though I do not understand why eliminating gravity (setting it to (0,0,0)) leads to a static system (without any acceleration from the capillary bridges).

Revision history for this message
Jérôme Duriez (jduriez) said :
#4

Then, we would need to see a MWE [*] demonstrating your problem. (You may also precise whether the initial capillaryBridge.py file was "working" or not on your system)

[*] https://www.yade-dem.org/wiki/Howtoask, point 3. in particular