CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   rhoSimplecFoam, solving rho (

j-blindi October 10, 2011 14:13

rhoSimplecFoam, solving rho

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

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

j-blindi October 13, 2011 04:46


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.


j-blindi October 13, 2011 05:57

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 October 14, 2011 03:01

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

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 October 14, 2011 04:54

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


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

All times are GMT -4. The time now is 13:12.