incompressible thermo type (in buoyantSimpleFoam)
 

Hi all,
I'm trying to compile the buoyantSimpleFoam with basicRhoThermo type to use the <incompressible> prop. to run simulation with water liquid (I saw this set up in chtMultiRegionLiquidFoam).
I tested the chtMultiRegionFoam with <incompressible> thermo type for liquid water and it work well, but it is transient. I want steady. So I tryied the chtMultiRegionSimpleFoam with <incompressible> thermo type and it crashes in hEqn because it divide rho/psi where psi = 0 !
If I use the hEqu of buoyantSimpleFoam it works but the temperature field is not solved !
Why the chtMultiRegionSimpleFoam with <incompressible> thermo type does not work? Should it work?
Did anyone meet this problem?
Thanks for any help, Andrea

Hi,
in the hEqn I changed the term "rho/psi" with "p" (they have same dimensions).
Now it works but I think the T is not well solved.
I attached the same model solved with buoyanPimpleFoam (bottom) and buoyantSimpleFoam modified (top).
Note that for both I'm using basicRhoThermo type <incompressible> and I added porous medium.

Thanks for any suggestion

Andrea

Hi Andrea,

You show here results for a case using buoyantPimpleFoam with the basicRhoThermo type <incompressible>. I also tried the same configuration, but convergence is terrible. There may be many reasons for bad convergence of course, but I am wondering if you also faced any convergence problems here. Have you made any changes to buoyantPimpleFoam?

About the buoyantSimpleFoam, the default thermo type is basicPsiThermo. What kind of changes have to be made to run here basicRhoThermo type <incompressible>?

Thanks for any answers.
Maria

Hi Maria,for pimple I did not find any strange problem. I used incompressible thermo type (check the multiRegionLiquidHeater tutorial).
You are right, for thesimple the default thermo type is the basicPsiThermo type. I changed it in the .C / .H code adding the basicRhoThermo.
But, as you can see from my results, I think this is not enough!
Let me know if you are going on with this,
Andrea

Hello OpenFOAMers,
I am trying to run a multiRegion case using chtMultiRegionFoam solver and using icoPolynomial equationOfState for fluid region and its not working.
I want to use chtMultiRegionFoam with <incompressible> thermo type for fluid region and to use icoPolynomial for temperatur dependent properties. Please let me know how I can do this. I mean where and what modification are required.

Thanks
Asad

Hello Asad,

it should be possible just to use icoPolynomial to set your fluid.

Helli Tobi,
Thanks for your immediate response. I have been using following thermotype for fluid region

and in the start simulation does not show any error but its just stops after some itterations--> courant number shoots up to several hundreds and simulation stops. I also checked my pressue BC they seemed quite reasonable.
I just cant figure out what is wrong with my case...

0 folder of my case...
 

please find attached a 0 folder of my case and please let me know if there is any problem BCs.
Thanks
Asad

I can not find any error in U,p and p_rgh.
Well I wanted to simulate a coat store like that one http://www.thermografie-team.de/news...htspeicher.jpg without success.

My problem was the p and p_rgh initialisation.
After a few modification it was posible to solve my problem with steady state conditions but not with transient solver.

Maybe you find any hints in that thread: http://www.cfd-online.com/Forums/ope...watertank.html

especially with the p and p_rgh initialisation