CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Problem of InterFoam with LES SpalartAllmarasIDDES

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 13, 2013, 10:27
Question Problem of InterFoam with LES SpalartAllmarasIDDES
  #1
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13
keepfit is on a distinguished road
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:

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:

0.tar.gz
Constant
system.tar.gz

Best,

David

Last edited by keepfit; May 13, 2013 at 11:05.
keepfit is offline   Reply With Quote

Old   May 13, 2013, 12:08
Default
  #2
Member
 
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16
timo_IHS is on a distinguished road
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?
timo_IHS is offline   Reply With Quote

Old   May 13, 2013, 13:53
Default
  #3
Senior Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 13
keepfit is on a distinguished road
Quote:
Originally Posted by timo_IHS View Post
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?
keepfit is offline   Reply With Quote

Old   August 29, 2013, 12:21
Default
  #4
Member
 
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13
vishal_s is on a distinguished road
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;
}


// ************************************************** *********************** //
vishal_s is offline   Reply With Quote

Reply

Tags
ddes, interfoam, spalart allmaras

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
rhoPimpleFoam LES problem with OF-2.2.0 halowine OpenFOAM Running, Solving & CFD 1 May 17, 2013 17:13
Parallel motorBike LES computing problem. xiexiehezuo OpenFOAM Running, Solving & CFD 1 August 22, 2012 17:44
Wmake problem interFoam solver feijooos OpenFOAM Running, Solving & CFD 4 December 8, 2008 12:01
LES problem Jose FLUENT 5 September 29, 2003 06:04
problem while looping over cells in LES Jimmy FLUENT 0 September 5, 2002 16:33


All times are GMT -4. The time now is 09:32.