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
Code:
/*--------------------------------*- 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
Code:
/*--------------------------------*- 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
Code:
/*--------------------------------*- 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;
}
}
// ************************************************************************* //