[guenter]

Hello,

I'm working on chtMultiRegionFoam and simulate the cooling of a fluid region, which is surrounded with a solid region. There are no inlets or outlets. If I use in thermophysicalProperties only constant values my simulation works. But if I want to use temperature dependent values
( I use:
transport polynomial;
thermo hPolynomial;
equationOfState icoPolynomial; )
then my simulation crashes. I found the following thread in this forum:
Introducing polynomial thermophysical properties

In this thread Praveen aksed:

This is quite an old post but I am running into similar issues with convergence in chtMultiRegionFoam with the icoPolynomial EOS model. I do have a closed liquid domain with no inlet or outlet boundaries. Is there a way to deal with such domains to account for property variations?

My question is more or less the same. How can I use temperature dependent values in closed regions and achieve convergence?

I tried to use inletOutlet boundary conditions for U and p_rgh. Then the simulation runs, but the cooling is to fast and unrealistic.

Thank you very much.