Pure and Pseudo-Pure fluid properties

Introduction

Nearly all the fluids modeling in CoolProp are based on Helmholtz energy formulations. This is a convenient construction of the equation of state because all the thermodynamic properties of interest can be obtained directly from partial derivatives of the Helmholtz energy.

It should be noted that the EOS are typically valid over the entire range of the fluid, from subcooled liquid to superheated vapor, to supercritical fluid.

Annoyingly, different authors have selected different sets of nomenclature for the Helmholtz energy. For consistency, the nomenclature of Lemmon will be used here. Also, some authors present results on a mole-basis or mass-basis, further complicating comparisons.

Thermodynamic properties of Fluid

In general, the EOS are based on non-dimensional terms \(\delta\) and \(\tau\), where these terms are defined by

\[ \begin{align}\begin{aligned}\delta=\rho/\rho_c\\\tau=T_c/T\end{aligned}\end{align} \]

where \(\rho_c\) and \(T_c\) are the critical density of the fluid if it is a pure fluid. For pseudo-pure mixtures, the critical point is typically not used as the reducing state point, and often the maximum condensing temperature on the saturation curve is used instead.

The non-dimensional Helmholtz energy of the fluid is given by

\[\alpha=\alpha^0+\alpha^r\]

where \(\alpha^0\) is the ideal-gas contribution to the Helmholtz energy, and \(\alpha^r\) is the residual Helmholtz energy contribution which accounts for non-ideal behavior. For a given set of \(\delta\) and \(\tau\), each of the terms \(\alpha^0\) and \(\alpha^r\) are known. The exact form of the Helmholtz energy terms is fluid dependent, but a relatively simple example is that of Nitrogen, which has the ideal-gas Helmholtz energy of

\[\alpha^0=\ln\delta+a_1\ln\tau+a_2+a_3\tau+a_4\tau^{-1}+a_5\tau^{-2}+a_6\tau^{-3}+a_7\ln[1-\exp(-a_8\tau)]\]

and the non-dimensional residual Helmholtz energy of

\[\alpha^r=\sum_{k=1}^{6}{N_k\delta^{i_k}\tau^{j_k}}+\sum_{k=7}^{32}{N_k\delta^{i_k}\tau^{j_k}\exp(-\delta^{l_k})}+\sum_{k=33}^{36}{N_k\delta^{i_k}\tau^{j_k}\exp(-\phi_k(\delta-1)^2-\beta_k(\tau-\gamma_k)^2)}\]

and all the terms other than \(\delta\) and \(\tau\) are fluid-dependent correlation parameters.

The other thermodynamic parameters can then be obtained through analytic derivatives of the Helmholtz energy terms. For instance, the pressure is given by

\[p=\rho RT\left[1+\delta\left(\frac{\partial \alpha^r}{\partial \delta}\right)_{\tau} \right]\]

and the specific internal energy by

\[\frac{u}{RT}=\tau \left[\left(\frac{\partial \alpha^0}{\partial \tau}\right)_{\delta}+ \left(\frac{\partial \alpha^r}{\partial \tau}\right)_{\delta} \right]\]

and the specific enthalpy by

\[\frac{h}{RT}=\tau \left[\left(\frac{\partial \alpha^0}{\partial \tau}\right)_{\delta}+ \left(\frac{\partial \alpha^r}{\partial \tau}\right)_{\delta} \right] +\delta\left(\frac{\partial \alpha^r}{\partial \delta}\right)_{\tau}+1\]

which can also be written as

\[\frac{h}{RT}=\frac{u}{RT}+\frac{p}{\rho RT}\]

The specific entropy is given by

\[\frac{s}{R}=\tau \left[\left(\frac{\partial \alpha^0}{\partial \tau}\right)_{\delta}+ \left(\frac{\partial \alpha^r}{\partial \tau}\right)_{\delta} \right]-\alpha^0-\alpha^r\]

and the specific heats at constant volume and constant pressure respectively are given by

\[ \begin{align}\begin{aligned}\frac{c_v}{R}=-\tau^2 \left[\left(\frac{\partial^2 \alpha^0}{\partial \tau^2}\right)_{\delta}+ \left(\frac{\partial^2 \alpha^r}{\partial \tau^2}\right)_{\delta} \right]\\\frac{c_p}{R}=\frac{c_v}{R}+\dfrac{\left[1+\delta\left(\frac{\partial \alpha^r}{\partial \delta}\right)_{\tau}-\delta\tau\left(\frac{\partial^2 \alpha^r}{\partial \delta\partial\tau}\right)\right]^2}{\left[1+2\delta\left(\frac{\partial \alpha^r}{\partial \delta}\right)_{\tau}+\delta^2\left(\frac{\partial^2 \alpha^r}{\partial \delta^2}\right)_{\tau}\right]}\end{aligned}\end{align} \]

The EOS is set up with temperature and density as the two independent properties, but often other inputs are known, most often temperature and pressure because they can be directly measured. As a result, if the density is desired for a known temperature and pressure, it can be obtained iteratively. The following algorithm is used to obtain a reasonable guess for the initial value for the iterative solver:

1. If the fluid is superheated, use a guess of ideal gas (\(\rho=p/(RT)\))

2. If the fluid is subcooled, use a guess of saturated liquid density

3. If the fluid is supercritical, use a guess of ideal gas (\(\rho=p/(RT)\))

4. No solution for density as a function of temperature and pressure if the fluid is two-phase

The documentation of the CoolProp.CoolProp module, or the CoolProp.State module are also available.

List of Fluids

Note

You can click on the fluid name to get more information about the fluid, or click on a bracketed reference to be taken to the reference for the fluid

All the fluids included in CoolProp

Name	EOS	\(c_{p0}\)	\(\lambda\)	\(\eta\)	melt	\(\sigma\)
1-Butene	[1]					[2]
Acetone	[3]					[2]
Air	[4]		[5]	[5]
Ammonia	[6]		[7]	[8]		[2]
Argon	[9]		[5]	[5]	[9]	[2]
Benzene	[10]		[11]	[12]		[2]
CarbonDioxide	[13]		[14]	[15]	[13]	[2]
CarbonMonoxide	[3]				[16]	[2]
CarbonylSulfide	[3]					[2]
CycloHexane	[17]			[18]	[19]	[2]
CycloPropane	[20]	[20]				[21]
Cyclopentane	[22]		[23]	[24]		[21]
D4	[25]					[21]
D5	[26]					[21]
D6	[27]					[21]
Deuterium	[28]					[2]
Dichloroethane	[25]
DiethylEther	[29]					[21]
DimethylCarbonate	[30]					[21]
DimethylEther	[31]			[32]		[2]
Ethane	[33]		[34]	[34]	[33]	[2]
Ethanol	[35]		[36]	[37]	[38]	[2]
EthylBenzene	[39]		[40]			[21]
Ethylene	[41]				[41]	[2]
EthyleneOxide	[25, 42, 43]
Fluorine	[44]	[44]			[44]	[2]
HFE143m	[45]
HeavyWater	[46]		[47]	[47]		[48]
Helium	[49]		[50]	[51]	[52]	[2]
Hydrogen	[53]		[54]	[55]	[52]	[2]
HydrogenChloride	[56]
HydrogenSulfide	[3]			[57]		[2]
IsoButane	[58]		[59]	[60]	[58]	[2]
IsoButene	[1]					[2]
Isohexane	[3]					[2]
Isopentane	[3]		[23]	[24]	[61]	[2]
Krypton	[3]				[62]	[2]
MD2M	[63]					[21]
MD3M	[26]					[21]
MD4M	[26]					[21]
MDM	[63]					[21]
MM	[25, 64]					[21]
Methane	[65]		[66]	[67]	[68]	[2]
Methanol	[69]		[70]	[71]	[72]	[2]
MethylLinoleate	[73]					[21]
MethylLinolenate	[73]
MethylOleate	[73]					[21]
MethylPalmitate	[73]					[21]
MethylStearate	[73]					[21]
Neon	[74]				[75]	[2]
Neopentane	[3]					[21]
Nitrogen	[76]		[5]	[5]	[76]	[2]
NitrousOxide	[3]					[2]
Novec649	[77]
OrthoDeuterium	[28]
OrthoHydrogen	[53]
Oxygen	[78, 79]		[5]	[5]	[80]	[2]
ParaDeuterium	[28]
ParaHydrogen	[53]		[54]		[80]	[2]
Propylene	[81]		[82]	[82]	[61]	[2]
Propyne	[20]	[20]				[2]
R11	[83]	[83]	[84]	[85]		[2]
R113	[86]					[2]
R114	[87]	[87]				[2]
R115	[88]
R116	[3]		[82]	[82]		[2]
R12	[86]		[84]	[85]		[2]
R123	[89]	[89]	[90]	[91]		[2]
R1233zd(E)	[92]			[93]		[94]
R1234yf	[95]		[96]	[93]		[2]
R1234ze(E)	[97]		[96]	[93]		[21]
R1234ze(Z)	[98]					[94]
R124	[99]		[82]	[93]		[2]
R1243zf	[98]
R125	[100]		[101]	[102]		[2]
R13	[87]	[87]	[82]	[82]		[2]
R134a	[103]		[84]	[82]		[2]
R13I1	[88]
R14	[87]	[87]	[82]	[82]		[2]
R141b	[3]		[82]	[82]		[2]
R142b	[3]		[82]	[82]		[2]
R143a	[104]		[84]	[85]		[21]
R152A	[105]		[106]	[93]		[2]
R161	[107]					[2]
R21	[87]	[87]				[2]
R218	[3]		[82]	[82]		[2]
R22	[108]		[84]	[93]		[2]
R227EA	[88]		[82]	[82]		[2]
R23	[109]		[110]	[110]		[2]
R236EA	[111]		[82]	[82]		[2]
R236FA	[112]		[82]	[82]		[2]
R245ca	[113]					[2]
R245fa	[114]		[82]	[93]		[2]
R32	[115]		[82]	[93]		[2]
R365MFC	[88]					[2]
R40	[29]					[21]
R404A	[116]		[117]	[118]		[119]
R407C	[116]		[117]	[118]		[119]
R41	[3]					[2]
R410A	[116]		[117]	[118]		[119]
R507A	[116]		[117]	[118]		[119]
RC318	[87]	[87]	[82]	[82]		[2]
SES36	[120]
SulfurDioxide	[121]					[2]
SulfurHexafluoride	[122]		[123]	[124]		[2]
Toluene	[3]		[125]	[126]		[2]
Water	[127]		[128]	[129]	[130]	[2]
Xenon	[3]					[2]
cis-2-Butene	[1]					[21]
m-Xylene	[39]		[40]	[131]		[21]
n-Butane	[58]		[132]	[133]	[58]	[2]
n-Decane	[3]		[134]	[135]		[2]
n-Dodecane	[136]		[137]	[137]		[2]
n-Heptane	[138]	[139]	[140]	[141]		[2]
n-Hexane	[142]		[143]	[144]		[2]
n-Nonane	[3]		[134]	[135]		[2]
n-Octane	[145]		[134]	[135]		[2]
n-Pentane	[146]		[23]	[67]	[61]	[2]
n-Propane	[147]		[148]	[149]	[61]	[2]
n-Undecane	[150]					[21]
o-Xylene	[39]		[40]	[151]		[21]
p-Xylene	[39]		[40]	[152]		[21]
trans-2-Butene	[1]					[21]

Ideal Curves

The so-called ideal curves demonstrate the good (or not!) extrapolative behavior of the equation of state. A few common curves are defined. They should ideally look similar to the figure shown below. Lemmon [100] provides a solid coverage of the ideal curves.

Ideal Curve

\[Z = 1\]

where \(Z\) can be given by

\[Z = \frac{pv}{RT}\]

Boyle Curve

Defined by

\[\left.\frac{\partial Z}{\partial v}\right|_{T} = 0\]

which can be expanded as (p can be expressed as a function of \(T\) and \(v\))

\[\left.\frac{\partial Z}{\partial v}\right|_{T} = \dfrac{p - \rho\left.\frac{\partial p}{\partial \rho}\right|_{T}}{RT}\]

Joule Inversion Curve

Defined by

\[\left.\frac{\partial Z}{\partial T}\right|_{v} = 0\]

which can be expanded as (p can be expressed as a function of \(T\) and \(v\))

\[\left.\frac{\partial Z}{\partial T}\right|_{v} = \frac{v}{R}\dfrac{T\left.\frac{\partial p}{\partial T}\right|_{v} - p}{T^2}\]

Joule-Thomson Curve

Defined by

\[\left.\frac{\partial Z}{\partial T}\right|_{p} = 0\]

which can be expanded as (v can be expressed as a function of \(T\) and \(p\))

\[\left.\frac{\partial Z}{\partial T}\right|_{p} = \frac{p}{R}\dfrac{T\left.\frac{\partial v}{\partial T}\right|_{p} - v}{T^2}\]

where

\[\left.\frac{\partial v}{\partial T}\right|_{p} = -\frac{\left.\frac{\partial \rho}{\partial T}\right|_{p}}{\rho^2}\]

Ideal Curves for Refrigerant R125

(Source code, png, .pdf)