Hi, experts,
Recently I try to couple cfdem inserting superquadric models. But the error always shows in my case like "ERROR on proc 0: Using default surface model for non-spherical particles! (../surface_model_default.h:144)" .
After I seeking for the code, it told that
"if(sidata.is_non_spherical)
error->one(FLERR,"Using default surface model for non-spherical particles!");"
I am a little confused about the meaning of the "sidata". Could anyone give some explanation and advice .
Thx so much,
Clyde
mschramm | Wed, 04/13/2022 - 18:43
sidata and surface models
Hello,
sidata is a data structure that contains data of a single interaction between two objects (atom to atom or atom to wall).
Basically, as LIGGGHTS loops through atom interactions, a new sidata is used to describe this interaction. Sidata is then sent to your pair models to apply forces.
surface models describe how an atom interacts with a surface. This is specified in your pair_style and wall/gran commands with the keyworld "surface".
When using superquadratics you must set the model to superquadric for both the pair_style and wall/gran commands...
pair_style gran model hertz tangential history rolling_friction epsd2 surface superquadric
and
fix granwalls all wall/gran model hertz tangential history rolling_friction epsd2 surface superquadric mesh n_meshes 1 meshes cad