[Gubevens]

Hey Guys,

I`m running a simulation with rhoSimpleFoam solver, LRR turbulence model and a mesh with 1.8 M tetraedrons.

During the first 100 Iterations my epsilons residuals are constant (like 0.3). After that it rises to one and falls again. It does not blow up.
See the attached picture.
In my eyes this is not ok, so my question is how can i fix that to get a convergence? Or is this behaviour not considerable?

info: the same case does converge with kEpsilon model after 4500 Iterations. CheckMesh Ok.

Thank you in advance

----------------------------------------------------------------------
Here is the fvSchemes Setup
----------------------------------------------------------------------
ddtSchemes
{
default steadyState;
}

gradSchemes
{
default cellLimited Gauss linear 1.0;
}

divSchemes
{
default bounded gauss upwind;;

div(phi,U) bounded Gauss upwind;
div(((rho*(thermo:mu|(p*thermosi)))*dev2(T(grad( U))))) Gauss linear;
div(phi,K) bounded Gauss linear;
div(phi,h) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,Ekp) bounded Gauss upwind;
div((rho*R)) Gauss linear;
div(phi,R) bounded Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}