here a test case for cohesive capillary viscous model .
In the script you find "easo" model, just change "easo" with "washino" to switch to Washino's model.
I set a liquid source on one atom and run for 1 timestep, so I have an initial superficial liquid on one particle. This allows to check the liquid migration after the impact.
The collision happens at timestep 1000000.
Check files dump985000.collision and dump1200000.collison, to see velocities before and after the collision. Moreover, look at the liquid on particle two, after the collision you will find some liquid on this particle, before the collision the superficial liquid was 0.
In another post I suggested I had some problem because the final velocity was higher than the initial velocity. This problem is not here anymore... probably I made a mistake in setting the model parameters.
By the way, the final velocity for easo cohesive model seems too low.
And... if someone have some suggestion on how to dump local information (forces), see line 77-78 of the script, please share :)
I remember Christoph Kloss suggested something about a history value... may be an hint for someone (not for me so far)
|particle information - copy in the same directory of the script||852 bytes|
|run liggghts < in.capillary.txt||2.49 KB|