[bingo641]

Dear Openfoam users,

I am doing the simulation of a diffuser with simplefoam and turbulence model. I have tried several times both by k omega and k omega SST, but no result makes sense.

In terms of the k omega model, the inlet velocity profile is wrong. In terms of the k omega SST model, I have got many different weird results (most of them has irregular vortex or singularity with larger velocity).

Here are my settings for simulation.

k: fixedvalue (a small value like 1*10^-12)

omega: omegaWallFunction

nut: calculated 0 (I also tried nutUSpaldingWallFunction)

boundary: top : wall
bottom: wall
inlet: patch
outlet: patch
frontAndBack: empty

fvschemes:
ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,k) bounded Gauss linearUpwind grad(k); (also tried //limitedLinear 1;//linearUpwind grad(k)
[/U][/I][/I] div(phi,omega) bounded Gauss linearUpwind grad(omega); // limitedLinear 1;//linearUpwind grad(omega);

div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}


// ************************************************** *********************** //

Could you please tell me what I should change? Thanks!