[keepfit]

Hi Foamers,

I did some simple 2D interFoam simulations to test and compare different turbulence models. So far the models such as dummy laminar, RAS KOmegaSST and LES SpalartAllmarasIDDES were tested to model stream flow over a slope.

There is no problem with laminar and ras KOmegaSST, however, as there is no examples of InterFoam LES SpalartAllmarasIDDES model, I tried to create such multiphase case, but the result looks odd, and the simulation blows up after 2 second.

InterFoam - LES SpalartAllmarasIDDES:



The results of Laminar and KOmegaSST look more reasonable:

Laminar:


and InterFoam KOmegaSST:



__________________________________________________

I tried different fvsolution and fvschemes, but could not figure out this issue.

So I wonder, is there any available cases related InterFoam LES SpalartAllmarasIDDES ?

fvSchemes:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default backward;
}

d2dt2Schemes
{
}

gradSchemes
{
    default         Gauss linear;

    grad(nuTilda)   cellLimited Gauss linear 1;
    grad(U)         cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss LUST unlimitedGrad(U);
    //div(phi,U)      Gauss linearUpwind unlimitedGrad(U);
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;

    div((nuEff*dev(T(grad(U))))) Gauss linear;

    div(rho*phi,U)  Gauss linearUpwind grad(U);
    div((muEff*dev(T(grad(U))))) Gauss linear;
    div(phi,B)      Gauss limitedLinear 1;
    div(B)          Gauss linear;


}

laplacianSchemes
{
    default         Gauss linear limited 0.33;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.33;
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha1;
}

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

fvSolution:

 /*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    pcorr
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-05;
            relTol          0;
            smoother        DICGaussSeidel;
            nPreSweeps      0;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration false;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
        tolerance       1e-10;
        relTol          0;
        maxIter         100;
    }

    
    p_rgh
    {
        solver          GAMG;
        tolerance       1e-07;
        relTol          0.01;
        smoother        DIC;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    p_rghFinal
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-07;
            relTol          0;
            nVcycles        2;
            smoother        DICGaussSeidel;
            nPreSweeps      2;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
        tolerance       1e-07;
        relTol          0;
        maxIter         20;
    }


    U
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-08;
        relTol           0;
    };

    UFinal
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-08;
        relTol           0;
    };

    k
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-07;
        relTol           0;
    };

    B
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-07;
        relTol           0;
    };


    nuTilda
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-07;
        relTol           0;
    };

    "(nuTilda)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0;
    };
}

PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 1;
}

PIMPLE
{
    momentumPredictor no;
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;

    nAlphaCorr      1;
    nAlphaSubCycles 3;
    cAlpha          1;

}

relaxationFactors
{
    default         0;
    p               0.3;
    U               0.7;
    nuTilda         0.7;

    fields
    {
    }
    equations
    {
    }
}


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

constant/LESProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      LESProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

LESModel            SpalartAllmarasIDDES;

delta               cubeRootVol;

printCoeffs         on;

laminarCoeffs
{
}

oneEqEddyCoeffs
{
    ck               0.07;
    ce               1.05;
}

dynOneEqEddyCoeffs
{
    ce               1.05;
    filter          simple;
}

locDynOneEqEddyCoeffs
{
    ce               1.05;
    filter          simple;
}

SmagorinskyCoeffs
{
    ce               1.05;
    ck               0.07;
}

Smagorinsky2Coeffs
{
    ce               1.05;
    ck               0.07;
    cD2              0.02;
}

spectEddyViscCoeffs
{
    ce               1.05;
    cB               8.22;
    cK1              0.83;
    cK2              1.03;
    cK3              4.75;
    cK4              2.55;
}

dynSmagorinskyCoeffs
{
    ce               1.05;
    filter          simple;
}

mixedSmagorinskyCoeffs
{
    ce               1.05;
    ck               0.07;
    filter          simple;
}

dynMixedSmagorinskyCoeffs
{
    ce               1.05;
    filter          simple;
}

LRRDiffStressCoeffs
{
    ce               1.05;
    ck               0.09;
    c1               1.8;
    c2               0.6;
}

DeardorffDiffStressCoeffs
{
    ce               1.05;
    ck               0.09;
    cm               4.13;
}

SpalartAllmarasCoeffs
{
    alphaNut         1.5;
    Cb1              0.1355;
    Cb2              0.622;
    Cw2              0.3;
    Cw3              2;
    Cv1              7.1;
    Cv2              5.0;
    CDES             0.65;
    ck               0.07;
}

SpalartAllmarasDDESCoeffs
{
    alphaNut         1.5;
    Cb1              0.1355;
    Cb2              0.622;
    Cw2              0.3;
    Cw3              2.0;
    Cv1              7.1;
    Cv2              5.0;
    CDES             0.65;
    ck               0.07;
}

SpalartAllmarasIDDESCoeffs
{
     delta  IDDESDelta;

     IDDESDeltaCoeffs
     {
         deltaCoeff       1;

         //according to src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C
         cw               0.15;
     }

     alphaNut         1.5;
     kappa            0.4187;
     Cb1              0.1355;
     Cb2              0.622;
     Cw2              0.3;
     Cw3              2.0;
     Cv1              7.1;
     Cv2              5.0;
     CDES             0.65;
     ck               0.07;

//These are the default values
//     sigmaNut        0.66666;
//     kappa           0.41;
//     Cb1             0.1355;
//     Cb2             0.622;
//     Cv1             7.1;
//     Cv2             5;
//     CDES            0.65;
//     ck              0.07;
//     Cw2             0.3;
//     Cw3             2;
//     fwStar          0.424;
//     cl              3.55;
//     ct              1.63;

}

cubeRootVolCoeffs
{
    deltaCoeff      1;
}

PrandtlCoeffs
{
    delta           cubeRootVol;
    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }
    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        maxDeltaRatio   1.1;
    }
    Cdelta           0.158;
}

vanDriestCoeffs
{
    delta           cubeRootVol;
    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }
    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        maxDeltaRatio   1.1;
    }
    Aplus            26;
    Cdelta           0.158;
}

smoothCoeffs
{
    delta           cubeRootVol;
    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }
    maxDeltaRatio   1.1;
}

kappa            0.4187;

wallFunctionCoeffs
{
    E                9;
}

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

Other files can be found in the attached files:

0.tar.gz
Constant
system.tar.gz

Best,

David

[timo_IHS]

One problem is, that you are using 2D for IDDES (LES).
It would be interesting to compare this simulation with a 3D simulation.
To which wall distance does your IDDES simulation run in RANS?
Can you compare the viscosity of your simulations?

[keepfit]

Originally Posted by timo_IHS

One problem is, that you are using 2D for IDDES (LES).
It would be interesting to compare this simulation with a 3D simulation.
To which wall distance does your IDDES simulation run in RANS?
Can you compare the viscosity of your simulations?

The frontAndback type is set "wall" in LES (boundary) other than "empty" in the case of 2D, so is it considered 3D?

As for wall distance, can you explain it for more detail?

[vishal_s]

I an new to Openfoam , I am trying to use Spalart-Allmaras model for my nozzle flow in interFoam. I incoporated nuTilda and nutildaFinal in the fvsolutions of the system file, but I a getting this error. can any one tell me, how to solve this error.


P { margin-bottom: 0.21cm; } Reading g
Calculating field g.h

time step continuity errors : sum local = 9.77995e-05, global = -9.77995e-05, cumulative = -9.77995e-05
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.94573e-11, No Iterations 194
time step continuity errors : sum local = 6.7998e-15, global = -4.64687e-17, cumulative = -9.77995e-05
Courant Number mean: 0.00278549 max: 0.285502

Starting time loop

Courant Number mean: 0.00193436 max: 0.198265
Interface Courant Number mean: 0 max: 0
deltaT = 0.000694444
Time = 0.000694444

MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220219 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220388 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220558 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220728 Min(alpha1) = 0 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0419123, No Iterations 85
time step continuity errors : sum local = 4.60931e-06, global = 9.84271e-09, cumulative = -9.77897e-05
DICPCG: Solving for p_rgh, Initial residual = 0.000734712, Final residual = 3.64035e-05, No Iterations 65
time step continuity errors : sum local = 1.06533e-05, global = 2.14005e-06, cumulative = -9.56496e-05
DICPCG: Solving for p_rgh, Initial residual = 0.00326193, Final residual = 8.82474e-08, No Iterations 123
time step continuity errors : sum local = 4.04865e-09, global = 5.87426e-11, cumulative = -9.56496e-05
smoothSolver: Solving for nuTilda, Initial residual = 0, Final residual = 0, No Iterations 0
--> FOAM Warning :
From function tmp<volScalarField> SpalartAllmaras::k() const
in file SpalartAllmaras/SpalartAllmaras.C at line 243
Turbulence kinetic energy not defined for Spalart-Allmaras model. Returning zero field
--> FOAM Warning :
From function tmp<volScalarField> SpalartAllmaras::k() const
in file SpalartAllmaras/SpalartAllmaras.C at line 243
Turbulence kinetic energy not defined for Spalart-Allmaras model. Returning zero field
ExecutionTime = 0.18 s ClockTime = 1 s

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

This is the what I incorporated in the SYSTEM/FVSOLUTIONS

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

solvers
{
pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

p_rghFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
}



nuTildacorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0;
}

nuTilda
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

nuTildaFinal
{
solver smoothSolver;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
}

}

PIMPLE
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 2;
}


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