[alanwelisson]

Hello everyone,

I am trying to run a case with a custom solver based on reactingFoam and I need to use the Peng–Robinson equation of state. This is my thermoType dictionary:

thermoType
{
    type            hePsiThermo;
    mixture         multiComponentMixture;
    transport       const;
    thermo          hConst;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Every time I try to run the case, I get this error:

--> FOAM FATAL IO ERROR: (openfoam-2212)
Unknown psiReactionThermo type 

Valid psiReactionThermo types :

23
(
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<eConst<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<multiComponentMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<multiComponentMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<multiComponentMixture<sutherland<janaf<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<pureMixture<const<eConst<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<pureMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<pureMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<pureMixture<sutherland<janaf<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<reactingMixture<const<eConst<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<reactingMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<reactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<reactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<singleStepReactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<singleStepReactingMixture<sutherland<janaf<perfectGas<specie>>,sensibleInternalEnergy>>>
hePsiThermo<veryInhomogeneousMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<veryInhomogeneousMixture<sutherland<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
hePsiThermo<veryInhomogeneousMixture<sutherland<janaf<perfectGas<specie>>,sensibleEnthalpy>>>
)

type         mixture                    transport   thermo  equationOfState  specie  energy                  

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       eConst  perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  multiComponentMixture      const       hConst  perfectGas       specie  sensibleEnthalpy        
hePsiThermo  multiComponentMixture      sutherland  janaf   perfectGas       specie  sensibleEnthalpy        
hePsiThermo  multiComponentMixture      sutherland  janaf   perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  pureMixture                const       eConst  perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  pureMixture                const       hConst  perfectGas       specie  sensibleEnthalpy        
hePsiThermo  pureMixture                sutherland  janaf   perfectGas       specie  sensibleEnthalpy        
hePsiThermo  pureMixture                sutherland  janaf   perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  reactingMixture            const       eConst  perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  reactingMixture            const       hConst  perfectGas       specie  sensibleEnthalpy        
hePsiThermo  reactingMixture            sutherland  janaf   perfectGas       specie  sensibleEnthalpy        
hePsiThermo  reactingMixture            sutherland  janaf   perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  singleStepReactingMixture  sutherland  janaf   perfectGas       specie  sensibleEnthalpy        
hePsiThermo  singleStepReactingMixture  sutherland  janaf   perfectGas       specie  sensibleInternalEnergy  
hePsiThermo  veryInhomogeneousMixture   const       hConst  perfectGas       specie  sensibleEnthalpy        
hePsiThermo  veryInhomogeneousMixture   sutherland  hConst  perfectGas       specie  sensibleEnthalpy        
hePsiThermo  veryInhomogeneousMixture   sutherland  janaf   perfectGas       specie  sensibleEnthalpy        


file: constant/thermophysicalProperties.thermoType at line 21 to 27.

From what I understand, the available combinations in my build only show perfectGas variants. My questions is: Is it possible to use PengRobinsonGas with a reactingFoam-based solver (i.e., hePsiThermo / psiReactionThermo), or is that EOS only supported for certain top-level thermo types?