k-epsilon model - k term in the U-equation
 

Hi,

I am having a close look at the k-epsilon model implemented for incompressible flows in Openfoam.

in the following piece of code:
tmp<fvVectorMatrix> kEpsilon::divDevReff(volVectorField& U) const
{
return
(
- fvm::laplacian(nuEff(), U)
- fvc::div(nuEff()*dev(fvc::grad(U)().T()))
);
}

The k term in the Boussinesq eqauation ( 2/3 * k * delta_ij ) does not appear.
Where can I find this term in the implementation?

Thanks for your help!
Gregory

You can find it in the file kEpsilon.C. I made a copy of it here:


I see that the minus sign is reversed ~? maybe it is multiplied with minus somewhere else !!!!!