equationOfState icoTabulated error

carl_r · October 23, 2020, 15:03

I'm trying to implement temperature-dependent density in chtMultiRegionFoam, in the openFoam foundation version V8.

In the user guide I found:
https://cfd.direct/openfoam/user-gui...36-2730007.1.5

Code:

7.1.5 Equation of state
icoTabulated
Tabulated data for an incompressible fluid using (T, ρ) value pairs, e.g.  
rho ( (200 1010) (400 980) );

I see this is a new feature in V8.

When I implement this in my model I get one of the below errors (depending on whether the regions is fluid or solid thermotype)

I read this that there are no valid thermophysical properties using equationOfState = icoTabulated. Is this correct? Is there a way to implement temperature-dependent density using icoTabulated.

The material is incompressible, it can be fluid or solid depending on the specific model.

Code:

--> FOAM FATAL ERROR:
Unknown solidThermo type
thermoType
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          eConst;
    equationOfState icoTabulated { rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ; };
    specie          specie;
    energy          sensibleInternalEnergy;
}


Valid solidThermo types are:

heSolidThermo  pureMixture  constAnIso   eConst       rhoConst  specie  sensibleInternalEnergy

heSolidThermo  pureMixture  constAnIso   ePolynomial  rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  constAnIso   ePower       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  constIso     eConst       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  constIso     ePolynomial  rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  constIso     ePower       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  exponential  eConst       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  exponential  ePolynomial  rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  exponential  ePower       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  polynomial   eConst       rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  polynomial   ePolynomial  rhoConst  specie  sensibleInternalEnergy
heSolidThermo  pureMixture  polynomial   ePower       rhoConst  specie  sensibleInternalEnergy


    From function static typename Table::iterator Foam::basicThermo::lookupThermo(const Foam::dictionary&, Table*, int, const char**, const Foam::word&) [with Thermo = Foam::solidThermo; Table = Foam::HashTable<Foam::autoPtr<Foam::solidThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>; typename Table::iterator = Foam::HashTable<Foam::autoPtr<Foam::solidThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>::iterator]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/basic/lnInclude/basicThermoTemplates.C at line 46.

FOAM exiting

Code:

--> FOAM FATAL ERROR:
Unknown rhoReactionThermo type
thermoType
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          eConst;
    equationOfState icoTabulated { rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ; };
    specie          specie;
    energy          sensibleInternalEnergy;
}


Valid rhoReactionThermo types are:

heRhoThermo  multiComponentMixture  WLF            eConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy

heRhoThermo  multiComponentMixture  WLF            eConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  WLF            eConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  WLF            hConst       adiabaticPerfectFluid     specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  WLF            hConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  WLF            hConst       rPolynomial               specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  WLF            hConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  WLF            hConst       rhoConst                  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  WLF            hConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          eConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       Boussinesq                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       adiabaticPerfectFluid     specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       perfectGas                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       rPolynomial               specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          hConst       rhoConst                  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          hConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          janaf        Boussinesq                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          janaf        Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          janaf        PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          janaf        PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          janaf        incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          janaf        incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  const          janaf        perfectGas                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  const          janaf        perfectGas                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  logPolynomial  hPolynomial  icoPolynomial             specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  logPolynomial  hPolynomial  icoPolynomial             specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  polynomial     hPolynomial  icoPolynomial             specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  polynomial     hPolynomial  icoPolynomial             specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     eConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     eConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     eConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     eConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     hConst       Boussinesq                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     hConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     hConst       PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     hConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     hConst       incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     hConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     hConst       perfectGas                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     hConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     janaf        Boussinesq                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     janaf        Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     janaf        PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     janaf        PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     janaf        incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     janaf        incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  multiComponentMixture  sutherland     janaf        perfectGas                specie  sensibleEnthalpy
heRhoThermo  multiComponentMixture  sutherland     janaf        perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            eConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            eConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            eConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            hConst       adiabaticPerfectFluid     specie  sensibleEnthalpy
heRhoThermo  pureMixture            WLF            hConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            hConst       rPolynomial               specie  sensibleEnthalpy
heRhoThermo  pureMixture            WLF            hConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  pureMixture            WLF            hConst       rhoConst                  specie  sensibleEnthalpy
heRhoThermo  pureMixture            WLF            hConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          eConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       Boussinesq                specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       adiabaticPerfectFluid     specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       adiabaticPerfectFluid     specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       perfectGas                specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       rPolynomial               specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       rPolynomial               specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          hConst       rhoConst                  specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          hConst       rhoConst                  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          janaf        Boussinesq                specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          janaf        Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          janaf        PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          janaf        PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          janaf        incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          janaf        incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            const          janaf        perfectGas                specie  sensibleEnthalpy
heRhoThermo  pureMixture            const          janaf        perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            logPolynomial  hPolynomial  icoPolynomial             specie  sensibleEnthalpy
heRhoThermo  pureMixture            logPolynomial  hPolynomial  icoPolynomial             specie  sensibleInternalEnergy
heRhoThermo  pureMixture            polynomial     hPolynomial  icoPolynomial             specie  sensibleEnthalpy
heRhoThermo  pureMixture            polynomial     hPolynomial  icoPolynomial             specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     eConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     eConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     eConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     eConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     hConst       Boussinesq                specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     hConst       Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     hConst       PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     hConst       PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     hConst       incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     hConst       incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     hConst       perfectGas                specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     hConst       perfectGas                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     janaf        Boussinesq                specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     janaf        Boussinesq                specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     janaf        PengRobinsonGas           specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     janaf        PengRobinsonGas           specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     janaf        incompressiblePerfectGas  specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     janaf        incompressiblePerfectGas  specie  sensibleInternalEnergy
heRhoThermo  pureMixture            sutherland     janaf        perfectGas                specie  sensibleEnthalpy
heRhoThermo  pureMixture            sutherland     janaf        perfectGas                specie  sensibleInternalEnergy


    From function static typename Table::iterator Foam::basicThermo::lookupThermo(const Foam::dictionary&, Table*, int, const char**, const Foam::word&) [with Thermo = Foam::rhoReactionThermo; Table = Foam::HashTable<Foam::autoPtr<Foam::rhoReactionThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>; typename Table::iterator = Foam::HashTable<Foam::autoPtr<Foam::rhoReactionThermo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::string::hash>::iterator]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/basic/lnInclude/basicThermoTemplates.C at line 46.

FOAM exiting

jherb · October 24, 2020, 19:44

This looks like you have some syntax errors in your thermophysicalProperties file. Look around for examples/tutorials using icoPolynomial for equation of state and start from there.

KKlee · February 8, 2022, 02:21

Hello Carl ,

Have you solved it? I think I am facing the same issue, I am really appreciate if you can help to present the right thermalProerties file.

mikulo · March 14, 2022, 07:05

Hey,

you can use only what is written on the list for heSolidThermo and heRhoThermo, say for e.g.

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

and for solid

thermoType
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo eConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}

your mistake is just simple a wrong combination..

for this... you have combined solid and liquid thermo properties dictionary, which is really wrong.

thermoType
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo eConst;
equationOfState icoTabulated { rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ; };
specie specie;
energy sensibleInternalEnergy;
}

change to,

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

at thermophysicalProperties dictionary..

...

mixture
{

equationOfState
{
rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ;
}

}

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries	NickG	OpenFOAM Installation	3	December 30, 2019 01:21
[blockMesh] blockMesh with double grading.	spwater	OpenFOAM Meshing & Mesh Conversion	92	January 12, 2019 10:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh	gschaider	OpenFOAM Community Contributions	300	October 29, 2014 19:00
OpenFOAM without MPI	kokizzu	OpenFOAM Installation	4	May 26, 2014 10:17
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 18:51

October 24, 2020, 19:44		#2
jherb Senior Member Joachim Herb Join Date: Sep 2010 Posts: 650 Rep Power: 22	This looks like you have some syntax errors in your thermophysicalProperties file. Look around for examples/tutorials using icoPolynomial for equation of state and start from there.

February 8, 2022, 02:21		#3
KKlee New Member lee wei Join Date: Dec 2020 Posts: 11 Rep Power: 6	Hello Carl , Have you solved it? I think I am facing the same issue, I am really appreciate if you can help to present the right thermalProerties file.

March 14, 2022, 07:05		#4
mikulo Member Join Date: Nov 2020 Posts: 53 Rep Power: 6	Hey, you can use only what is written on the list for heSolidThermo and heRhoThermo, say for e.g. thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } and for solid thermoType { type heSolidThermo; mixture pureMixture; transport constIso; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } your mistake is just simple a wrong combination.. for this... you have combined solid and liquid thermo properties dictionary, which is really wrong. thermoType { type heSolidThermo; mixture pureMixture; transport constIso; thermo eConst; equationOfState icoTabulated { rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ; }; specie specie; energy sensibleInternalEnergy; } change to, thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState icoTabulated; specie specie; energy sensibleInternalEnergy; } at thermophysicalProperties dictionary.. ... mixture { equationOfState { rho ( ( 0 1010 ) ( 350 1000 ) ( 800 980 ) ) ; } }