nonorthogonality vs skewness
hi,
i was wondering whether someone could provide some insight regarding the difference between nonorthogonality and skewness in terms of corrections available in the OF fvSchemes file. From what i had read, nonorthogonality 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 nonorthogonality when it seems an nonorthogonal 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 nonorthogonality 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 OF2.1.1 to do steadystate 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 
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 pv coupling will address it as you go along. nOrthoCorrectors is mainly useful where you need extra loops to correct for nonorthogonality in a single time step. 2) the nonortho corrections used in OF is the overrelaxed method as per HRV's thesis 3) As for the nonorthogonality vs skewness: i) nonorthogonality is mainly related to the nonparallel 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 nonorthogonality ... ii) how much nonorthogonal correction to apply i.e. is this decision purely based on the tradeoff 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 
Thank you for sharing

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:
Quote:
cheer and best regards jonathan 
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 (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (3.43443e20 4.80201e21 5.39969e15) OK. Max cell openness = 4.90788e16 OK. Max aspect ratio = 739.181 OK. Minimum face area = 2.55578e06. Maximum face area = 10.688. Face area magnitudes OK. Min volume = 1.6488e08. Max volume = 0.0751724. Total volume = 274.745. Cell volumes OK. Mesh nonorthogonality Max: 79.9874 average: 11.9844 *Number of severely nonorthogonal faces: 715. Nonorthogonality check OK. <<Writing 715 nonorthogonal 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 ? 
All times are GMT 4. The time now is 13:51. 