[rezvani]

hello. I want to use the Tait equation of state in compressibleInterFoam. in some documents, they said this is New compressibility models in OpenFOAM® but in help or tutorial, I can't see any guide to use it.
Thank You

[cbarring]

The Tait equation is the adiabaticPerfectFluid EOS. I am using the code below within the thermophysicalProperties.water file in compressibleInterFoam in OF7 (In OF8 & 9 I get a bunch of errors that I am trying to figure out).

I did a curve fit to isentropic liquid water data to get the gamma and B values. These are relevant for my particular project. In the open literature you will find a lot of different values for gamma and B. You will need to figure out the values that fit the fluid properties for your particular situation. This will require some (non-linear) curve fitting. This is really important or you may get very strange numbers. The values below for gamma and B may be totally wrong for your situation.

thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo eConst;
equationOfState adiabaticPerfectFluid;
specie specie;
energy sensibleInternalEnergy;
}

mixture
{
specie
{
nMoles 1;
molWeight 18.01528;
}

equationOfState
{
// water
rho0 1000; // approx
p0 100000;
gamma 7.15;
B 304913000;
}

thermodynamics
{
Cv 4186;
Hf 0;
}
transport
{
mu 3.645e-4;
Pr 6.99;
}
}

[Beatcfd]

Originally Posted by cbarring

The Tait equation is the adiabaticPerfectFluid EOS. I am using the code below within the thermophysicalProperties.water file in compressibleInterFoam in OF7 (In OF8 & 9 I get a bunch of errors that I am trying to figure out).

I did a curve fit to isentropic liquid water data to get the gamma and B values. These are relevant for my particular project. In the open literature you will find a lot of different values for gamma and B. You will need to figure out the values that fit the fluid properties for your particular situation. This will require some (non-linear) curve fitting. This is really important or you may get very strange numbers. The values below for gamma and B may be totally wrong for your situation.

thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo eConst;
equationOfState adiabaticPerfectFluid;
specie specie;
energy sensibleInternalEnergy;
}

mixture
{
specie
{
nMoles 1;
molWeight 18.01528;
}

equationOfState
{
// water
rho0 1000; // approx
p0 100000;
gamma 7.15;
B 304913000;
}

thermodynamics
{
Cv 4186;
Hf 0;
}
transport
{
mu 3.645e-4;
Pr 6.99;
}
}

While I'm applying your set to my case, it comes error that：“request for volScalarField thermo: psi from objectRegistry region0 failed”.

Do you know what's wrong or maybe my settings wrong?

[Beatcfd]

Originally Posted by Beatcfd

While I'm applying your set to my case, it comes error that：“request for volScalarField thermo: psi from objectRegistry region0 failed”.

Do you know what's wrong or maybe my settings wrong?

Or if possible，could you upload your "0" case files, I'm wondering if these initial files' setting may be somewhere wrong.
Thank you！！！