Nonorthogonal correction in SimpleFOAM
I am using a nonorthogonal mesh in the calculation. I tried to put the nonorthogonal as 1 to add some correction. But when I read the code of simpleFoam. I can not understand what difference it can make. It seems the correction just solve the matrix one more time, with little change in the matrix. The change of the matrix is only from
pEqn.setReference(pRefCell, pRefValue) So I cannot understand what is going on with the nonorthogonal correction. Can any one explain it? Thanks a lot. Best, 
Nonorthogonal correction is not recommended for steadystate solvers as it is not needed to have a converged solution every timestep. Just put it to 0 and set limiters in your fvSchemes file in order to handle a nonorthogonal mesh.

Hi Cedric,
Thanks for your reply. Do you mean I only need to apply the correction in the laplacian terms. For example laplacian (phi, *) Gauss linear corrected Thanks a lot. Best, 
If your max orthogonality is < 60, you can indeed use "Gauss linear corrected" for your laplacian scheme, if higher you can use "Gauss linear limited 0.333". Your snGradScheme has to be edited accordingly (either corrected or limited 0.333). You can also limit the other schemes.
For example: Code:
gradSchemes 
Thanks a lot, Cedric.
In your example, you used gradSchemes { default cellLimited Gauss linear 1; } This is my first time to know there is a cellLimited scheme. I am curious about the difference between cellLimited and other limited schemes? If possible, could you also explain what is the function of the numbers following the limited scheme. Thanks again. Best, 
Quote:
e.g. a cell limited Gauss scheme, grad(p) cellLimited Gauss linear 1; for more details you can refer http://www.openfoam.org/docs/user/fvSchemes.php (4.4.3 Gradient schemes) 
This very nice document illustrates a lot about the choice of schemes, underrelaxation, pre conditioners for OpenFOAM.
Special treatment for nonorthogonal meshes. very well said "Who owns the mesh, owns the solution!!" http://www.dicat.unige.it/guerrero/o...sandtricks.pdf 
All times are GMT 4. The time now is 08:02. 