CFD中文网

酱

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多的液体密度的，是我理解错了吗？到底是怎么算这个表的？

酱

lz解决了吗？我也遇到同样的问题了，高压气流冲击真空环境，使用rhoCentralFoam，算了一段时间之后就出现negative initial temperature t0，改了好久，可以推迟出现的时间，但总会出现，不知道怎么弄了

CFD中文网

酱爆鸡丁

帖子