Reference Temperature

[cwcatkins]

Hi,

I am currently doing a write up of thermodynamic models and have a query regarding the model in Mutation++. There appears to be an inconsitency between the code and the Mutation++ write-up in the 2014 paper. In Eq. 12 there is no subtraction of the energy at the reference temperature. However, in the code there seems to be:
LOOP(op(h[i], mp_hform[i] - mp_part_sst[i]))

This also links to the discussion on negative energies as the formulation in the paper would always give a positive energy/enthalpy for N2 and O2. However, running mppequil with air5 at 250K would gives a negative enthalpy/energy.

Please could you clarify what formulation is used?

Thanks,

Chay