December 8, 2014, 00:18
|
Implementing kEpsilon turbulance mode with the airfoil example.....
|
#1
|
Senior Member
CFD
Join Date: Nov 2010
Location: United States
Posts: 243
Rep Power: 17
|
Hey Foamers,
I have implemented kOmega turbulence model in the airfoil2d example. However, when I switched to kepsilon turbulence model following all required steps, simpleFoam was unable to run. In the 0 folder, I put k and epsilon folders, and I deleted the file nuTilda and that is based on what I have read on the CFD forum.
My k is:
Quote:
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0.01;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.01;
}
outlet
{
type zeroGradient;
}
walls
{
type kqRWallFunction;
value uniform 0.01;
}
frontAnBack
{
type empty;
}
}
|
My epsilon file is
Quote:
dimensions [ 0 2 -3 0 0 0 0 ];
internalField uniform 4.85;
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type zeroGradient;
}
walls
{
type epsilonWallFunction;
value $internalField;
}
frontAndBack
{
type empty;
}
}
|
fscheme is
Quote:
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}
|
fsolution is
Quote:
solvers
{
p GAMG
{
tolerance 1e-16;
relTol 0.01;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
minIter 0;
maxIter 1000;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
scaleCorrection true;
directSolveCoarsest false;
};
Urel smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};
k smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};
epsilon smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};
omega smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
Urel 0.6;
k 0.5;
epsilon 0.5;
omega 0.5;
R 0.5;
nuTilda 0.5;
}
}
|
Can somebody explain why it is not working?
Thank you....
|
|
|