| ch_jaehnel |
July 11, 2018 08:26 |
CubeRootVol in IDDES for interFoam
1 Attachment(s)
Hi Foamers,
I have a problem with IDDES in interFoam. Since I moved to OF5.0 I get this error:
Code:
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:
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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:
Code:
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:
Code:
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
Code:
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 |
