[zaripov]

I'm solving simple case with flow in сoaxial circular tube, and found one strange problem.

The problem is that rho field is constant and equals 1 everywhere except walls. Although using equation of state density must be near 1.12.

Solver is rhoSimplecFoam.

fvSolution:

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 20;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    "(U|h|k|epsilon)"
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 20;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    rhoMin          rhoMin [1 -3 0 0 0] 0.1;
    rhoMax          rhoMax [1 -3 0 0 0] 1.0;
    transonic       no;
}

relaxationFactors
{
    fields
    {
        p               1;
        rho             1;
    }
    equations
    {
        U               0.9;
        h               0.8;
        k               0.9;
        epsilon         0.9;
    }
}

thermophysicalProperties:

thermoType      hPsiThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>;

mixture
{
    specie
    {
        nMoles      1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1007;
        Hf          0;
    }
    transport
    {
        As          1.4792e-06;
        Ts          116;
    }
}

What else should I show?

P. S. Sorry for my English.

[nimasam]

well in fvSolution you defined the min and max of density , increase rhoMax, forexample put it 2

[zaripov]

Originally Posted by nimasam

well in fvSolution you defined the min and max of density , increase rhoMax, forexample put it 2

Yes, I see it this night. Stupid mistake.

Now everything is ok.