[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

[olesen]

Originally Posted by bastil

Dear Solver-Experts,

I am running a case in rhoSimpleFoam and I get exploding values for temperature and especially density after few iterations

Have you tried severely underrelaxing rho? For example, a value of 0.05 might get you over the rough spots.

[bastil]

Originally Posted by olesen

Have you tried severely underrelaxing rho? For example, a value of 0.05 might get you over the rough spots.

Thanks Mark. Yes I do this. I still have problems in some critical areas. I am also wondering about BCs at walls. We use for walls:
fixedValue, 1e-20 for k and epsilon (for LaunderSharma-k-epsilon)
zeroGradient for temperature
I am not sure for temperature... The only similar tutorial case is NacaAirfoil from turbFoam. This uses zeroGradient for temperature but also for k and epsilon which is obviously wrong?
Any other ideas? Thanks.

Regards BastiL