CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   rhoSimplecFoam, solving rho (https://www.cfd-online.com/Forums/openfoam-solving/93271-rhosimplecfoam-solving-rho.html)

 j-blindi October 10, 2011 14:13

rhoSimplecFoam, solving rho

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 October 13, 2011 04:46

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 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

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

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

Best
Jason

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

 All times are GMT -4. The time now is 20:33.