[jet_engine]

Hi! I'm running a kwSST simpleFoam simulation with a cylinder and a plate and apparently I'm having trouble with the convergence of the different magnitudes. The meshing process has been done with snappyHexMesh (all 3 steps run succesfully). The geometry of the case and a previous "lower-level" mesh can be seen in the attachments, as well as the plot of the residuals.

Does anybody have any idea why the magnitudes are not converging? Maybe a fvSchemes problem, a Boundary Conditions problem, the need to run potentialFoam before to initialise the field...? The different relevant files are attached below as well.

Thanks in advance!

U:

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value      uniform (26 0 0);
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (26 0 0);
    }
    lateralWall
    {
        type            symmetry;
    }
    upperWall
    {
        type            symmetry;
    }
    lowerWall
    {
        type            symmetry;
    }
    Cylinder
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Plate_round
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

p:

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    upperWall
    {
        type            symmetry;
    }
    lowerWall
    {
        type            symmetry;
    }
    lateralWall
    {
        type            symmetry;
    }
    Cylinder
    {
        type            zeroGradient;
    }
    Plate_round
    {
        type            zeroGradient;
    }
}

k:

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 0.02535;
    }
    outlet
    {
        type            zeroGradient;
    }
    upperWall
    {
        type            symmetry;
    }
    lowerWall
    {
        type            symmetry;
    }
    lateralWall
    {
        type            symmetry;
    }
    Cylinder
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
    Plate_round
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
}

omega:

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

internalField   uniform 162.5;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 162.5;
    }
    outlet
    {
        type            zeroGradient;
    }
    upperWall
    {
        type            symmetry;
    }
    lowerWall
    {
        type            symmetry;
    }
    lateralWall
    {
        type            symmetry;
    }
    Cylinder
    {
        type            omegaWallFunction;
        value           uniform 162.5;
    }
    Plate_round
    {
        type            omegaWallFunction;
        value           uniform 162.5;
    }
}

nut:

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            calculated;
        value           uniform 5e-4;
    }
    outlet
    {
        type            zeroGradient;
    }
    upperWall
    {
        type            symmetry;
    }
    lowerWall
    {
        type            symmetry;
    }
    lateralWall
    {
        type            symmetry;
    }
    Cylinder
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }
    Plate_round
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }
}

fvSchemes:

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    grad(U)         cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss linearUpwindV grad(U);
    div(phi,k)      bounded Gauss upwind;
    div(phi,omega)  bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

fvSolution:

solvers
{
    p
    {
        solver           GAMG;
        tolerance        1e-8;
        relTol           0.01;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;
    }

    Phi
    {
        $p;
    }

    U
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    }

    k
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    }

    omega
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 5;
    consistent yes;

    residualControl
    {
        p               1e-7;
        U               1e-7;
        "(k|epsilon|omega)" 1e-7;
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
    equations
    {
        U               0.9;
        k               0.7;
        omega           0.7;
    }
}

cache
{
    grad(U);
}