nonorthogonal corrections, when and where to apply?
 

Hello!
I found that it is necessary to apply nonorthogonal corrections when working with skew meshes. What exactly are the necessary steps?

Do I merely add a loop
around each "solve(...)"?
The Programmer's guide mentions that snGradCorrection is also needed on nonorthogonal meshes. But snGradCorrection isn't called in any of the solvers.

What if I solve several equations in a segregated solver?
Do I run the loop on each individual equation:
Or do I run the loop over the whole sequence:
In the example: http://openfoamwiki.net/index.php/Ho...ation_to_solve
TEqn.solve() is added *outside* the non-orthogonal corrections loop. Why? Is TEqn not affected by non-orthogonality?

If you know a good explanation for understanding the details of non-orthogonal corrections, please tell.

Thank you
Andrej

I am not an expert but I think the non-orthogonal corrections are applied as implicit and explicit term in the discretization of the equation to be solved. Look into Jasaak's thesis to see how it is done. For example the poison pressure equation solution has a laplacian term which can benefit from non-orthogonal correction. So whenever the loop is called the explicit correction term is modified thereby refining the solution. It is all automatic so you don't have to add a snGradCorrection term yourself, just specify that correction should be applied in the discretized equation in fvSchemes.
Hope it helps.

Daniel, thank you. That answers my first question.

Now what if I have to solve a sequence of equations?

And why is TEqn not corrected in the example from the Wiki?

Thank you
Andrej

I think solution of the pressure equation with non-orthogonal correction should finish first before solving the transport equation. Once you know the corrected velocity field then you can solve transport equation of T. I suppose you can apply the correction to the diffusion term of T as well ( any implicit laplacian or snGrad term on a grid with skewed cells will benefit from multiple corrections ). So in case of T , i guess the correction is only applied once. It maybe that the pressure equation benefits more from multiple corrections depending on the skeweness of the grid.
Hope it helps
Wiki link: http://www.cfd-online.com/Wiki/Diffusion_term

Hi,
This is an old thread, but if someone comes across it due to their interest in non-orthogonality, have a look at the following videos which describe it quite well:
https://www.youtube.com/watch?v=Dg2PwX6yxYY
https://www.youtube.com/watch?v=yU7r8mYK3bs