github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-02-09 05:15:45 -05:00
Code Issues Packages Projects Releases Wiki Activity

HEOS can call transport properties

Browse Source 
Signed-off-by: Ian bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian bell
2014-05-19 22:39:04 +02:00
parent 69f76386b1
commit 41dacf83c6
2 changed files with 25 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
View File

@@ -23,6 +23,7 @@
#include "MatrixMath.h"
#include "VLERoutines.h"
#include "FlashRoutines.h"
#include "TransportRoutines.h"
namespace CoolProp {
@@ -127,6 +128,29 @@ long double HelmholtzEOSMixtureBackend::calc_surface_tension(void)
throw NotImplementedError(format("surface tension not implemented for mixtures"));
}
}
long double HelmholtzEOSMixtureBackend::calc_viscosity(void)
{
if (is_pure_or_pseudopure)
{
// Dilute part
long double eta_dilute = TransportRoutines::dilute_gas_viscosity(*this);
// Residual part
long double B_eta_initial = TransportRoutines::initial_density_dependence_viscosity_term(*this);
long double rho = rhomolar();
long double initial_part = eta_dilute*B_eta_initial*rhomolar();
long double delta_eta_h = TransportRoutines::modified_Batschinski_Hildebrand_viscosity_term(*this);
long double eta_residual = initial_part + delta_eta_h;
// Critical part
long double eta_critical = 0;
return eta_dilute + eta_residual + eta_critical;
}
else
{
throw NotImplementedError(format("viscosity not implemented for mixtures"));
}
}
long double HelmholtzEOSMixtureBackend::calc_Ttriple(void)
{
double summer = 0;

1
src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.h
View File

@@ -117,6 +117,7 @@ public:
long double calc_d2alpha0_dTau2(void);
long double calc_surface_tension(void);
long double calc_viscosity(void);
long double calc_Tmax(void);
long double calc_pmax(void);