github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-04-23 03:00:17 -04:00
Code Issues Packages Projects Releases Wiki Activity

Implement PIP for REFPROP; closes #658

Browse Source
This commit is contained in:
Ian Bell
2015-05-08 17:50:41 -06:00
parent ddc0fb3fcd
commit a4981443a5
2 changed files with 17 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/Backends/REFPROP/REFPROPMixtureBackend.cpp
View File

@@ -522,6 +522,21 @@ CoolPropDbl REFPROPMixtureBackend::calc_melting_line(int param, int given, CoolP
}
}
CoolPropDbl REFPROPMixtureBackend::calc_PIP(void)
{
// Calculate the PIP factor of Venkatharathnam and Oellrich, "Identification of the phase of a fluid using
// partial derivatives of pressure, volume,and temperature without reference to saturation properties:
// Applications in phase equilibria calculations"
double t = _T, rho = _rhomolar/1000.0, // mol/dm^3
x = 0,p = 0,e = 0,h = 0,s = 0,cv = 0,cp = 0,w = 0,Z = 0,hjt = 0,A = 0,G = 0,
xkappa = 0,beta = 0,dPdrho = 0,d2PdD2 = 0,dPT = 0,drhodT = 0,drhodP = 0,
d2PT2 = 0,d2PdTD = 0,spare3 = 0,spare4 = 0;
//subroutine THERM2 (t,rho,x,p,e,h,s,cv,cp,w,Z,hjt,A,G,
// & xkappa,beta,dPdrho,d2PdD2,dPT,drhodT,drhodP,
// & d2PT2,d2PdTD,spare3,spare4);
THERM2dll(&t,&rho,&(mole_fractions[0]),&p,&e,&h,&s,&cv,&cp,&w,&Z,&hjt,&A,&G, &xkappa,&beta,&dPdrho,&d2PdD2,&dPT,&drhodT,&drhodP,&d2PT2,&d2PdTD,&spare3,&spare4);
return 2-rho*(d2PdTD/dPT- d2PdD2/dPdrho);
};
CoolPropDbl REFPROPMixtureBackend::calc_viscosity(void)
{

2
src/Backends/REFPROP/REFPROPMixtureBackend.h
View File

@@ -83,6 +83,8 @@ public:
/// Returns true if REFPROP is supported on this platform
bool REFPROP_supported(void);
CoolPropDbl calc_PIP(void);
CoolPropDbl calc_cpmolar_idealgas(void);
/// Set the fluids in REFPROP DLL by calling the SETUPdll function