lfgmarc 
July 14, 2012 20:25 
A question about UEqn sonicFoam
Hi, I have a question about the UEqn implemented into sonicFoam, the code of the UEqn is:
Code:
fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence>divDevRhoReff(U)
);
My question is about the term that is computed from the turbulence model, i.e if Kepsilon model is selected, this term results:
Code:
tmp<fvVectorMatrix> kEpsilon::divDevRhoReff(volVectorField& U) const
{
return
(
 fvm::laplacian(muEff(), U)
 fvc::div(muEff()*dev2(fvc::grad(U)().T()))
);
}
Initially I thought that it corresponded to
but after examining the definition of dev2 in src/OpenFOAM/primitives/Tensor/TensorI.H
I found the following:
is the rateofstrain tensor given by
but according to the definition of DivDevRhoReff
And manipulating the expression (1) I conclude that the term
Code:
+ turbulence>divDevRhoReff(U)
does not correspond to (1).
If someone can help me with my confusion about this term, I was very thankful
