[alundilong]

Here is the function return effective stress tensor term in momentum equation from GenEddyVis.C (OpenFOAM-2.1.1).

tmp<fvVectorMatrix> GenEddyVisc::divDevRhoBeff(volVectorField& U) const
{
return
(
- fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(T(fvc::grad(U))))
);
}

my question is
where is the term for the trace of the subgrid scale stresses?
It should looks like this or others..., anyway, k them is missing, and I didn't see anywhere account for it.

tmp<fvVectorMatrix> GenEddyVisc::divDevRhoBeff(volVectorField& U) const
{
return
(
- fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(T(fvc::grad(U))))
+ 2.0/3.0*fvc::grad(rho()*k())
);
}

My derivation is attached, please give me some suggestions?