Introducing viscous forces in a liquid bridge

Hi,

1.
If you have a closed form expression, then it is quite easy to add the force to each body from Python (then from C++ if the extra cost of this operation -- always greater in Python than in C++ -- is a problem).

E.g.
O.forces.addF(i.id1,- yourForceVectorHere) # i is an interaction
O.forces.addF(i.id2,yourForceVectorHere)

2.
Interaction forces proportional to relative velocity already exist in YADE, and even with extra fluid-induced forces: see [*] and related classes. See examples at [**]

3.
Probably many other ideas: YADE being open source, you can implement everything you want ! (typically with new Engine classes)

Maybe you can try to give more details if you want more answers.

[*] https://yade-dem.org/doc/yade.wrapper.html?highlight=viscel#yade.wrapper.ViscElCapMat
[**] https://gitlab.com/yade-dev/trunk/tree/master/examples/capillary

Hello Jerome,

Thanks for the feedback.
ViscElCapMat would indeed be the closest to what I am looking for but I can not find the theoretical details of the law, in particular the implementation of the viscous contributions.
In my case and according to theory, the normal viscous force is proportional to the relative translational velocity between two particles with a damping coefficient which is a function of the separation distance between the two particles.
The tangential viscous force is proportional to both the relative translational and rotational velocities with two damping coefficients which are again functions of the separation distance between the particles (Washino et al, Powder Technology 302 (2016) 100-107).
So this sounds more tricky....
??
My best wishes

Fr.

As for more details about ViscElCapMat, you might found some in https://link.springer.com/article/10.1007%2Fs10035-014-0527-z (if I'm not mistaken, that's the paper related to this code). Source code is at [*]

Maybe you do need a new consitutive law. In this case, C++ coding new Material, Ip2 and Law2 functors would be the (only) way to go. You may take inspiration from existing contact models ([*] for instance). The procedure also had been discussed somewhere in the mailing list.

Anyhow, significant understanding of YADE internals and C++ would help / be gained in the process..

[* https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/ViscoelasticCapillarPM.hpp and .cpp

Hello Jerome,

Yes, ViscElCapMat is more precisely a viscoelastic model of contact (very usual) added to a capillary force so this is not what I wanted to model: adding the viscous contribution of the liquid bridge between two particles. So I probably would need a new constitutive law as you suggest...

Thanks,

Fr.

Hi,
There is probably no need to get both capillary force and viscous force from a single model.
The viscous (lubrication) forces are given by Law2_ScGeom_ImplicitLubricationPhys.
Bruno

Hi Bruno,

That's a very interesting option; the previous theory of viscous contribution of bridges (evoked above) was derived
indeed from Reynolds equation of lubrication.
Where can I get all the theoretical details of Law2_ScGeom_ImplicitLubricationPhys and references ??
I don't see viscosity as an input parameter ??
Aren't there any technical difficulties with YADE dealing with two different interaction laws for the same
 couple of particles ??

All the best,

Fr.

>Where can I get all the theoretical details of Law2_ScGeom_ImplicitLubricationPhys and references ??

There is paper in preparation, but overall lubrication is lubrication, there are plenty papers about it.
See for instance https://link.springer.com/article/10.1007/s10035-015-0560-6

>I don't see viscosity as an input parameter ??

You should.

> two different interaction laws

Capillary forces do not come from an interaction law actually, hence no issue in principle.

Bruno

Thanks Anton Gladky, that solved my question.

See maybe Ip2_FrictMat_FrictMat_LubricationPhys.eta for the fluid viscosity in Lubrication classes