WJPeace1 |
July 12, 2023 12:44 |
Issue with implementing kOmegaSST model in steady state stirred tank simulation
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:
Code:
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;
}
}
// ************************************************************************* //
Code:
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;
}
}
// ************************************************************************* //
Code:
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;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Code:
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;
}
}
// ************************************************************************* //
Code:
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!
|