CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Deviatoric Part of the reynolds tress tensor (

make September 22, 2008 10:55

In OpenFOAM-1.4.1 the reynolds
In OpenFOAM-1.4.1 the reynolds stress tensor in the kEpsilon model has been calculated using
(2.0/3.0)*I)*k_ -(nut_)*2*dev(symm(fvc::grad(U_)))
which seems plausible to me.

Why is it in the RNGkEpsilon
((2.0/3.0)*I)*k_ - nut_*2*symm(fvc::grad(U_))
Where is here the deviatoric part? Isn't that wrong? Shouldn't it be
((2.0/3.0)*I)*k_ - nut_*2*dev(symm(fvc::grad(U_)))

Also for the divR term it is calculated using
- fvm::laplacian(nuEff(), U)
- fvc::div(nuEff()*dev2(fvc::grad(U)().T())
Shouldn't it be dev2?


pbo September 22, 2008 12:11

Matthias, it seems you're l

it seems you're looking at the incompressible models here, for which

dev(symm(fvc::grad(U_))) = symm(fvc::grad(U_))


tr(symm(fvc::grad(U_))) = div(U) = 0


make September 22, 2008 14:02

yep I was too small minded and
yep I was too small minded and didn't catch that point.
Nevertheless, no I can put my mind at rest and perhaps others too.
Patrick, thank's a lot for your quick reply.


All times are GMT -4. The time now is 07:31.