CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Mesh skewness correction

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

Like Tree1Likes
  • 1 Post By joegi.geo

Reply
 
LinkBack Thread Tools Display Modes
Old   November 21, 2012, 15:44
Default Mesh skewness correction
  #1
New Member
 
Jakub Pola
Join Date: Feb 2011
Posts: 17
Rep Power: 5
stainboy is on a distinguished road
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.
stainboy is offline   Reply With Quote

Old   October 23, 2013, 14:07
Default
  #2
Member
 
Malik
Join Date: Dec 2012
Location: Austin, USA
Posts: 52
Rep Power: 3
malaboss is on a distinguished road
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...
malaboss is offline   Reply With Quote

Old   October 23, 2013, 14:29
Default
  #3
New Member
 
Jakub Pola
Join Date: Feb 2011
Posts: 17
Rep Power: 5
stainboy is on a distinguished road
Hey,

Neither do I understand. It would be great to hear little bit about that. Does anyone have some article.

Jakub.
stainboy is offline   Reply With Quote

Old   October 24, 2013, 09:47
Default
  #4
Member
 
joegi
Join Date: Nov 2009
Location: genoa
Posts: 64
Rep Power: 6
joegi.geo is on a distinguished road
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.
joegi.geo is offline   Reply With Quote

Old   November 27, 2013, 10:09
Default
  #5
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 918
Rep Power: 14
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by joegi.geo View Post
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 07:26.
RodriguezFatz is offline   Reply With Quote

Reply

Tags
correction, discretisation, fvshemes, skewness

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] Unstructure Meshing Around Imported Plot3D Structured Mesh ICEM kawamatt2 ANSYS Meshing & Geometry 17 December 20, 2011 11:45
engrid -> save as .stl with boundarie codes Zymon enGrid 31 August 29, 2011 13:40
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 11:55
Severe nonorthogonality and severe skewness problem qtian OpenFOAM Running, Solving & CFD 2 January 22, 2008 18:47
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10


All times are GMT -4. The time now is 12:38.