|
[Sponsors] |
nonorthogonal corrections, when and where to apply? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 7, 2012, 06:00 |
nonorthogonal corrections, when and where to apply?
|
#1 |
New Member
Andrej
Join Date: Feb 2012
Posts: 2
Rep Power: 0 |
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 Code:
for (int nonOrth=0; nonOrth<nNonOrthCorr; nonOrth++){ 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: Code:
for (int nonOrth=0; nonOrth<nNonOrthCorr; nonOrth++){ Eqn1.solve(); } for (int nonOrth=0; nonOrth<nNonOrthCorr; nonOrth++){ Eqn2.solve(); } Code:
for (int nonOrth=0; nonOrth<nNonOrthCorr; nonOrth++){ Eqn1.solve(); Eqn2.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 |
|
February 7, 2012, 09:13 |
|
#2 |
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 16 |
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. |
|
February 7, 2012, 10:45 |
|
#3 |
New Member
Andrej
Join Date: Feb 2012
Posts: 2
Rep Power: 0 |
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 |
|
February 7, 2012, 16:36 |
|
#4 | |
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 16 |
Quote:
Hope it helps Wiki link: http://www.cfd-online.com/Wiki/Diffusion_term |
||
May 5, 2021, 13:29 |
|
#5 |
New Member
Gerhard
Join Date: Mar 2017
Posts: 26
Rep Power: 9 |
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 |
|
|
|