[Jonathan]

Hi everyone,

I am solving a turbine case in OF2.1.1 using SRFSimpleFoam using kOmegaSST. I am having some bounding and convergence issues but am confused about a particular term in the divSchemes subdictionary in fvSchemes:

Does anyone know which equation:

div((nuEff*dev(T(grad(Urel))))) Gauss upwind; //linear;

comes from and why i cannot using Gauss upwind to discretize it - I get fatal read beyond EOF error if i dont use linear with it.

thanks very much in advance
jonathan

[immortality]

Hi Jonathan
I had a problem like this with tauMC in rhoCentralFoam and Bruno said me that its a tensor field not a scalar or vector field that can use upwind schemes.
if you can get more information or others its nice to share.

[Jonathan]

hi,

so the term comes from the turbulence model (kOmegaSST):

From kOmegaSST.H

        //- Return the source term for the momentum equation
        virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;

From kOmegaSST.C

tmp<fvVectorMatrix> kOmegaSST::divDevReff(volVectorField& U) const
{
    return
    (
      - fvm::laplacian(nuEff(), U)
      -  fvc::div(nuEff()*dev(T(fvc::grad(U))))
    );
}

As immortality says, this might be a tensor field, and so the standard divSchemes as per the UG do not apply, but i am still not 100% sure.

[wyldckat]

Greetings to all!

@Jonathan: The post Ehsan was talking about is this one: http://www.cfd-online.com/Forums/ope...tml#post429375 - post #4.
And yes, this is basically the same pattern of mathematics going on in this case.

If you run the following command, you'll see that all tutorials in OpenFOAM use "Gauss linear" for this particular kind of div:

find $FOAM_TUTORIALS -name fvSchemes | xargs grep "T(grad(U"

Best regards,
Bruno

[Jonathan]

Originally Posted by wyldckat

Greetings to all!

@Jonathan: The post Ehsan was talking about is this one: http://www.cfd-online.com/Forums/ope...tml#post429375 - post #4.
And yes, this is basically the same pattern of mathematics going on in this case.

If you run the following command, you'll see that all tutorials in OpenFOAM use "Gauss linear" for this particular kind of div:

find $FOAM_TUTORIALS -name fvSchemes | xargs grep "T(grad(U"

Best regards,
Bruno

Hi Bruno,

Thanks very much - your help is much appreciated.

Best regards
Jon