github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-24 05:18:04 -05:00
Code Issues Packages Projects Releases Wiki Activity

Added highly approximate methods for conductivity and viscosity for mixtures; Closes #307

Browse Source 
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2015-02-02 14:01:12 -05:00
parent cf4a7852e0
commit 23912165fc
1 changed files with 17 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -20,6 +20,7 @@
//#include "CoolProp.h"
#include "HelmholtzEOSMixtureBackend.h"
#include "HelmholtzEOSBackend.h"
#include "Fluids/FluidLibrary.h"
#include "Solvers.h"
#include "MatrixMath.h"
@@ -445,7 +446,14 @@ long double HelmholtzEOSMixtureBackend::calc_viscosity(void)
}
else
{
throw NotImplementedError(format("viscosity not implemented for mixtures"));
set_warning_string("Mixture model for viscosity is highly approximate");
long double summer = 0;
for (std::size_t i = 0; i < mole_fractions.size(); ++i){
shared_ptr<HelmholtzEOSBackend> HEOS(new HelmholtzEOSBackend(components[i]));
HEOS->update(DmolarT_INPUTS, _rhomolar, _T);
summer += mole_fractions[i]*log(HEOS->viscosity());
}
return exp(summer);
}
}
long double HelmholtzEOSMixtureBackend::calc_conductivity_background(void)
@@ -543,7 +551,14 @@ long double HelmholtzEOSMixtureBackend::calc_conductivity(void)
}
else
{
throw NotImplementedError(format("viscosity not implemented for mixtures"));
set_warning_string("Mixture model for conductivity is highly approximate");
long double summer = 0;
for (std::size_t i = 0; i < mole_fractions.size(); ++i){
shared_ptr<HelmholtzEOSBackend> HEOS(new HelmholtzEOSBackend(components[i]));
HEOS->update(DmolarT_INPUTS, _rhomolar, _T);
summer += mole_fractions[i]*HEOS->conductivity();
}
return summer;
}
}
long double HelmholtzEOSMixtureBackend::calc_Ttriple(void)