[palarcon]

Hi all!

I have been asigned the task of performing a basic simulation with the idea of simulating a piston.
To do so, the only information that I received was the dimension of the design domain and the description of the Peng Robinson Model, as can be seen in the following code snippet

mixture	// air at room temperature (293 K)
{
    specie
    {
		nMoles		1;		// reference value for this case
        molWeight   52.042;
    }
	equationOfState
	{
		Tc  		351.26;
		Vc  		0.00235849056603773584905660377358;
		Pc  		5.782e6;
		omega		0.2769;
	}
    thermodynamics
    {
        CpCoeffs<8> ( 979.08 0 0 0 0 0 0 0 );
        Hf          0;
        Sf          0;
    }
    transport
    {
        muCoeffs<8>		( -1.975603342276E-06 5.405129380395E-08 -1.668826016976E-11 0 0 0 0 0 );
        kappaCoeffs<8>	( -8.796126980000E-03 1.098461600000E-04 -2.345618740000E-07 3.75538549E-10 0 0 0 0 );

    }
}

After a short research I came to the conclusion that I can either use the solvers rhoCentralFoam or rhoPimpleFoam (by the way, I'm using OpenFOAM 8) and I decided to start with rhoCentralFoam as a base, considering that the problem was a good start for what I wanted to simulate.

Problem is when I configure the thermoType dictionary I receive the error Unknown psiThermo type. Does anybody has an idea how to solve it?

The thermoType dictionary is as follows:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;	// sensibleEnthalpy - sensibleInternalEnergy 
}

Thank you all in advance