CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [snappyHexMesh] Too High Skewness (https://www.cfd-online.com/Forums/openfoam-meshing/226716-too-high-skewness.html)

Mateus Elias May 6, 2020 02:33

Too High Skewness
 
I create a snappyHexMesh and, when I run the checkMesh to confirm it, in the Time 2 an 3 occour an error in the Max. skewness = 5.56.

I would like to know if the accuracy of the simulation would be affected by this high skewness? What parameters I could change to catch a better result?

Kind regards.

HPE May 6, 2020 07:31

google is your friend - discussed in this forum and elsewhere many times - please do consider to search "skewness" and "OpenFOAM".

Mateus Elias May 14, 2020 21:01

Thank you for your wish to help.
But your answer was to general. I was expecting something more specific.
But thank you anyway.

HPE May 15, 2020 16:42

Hi,

I know.. But please do consider to isolate and focus the questions above which were very broad. I believe, most people tend to shy away from answering a question which could be answered with pages of information.

A quick search revealed the following, for instance (link):

"Mesh skewness reduces the diagonal dominance of the discrete Poisson operator which slows down convergence. ADI type preconditioners also become less effective.

Convergence rate aside, the accuracy might also drop when mesh is highly skewed. This results mainly from the way in which the face centered (in FV schemes) pressure gradients are computed using cell-centered pressure values. Usually a second order central-difference approximation is used and the accuracy might drop to first order for very high skewness (or even stretching). It is perhaps more appropriate to compute pressure gradients using multi-dimensional stencils (which could be non-symmetric unlike the central difference operators). I am not sure about the effects of directional biasing in this way of computing pressure gradients (since there is no physical basis for biased stencils here unlike upwinding for advective operators). However, the matrix corresponding the discrete Laplacian operator would still be symmetric.

Some the grid layout issues have been discussed in the introduction of the paper by Zang et al. (JCP, 1994 or 1995). You can find a concise discussion on the problems/advantages of the grid system (variable layout) used for incompressible flows. You might find some useful discussions in the papers they refer to as well."

Hope this helps for a good start.


All times are GMT -4. The time now is 11:33.