Contact Models for superquadric particles

Submitted by Sherif on Sat, 05/15/2021 - 23:00

I want to ask about using different contact models for superquadric particles. I can see only the following contact models are available in the "style_contact_model.whitelist" file for superquadric particles:

GRAN_MODEL(HOOKE, TANGENTIAL_HISTORY, COHESION_OFF, ROLLING_OFF, SURFACE_SUPERQUADRIC)
GRAN_MODEL(HOOKE, TANGENTIAL_NO_HISTORY, COHESION_OFF, ROLLING_OFF, SURFACE_SUPERQUADRIC)
GRAN_MODEL(HOOKE, TANGENTIAL_NO_HISTORY, COHESION_OFF, ROLLING_CDT, SURFACE_SUPERQUADRIC)
GRAN_MODEL(HOOKE, TANGENTIAL_NO_HISTORY, COHESION_SJKR, ROLLING_OFF, SURFACE_SUPERQUADRIC)

I have added different combinations of models in this file and recompiled liggghts so that the new models are considered optimized when I run the model (the error that the contact model is unoptimized does not pop up anymore).

Are these new added models really applied to my superquadric particles interactions or the resulted interactions are fake? I mean in other words, does making any combination of contact models valid? If yes why the developers didn't add all possible combinations to the whitelist?

mschramm | Tue, 05/18/2021 - 03:38

Hello,
Rest assured that when you see the "unoptimized" warning message that the contact model combination you selected is being used.
It simply has not been "optimized". I am not 100% on why they did not make it so all of the combinations were optimized but it probably either deals with keeping the executable small or keeping the memory low.

Sherif | Tue, 05/18/2021 - 13:06

Thanks Mathew, I really appreciate your responds to different questions.