[alextau]

Hello,

I am using compressibleInterFoam solver and trying to apply peng robinson eqaution of state on methane gas.

The followng code in thermoPhysicalProperties.methane:

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

mixture
{
specie
{
molWeight 16.04;
}
thermodynamics
{
Cp 3300;
Hf 0;
}
transport
{
mu 1.39e-05;
Pr 0.8;
}

equationOfState
{
Tc 190.4;
Vc 99.2E-03;
Pc 4.6E+6;
Zc 0.85;
omega 0.011;
}
}


and the output is :


31]
[31]
[31] file: IOstream.thermoType at line 0.
[31]
[31] From static typename Table::iterator Foam::basicThermo::lookupThermo(const Foam::dictionary&, Table*, std::initializer_list<const char*>, const Foam::word&) [with Thermo = Foam::rhoThermo; Table = Foam::HashTable<Foam::autoPtr<Foam::rhoThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>; typename Table::iterator = Foam::HashTable<Foam::autoPtr<Foam::rhoThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>::iterator]
[31] in file lnInclude/basicThermoTemplates.C at line 81.
[31]
FOAM parallel run exiting
[31]
ibleInternalEnergy
heRhoThermo reactingMixture const eConst rhoConst specie sensibleInternalEnergy
heRhoThermo reactingMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst perfectFluid specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst rhoConst specie sensibleEnthalpy
heRhoThermo reactingMixture polynomial hPolynomial icoPolynomial specie sensibleEnthalpy
heRhoThermo reactingMixture polynomial hPolynomial icoPolynomial specie sensibleInternalEnergy
heRhoThermo reactingMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture sutherland janaf incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo reactingMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo singleStepReactingMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo singleStepReactingMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo veryInhomogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy
[27]
[27]
[27] file: IOstream.thermoType at line 0.
[27]
[27] From static typename Table::iterator Foam::basicThermo::lookupThermo(const Foam::dictionary&, Table*, std::initializer_list<const char*>, const Foam::word&) [with Thermo = Foam::rhoThermo; Table = Foam::HashTable<Foam::autoPtr<Foam::rhoThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>; typename Table::iterator = Foam::HashTable<Foam::autoPtr<Foam::rhoThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>::iterator]
[27] in file lnInclude/basicThermoTemplates.C at line 81.
[27]
FOAM parallel run exiting




Well, it seems that PengRobinson is not even an option here. I checked in src that it is there. What am I doing wrong?

Please help,
Alex

[dasith0001]

Hi,

What OF version are you running ?

Anything above OF 8 should work, I run on OF10 and the dynamicCode set it up for you and it works fine.

Cheers,
Dasith