Problem about fvMatrix in OpenFOAM, over volume or not
 

Dear All foamers:

I have a question about the equations in OpenFOAM. When we set those equations, like UEqn and PEqn, is every term over by cell volume or not? Like in UEqn, the fvm::div term is value integral the cell volume, but the grad(p) term is divided by cell volume. That really confused me,

Can anybody give me a clear explanation about this issue? Or where can I find the answer about this??

thank you all in advance.

Problem about fvMatrix in OpenFOAM, over volume or not
 

for the source term in UEqn should follow similar expression of other ddt or div terms , for example fvm::ddt(rho, U) terms already includes mesh.V() as (see in src finiteVolume code for schemes)

fvm.source() = rDeltaT*mesh().V()*
(
coefft0*vf.oldTime().internalField()
- coefft00*vf.oldTime().oldTime().internalField()
);

for updating any source terms following fvm.source() function just update the function that does not need to divide by V. See the example in OpenFoam >src lagrangian> intermediate> KinematicCloudI.H for SU(U) source term.