small confusion on the div() term.
 

Hii Foamers,
I am modifying one of the turbulence models. In that process I have to introduce a new term in the source code, D/Dt(k_/omega_) which can be written as d/dt(k_/omega_)+ fvc::div(phi_,k_/omega_) for incompressible flows. Also considering steady state, we have:

D/Dt(k_/omega_) = fvc::div(phi_,k_/omega_)

Now my doubt is::in openfoam, the RHS of equation is solved based on FVM so, when these terms are discretized, is the cell-volume already included in it or do we need to specify it in the source code itself?

Or to put it in another words......
what I want is "D/Dt(k/omega)" and not "integral of D/Dt(k/omega) dV".
i.e. when I assign this new Field,
volScalarField DDtkw = fvc::div(phi_,k_/omega_)
what I get is DDtkw = integral of div(phi, k/omega) dV
Did I get it right here?

Nobody with insights in this, really??
I had even posted this problem in another location, yet no response in there too....

http://www.cfd-online.com/Forums/ope...nce-terms.html

Maybe the question is rather clearer in the above mentioned forum.
Anyone???? help!!