[Tobi]

Hello all,

I am working on a coffee - cht - case. Therefor I am using three regions:

- air
- coffee (water)
- cup (solid/ceramic)

Meshing with sHM is not a problem but after initializing (constructors of the code) I will get a negative Co-Number of my water region.

To simulate water I first used the thermodynamicProperties of the chtMultiRegionFoam tutorial - liquid heater. Well I don't want have a constant density so I decided to use the polynomial function.


With the tutorial properties I get no a negative number
With the polynomial functions I get negative co - numbers


Therefor I think that my polynomial coefficients are wrong because Co is calculated with U, dt, dx and mu. The first three variables are the same only the mu value is calculated with the polynomial functions. I checked the functions with gnu plot and open office and there is everything okay.

Additionally I checked the sign of the coefficient values but there is no error too.

But it should be the reason, isn ´t it?
I will try to use just one coefficient - the first one. This should be the same as constant values then, isn´t it?

Does someone know that problem?
I think I am on the right way and will you keep posted. But the strange thing is that gnu plot and open office is calculating everything correct.

The only thing that could be is that I get unphysical values for the temperature so that the polynomial function is not guilty anymore.

Regards
Tobi

[Tobi]

Hi all,

density polynomial was wrong. I missed a decade.

0.05 --> 0.005

Regards
Tobi