|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
CK
Join Date: Apr 2011
Posts: 1
Rep Power: 0 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
zhiwei liu
Join Date: Mar 2009
Posts: 22
Rep Power: 17 ![]() |
Hi,
according to the Gauss' divergence theorem, the gradient at node p is calculated by integrate phi at the closed surface around the node. |
|
![]() |
![]() |
![]() |
Tags |
gaussgrad |
|
|