[shiro44]

Hi guys, I am new to OpenFoam.
I was running a simulation for a steady flow using SimpleFoam.
The geometry is quite simple, two finite vertical walls and a costant horizontal mean wind upon them.
The problem is that k oscillates a lot vertically and its values are like, in general, 2 orders of magnitude higher than expected. Also it's converged after 9000 iterations.


The fvSchemes is:

 

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div(phi,omega)  bounded Gauss limitedLinear 1;
    div(phi,v2)     bounded Gauss limitedLinear 1;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div(nonlinearStress) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

The fvSolution is:

solvers
{
    U
    {
            solver          PBiCG;
            preconditioner  DILU;
            tolerance       1e-6;
            relTol          0;
    }
    "(p|p_rgh)"
    {
            solver          PCG;
            preconditioner  DIC;
            tolerance       1e-6;
            relTol          0;
    }


    "(U|k|epsilon|omega|f|v2)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0.1;
    }
}

SIMPLE
{
        pRefCell  0;
        pRefValue 0;
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-3;
        U               1e-4;
        "(k|epsilon|omega|f|v2)" 1e-3;
    }

}

relaxationFactors
{
    equations
    {
        U               0.75
    }
}

Any suggestions?