[giammy92]

but my previous are unbelieveble! I have changed my mixing leght, and increasing it until 100, i have achieved the same results for Cl!! What means this?

[giammy92]

My current dat are:

K

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

internalField uniform 0.0025847;

boundaryField
{
TRAILING_EDGE
{
type fixedValue;
value uniform 1e-12;
}
SUCTION_SIDE
{
type fixedValue;
value uniform 1e-12;
}
FAR_FIELD
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.001;
value uniform 0.0025847;
}
PRESSURE_SIDE
{
type fixedValue;
value uniform 1e-12;
}
frontAndBackPlanes
{
type empty;
}
}

NUT

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

internalField uniform 0;

boundaryField
{
TRAILING_EDGE
{
type calculated;
value uniform 0;
}
SUCTION_SIDE
{
type calculated;
value uniform 0;
}
FAR_FIELD
{
type calculated;
value uniform 0;
}
PRESSURE_SIDE
{
type calculated;
value uniform 0;
}
frontAndBackPlanes
{
type empty;
}
}

OMEGA

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

internalField uniform 1769.45;

boundaryField
{
TRAILING_EDGE
{
type omegaWallFunction;
value uniform 1769.45;
}
SUCTION_SIDE
{
type omegaWallFunction;
value uniform 1769.45;
}
FAR_FIELD
{
type turbulentMixingLengthFrequencyInlet;
mixingLength 0.4;
value uniform 1769.45;
}
PRESSURE_SIDE
{
type omegaWallFunction;
value uniform 1769.45;
}
frontAndBack
{
type empty;
}
}

P

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

internalField uniform 0;

boundaryField
{
TRAILING_EDGE
{
type zeroGradient;
}
SUCTION_SIDE
{
type zeroGradient;
}
FAR_FIELD
{
type freestreamPressure;
}
PRESSURE_SIDE
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}
}

U

imensions [0 1 -1 0 0 0 0];

internalField uniform (40.447 9.338 0); // (41.511*cos13 41.511*sin13)

boundaryField
{
TRAILING_EDGE
{
type fixedValue;
value uniform (0 0 0);
}
SUCTION_SIDE
{
type fixedValue;
value uniform (0 0 0);
}
FAR_FIELD
{
type freestream;
freestreamValue uniform (40.447 9.338 0);
}
PRESSURE_SIDE
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBackPlanes
{
type empty;
}
}

[giammy92]

and my now i'm trying with this fvschemes and fvsolution:

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss linear;
div(phi,k) bounded Gauss limitedLinear 1;
div(phi,omega) 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;
}



solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.01;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-08;
relTol 0.1;
}

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}

}

SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;

residualControl
{
p 1e-5;
U 1e-5;
k 1e-5;
omega 1e-5;
}
}

potentialFlow
{
nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.4;
k 0.4;
omega 0.4;
}
}

cache
{
grad(U);
}

[giammy92]

Do you think that for omega's bc at wall i have to use Menter's recommended wall BC, i. e. omegawall=60*nu/(beta*y^2), with nu=kinematic viscosity at the wall, beta=0.075 and y=normal distance between the first fluid node and the nearest wall ? But how i can calculate y and nu at wall? Is the same of nu calculated with mu/rho?
Or if i want use omegawallfunction i have before implement the code with:
omega[faceCellI] = omegaLog;

i refer to this discussion http://www.cfd-online.com/Forums/ope...estigated.html