[kurne]

Dear All

I am using OpenFOAM 1.5-dev and my case is of turbo machinery and using MRFSimpleFoam solver.I have given following BC
Inlet

U            flowRateInletVelocity
P             zeroGradient
k             fixedValue;
              value  uniform 1.057;
epsilon     fixedValue;
              uniform 11290;

Outlet

U            zeroGradient;
P             fixedMean
k             zeroGradient;;
epsilon     fzeroGradient;

fvschemes

ddtSchemes
{
    default steadyState;
}

gradSchemes
{
    default                                 Gauss linear;
    grad(p)                                Gauss linear;
    grad(U)                                Gauss linear;
}

divSchemes
{
    default                                 none;
    div(phi,U)                             Gauss upwind;
    div(phi,k)                             Gauss upwind;
    div(phi,epsilon)                      Gauss upwind;
    div((nuEff*dev(grad(U).T())))   Gauss linear;
}

laplacianSchemes
{
    default                                  none;
    laplacian(nuEff,U)                   Gauss linear corrected;
    laplacian((1|A(U)),p)               Gauss linear corrected;
    laplacian(DkEff,k)                   Gauss linear corrected;
    laplacian(DepsilonEff,epsilon)    Gauss linear corrected;
    laplacian(1,p)                         Gauss linear corrected;
}

interpolationSchemes
{
    default                                   linear;
    interpolate(U)                         linear;
}

snGradSchemes
{
    default                                   corrected;
}

fluxRequired
{
    default                                   no;
    p;
}

fvSolution

{
      p GAMG
    {
        tolerance                            1e-8;
        relTol                                 0.05;

        smoother                            GaussSeidel;

        cacheAgglomeration              true;

        nCellsInCoarsestLevel           20;
        agglomerator                      faceAreaPair;
        mergeLevels                       1;
    };

    U smoothSolver
    {
        smoother                           GaussSeidel;
        nSweeps                            2;
        tolerance                           1e-7;
        relTol                                0.01;
    };

    k smoothSolver
    {
        smoother                          GaussSeidel;
        nSweeps                           2;
        tolerance                          1e-7;
        relTol                               0.1;
    };

    epsilon smoothSolver
    {
        smoother                          GaussSeidel;
        nSweeps                           2;
        tolerance                         1e-7;
        relTol                               0.1;
    };
}

SIMPLE
{
    nNonOrthogonalCorrectors      2;
    pRefCell                              0;
    pRefValue                            0;
}

relaxationFactors
{
    p                                       0.5;
    U                                      0.5;
    k                                      0.3;
    epsilon                              0.3;
}

I am facing the problem of convergence of solver.I am attaching the residual plots and log file.I have done number of changes in it but didn't get the convergence.Will anybody help me.

I am waiting.Thanks In Advance.

[kurne]

Dear All

Yet i didn't get rid out of it.I am need of help from you people.

I am waiting.

Thanks In Advance.

[betakv]

Try to use inletOutlet boundary condition for velocity outlet instead zeroGradient