[asalar.14]

Hi everyone
I am not able to find the best fvSchemes and fvSolution for my case. need your advise for best setting



The case is a 3D channel (Length= 1.52, Height= 0.5, width- 0.56) running with interFoam RANS kOmegaSST model having a circular cylinder in vicinity of the upstream boundary.
Mesh non-orthogonality Max: 63.8068 average: 5.53946


I am checking momentum changes in two zones ( 1) the wake zone =4cm after the cylinder and 2) after the wake zone= 12cm after the cylinder)

my optimum mesh size is (0.013 × 0.01 × 0.008)
my deltaT is set to be 0.01 with maxCo=0.8


My fvSchemes is:


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2212 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default cellMDLimited Gauss linear 0.5;
grad(U) cellMDLimited Gauss linear 0.5;
grad(p) Gauss linear;
}

divSchemes
{
div(phi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(rhoPhi,U) Gauss upwind;
div(phi,omega) Gauss linearUpwind default;
div(phi,k) Gauss linearUpwind default;
div(phirb,alpha) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear limited 0.777;;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default limited 0.777;
}
wallDist
{
method meshWave;
}



// ************************************************** *********************** //



and my fvSolution is:


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2212 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
"alpha.water.*"
{
nAlphaCorr 1;
nAlphaSubCycles 3;
cAlpha 1;

MULESCorr yes;
nLimiterIter 3;

solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;
}

p_rgh
{
solver GAMG;
tolerance 1e-04;
relTol 0.01;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 1000;
agglomerator faceAreaPair;
mergeLevels 1;
}

p_rghFinal
{
$p_rgh;
tolerance 1e-06;
relTol 0;
}

"pcorr.*"
{
$p_rghFinal;
tolerance 0.001;
}

"U.*"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
nSweeps 1;
minIter 3;
}

"(k|omega|nuTilda).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-08;
relTol 0;
}
}

PIMPLE
{
momentumPredictor yes;
nCorrectors 3;
nOuterCorrectors 1;
nNonOrthogonalCorrectors 1;
}


relaxationFactors
{
fields
{

"p." 0.7;
}

equations
{
"p.*" 0.7;
"U.*" 0.7;
"k.*" 0.7;
"omega.*" 0.7;
".*" 1;
}
}


// ************************************************** *********************** //





Now what needs to be changed or omitted for better accuracy, faster convergence and stability?