[j-blindi]

Hi,

I've got a problem with rhoSimplecFoam.
I got a case which is quite similar to the squareBend in the tutorials but with symmetryPlanes on each side.

My problem is, that although I copied all BCs, the solver won't solve rho. It says
"rho max/min : 1 1"
in each time step.
But when I run squareBend it solves rho perfectly.

Any ideas where my problem could be?

Thanks a lot.

Best regards
Jason

Edit:
And I can't edit rho in any way, too. If I create a rho-BC the solver ignores it.

[j-blindi]

Hi,

I now know where the problem is.
The temperature in my case is quite low, so the initial density would be >1. But somehow the solver doesn't solve densities above 1.
I tried different temperatures and it all works fine with densities <1 but as soon as I get values above 1, the value is fixed at 1.

Although I know the problem, I don't know the solution. Have you got any ideas? I would be very tankful.

Best
Jason

[j-blindi]

I now switched to rhoPimpleFoam because it solves the density correctly. It's not the solution for the problem but it's ok for now....

[cosimobianchini]

Maybe you forgot changing the rhoMax entry in the SIMPLE subDict of the fvSolution dictionary

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

If you will put reasonable values for your case I do not doubt that the solution will show the actual density values.

[j-blindi]

Thank you very very much, I think that's the solution .

Best
Jason

PS: rhoPimpleFoam was not the right solver....