specific heat for mixtures

[alexDfluid]

Hello guys,

I would like to know if it is possible to calculate the specific heats c_v and c_p for mixtures in refprop using the python interface?
For example:
Propane - Nitrogen mixture
mixture ratio z [0.5,0.5]
temperature T
pressure P

I use TPFLSH which gives reasonable results for the single phase region.
But the documentation mentions, that cp is not defined for 2-phase states and Cv for 2-phase states is not calculated.
As far as I've seen Cv2phase is not applicable to mixtures, but only for pure substances?

Any possibility to get the specific heats for mixtures?

best,
Alex

[ianhbell]

c_p is not defined for two-phase states for pure fluids. Thinking a little bit about your question, I can agree that one should be able to define a c_p for a two-phase mixture because dh/dT|P is not ill-defined for a two-phase VLE mixture. But REFPROP does not currently support such a calculation. How you could do this is to take your inputs from REFPROP and plug them into CoolProp, which has the first_partial_deriv function that should in theory give you what you are asking for (I think).

[ianhbell]

Now that I say that, there is also this function: https://refprop-docs.readthedocs.io/en/latest/DLL/legacy.html#f/_/DHD1dll

[ianhbell]

So you would calculate the density as usual, and then plug that density into DHD1dll function

[alexDfluid]

Thank you for the very fast reply.
I tried to use DHD1dll in the two phase region where 0 <= q <=1 for the Propane/Nitrogen system at 280K and 50bar.
Unfortunately, there seem to be a discontinuity.

Any further suggestions how to get c_v and c_p in a two phase mixture?

[ianhbell]

Hmm, I am not entirely sure the very idea of a two-phase c_p is well-posed for a mixture.

[nist-aharvey]

Cp and Cv are both defined for a mixture in the 2-phase region as just their normal definitions, taken with constant overall composition. So for example Cp is (dh/dT)_p. [I think 2-phase Cp may become undefined if you are at an azeotrope, but otherwise should be a well-defined quantity.]
However, I don't think there is any straightforward way to calculate it; because vaporization and composition changes in the phases happen when heat is added it is not anything simple like a combination of (dh/dT)_p of the individual phases.
So my only suggestion would be to evaluate the derivative as a finite difference, being careful that you keep the appropriate variables constant when you increment the temperature.
To be safe, I would also apply the flag to set the reference state to "each pure component" since sometimes mixture reference states can cause unexpected results when compositions are changing.

[alexDfluid]

Thank you for your response.
I was hoping for a more simple solution for the calculation of the specific heats, but it seems that I will have to apply your suggestion.