[sfigato]

Hi Foamers,

I am doing a simulation of an hydraulic orifice by using simpleFoam. I got many problems about convergence. After several trials, I got a good solution under the physical point of view.

Nevertheless, I still have some doubts regarding the convergence.
Here there are my fv Schemes and my fvSolution files:

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
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 PCG;
preconditioner DIC;
tolerance 1e-12;
relTol 0.1;
minIter 1;
//maxIter 100;
}*/

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

/*U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 0.1;
minIter 1;
//maxIter 100;
}*/

U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-12;
relTol 0.1;
nSweeps 1;
minIter 1;
//maxIter 100;
}

/*k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 0.1;
minIter 1;
//maxIter 100;
}*/

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-12;
relTol 0.1;
nSweeps 1;
minIter 1;
//maxIter 100;
}

/*omega
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-12;
relTol 0.1;
minIter 1;
//maxIter 100;
}*/

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-12;
relTol 0.1;
nSweeps 1;
minIter 1;
//maxIter 100;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;

}

potentialFlow
{
nNonOrthogonalCorrectors 10;
}

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

cache
{
grad(U);
}

As attachments there are all the data regarding convergence!


Although all the variables seems to converge (except for Ux), some of them still oscillate. Must I worry about it??

Additionally, I forced the convergence by edit the tolerance, unless all the variables would have decreased until reach very (VERY!!) low values!

Anyone can tell me if my solution converge or not (and if does it reach a steady state)?

Moreover, what exactly means the continuity graph (plotted by using pyFoam)?

Thanks in advance
Best Regards
Marco