[hamsteri15]

I'm comparing a couple of sgs models for turbulent pipe flow and so far I have been able to run Smagorinsky and oneEqEddy properly. For my last sgs model I would like to run mixedSmagorinsky but my fvschemes produce a long error:

--> FOAM FATAL IO ERROR:
keyword div((dev((((surfaceSum((magSf*interpolate(sqr(U))) )|surfaceSum(magSf))-sqr((surfaceSum((magSf*interpolate(U)))|surfaceSum (magSf))))+(((0.666667)*((0.5*((surfaceSum((magSf* interpolate(magSqr(U))))|surfaceSum(magSf))-magSqr((surfaceSum((magSf*interpolate(U)))|surface Sum(magSf)))))+((((2*ck)|ce)*sqr(delta))*magSqr(de v(symm(grad(U)))))))-(nuSgs*twoSymm(grad(U))))))+(-nuEff*dev(twoSymm(grad(U)))))) is undefined in dictionary "/home/ubuntu/OpenFOAM/ubuntu-2.2.2/run/Kandi/putki/system/fvSchemes.divSchemes"

my divschemes:
divSchemes
{
default none;
div(phi,U) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div((nu*dev(T(grad(U))))) Gauss linear;
div((grad(U)&U)) Gauss linear;
div(phi,k) Gauss limitedLinear 1; //*
div(phi,B) Gauss limitedLinear 1; //*
div(B) Gauss linear; //*
}

[vut]

Dear hamsteri15,

I had the same error running mixed Smagorinsky model.

Have you got the solution?

Regards,

vut

[msuaeronautics]

The code is just telling you that there is a term hanging around that you have not told it how to discretize.

As a remedy, I would just place Gauss linear in the default category.

In other words, instead of:

default none;

Place:

default Gauss linear;

That will discretize this remaining term. Now, as to what that term means- I do not know, because I do not know much about the implementation.

[vut]

I tried almost the same way: paste these formula and use Gauss linear.

Thanks anyway. It seems clearer now for fvScheme for me

Nice weekend,

vut