sonicFoam: Supported energy type is e, thermodynamics package provides h

mm66 · March 3, 2021, 17:58

Dear Foamers,

I have been working with sonicFoam to get the pressure drop of a sonic flow. The default settings uses perfectGas as the equation of state:

Code:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Consequently, the pressure drop values were very low. Therefore, I decided to switch to Peng-Robinson and according to the allowed combinations changed the thermoType to:

Code:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

However, it returns the following error:

Code:

--> FOAM FATAL ERROR: 
Supported energy type is e, thermodynamics package provides h

I tried with sensibleInternalEnergy but it also resulted in the following error:

Code:

unknown psiThermo type

Any ideas? I really appreciate any input here.

Thanks,
MJ

jnanabrao · April 2, 2021, 07:12

It appears the sonicFoam energy equation is written for 'e', which should be internal energy. The rhoPimpleFoam solver's energy equation seems to be written for 'he' which I presume is enthalpy. I am aware that in the PIMPLE dictionary(in fvSolutions) for rhoPimpleFoam, in newer versions of OpenFoam(>6), rhoPimpleFoam can work equivalent to sonicFoam with a 'transonic' flag. Maybe this is what will work for you. Some newer OpenFoam versions have deprecated sonicFoam completely.

March 3, 2021, 17:58	sonicFoam: Supported energy type is e, thermodynamics package provides h	#1
mm66 Member Join Date: Mar 2019 Posts: 81 Rep Power: 7	Dear Foamers, I have been working with sonicFoam to get the pressure drop of a sonic flow. The default settings uses perfectGas as the equation of state: Code: thermoType { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Consequently, the pressure drop values were very low. Therefore, I decided to switch to Peng-Robinson and according to the allowed combinations changed the thermoType to: Code: thermoType { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState PengRobinsonGas; specie specie; energy sensibleEnthalpy; } However, it returns the following error: Code: --> FOAM FATAL ERROR: Supported energy type is e, thermodynamics package provides h I tried with sensibleInternalEnergy but it also resulted in the following error: Code: unknown psiThermo type Any ideas? I really appreciate any input here. Thanks, MJ

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
time step continuity error increases with time_SRFSimplefoam	mostafa kamal	OpenFOAM Running, Solving & CFD	7	October 2, 2019 03:00
Compression instead of expansion	EnricoDeFilippi	OpenFOAM Running, Solving & CFD	1	October 8, 2018 11:19
multiphaseEulerFoam/tank discharge	ves	OpenFOAM Running, Solving & CFD	8	December 31, 2017 15:59
LES supersonic free jet	martyn88	OpenFOAM	22	April 17, 2015 07:00
Pressure instability with rhoSimpleFoam	daniel_mills	OpenFOAM Running, Solving & CFD	44	February 17, 2011 18:08

April 2, 2021, 07:12		#2
jnanabrao Member Jnana Bhaskar Rao Join Date: Mar 2020 Posts: 46 Rep Power: 6	It appears the sonicFoam energy equation is written for 'e', which should be internal energy. The rhoPimpleFoam solver's energy equation seems to be written for 'he' which I presume is enthalpy. I am aware that in the PIMPLE dictionary(in fvSolutions) for rhoPimpleFoam, in newer versions of OpenFoam(>6), rhoPimpleFoam can work equivalent to sonicFoam with a 'transonic' flag. Maybe this is what will work for you. Some newer OpenFoam versions have deprecated sonicFoam completely.