gauss grad no correction for coupled boundary??!!
In the code for gauss gradient (gaussGrad.c), the grad calcuation at the boundary is done as:
const vectorField n
gGrad.boundaryField()[patchi] += n *
- (n & gGrad.boundaryField()[patchi])
So is there no correction for coupled patches?? I am getting embarrassingly wrong values of gradient at cyclic patch.
Second question: In the same code, calculation of grad has been done for three sets of mesh entities:
1. On face centers on patch
2. On cell centers adjacent to patch
3. On all inside cell centers
Is my understanding correct?
|All times are GMT -4. The time now is 20:22.|