[bastil]

Dear Solver-Experts,

I am running a case in rhoSimpleFoam and I get exploding values for temperature and especially density after few iterations, but no divergence. This even happens for turbulence turned off. I have to blame a few cells in my prism layer (low-Re, layer next to wall) for this.
I am trying to bound rho and h in some way, but without success so far. I tried:

interpolationSchemes
{
      interpolate(rho) limitedLimitedLinear phi 1.0 1.0 1.25;

...
}

to bound rho between 1.0 and 1.25 but this does not work at all. I get higer densities after 4 iterations... How can these values be bounded? Thanks.

Regards