Simulation Box With Dryer and Wetter Regions of Cohesive Particles

Submitted by adguerra on Thu, 08/13/2020 - 18:00

Hello All,

I would like to study a situation where there are multiple groups of cohesive particles with different initial surface liquid contents in the same simulation box. To attempt this I have used
"pair_style gran model hertz tangential history cohesion easo/capillary/viscous"
and I tried to change the "surfaceLiquidContentInitial" to a "peratomtype" so that I could define multiple types of atoms with different surface liquid contents. This failed however, and I got the error
"Fix property/global (id m8): scalar style required for fix property/global variable surfaceLiquidContentInitial for usage with cohesion_model easo/capillary/viscous."
I tried with washino/capillary/viscous as well and got the same error.

I was able to run a simulation where some particles were totally dry and some particles were wet by using a hybrid pair style:
"pair_style hybrid gran model hertz tangential history cohesion easo/capillary/viscous &
bubble model hertz tangential off"
Such that there were some regions with cohesion and some with no cohesion. But I would like to implement different regions with different amounts of cohesion. Is this possible?

Thank you!!