[yhaomin2007]

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.

[Nazmul Khan]

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.