[pvpnrao]

I have the .stl file of a closed geometry (a simple cylinder) which I superimposed on a Cartesian hex mesh (a bigger cube) and ran snappyHexMesh. Everything works fine.

However, snappyhexmesh deletes the mesh from the region lying inside .stl file geometry.

Is there some option in snappyhexmeshdict that I can switch-off to avoid it from deleting that mesh?

I need to retain mesh on both inside and outside of the .stl surface

[Antimony]

Hi,

You would need to specify and set cellZoneInside to inside, among others.

See page 31 for more details: https://openfoamwiki.net/images/f/f0...SlidesOFW7.pdf

Cheers,
Antimony

[KaLium]

Do you want different cellzone to cylinder? Are you planning to solve it as solid?

For quite similar problem, my code was something like this:

Surface1
{
level (2 2);
faceType baffle;
cellZone fluid1;
cellZoneInside solid1;
}

locationsInMesh
((-5.0, 0.0, 0.0) fluid1)
((-1.0, 0.0, 0.0) solid1)

[pvpnrao]

Antimony ... KaLium

Thank you very much for your responses.

KaLium: Yes, I want to model the internal mesh as solid. I will try your suggestion.
If is it not too much to ask. I am also struggling to reference the solid cellzone in the solver. How do you realize a solid within the fluid region in the solver, I believe this is done my adding a forcing term in on the right hand side using fvOptions which will only work for on solid cellzone...but I have not been successful in getting it to work ...
can you please give me some insights on how to get it working
Do you have any sample code ?

how is the solid-liquid interface treated in openfoam

Thanks in advance