|
[Sponsors] |
Problem of InterFoam with LES SpalartAllmarasIDDES |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 13, 2013, 10:27 |
Problem of InterFoam with LES SpalartAllmarasIDDES
|
#1 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14 |
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 { } } // ************************************************************************* // 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; } // ************************************************************************* // 0.tar.gz Constant system.tar.gz Best, David Last edited by keepfit; May 13, 2013 at 11:05. |
|
May 13, 2013, 12:08 |
|
#2 |
Member
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16 |
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? |
|
May 13, 2013, 13:53 |
|
#3 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14 |
Quote:
As for wall distance, can you explain it for more detail? |
||
August 29, 2013, 12:21 |
|
#4 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
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; } // ************************************************** *********************** // |
|
Tags |
ddes, interfoam, spalart allmaras |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |