|April 13, 2011, 23:44||
Join Date: Apr 2011
Posts: 1Rep Power: 0
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.
|August 24, 2011, 23:02||
Join Date: Mar 2009
Posts: 22Rep Power: 8
according to the Gauss' divergence theorem,
the gradient at node p is calculated by integrate phi at the closed surface around the node.