CFD Online Discussion Forums - differences in OF 2.2

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - differences in OF 2.2 (https://www.cfd-online.com/Forums/openfoam-solving/115108-differences-2-2-a.html)

differences in OF 2.2

i want to run a case i made in OF 2.1.1 into OF2.2 but it doesn't know thermophysical model.

It might be helpful to post the name of the model you want to use and the error message you get in OF 2.2 so we at least can check what is wrong. I'm guessing the model has gotten a different name in OF 2.2 for some reason.

Greetings to all!

OpenFOAM 2.2 is now easier to use! ;) Check the section "Dictionary-based Thermodynamics Selection" on the release notes: http://www.openfoam.org/version2.2.0....php#x10-40000

Best regards,
Bruno

this is the error:

Code:

PIMPLE: no residual control data found. Calculations will employ 7 corrector loops



Reading thermophysical properties



Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>

[0] [1]

[1]

[1] --> FOAM FATAL ERROR:

[1] Unknown psiThermo type hPsiThermo<pureMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>



Valid psiThermo types are:



34

(

hePsiThermo<homogeneousMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<homogeneousMixture<sutherland<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<homogeneousMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<inhomogeneousMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<inhomogeneousMixture<sutherland<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<inhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<multiComponentMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<multiComponentMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<pureMixture<const<eConst<perfectGas<specie>>,sensibleInternalEnergy>>>

hePsiThermo<pureMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<pureMixture<const<hConst<perfectGas<specie>>,sensibleInternalEnergy>>>

hePsiThermo<pureMixture<sutherland<eConst<perfectGas<specie>>,sensibleInternalEnergy>>>

hePsiThermo<pureMixture<sutherland<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<pureMixture<sutherland<hConst<perfectGas<specie>>,sensibleInternalEnergy>>>

hePsiThermo<pureMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<pureMixture<sutherland<janaf<perfectGas<specie>>,sensibleInternalEnergy>>>

hePsiThermo<reactingMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<reactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<singleStepReactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<veryInhomogeneousMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<veryInhomogeneousMixture<sutherland<hConst<perfectGas<specie>>,sensibleEnthalpy>>>

hePsiThermo<veryInhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>

heheuPsiThermo<egrMixture<const<hConst<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<egrMixture<sutherland<janaf<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<egrMixture<sutherland<janaf<perfectGas<specie>>,absoluteInternalEnergy>>>

heheuPsiThermo<homogeneousMixture<const<hConst<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<homogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<homogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteInternalEnergy>>>

heheuPsiThermo<inhomogeneousMixture<const<hConst<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<inhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<inhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteInternalEnergy>>>

heheuPsiThermo<veryInhomogeneousMixture<const<hConst<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<veryInhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteEnthalpy>>>

heheuPsiThermo<veryInhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,absoluteInternalEnergy>>>

)



[1]

[1]

[1]     From function psiThermo::New

[1]     in file lnInclude/basicThermoTemplates.C at line 123.

[1]

FOAM parallel run exiting

[1]

whats the difference between sensibleEnthalpy and sensibleInternalEnergy?
is former more suitable for gases and later for liquids?is there any more?

Hi Ehsan,

According to the Doxygen documentation:

sensibleEnthalpy: http://foam.sourceforge.net/docs/cpp/a01976.html

Quote:

Thermodynamics mapping class to expose the sensible enthalpy function as the standard enthalpy function h(T).
sensibleInternalEnergy: - http://foam.sourceforge.net/docs/cpp/a01977.html

Quote:

Thermodynamics mapping class to expose the sensible internal energy function as the standard internal energy function e(T).

Comparing the tutorial "compressible/rhoPimpleFoam/ras/angledDuct/constant/thermophysicalProperties" on each version:

OpenFOAM 2.1.x: https://github.com/OpenFOAM/OpenFOAM...icalProperties
OpenFOAM 2.2.x: https://github.com/OpenFOAM/OpenFOAM...icalProperties

Interesting... it looks like the selection of "h" or "e" depends on a run time option: http://www.openfoam.org/version2.2.0/thermophysical.php

Strange... something doesn't add up... I'll keep looking at the code and let you know if I can figure how exactly "h" or "e" are chosen!?

edit: :eek: Now I get it... it's chosen in the option "energy":

Quote:

Originally Posted by http://www.openfoam.org/version2.2.0/thermophysical.php

The syntax include the energy keyword, in which the user specifies the form of energy to be used in the solution, e.g. sensibleEnthalpy, sensibleInternalEnergy, absoluteEnthalpy.

At first, it looked like "h" and "e" were selected directly from the command line... but it's this option "energy" that helps choose "h" or "e".

Best regards,
Bruno

thanks dear Bruno.at which problems we use absoluteEnthalpy and in which sensibleEnthalpy?

Hi Ehsan,

I've checked the tutorials for ideas... the following command gathers information about the tutorials that are using "absoluteEnthalpy":

Code:

grep -R 'absoluteEnthalpy' $FOAM_TUTORIALS

Only the following ones appeared:

Code:

combustion/engineFoam/kivaTest/constant/thermophysicalProperties:    energy          absoluteEnthalpy;

combustion/XiFoam/ras/moriyoshiHomogeneous/constant/thermophysicalProperties:    energy          absoluteEnthalpy;

combustion/PDRFoam/flamePropagationWithObstacles/constant/thermophysicalProperties:    energy          absoluteEnthalpy;

After checking all other tutorials that use "sensibleEnthalpy", it seems that "absoluteEnthalpy" should be used along with the type "heheuPsiThermo", which seems to be specific to the "reactionThermo" library.

I'm not familiar with the specific details for each, but the tutorials seem to indicate that "sensibleEnthalpy" is the default option...
Although for this logic, it would mean that "absoluteInternalEnergy" is rarely used... there is only one reference to it in the tutorial "combustion/engineFoam/kivaTest".

I can only guess that the specific reasons depends on the chemical formulations behind it all.

Best regards,
Bruno

Dear wyldchat:

I want use the thermo type as the following:

Code:

 thermoType

{

type hePsiThermo;

mixture reactingMixture;

transport sutherland;

thermo janaf;

energy absoluteEnthalpy;

equationOfState perfectGas;

specie specie;

}

but class hePsiThermo only have "sensibleEnthalpy" insteadingof "absoluteEnthalpy".In 2.2.1,only when I use heheuPsiThermo,can I use "absoluteEnthalpy".I think my defined termotype is not available in2.2.1.I According to words "if the particular combination of models is not available in the list, that combination can be compiled into a library or the application (as always). ", I have to according to class heheuPsiThermo to add "absoluteEnthalpy" in class hePsiThermo .However,when I run my case,the error showed me that there no change with "Valid psiThermo types ".In other word ,my modifies did not work.the file I modified is "thermophysicalModels/reactionThermo/psiReactionThermos.C"I just add the following code:

Code:

makeReactionThermo

(

PsiThermo,

PsiReactionThermo,

hePsiThermo,

reactingMixture,

sutherland,

absoluteEnthalpy,

janafThermo,

perfectGas,

specie;

)

and #include "absoluteEntalpy.H" to psiReactionThermos.C.As you know ,it did not work.I indeed have some troubles in how to change codes to combinate new thermo type.could you tell me how to modify psiReactionThermos.C to make my difined thermo type :

Code:

{

type hePsiThermo;

mixture reactingMixture;

transport sutherland;

thermo janaf;

energy absoluteEnthalpy;

equationOfState perfectGas;

specie specie;

}

to be useful?
Or can I build myThermo library that include thermotype what I want to use?Please someone could give some advise!

Please!here is my new thermotype library myThermo.zip

Greetings zqlhzx,

Sorry, unfortunately I haven't had the time to look into this. And I saw just now that you've asked about this on several other threads and have gotten a few answers, so I'll compile the list of threads you've gotten answers so far on your latest thread: http://www.cfd-online.com/Forums/ope...in2-2-1-a.html

Best regards,
Bruno

Dear wyldckat:
Thank you for your reply!If you have some times,please take a look on my problem!Thanks in advance!

hi
I want to use homogeneousmixture in engineFoam > kivaTest but it doesn't work, can you guide me please?