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

nonorthogonality vs skewness

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

Like Tree31Likes
  • 1 Post By Jonathan
  • 16 Post By Jonathan
  • 13 Post By Jonathan
  • 1 Post By miladrakhsha

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 9, 2013, 09:46
Default nonorthogonality vs skewness
  #1
Senior Member
 
Join Date: Mar 2010
Posts: 175
Rep Power: 17
Jonathan is on a distinguished road
hi,

i was wondering whether someone could provide some insight regarding the difference between non-orthogonality and skewness in terms of corrections available in the OF fvSchemes file.

From what i had read, non-orthogonality affects only the diffusion terms and can be fixed using 'corrected' schemes in fvSchemes. What i dont understand is how the convection terms 'seem' to remain unaffected by nonorthogonality, or why there is a differentiation made between skewness and non-orthogonality when it seems an non-orthogonal mesh will probably be 'skew' as well.

Can anyone point me in the right direction for some answers or provide a brief explanation as to how skewness and non-orthogonality differ and why they affect only the diffusion and convection terms respectively, and whether these should be dealt with in the fvSchemes file or by running additional nonOrtho corrector loops in the fvSolution SIMPLE control dictionary.

I am using a modified version of SRFSimpleFoam in OF-2.1.1 to do steady-state simulations of a turbine rotor. checkMesh gives me some warnings about nonOrtho faces, but mesh skewness is OK. I am using linear interpolations for the grad and laplacian terms and upwind discretisation for the convective terms.

many thanks in advance
kind regards
Jonathan
Ramzy1990 likes this.
Jonathan is offline   Reply With Quote

Old   June 12, 2013, 08:41
Default
  #2
Senior Member
 
Join Date: Mar 2010
Posts: 175
Rep Power: 17
Jonathan is on a distinguished road
so i have some answers for myself - i'll just share them in case someone else finds them usefull ...

1) Unless you are running an unsteady simulation, you dont really need to include additional nOrthoCorrector loops in your simulation setup, as the iterative nature of the p-v coupling will address it as you go along. nOrthoCorrectors is mainly useful where you need extra loops to correct for non-orthogonality in a single time step.

2) the nonortho corrections used in OF is the over-relaxed method as per HRV's thesis

3) As for the nonorthogonality vs skewness:

i) non-orthogonality is mainly related to the non-parallel of the face normal vector, while

ii)skewness is more related to the fact that distorted meshes result in the intersection of the vector between the adjacent nodes not being at the face centroid

iii) so a nonorthogonal mesh may not necessarily be skew

What i am still trying to figure out is:

i) why convective terms do not need to be corrected for as the diffusion terms are? Seems strange since both terms use a face normal vector to calculate the face fluxes and therefore should be affected by non-orthogonality ...

ii) how much non-orthogonal correction to apply i.e. is this decision purely based on the trade-off between stability (nonortho correction reducing stability as reduces the diagonal dominance of the system) and accuracy?

Many thanks for any help in advance,
regards
jon
Jonathan is offline   Reply With Quote

Old   April 24, 2014, 20:53
Default
  #3
New Member
 
Francis
Join Date: Mar 2014
Posts: 1
Rep Power: 0
Mazart is on a distinguished road
Thank you for sharing
Mazart is offline   Reply With Quote

Old   April 25, 2014, 04:15
Default
  #4
Senior Member
 
Join Date: Mar 2010
Posts: 175
Rep Power: 17
Jonathan is on a distinguished road
Hi - no problem!

Since it seems a couple people are reading this post, i will update some of the things i did not yet understand since last time:

Quote:
i) why convective terms do not need to be corrected for as the diffusion terms are? Seems strange since both terms use a face normal vector to calculate the face fluxes and therefore should be affected by non-orthogonality ...
This is because the correction applies to the gradient term only, and therefore, since diffusion only relies on the this term, the convective terms are not affected by it.

Quote:
ii) how much non-orthogonal correction to apply i.e. is this decision purely based on the trade-off between stability (nonortho correction reducing stability as reduces the diagonal dominance of the system) and accuracy?
Yes, this only depends on the trade-off between stability and accuracy. Since correction reduces diagonal dominance of the coefficient matrices, it comes down to that only.

cheer and best regards
jonathan
bmikuz, arvindpj, Metti and 10 others like this.
Jonathan is offline   Reply With Quote

Old   September 1, 2014, 00:39
Default
  #5
New Member
 
miladrakhsha
Join Date: Aug 2012
Posts: 29
Rep Power: 14
miladrakhsha is on a distinguished road
I was wandering if somebody can give some hint on this one.

Checking geometry...
Overall domain bounding box (-50 -50.0901 -0.035) (50.29 50 0)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.43443e-20 4.80201e-21 5.39969e-15) OK.
Max cell openness = 4.90788e-16 OK.
Max aspect ratio = 739.181 OK.
Minimum face area = 2.55578e-06. Maximum face area = 10.688. Face area magnitudes OK.
Min volume = 1.6488e-08. Max volume = 0.0751724. Total volume = 274.745. Cell volumes OK.
Mesh non-orthogonality Max: 79.9874 average: 11.9844
*Number of severely non-orthogonal faces: 715.
Non-orthogonality check OK.
<<Writing 715 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.47068 OK.
Coupled point location match (average 0) OK.

Mesh OK.


Does this high skewness affect the accuracy of my solution? I know that maximum admissible skewness is defined as 4. I just wanted to know can the skewness be the source of inaccuracy or not. And also how one can deal with the skewed cells in OpenFOAM and is there any remedy for this problem ?
elham usefi likes this.
miladrakhsha is offline   Reply With Quote

Reply

Tags
correction, nonrothogonal, skewness

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Netgen] Import netgen mesh to OpenFOAM hsieh OpenFOAM Meshing & Mesh Conversion 32 September 13, 2011 05:50
TGrid High Skewness khimkhim ANSYS Meshing & Geometry 1 February 17, 2010 23:36
problem with skewness jjchristophe FLUENT 1 January 14, 2010 21:22
Severe nonorthogonality and severe skewness problem qtian OpenFOAM Running, Solving & CFD 2 January 22, 2008 18:47
Skewness Taskin FLUENT 5 April 6, 2005 15:21


All times are GMT -4. The time now is 02:14.