github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-08 21:48:11 -05:00
Code Issues Packages Projects Releases Wiki Activity

R23 conductivity hardcoded

Browse Source 
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-05-27 10:00:28 +02:00
parent b3c2c03192
commit c2a1fa7317
9 changed files with 57 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
dev/fluids/R125.json
View File

@@ -336,17 +336,17 @@
"residual": {
"B": [
-0.0072941,
0.0110497,
0.0416339,
-0.0289236,
-0.0311487,
0.01105,
0.041634,
-0.028924,
-0.031149,
0.0278399,
0.0112682,
-0.01211,
-0.00138322,
0.00211196
],
"T_reducing": 339.17,
"T_reducing": 339.173,
"T_reducing_units": "K",
"d": [
1,

4
dev/fluids/R23.json
View File

@@ -273,6 +273,10 @@
],
"NAME": "R23",
"TRANSPORT": {
"conductivity": {
"BibTeX": "Shan-ASHRAE-2000",
"hardcoded": "R23"
},
"viscosity": {
"BibTeX": "Shan-ASHRAE-2000",
"hardcoded": "R23"

12
doc/transport_table/table.tex
View File

@@ -48,20 +48,20 @@ For a description of the conversion between $\mathfrak{S}$ and $\Omega$, see Ves
\begin{tabular}{cccp{3in}p{8in}p{3in}}
\hline\hline
Fluid & Data & Reference & Dilute & Residual & Critical \\ \hline
Nitrogen, Argon, Oxygen, Air & & Lemmon 2004 & & &\\
CO2 & & Vesovic JPCRD 1990 \\
Ammonia & & Tufeu 1984 \\
Ethane & & Friend JPCRD 1991 \\
Helium & & Hands Cryo 1981\\
R123 & & Laesecke IJR 1996 \\
R23 & & Shan ASHRAE 2000 \\
Water & & Huber JPCRD 2012\\
\hline\hline STANDARD BELOW THIS LINE \\ \hline\hline
R152A & & Krauss IJT 1996 \\
R404A, R407C, R401A, R507A & & Geller IJT 2001 \\
\hline\hline DONE BELOW THIS LINE\\ \hline\hline
Water & & Huber JPCRD 2012\\
Nitrogen, Argon, Oxygen, Air & & Lemmon 2004 & & &\\
Ethane & & Friend JPCRD 1991 \\
CO2 & & Vesovic JPCRD 1990 \\
R123 & & Laesecke IJR 1996 \\
R152A & & Krauss IJT 1996 \\
R1234ze(E) & & Perkins JCED 2011 \\
Propane & & Marsh JCED 2002\\
Dodecane & & Huber EF 2004 \\

1
include/CoolPropFluid.h
View File

@@ -197,6 +197,7 @@ public:
};
enum ConductivityDiluteEnum {
CONDUCTIVITY_HARDCODED_WATER,
CONDUCTIVITY_HARDCODED_R23,
CONDUCTIVITY_NOT_HARDCODED
};
ViscosityDiluteVariables viscosity_dilute;

3
src/Backends/Helmholtz/Fluids/FluidLibrary.h
View File

@@ -613,6 +613,9 @@ protected:
if (!target.compare("Water")){
fluid.transport.hardcoded_conductivity = CoolProp::TransportPropertyData::CONDUCTIVITY_HARDCODED_WATER; return;
}
else if (!target.compare("R23")){
fluid.transport.hardcoded_conductivity = CoolProp::TransportPropertyData::CONDUCTIVITY_HARDCODED_R23; return;
}
else{
throw ValueError(format("hardcoded residual conductivity term [%s] is not understood for fluid %s",target.c_str(), fluid.name.c_str()));
}

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

@@ -221,6 +221,8 @@ long double HelmholtzEOSMixtureBackend::calc_conductivity(void)
{
case CoolProp::TransportPropertyData::CONDUCTIVITY_HARDCODED_WATER:
return TransportRoutines::conductivity_hardcoded_water(*this);
case CoolProp::TransportPropertyData::CONDUCTIVITY_HARDCODED_R23:
return TransportRoutines::conductivity_hardcoded_R23(*this);
default:
throw ValueError(format("hardcoded viscosity type [%d] is invalid for fluid %s", components[0]->transport.hardcoded_conductivity, name().c_str()));
}

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

@@ -683,5 +683,31 @@ long double TransportRoutines::conductivity_hardcoded_water(HelmholtzEOSMixtureB
return (lambdabar_0*lambdabar_1+lambdabar_2)*lambdastar;
}
long double TransportRoutines::conductivity_hardcoded_R23(HelmholtzEOSMixtureBackend &HEOS){
double B1 = -2.5370, // [mW/m/K]
B2 = 0.05366, // [mW/m/K^2]
C1 = 0.94215, // [-]
C2 = 0.14914, // [mW/m/K^2]
DeltaGstar = 2508.58, //[J/mol]
rhoL = 68.345, // [mol/dm^3] = [mol/L]
rhocbar = 7.5114, // [mol/dm^3]
DELTAlambda_max = 25, //[mW/m/K]
Ru = 8.31451, // [J/mol/K]
Tc = 299.2793, //[K]
T = HEOS.T(); //[K]
double lambda_DG = B1 + B2*T;
double rhobar = HEOS.rhomolar()/1000; // [mol/L]
double lambda_L = C2*(rhoL*rhoL)/(rhoL-rhobar)*sqrt(T)*exp(rhobar/(rhoL-rhobar)*DeltaGstar/(Ru*T));
double chi = rhobar - rhocbar;
double tau = T - Tc;
double DELTAlambda_c = 4*DELTAlambda_max/((exp(chi)+exp(-chi))*(exp(tau)+exp(-tau)));
return (pow((rhoL-rhobar)/rhoL,C1)*lambda_DG+pow(rhobar/rhoL,C1)*lambda_L+DELTAlambda_c)/1e3;
}
}; /* namespace CoolProp */

3
src/Backends/Helmholtz/TransportRoutines.h
View File

@@ -167,7 +167,8 @@ public:
static long double conductivity_residual_polynomial_and_exponential(HelmholtzEOSMixtureBackend &HEOS);
static long double conductivity_hardcoded_water(HelmholtzEOSMixtureBackend &HEOS);
static long double conductivity_hardcoded_R23(HelmholtzEOSMixtureBackend &HEOS);
}; /* class TransportRoutines */

10
src/Tests/CoolProp-Tests.cpp
View File

@@ -121,7 +121,7 @@ vel("R123", "T", 265, "Dmass", 1545.8, "V", 627.1e-6, 1e-3),
vel("R123", "T", 265, "Dmass", 1.614, "V", 9.534e-6, 1e-3),
vel("R123", "T", 415, "Dmass", 1079.4, "V", 121.3e-6, 1e-3),
vel("R123", "T", 415, "Dmass", 118.9, "V", 15.82e-6, 1e-3),
//
// Krauss, IJT, 1996
vel("R152A", "T", 242, "Dmass", 1025.5, "V", 347.3e-6, 1e-3),
vel("R152A", "T", 242, "Dmass", 2.4868, "V", 8.174e-6, 1e-3),
@@ -314,7 +314,7 @@ vel("ParaHydrogen", "T", 18, "Dmass", 1e-13, "L", 13.643e-3, 1e-4),
vel("ParaHydrogen", "T", 18, "Dmass", 75, "L", 100.52e-3, 1e-4),*/
// Some of these don't work
vel("R125", "T", 341, "Dmass", 600, "L", 0.056564297849487, 2e-4),
vel("R125", "T", 341, "Dmass", 600, "L", 0.0565642978494, 2e-4),
vel("R125", "T", 200, "Dmass", 1e-13, "L", 0.007036843623086, 2e-4),
vel("IsoButane", "T", 390, "Dmass", 387.09520158645068, "L", 0.063039, 1e-4),
vel("IsoButane", "T", 390, "Dmass", 85.76703973869482, "L", 0.036603, 2e-4),
@@ -405,6 +405,12 @@ vel("Water", "T", 647.35, "Dmass", 322, "L", 1443.75556e-3, 1e-6),
vel("Water", "T", 647.35, "Dmass", 372, "L", 650.319402e-3, 1e-6),
vel("Water", "T", 647.35, "Dmass", 422, "L", 448.883487e-3, 1e-6),
vel("Water", "T", 647.35, "Dmass", 750, "L", 600.961346e-3, 1e-6),
// From Shan, ASHRAE, 2000
vel("R23", "T", 180, "Dmolar", 21097, "V", 143.19e-3, 1e-4),
vel("R23", "T", 420, "Dmolar", 7564, "V", 50.19e-3, 1e-4),
vel("R23", "T", 370, "Dmolar", 32.62, "V", 17.455e-3, 1e-4),
};
TEST_CASE_METHOD(TransportValidationFixture, "Compare thermal conductivities against published data", "[conductivity]")