[WJPeace1]

Hi all,

I'm trying to simulate a stirred tank using simpleFoam. I've managed to get it working with the k-epsilon model, but have now tried to switch it to the k-omega SST model as that can more accurately capture the flow dynamics within a stirred tank. I have had issues with my solution not converging, particularly noticeable in the bounding k and omega values.
My initial conditions are shown below:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    top
    {
        type            slip;
    }
	
    walls
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
	
	baffles
	{
		type			kqRWallFunction;
		value			uniform 0;
	}

    impeller
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
}


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

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    top
    {
        type            slip;
    }
	
    walls
    {
        type            nutLowReWallFunction;
        value           uniform 0;
    }
	
	baffles
	{
		type			nutLowReWallFunction;
		value			uniform 0;
	}

    impeller
    {
        type            nutLowReWallFunction;
        value           uniform 0;
    }
}



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

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

dimensions          [0 2 -2 0 0 0 0];

internalField       uniform 0;

boundaryField
{
    top
    {
        type            zeroGradient;
    }
	
    walls
    {
        type            zeroGradient;
    }
	
	baffles
	{
		type			zeroGradient;
	}

    impeller
    {
        type            zeroGradient;
    }
}



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

oamFile
{
    version     2.0;
    format      binary;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
    top
    {
        type            slip;
    }
	
    walls
    {
        type            noSlip;
    }
	
	baffles
	{
		type			noSlip;
	}

    impeller
    {
        type            noSlip;
    }
}


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

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 2.659;

boundaryField
{
    top
    {
        type            slip;
    }
	
    walls
    {
        type            omegaWallFunction;
        value           $internalField;
    }
	
	baffles
	{
		type			omegaWallFunction;
		value			$internalField;
	}

    impeller
    {
        type            omegaWallFunction;
        value           $internalField;
    }
}


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

Any help would be greatly appreciated!