CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Covergence Problem with simpleFoam (http://www.cfd-online.com/Forums/openfoam-solving/112743-covergence-problem-simplefoam.html)

sfigato February 4, 2013 12:11

Covergence Problem with simpleFoam
 
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 February 6, 2013 08:31

Hi Foamers,

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

Regards
Marco


All times are GMT -4. The time now is 01:09.