[SHY22]

Hi everyone!
I m a newer with the openfoam and I want to simulate the Bartolomej subcooled boiling using reactingTwoPhaseEulerFoam.
There are so many questions with me.
first of all:How to set the thermophysicalProperties？
According to a literature by Qin Hao: set the physical properties of the gas to be saturated, and the physical properties of the fluid change with temperature. I don't know how this should be set up, can you please give me some advice or any information to refer to?

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object thermophysicalProperties.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;//perfectFluid;//
specie specie;
energy sensibleEnthalpy;//sensibleInternalEnergy;
}

mixture
{
specie
{
molWeight 18.0153;
}
equationOfState
{
R 3000;
rho0 868;
rho 868;
//rhoCoeffs<8> ( 1513.93 -1.3688 0 0 0 0 0 0 ); // rho = a + bT

}
thermodynamics
{
Sf 0;
Hf 0;
// CpCoeffs<8> ( 654.38 8.1036 0 0 0 0 0 0 )
Cp 4954.1;
// Cv 3136.7;
Tref 530.59;
Href 1122200;
//Eref 1116500;
}
transport
{

mu 0.00010293;
Pr 0.8377;
}
}


FoamFile
{
version 2.0;
format ascii;
class dictionary;
object thermophysicalProperties.gas;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;// sensibleInternalEnergy;
}

mixture
{
specie
{
molWeight 18.0153;
}
equationOfState
{
R 1695;
rho 22.697;
}
thermodynamics
{
Hf 0;
Cp 4227;
// Cv 2532;
Tref 530.59;
Href 2797900;
//Eref 2599700;
}
transport
{
mu 1.771e-5;
Pr 1.423;
}
}


// ******************************