[Jiricbeng]

Hi there,

please help me with the following:
I wish to calculate 3D geometry as an incompressible k-epsilon model with boundary conditions - velocity inlet and pressure outlet (which is likely the most humble task).
Below you can see the individual files. I would be very grateful for any help. I think the problem is especially in settings of k and epsilon..
Thank you very much


wall-4 is inlet
wall-5 is outlet.


The folder 0/p:


internalField uniform 0;
boundaryField
{
//- Set patchGroups for constraint patches
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
wall-4
{
type zeroGradient;
}
wall-5
{
type fixedValue;
value uniform 0;
}
wall-3
{
type zeroGradient;
}
wall-6
{
type zeroGradient;
}
wall-7
{
type zeroGradient;
}
}
***************



The folder 0/U:

internalField uniform (0 0 0);
boundaryField
{
//- Set patchGroups for constraint patches
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
wall-4
{
type fixedValue;
value uniform (0 0 1);
}
wall-5
{
type zeroGradient;
}
wall-3
{
type zeroGradient;
}
wall-6
{
type zeroGradient;
}
wall-7
{
type zeroGradient;
}
}
**********************




The folder 0/epsilon:


dimensions [0 2 -3 0 0 0 0];

internalField uniform 8.87;
boundaryField
{
wall-3
{
type zeroGradient;
value uniform 8.87;
}
wall-4
{
type freestream;
freestreamValue uniform 8.87;
}
wall-5
{
type zeroGradient;
}
wall-6
{
type zeroGradient;
value uniform 8.87;
}
wall-7
{
type zeroGradient;
value uniform 8.87;
}
}
*******************



The folder 0/k:


dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.0158;
boundaryField
{
wall-3
{
type zeroGradient;
value uniform 0.0158;
}
wall-4
{
type freestream;
freestreamValue uniform 0.0158;
}
wall-5
{
type zeroGradient;
}
wall-6
{
type zeroGradient;
value uniform 0.0158;
}
wall-7
{
type zeroGradient;
value uniform 0.0158;
}
}



The folder 0/nut:

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0.14;
boundaryField
{
wall-4
{
type calculated;
value $internalField;
}
wall-5
{
type calculated;
value $internalField;
}
wall-3
{
type nutkWallFunction;
value $internalField;
}
wall-6
{
type nutkWallFunction;
value $internalField;
}
wall-7
{
type nutkWallFunction;
value $internalField;
}
}
****************




The folder system/controlDict:

libs
(
"libOpenFOAM.so"
"libincompressibleTurbulenceModel.so"
"libincompressibleRASModels.so"
);
application simpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 150;

deltaT 1;

writeControl timeStep;

writeInterval 10;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable true;
**************************




The folder system/fvSchemes:

ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss limitedLinearV 1;
div(phi,k) bounded Gauss limitedLinear 1;
div(phi,epsilon) bounded Gauss limitedLinear 1;
div(phi,R) bounded Gauss limitedLinear 1;
div(R) Gauss linear;
div(phi,nuTilda) bounded Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}
*********************





The folder system/fvSolution:

solvers
{
p
{
solver GAMG;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
tolerance 1e-06;
relTol 0.05;
}
pFinal
{
$p;
tolerance 1e-06;
relTol 0;
}

"(U|k|epsilon)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-05;
relTol 0.1;
}

"(U|k|epsilon)Final"
{
$U;
tolerance 1e-05;
relTol 0;
}
nuTilda
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.1;
}
}
PIMPLE
{
nOuterCorrectors 4;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p 1e-2;
U 1e-3;
"(k|epsilon)" 1e-3;
}
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
"epsilon.*" 0.7;
}
}
cache
{
grad(U);
}