[pino@JAEA]

I have a question for gaussGrad.

00106 forAll(mesh.boundary()[patchi], facei)
00107 {
00108 igGrad[pFaceCells[facei]] += pSf[facei]*pssf[facei];
00109 }
00110 }
00111
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.

[lzw2003]

Hi,

according to the Gauss' divergence theorem,
,
the gradient at node p is calculated by integrate phi at the closed surface around the node.