[sharonyue]

Hi guys,

I found this in Ueqn:

volTensorField Rc1
        (
            "Rc1",
            (((2.0/3.0)*I)*nuEff1)*tr(gradU1T) - nuEff1*gradU1T
        );

but in openfoam 2.0.x its like:

volTensorField Rca
        (
            "Rca",
            ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
        );

mathematically, there should be this sqr(Ct)*k in the code, see the last term:


I also found this thread:
http://www.cfd-online.com/Forums/ope...eulerfoam.html

Any ideas?

[sharonyue]

also :
from math equations;we know:



so if I translate it into code: does tr(grad(U))*I equal fvc::div(U)*I?

Thanks!!