Salome Viscous Layers
 

Anyone with experience. I am doing a simple example, 3D cylinder with a viscous layer (Total thickness: 1e-5 m, Number of layers: 10-20, Stretch factor: 1.2 (diameter of cylinder 0.03 m)). When I compute all is ok, but when I check the mesh with OpenFOAM, the mesh has plenty of mistakes and is not proper for simulation...even if I refine with triangles the surface.
Any Suggestion

Could you share the output of checkMesh?

PS: Guide: How to ask a question on the forums

I think is information useless but Sure! I am sharing the best mesh I could generate. However is not the the one I want, is not really refined near the wall (first layer 1e-6 m, then 5 layers, Stretch factor: 3.) I would like to add more layers but errors start to apear, like:

***Max cell openness
***Max aspect ratio
***Max skewness, etc etc.

In the present case It is just high number of severely non-orthogonal faces.

Thanks for ur time

In contrast, providing checkMesh output is critical. Plus, you should explain your case thoroughly without the reader asking for it.

Anyways.

As can be seen from the output:

Your grid contains tetrahedra and prisms (and a bit pyramids). With this tetrahedral topology, it is natural that checkMesh measures these quantities. It does not necessarily mean that your mesh is bad, but checkMesh's meshQuality metrics (and schemes) are mostly tuned for hexahedra.

I recommend you to use hexahedral mesh capabilities of OpenFOAM, or I speculate this, you can use OpenFOAM-extend for simulations using tetrahedral mesh.

Kind regards

My intention was not the reader asking, just someone that have tried to refine (with salome) near the wall. Sorry if I did not stated in a proper way the problem.

Do u mean that if I get "Failed n mesh checks" i can still run my simulation with good results?

If i would have to mesh a cylinder, definitively I would do it with blockMesh, but my geometry is more complicated (I already tried with blockMesh and is not easy to get good results when trying join different blocks (boundary layer in the order of 1e-6 m).
What I did is starting to work in a cylinder in salome (if i can not do it in a cylinder, definitively I would never do it in my geometry).

Is OpenFOAM-extend better when using tetrahedral meshes? I have some things (BC, post-proc and solvers) programed by myself and i got problems of compatibility every few months with of-dev (or from of-3 to of-5). I do not want to imagine if I jump to OpenFOAM-extend. Any recommendation or forum?

Thanks a lot

Yes, in principle, you can. Tetrahedra is commonly being used in industry. Though I speculate that OF schemes are not appropriate to compute on tetrahedra. Might someone else give some more-educated opinion, don't know.

IMHO, obtain STL surface mesh of your complicated geometry from SALOME (it provides good quality STLs as opposed to CAD software. e.g. SolidWorks etc.), then use snappyHexMesh to mesh it.

It was a suggestion only. OF-extend, AFAIK, has schemes more suitable to tetrahedra-meshed sims (for FSI). You can check forum for further information.

Hope these help.

About snappyHexMesh, I am afraid of having problems as with blockMesh (refineWallLayer) in a complex geometry. Reading the following link (https://cfd.direct/openfoam/user-gui...26-1980005.4.7)

4 if the validation criteria can be satisfied, insert mesh layers;
5 the mesh is checked again; if the checks fail, layers are removed

And seen what people experienced with boundary layers and snappyHexMesh (https://www.cfd-online.com/Forums/op...pyhexmesh.html)

I think i am right...

Do u have any experience using snappyHexMesh for boundary layers in the order of 1e-6 m (y+ < 1)?

Finally, I do not see the restriction u said about tetrahedral mesh

"By default OpenFOAM defines a mesh of arbitrary polyhedral cells in 3-D, bounded by arbitrary polygonal faces, i.e. the cells can have an unlimited number of faces where, for each face, there is no limit on the number of edges nor any restriction on its alignment. A mesh with this general structure is known in OpenFOAM as a polyMesh." (https://cfd.direct/openfoam/user-gui...h-description/)
Am i missing something?

No, I'm afraid.

OpenFOAM mesh topology capabilities allow polyhedra. My suspicion is however whether schemes you will use on tetrahedra were derived and implemented with the consideration of tetrahedra? For example, you may very likely have issues if you use Gauss linear for gradient terms within such a domain (but you can use leastSquares to ensure second-order accuracy independent of the mesh topology).

As I said, I only speculated about this, and encourage you to be more cautious when you will use tetrahedra.

Hope you resolve this matter.

PS: It would be very nice of you if you would share your experience herein in the future.

Generally, hex-meshes are better than tet-meshes, but OpenFOAM and its schemes can deal with both.

It is a very simple utility and does not always provide good layers, but it is worth looking at. That is, the refineWallLayer utility of OpenFOAM. I have used it sometimes to get through the viscous sublayer, i.e. to get y+ ~ 1.

Say you want to add five layers to a wall called pipeWall, it would look something like:

In a bash script you could do it more elegantly with something like