github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-02-12 06:44:59 -05:00
Code Issues Packages Projects Releases Wiki Activity

ECS ESRR need to be updated after conformal state solver - closes #337

Browse Source 
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-12-10 17:00:58 -05:00
parent 6b028024e9
commit 28c7b1dccd
1 changed files with 10 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/Backends/Helmholtz/TransportRoutines.cpp
View File

@@ -1035,8 +1035,12 @@ long double TransportRoutines::viscosity_ECS(HelmholtzEOSMixtureBackend &HEOS, H
conformal_state_solver(HEOS, HEOS_Reference, T0, rhomolar0);
// Update the reference fluid with the conformal state
// Update the reference fluid with the updated conformal state
HEOS_Reference.update(DmolarT_INPUTS, rhomolar0*psi, T0);
// Recalculate ESRR
f = HEOS.T()/T0;
h = rhomolar0/HEOS.rhomolar(); // Must be the ratio of MOLAR densities!!
// **********************
// Remaining calculations
@@ -1046,7 +1050,7 @@ long double TransportRoutines::viscosity_ECS(HelmholtzEOSMixtureBackend &HEOS, H
long double eta_resid = HEOS_Reference.calc_viscosity_background();
// The F factor
long double F_eta = sqrt(f)*pow(h, static_cast<long double>(-2.0/3.0))*sqrt(M/M0);
long double F_eta = sqrt(f)*pow(h, -2.0L/3.0L)*sqrt(M/M0);
// The total viscosity considering the contributions of the fluid of interest and the reference fluid [Pa-s]
long double eta = eta_dilute + eta_resid*F_eta;
@@ -1118,6 +1122,10 @@ long double TransportRoutines::conductivity_ECS(HelmholtzEOSMixtureBackend &HEOS
// Update the reference fluid with the conformal state
HEOS_Reference.update(DmolarT_INPUTS, rhomolar0*psi, T0);
// Recalculate ESRR
f = HEOS.T()/T0;
h = rhomolar0/HEOS.rhomolar(); // Must be the ratio of MOLAR densities!!
// The reference fluid's contribution to the conductivity [W/m/K]
long double lambda_resid = HEOS_Reference.calc_conductivity_background();