[namit91]

Hi everybody,

I am am new to OpenFoam and am working on simulating a mixing tank using SRFSimpleFoam. The agitator rotates at a speed of 200 rpm. The geometry is closed (surfaceCheck) and the the mesh is also OK (checkMesh).The maximum orthogonality is 59. The solver uses komegaSST model.

I've tried multiple Urel boundary conditions. However, my omega continuously diverges. Can someone please help?

smoothSolver: Solving for Urelx, Initial residual = 0.19671055, Final residual = 0.012671563, No Iterations 3
smoothSolver: Solving for Urely, Initial residual = 0.21467715, Final residual = 0.016775008, No Iterations 3
smoothSolver: Solving for Urelz, Initial residual = 0.34396204, Final residual = 0.020746431, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.11815782, Final residual = 0.0011343767, No Iterations 27
time step continuity errors : sum local = 0.17293179, global = 2.2537818e-017, cumulative = -3.368218e-014
smoothSolver: Solving for omega, Initial residual = 0.0025648008, Final residual = 8.5676094e-005, No Iterations 2
bounding omega, min: -153.83346 max: 265415.18 average: 5884.1895
smoothSolver: Solving for k, Initial residual = 0.16180719, Final residual = 0.012939662, No Iterations 2
ExecutionTime = 374.456 s ClockTime = 374 s


The p, Urel and omega values are below.

0/p

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

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    body
    {
        type            zeroGradient;
    }
    dish
    {
        type            zeroGradient;
    }
    agitator
    {
        type            zeroGradient;
    }
}


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

0/Urel

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

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            SRFVelocity;
        inletValue      uniform (0 0 0);
        relative        no;
        value           uniform (0 0 0);
    }
	body
    {
        type            SRFVelocity;
        inletValue      uniform (0 0 0);
        relative        no;
        value           uniform (0 0 0);
    }
	dish
    {
		type            SRFVelocity;
        inletValue      uniform (0 0 0);
        relative        no;
        value           uniform (0 0 0);
    }
	agitator
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}


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

0/omega

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

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

internalField   uniform 2.6;

boundaryField
{
    inlet
    {
        type            omegaWallFunction;
		value			$internalField;
    }
    body
    {
        type            omegaWallFunction;
		value			$internalField;
    }
    dish
    {
        type            omegaWallFunction;
		value			$internalField;
    }
    agitator
    {
        type            omegaWallFunction;
		value			$internalField;
    }
	
}


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

[namit91]

Hi,

I found something that works and gives a solution.

I just changed the bc SRFVelocity to SRFWallvelocity.