[fisch]

Hi Foamers,

does anybody know if div(phi,U) in an incompressible solver (like simpleFoam) is calculating the "full" divergence term of phi*U or is it directely introducing the incompressibility conditions that the calculation would be reduced to the half.

For a better understanding what i mean i will write it analytical for 2D motion (i will leave out the density):

Option a) div(U,U) = u*( du/dx+dv/dy) + u*du/dx + v*du/dy
Option b) div(U,U) = u*du/dx + v*du/dy (because the term in the brackets satisfies div(U)=0)

Which option is openFoam using?

Thanks a lot,
rupert