CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Problem of InterFoam with LES SpalartAllmarasIDDES (http://www.cfd-online.com/Forums/openfoam/117671-problem-interfoam-les-spalartallmarasiddes.html)

keepfit May 13, 2013 09:27

Problem of InterFoam with LES SpalartAllmarasIDDES
 
2 Attachment(s)
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:
http://24.media.tumblr.com/6e4795e32...saho3_1280.png


The results of Laminar and KOmegaSST look more reasonable:

Laminar:
http://24.media.tumblr.com/a5dad8a04...saho2_1280.png

and InterFoam KOmegaSST:
http://25.media.tumblr.com/5a136da32...saho1_1280.png


__________________________________________________

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:

Code:

/*--------------------------------*- 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
:

Code:

/*--------------------------------*- 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

Code:

/*--------------------------------*- 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:

Attachment 21691
Constant
Attachment 21692

Best,

David

timo_IHS May 13, 2013 11:08

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 May 13, 2013 12:53

Quote:

Originally Posted by timo_IHS (Post 427225)
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 August 29, 2013 11:21

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


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


All times are GMT -4. The time now is 20:42.