Hello,
I am having various problems while implementing SPH cases in LIGGGHTS. Firstly, I am not able to include .stl geometries into the fix wall/region/sph command as it does not accept a fix mesh/surface. I then converted the .stl into an ASCII VTK format and tried to introduce it as a region, which was successful. However, implementing this region in fix wall/region/sph is causing problems when it is an open geometry, such as a cylinder with one face open.
Secondly, the lattice and create_atoms commands restrict freedom of particle generation. For atom_style sph, is there any way I can introduce a factory with fix insert/rate/region or any other option ?
If someone can elaborate on the relation between h, lattice distance, and distance r0 used in the fix wall/region/sph command, it will be useful. The following is the most popular error I got:
Particle <> with the Coordinates x= , y= , z= is on or inside fix wall/region/sph surface.