[callumso]

Hi dear foamers,

I am running simulations for a complex geometry, mesh size around 5 million. I use k-omega SST model and wall function. From the residuals it seems all the equations converge all right. But omega is just always extremely large and goes to infinity in the orifice where mixing of flow occurs and is of interest. Consequently k goes to almost zero in this region where turbulence should be really high. Strangely even omega equation seems to converge.

Have anyone encountered something like this? below is the boundary condition I used for k and omega:

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

internalField   uniform 0.00027;

boundaryField
{
    PILOT
    {
        type            compressible::kqRWallFunction;
        value           uniform 0.00027;
    }
    PLENUM
    {
        type            zeroGradient;
    }
    INLET
    {
        type            fixedValue;
        value           uniform 0.00027;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    DUCT
    {
        type            compressible::kqRWallFunction;
        value           uniform 0.00027;
    }
    VANES
    {
        type            compressible::kqRWallFunction;
        value           uniform 0.00027;
    }
    PASSAGE
    {
        type            compressible::kqRWallFunction;
        value           uniform 0.00027;
    }
    BELLMOUTH
    {
        type            compressible::kqRWallFunction;
        value           uniform 0.00027;
    }
}

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

internalField   uniform 2.6;

boundaryField
{
    PILOT
    {
        type            compressible::omegaWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.00027;
    }
    PLENUM
    {
        type            zeroGradient;
    }
    INLET
    {
        type            fixedValue;
        value           uniform 0.1813;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    DUCT
    {
        type            compressible::omegaWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.00027;
    }
    VANES
    {
        type            compressible::omegaWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.00027;
    }
    PASSAGE
    {
        type            compressible::omegaWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.00027;
    }
    BELLMOUTH
    {
        type            compressible::omegaWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.00027;
    }
}

The numerical schems look like below:

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phid,p)     Gauss limitedLinear 1;
    div(phi,K)      Gauss linear;
    div(phi,h)      Gauss upwind;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,epsilon) Gauss limitedLinear 1;
    div(phi,R)      Gauss upwind;
    div(phi,omega)  Gauss limitedLinear 1;
    div((rho*R))    Gauss linear;
    div(R)          Gauss linear;
    div(U)          Gauss linear;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         none;
    laplacian(muEff,U) Gauss linear corrected;
    laplacian(mut,U) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DomegaEff,omega) Gauss linear corrected;
    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
    laplacian(alphaEff,h) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-08;
        relTol          0.01;
    }

    pFinal
    {
        $p;
        tolerance       1e-08;
        relTol          0;
    }

    "(rho|U|h|k|epsilon|omega)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0.01;
    }

    "(rho|U|h|k|epsilon|omega)Final"
    {
        $U;
        tolerance       1e-08;
        relTol          0;
    }

}

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors 4;
    nCorrectors     2;
    nNonOrthogonalCorrectors 1;
    rhoMin          rhoMin [ 1 -3 0 0 0 ] 0.5;
    rhoMax          rhoMax [ 1 -3 0 0 0 ] 2.0;

    residualControl
    {
        "(U|k|epsilon)"
        {
            relTol          0;
            tolerance       0.000001;
        }
    }
}

relaxationFactors
{
    fields
    {
        "p.*"           0.7;
        "rho.*"         0.9;
    }
    equations
    {
        "(U|h|k|epsilon|omega).*" 0.7;
    }
}

I fiddled with the numerical scheme a bit (e.g. use upwind scheme, use very low relaxation factors, etc). But nothing changed. I am still getting this weird infinite omega. Below is the residual trace:

Courant Number mean: 0.002384836 max: 10.06362
deltaT = 4.105397e-05
Time = 0.001018885

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.009422393, Final residual = 1.336213e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.005013185, Final residual = 1.673009e-06, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.004954721, Final residual = 1.658077e-06, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.008704607, Final residual = 1.068494e-08, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.02191301, Final residual = 8.992473e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 0.0001136703, Final residual = 1.065346e-06, No Iterations 17
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.221163e-06, global = -8.451316e-07, cumulative = -0.0001596296
rho max/min : 1.193396 1.160766
DICPCG: Solving for p, Initial residual = 0.006488465, Final residual = 2.687909e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 3.438664e-05, Final residual = 3.410989e-07, No Iterations 36
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.218193e-06, global = -8.500729e-07, cumulative = -0.0001604797
rho max/min : 1.193396 1.160765
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.006534682, Final residual = 9.631645e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.004523631, Final residual = 6.797623e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.004480604, Final residual = 6.662678e-07, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.009199319, Final residual = 1.208279e-08, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.008744652, Final residual = 3.671174e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 4.807335e-05, Final residual = 4.056757e-07, No Iterations 17
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.218818e-06, global = -8.48645e-07, cumulative = -0.0001613283
rho max/min : 1.193396 1.160793
DICPCG: Solving for p, Initial residual = 0.002611529, Final residual = 1.103145e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 1.45108e-05, Final residual = 1.270581e-07, No Iterations 42
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.217819e-06, global = -8.509643e-07, cumulative = -0.0001621793
rho max/min : 1.193396 1.160792
PIMPLE: iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 0.004572449, Final residual = 6.991547e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.003637689, Final residual = 2.368283e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.003609003, Final residual = 2.530072e-07, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.005076553, Final residual = 6.463705e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.004023061, Final residual = 1.699759e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 2.310766e-05, Final residual = 2.116908e-07, No Iterations 17
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.217778e-06, global = -8.498158e-07, cumulative = -0.0001630291
rho max/min : 1.193396 1.160799
DICPCG: Solving for p, Initial residual = 0.001204748, Final residual = 5.135975e-06, No Iterations 4
DICPCG: Solving for p, Initial residual = 6.987517e-06, Final residual = 6.530037e-08, No Iterations 41
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.217191e-06, global = -8.507411e-07, cumulative = -0.0001638799
rho max/min : 1.193396 1.160799
PIMPLE: iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 0.003182333, Final residual = 5.051115e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00280553, Final residual = 2.892812e-09, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.002785603, Final residual = 2.954405e-09, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.00268957, Final residual = 3.386828e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.00198696, Final residual = 8.607319e-06, No Iterations 4
DICPCG: Solving for p, Initial residual = 1.217912e-05, Final residual = 1.181804e-07, No Iterations 20
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.217629e-06, global = -8.508334e-07, cumulative = -0.0001647307
rho max/min : 1.193396 1.160801
DICPCG: Solving for p, Initial residual = 0.0005952442, Final residual = 2.607833e-06, No Iterations 4
DICPCG: Solving for p, Initial residual = 3.687937e-06, Final residual = 9.698358e-09, No Iterations 67
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.217111e-06, global = -8.50822e-07, cumulative = -0.0001655815
rho max/min : 1.193396 1.160801
DILUPBiCG: Solving for omega, Initial residual = 0.1111112, Final residual = 5.200637e-19, No Iterations 1
bounding omega, min: -1.173741e+38 max: 2.166221e+50 average: 8.778702e+43
DILUPBiCG: Solving for k, Initial residual = 0.048428, Final residual = 2.335364e-18, No Iterations 1
bounding k, min: -5.711563e-06 max: 0.0734092 average: 0.0005477621
PIMPLE: not converged within 4 iterations
ExecutionTime = 185.73 s ClockTime = 187 s

I went through some related threads but unfortunately none of the advice really helps.

I run the same case on Fluent and there is no problem at all. I tried a finer mesh, 12million cells. But it still doesn't work on OpenFOAM. Can anyone help me out with this? Thank you in advance.

Callum