I have a question for gaussGrad.
00106 forAll(mesh.boundary()[patchi], facei)
00108 igGrad[pFaceCells[facei]] += pSf[facei]*pssf[facei];
00112 igGrad /= mesh.V();
In the 112 line, why igGrad is divided by mesh.V().
In the case of FVM, I understand that volume integration of grad(p) is sum of surface value of p. So, I can not understand the meaning of 112 line.
Please comment or advice for me.
Thanks in advance.
according to the Gauss' divergence theorem,
the gradient at node p is calculated by integrate phi at the closed surface around the node.
|All times are GMT -4. The time now is 21:32.|