关于multiphaseEulerFoam中的物性问题
-
OpenFOAM v9自带算例wall boiling里面,管内沸腾,介质应该是R12,他里面的物性是通过tabulated的方式定义的,
thermoType { type heRhoThermo; mixture pureMixture; transport tabulated; thermo hTabulated; equationOfState rhoTabulated; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 120.914; } #include "$FOAM_TUTORIALS/resources/thermoData/wallBoiling-liquid" }
打开这个wallBoiling-liquid文件,
equationOfState { rho { low ( 2200000 325 ); high ( 3200000 425 ); values 26 26 ( ( 1214.612 1196.5038 1177.4967 1157.4193 1136.0403 1113.0341 1087.9176 1050.20560434 1002.48588455 954.766164748 907.04644495 859.326725151 811.607005353 763.887285555 716.167565757 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 ) ( 1215.0065 1196.9395 1177.9822 1157.9662 1136.6646 1113.7597 1088.7819 1055.45565866 1007.07208321 958.688507754 910.304932301 861.921356848 813.537781395 765.154205943 716.77063049 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 672.81766 ) ...
第一行数据应该是2.2MPa下,从325K 到425K 等分26份,算密度。
我用NIST查了一下R12的物性,2.2MPa,325K
的时候的确是1214.6,但是350.94K的时候就饱和了,密度是1074.7,我不知道wallBoiling-liquid这个表里面怎么会有600多的液体密度的,是我理解错了吗?到底是怎么算这个表的?