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 prediction of critical density based on Tc/pc ratio as suggested by A. Kazakov

Browse Source
This commit is contained in:
Ian Bell
2016-02-09 19:38:59 -07:00
parent ef400b92d4
commit 83602f6b25
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/Backends/Cubics/CubicBackend.cpp
View File

@@ -17,10 +17,11 @@ void CoolProp::AbstractCubicBackend::get_linear_reducing_parameters(double &rhom
T_r = 0;
double v_r = 0;
const std::vector<double> Tc = cubic->get_Tc(), pc = cubic->get_pc();
std::vector<double> rhoc(2,10139.128); rhoc[1] = 6856.886685;
for (std::size_t i = 0; i < mole_fractions.size(); ++i){
T_r += mole_fractions[i]*Tc[i];
v_r += mole_fractions[i]/rhoc[i];
// Curve fit from all the pure fluids in CoolProp (thanks to recommendation of A. Kazakov)
double v_c_Lmol = 2.14107171795*(Tc[i]/pc[i]*1000)+0.00773144012514; // [L/mol]
v_r += mole_fractions[i]*v_c_Lmol/1000.0;
}
rhomolar_r = 1/v_r;
}