Caso model Change

Submitted by marceloschwed on Mon, 12/05/2016 - 14:10

Hello everyone,

In my work, I'm using the Easo model for capillary bridges. The model works good for very wet systems. However, I'm working with a system that the particles are held together by liquid bridges just at their contact point. This means that bridge breakage criteria works differently, since the bridges can't really slide due to particle rotation.
My idea was to calculate the initial contact points of each bridge, update them each time step by means of linear and angular velocity, and compare the new distance with the rupture distance, in order to update the contact flag.
One of the problems that I'm having is that I need arrays arrays to allocate the contact points data, and I can't use really use pointers from other classes (example atom) because the model is written in a header file (.h).
I'm experiencing a lot of memory management issues as I can't dynamically allocate and deallocate the memory on header files.

I would like to know if you can help me somehow with this problem.

Thank you very much in advance.

Marcelo Schwed

ckloss's picture

ckloss | Mon, 12/12/2016 - 21:10

Hi Marcelo,

first of all, I disagree to the statement that the Easo model is a model for _very_ wet systems, it is designed for liquid bridges between the contact points, it does not work for very wet (slurry) systems

If you want to implement a different model, you can store any contact values (also initial contact point) using the history value system which is used in the Easo model or by tangential model history already  -you don't need to allocate any arrays manually

Best wishes and happy coding,
Christoph