[tareqkh]

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:

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

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

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

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....