[sfigato]

Hi Foamers,

we are simulating a valve with simpleFoam! We have carried out the first simulation by using a first order (upwind) numerical schems for divergence. So, we have switched it with a second order one (Gauss Linear). Here appears many problems with convergence!!

Anyone knows how set properly fvSchemes in order to have a seocnd order accurancy??

Here is our fvSchemes and fvSolutiuon:


ddtSchemes
{
default steadyState;
}

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

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV cellMDLimited Gauss linear 1;
div(phi,k) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,omega) Gauss linearUpwind cellLimited Gauss linear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear limited 0.5;
laplacian((1|A(U)),p) Gauss linear limited 0.5;
laplacian(DkEff,k) Gauss linear limited 0.5;
laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
laplacian(DomegaEff,omega) Gauss linear limited 0.5;
laplacian(DREff,R) Gauss linear limited 0.5;
laplacian(DnuTildaEff,nuTilda) Gauss linear limited 0.5;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}



__________________________________________________ ______-

solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.05;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration true;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;

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

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0.01;
nSweeps 1;
}

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0.01;
nSweeps 1;
}
}

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


relaxationFactors
{
default 0;
p 0.3;
U 0.7;
k 0.7;
omega 0.7;
nut 0.7;
}
}


Thx

Regards
Marco

[sfigato]

Hi Foamers,

no one has an idea about which problem could be???
Pleaseee answer!!

Regards
Marco