[tdog]

Hi guys
I hope you can help me. I am running a case in rhoSimplecFoam in OF 2.0.1, I copied the fvSolution from the example case squareBend for the solver, it looks like this:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p0
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.01;
}

p
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}

U0
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

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

U
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}

h0
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

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

h
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}

k0
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

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

k
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}

epsilon0
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
}

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

epsilon
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
rhoMin rhoMin [1 -3 0 0 0] 0.1;
rhoMax rhoMax [1 -3 0 0 0] 1.0;
transonic yes;
}

relaxationFactors
{
p 1;
rho 1; //0.1;
U 0.9;
h 0.95;
k 0.9;
epsilon 0.9;
}

relaxationFactors0
{
p 0.3;
rho 0.1;
U 0.7;
h 0.7;
k 0.7;
epsilon 0.7;
}

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


I am just wondering what effect it has when defining multiple linear solvers for the key parameters? E.g. the velocity field has three different solvers and the pressure has two. Furthermore there's also two sets of relaxation factor, but which one does OF use, both of them?

My results are running very fast with this setup, but they do not seem to converge...

[tdog]

Shoutout to anybody who might know?

[tdog]

Anyone knows?