[usermch2]

Hi,
I am using rhoSimpleFoam to simulate the flow between two coencentric cylinders with the outer one rotating. My simulation for constant viscosity went well.

{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
 specie
    {
        nmoles          1;
        molWeight       18;
    }
    thermodynamics
    {
        Cp              4180;
        Hf              0;
        Cv              4150;
    }
    equationOfState
    {
       rho             1000;
    }
    transport
    {
        mu              0.1;
        Pr              7.56;
    }
}

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

When I set the polynomial values for viscosity my simulaion is not converging. I tried setting muCoeffs<8> (0.1 0 0 0 0 0 0 0) to see if it will converge but didn't.

My thermophysicalProperties file:

thermoType
{   
    type            heRhoThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState icoPolynomial;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
specie
    {
        nmoles          1;
        molWeight       18;
    }
    equationOfState
    {
        rhoCoeffs<8>    ( 1000 0 0 0 0 0 0 0 );
    }
    thermodynamics
    {
         Hf              0;   
         Sf              0;
         CpCoeffs<8>     ( 4200 0 0 0 0 0 0 0 );
    }
    transport
    {
         muCoeffs<8>     ( 0.275  0.0080882353  0.0000108131  0.0000000064 0 0 0 0 );
         kappaCoeffs<8>  ( 0.598 0 0 0 0 0 0 0 );
    }
}

T file:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 350;

boundaryField
{
    innercylinder
    {
        type            fixedValue;
        value           uniform 340;
    }
    outercylinder
    {
        type            fixedValue;
        value           uniform 350;
    }
    top
    {
		type empty;
	}
    bottom
    {
		type empty;
	}
}

U file:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{

    innercylinder
    {
          type noSlip;
    }

    outercylinder
    {

        type            rotatingWallVelocity;
        origin          ( 0 0 0 );
        axis            ( 0 0 1 );
        omega           4; // rad/s

    }

    top
    {
		type empty;
	}
    bottom
    {
		type empty;
	}

}

fvsolution:

solvers
{
    p
    {
        solver          GAMG;
        smoother        FDIC;
	maxIter		10;
        tolerance       1e-6;
        relTol          0.0001;
    }

    "(U|e|k|epsilon)"
    {
        solver          GAMG;
        smoother        DILU;  
cacheAgglomeration      yes;
nCellsInCoarsestLevel   60;
nPreSweeps              0;
maxPreSweeps            2;
nPostSweeps             2;
maxPostSweeps           2;
nFinestSweeps           2;
	maxIter 	50;
        tolerance       1e-8;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    pMin      1; 
    pMax      1.5;
	pRefCell 0;
    transonic       no;
    consistent      no;
	pRefValue 1e5;
    residualControl
    {
        p               1e-4;
        U               1e-6;
	e		1e-6;
        "(k|omega|epsilon|h)" 1e-8;
    }
}

relaxationFactors
{
    fields
    {
        p               0.8;
    }
    equations
    {
        U               0.6;
        e               0.8;
        k               0.8;
        epsilon         0.8;
    }
}

Could someone help me or give hint? I really woul appreciate it.