# Mesh skewness correction

 Register Blogs Members List Search Today's Posts Mark Forums Read

 November 21, 2012, 16:44 Mesh skewness correction #1 New Member   Jakub Pola Join Date: Feb 2011 Posts: 17 Rep Power: 6 Dear Foamers, I have a question related to fvSchemes and interpolation schemes. I have read that mesh skewness affecs interpolation calculation. Is there a discretisation scheme which may help correcting this issue. I'm curious about skewLinear (Linear with skewness correction) from Numerical schemes. Does it somehow help with the calculations on skewed mesh? Does skewCorrected is the same scheme in OF 2.1.1? How to use it? Because wehn I want to put it in fvSchemes Code: ```interpolationSchemes { default skewCorrected; }``` OF throws an error "Discretisation scheme not specified". Thank you in advance for your help. Jakub.

 October 23, 2013, 14:07 #2 Member   Malik Join Date: Dec 2012 Location: Austin, USA Posts: 52 Rep Power: 4 Hi, I have almost the same issue : how to use the skewness correction in OpenFOAM. To me the skewness is critical especially when we do an interpolation. Hence it has to be mentionned in the interpolation schemes. Plus, it may be used in the div schemes of non surfaceField like U : when we integrate over a cell, we will need to interpolate the field from the center to the surface. I guess for a surfaceField, it is useless. If then we use : Code: ```interpolationSchemes { default linear skewCorrected; } divSchemes { div(U) Gauss linear skewCorrected; }``` it will run. But I compared the result with the result simulating Taylor Vertices on an unstructured grid using Code: ```interpolationSchemes { default linear; } divSchemes { div(U) Gauss linear; }``` And the resulting U and p are exactly the same. Then I think I don't really understand how this skewness correction work...

 October 23, 2013, 14:29 #3 New Member   Jakub Pola Join Date: Feb 2011 Posts: 17 Rep Power: 6 Hey, Neither do I understand. It would be great to hear little bit about that. Does anyone have some article. Jakub.

 October 24, 2013, 09:47 #4 Member   joegi Join Date: Nov 2009 Location: genoa Posts: 67 Rep Power: 7 Hi, Let me put my two cents. Mesh skeness will reduce the order of face integration. Basically it is the deviation from the face center to the actual face integration point (the intersection between the vector joining two neighboring CV and the common face), and this is already taken care of in OF. You can refer to Hrvoje Jasak phd thesis. If you have highly skewed cells you can use least squares for gradients calculation in order to reduce the numerical diffusion. Maybe increasing the number of corrector steps will help as well. jg RodriguezFatz likes this.

November 27, 2013, 11:09
#5
Senior Member

Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
Quote:
 Originally Posted by joegi.geo Mesh skeness will reduce the order of face integration. Basically it is the deviation from the face center to the actual face integration point (the intersection between the vector joining two neighboring CV and the common face), and this is already taken care of in OF. You can refer to Hrvoje Jasak phd thesis.
Hi, could you tell me, where in the code and where in the thesis I can see how the skewness correction is applied in OF?

Edit:
Ok, so after some more reading and searching I found the following:
Mesh weights (also used for central differencing in OF) for the cell-to-face interpolation are calculated using this formula:
http://openfoamwiki.net/index.php/Op...on_%28class%29
In the source code one can find the calculation in the "surfaceInterpolation.C" file, function is called "makeWeights()".
Now, obviously this is not just the standard linear weighting, as it is not only the distance of cell centers to the cell face that counts, but also the projection of the vectors on the face normal.
But also in the thesis I can not find any derivation of that fomula, I just found some words about the diffusion correction. Unfortunately, I could not adapt the lines to the convective / face value calculation.
__________________
The skeleton ran out of shampoo in the shower.

Last edited by RodriguezFatz; November 28, 2013 at 08:26.

 Tags correction, discretisation, fvshemes, skewness

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post kawamatt2 ANSYS Meshing & Geometry 17 December 20, 2011 12:45 Zymon enGrid 31 August 29, 2011 13:40 SSL FLUENT 2 January 26, 2008 12:55 qtian OpenFOAM Running, Solving & CFD 2 January 22, 2008 19:47 Joe CFX 2 March 26, 2007 18:10

All times are GMT -4. The time now is 10:15.