[prodo]

Hello,

I made a mesh by cfMesh, like attached picture.

and then solve this case by icoFoam

initial conditon is v=0.1m/s

but, you can see there is non-physical result on the mesh grading region.

so, how can i get proper result in this case?

this is my fvSolution, fvSchemes file.

thank you

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

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

divSchemes
{
    default                 Gauss upwind;
    div(phi,U)                  Gauss linearUpwindV Gauss linear;
    div((nuEff*dev(grad(U).T())))     Gauss linear;

}

laplacianSchemes
{
    default         Gauss linear orthogonal;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default     corrected;
 }

fluxRequired
{
    default         no;
    p               ;
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver           GAMG;
        tolerance        1e-7;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;

  /*    solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0;             */
    }

    U
    {


solver    smoothSolver;
        smoother    GaussSeidel;
        tolerance    1e-8;
        relTol        0;
        nSweeps        1;
/*        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0;            */
    }
}

PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}
relaxationFactors
{
    p               0.3;
    U               0.7;
    k               0.7;
    epsilon         0.7;
    omega            0.7;
}

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

[franjo_j]

Hello,

What happens when you replace linearUpwindV or the momentum equation with vanLeer, Gamma or limitedLinear?
In addition, what happens when you increase the number of sweeps in solver?

Regards,

Franjo