[ch_jaehnel]

Hi Foamers,


I have a problem with IDDES in interFoam. Since I moved to OF5.0 I get this error:

Attempt to cast type cubeRootVol to type IDDESDelta
    From function To &Foam::refCast(From &) [with To = Foam::LESModels::IDDESDelta, From = Foam::LESdelta] in file /sw/taurus/applications/OpenFOAM/OpenFOAM-5.0/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

This is the Outputfile:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 5.0
Exec   : /sw/taurus/applications/OpenFOAM/OpenFOAM-5.0/platforms/linux64IccDPInt32Opt/bin/interFoam -parallel
Date   : Jul 11 2018
Time   : 12:40:28
Host   : "taurusi6340.taurus.hrsk.tu-dresden.de"
PID    : 29205
I/O    : uncollated
Case   : /lustre/scratch2/s2665038/DMF_LES
nProcs : 24
Slaves : 
23
(
"taurusi6340.taurus.hrsk.tu-dresden.de.29206"
"taurusi6340.taurus.hrsk.tu-dresden.de.29207"
"taurusi6340.taurus.hrsk.tu-dresden.de.29208"
"taurusi6340.taurus.hrsk.tu-dresden.de.29209"
"taurusi6340.taurus.hrsk.tu-dresden.de.29210"
"taurusi6340.taurus.hrsk.tu-dresden.de.29211"
"taurusi6340.taurus.hrsk.tu-dresden.de.29212"
"taurusi6340.taurus.hrsk.tu-dresden.de.29213"
"taurusi6340.taurus.hrsk.tu-dresden.de.29214"
"taurusi6340.taurus.hrsk.tu-dresden.de.29215"
"taurusi6340.taurus.hrsk.tu-dresden.de.29216"
"taurusi6340.taurus.hrsk.tu-dresden.de.29217"
"taurusi6340.taurus.hrsk.tu-dresden.de.29218"
"taurusi6340.taurus.hrsk.tu-dresden.de.29219"
"taurusi6340.taurus.hrsk.tu-dresden.de.29220"
"taurusi6340.taurus.hrsk.tu-dresden.de.29221"
"taurusi6340.taurus.hrsk.tu-dresden.de.29222"
"taurusi6340.taurus.hrsk.tu-dresden.de.29223"
"taurusi6340.taurus.hrsk.tu-dresden.de.29224"
"taurusi6340.taurus.hrsk.tu-dresden.de.29225"
"taurusi6340.taurus.hrsk.tu-dresden.de.29226"
"taurusi6340.taurus.hrsk.tu-dresden.de.29227"
"taurusi6340.taurus.hrsk.tu-dresden.de.29228"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type LES
Selecting LES turbulence model SpalartAllmarasIDDES
Selecting LES delta type cubeRootVol
Selecting patchDistMethod meshWave
LES
{
    LESModel        SpalartAllmarasIDDES;
    delta           cubeRootVol;
    printCoeffs     on;
    turbulence      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;
    }
    SpalartAllmarasIDDESCoeffs
    {
        delta           IDDESDelta;
        IDDESDeltaCoeffs
        {
            deltaCoeff      1;
            cw              0.15;
        }
    }
    IDDESDeltaCoeffs
    {
        cw              0.15;
        deltaCoeff      1;
    }
    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;
    }
    Ce              1.048;
    sigmaNut        0.66666;
    Cb1             0.1355;
    Cb2             0.622;
    Cw2             0.3;
    Cw3             2;
    Cv1             7.1;
    Cs              0.3;
    CDES            0.65;
    ck              0.07;
}

The same error occurs for other delta versions. Here are my turbulence Properties:

simulationType  LES;

LES
{
LESModel            SpalartAllmarasIDDES;

delta               cubeRootVol;

printCoeffs         on;

turbulence            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;
}

SpalartAllmarasIDDESCoeffs
{
    delta IDDESDelta;
    IDDESDeltaCoeffs
    {
        deltaCoeff 1;// 1; 2 fuer hexMesh ?
        cw 0.15;
    }
}

IDDESDeltaCoeffs
{
    cw 0.15;
    deltaCoeff      1;// 1; 2 fuer hexMesh ?

}
cubeRootVolCoeffs
{
    deltaCoeff      1;// 1; 2 fuer hexMesh ?
}

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;
}

}

fvSchmes:

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default             none;

    div(rhoPhi,U)       Gauss linearUpwind grad(U);
    div(phi,alpha)      Gauss vanLeer;
    div(phirb,alpha)    Gauss linear;

    "div\(phi,(k|omega)\)"      Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

fvSolution

solvers
{
    "alpha.water.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 1;
        cAlpha          1;

        MULESCorr       yes;
        nLimiterIter    3;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          0;
    }

    "pcorr.*"
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-5;
            relTol          0;
            smoother        GaussSeidel;
        }
        tolerance       1e-5;
        relTol          0;
        maxIter         50;
    }

    p_rgh
    {
        solver           GAMG;
        tolerance        5e-9;
        relTol           0.01;

        smoother         GaussSeidel;



        maxIter          50;
    };

    p_rghFinal
    {
        $p_rgh;
        tolerance       5e-9;
        relTol          0;
    }

    "(U|k|omega).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        nSweeps         1;
        tolerance       1e-6;
        relTol          0.1;
    };
}

PIMPLE
{
    momentumPredictor no;
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
    equations
    {
        ".*" 1;
    }
}

Since I didnt find any solutions or similar problems like this yet, I am asking you. What could this mean? How can I get around this problem? Thank you for your help!



Is someone of you using IDDES in OpenFoam 5.0? Unfortunately I didn't find any tutorials for DES.



Christian

[S123]

Edited: Just found the solution to this problem!! Check below.

How to use DESModelRegions function object


Hello,

I am also facing the exact same problem. I was wondering whether you've found a solution to this or not?

Many thanks.