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

Added Nitrogen,Argon,Oxygen,Air viscosity

Browse Source 
Signed-off-by: Ian bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian bell
2014-05-21 16:19:10 +02:00
parent 646aa2f18f
commit b25b4e31d3
6 changed files with 404 additions and 36 deletions
Download Patch File Download Diff File Expand all files Collapse all files

94
dev/fluids/Air.json
View File

@@ -296,5 +296,97 @@
"rhoVtriple_units": "mol/m^3"
}
],
"NAME": "Air"
"NAME": "Air",
"TRANSPORT": {
"BibTeX": "Lemmon-IJT-2004",
"epsilon_over_k": 103.3,
"epsilon_over_k_units": "K",
"sigma_eta": 3.6e-10,
"sigma_eta_units": "m",
"viscosity": {
"dilute": {
"C": 2.66958e-08,
"a": [
0.431,
-0.4623,
0.08406,
0.005341,
-0.00331
],
"molar_mass": 0.0289586,
"molar_mass_units": "kg/mol",
"t": [
0,
1,
2,
3,
4
],
"type": "collision_integral"
},
"higher_order": {
"T_reduce": 132.6312,
"T_reduce_units": "K",
"a": [
1.072e-05,
1.122e-06,
2.019e-09,
-8.876e-06,
2.916e-08
],
"d1": [
1,
4,
9,
1,
8
],
"d2": [
1
],
"f": [
0
],
"g": [
1
],
"gamma": [
0,
0,
0,
-1,
-1
],
"h": [
0
],
"l": [
0,
0,
0,
1,
1
],
"p": [
1
],
"q": [
0
],
"rhomolar_reduce": 10447.7,
"rhomolar_reduce_units": "mol/m^3",
"t1": [
0.2,
0.05,
2.4,
0.6,
3.6
],
"t2": [
0
],
"type": "modified_Batschinski_Hildebrand"
}
}
}
}

99
dev/fluids/Argon.json
View File

@@ -382,5 +382,102 @@
"rhoVtriple_units": "mol/m^3"
}
],
"NAME": "Argon"
"NAME": "Argon",
"TRANSPORT": {
"BibTeX": "Lemmon-IJT-2004",
"epsilon_over_k": 143.2,
"epsilon_over_k_units": "K",
"sigma_eta": 3.35e-10,
"sigma_eta_units": "m",
"viscosity": {
"dilute": {
"C": 2.66958e-08,
"a": [
0.431,
-0.4623,
0.08406,
0.005341,
-0.00331
],
"molar_mass": 0.039948,
"molar_mass_units": "kg/mol",
"t": [
0,
1,
2,
3,
4
],
"type": "collision_integral"
},
"higher_order": {
"T_reduce": 150.687,
"T_reduce_units": "K",
"a": [
1.219e-05,
1.399e-05,
5.027e-09,
-1.893e-05,
-6.698e-06,
-3.827e-06
],
"d1": [
1,
2,
10,
5,
1,
2
],
"d2": [
1
],
"f": [
0
],
"g": [
1
],
"gamma": [
0,
0,
0,
-1,
-1,
-1
],
"h": [
0
],
"l": [
0,
0,
0,
2,
4,
4
],
"p": [
1
],
"q": [
0
],
"rhomolar_reduce": 13407.43,
"rhomolar_reduce_units": "mol/m^3",
"t1": [
0.42,
0.0,
0.95,
0.5,
0.9,
0.8
],
"t2": [
0
],
"type": "modified_Batschinski_Hildebrand"
}
}
}
}

94
dev/fluids/Nitrogen.json
View File

@@ -385,5 +385,97 @@
"rhoVtriple_units": "mol/m^3"
}
],
"NAME": "Nitrogen"
"NAME": "Nitrogen",
"TRANSPORT": {
"BibTeX": "Lemmon-IJT-2004",
"epsilon_over_k": 98.94,
"epsilon_over_k_units": "K",
"sigma_eta": 3.656e-10,
"sigma_eta_units": "m",
"viscosity": {
"dilute": {
"C": 2.66958e-08,
"a": [
0.431,
-0.4623,
0.08406,
0.005341,
-0.00331
],
"molar_mass": 0.02801348,
"molar_mass_units": "kg/mol",
"t": [
0,
1,
2,
3,
4
],
"type": "collision_integral"
},
"higher_order": {
"T_reduce": 126.192,
"T_reduce_units": "K",
"a": [
1.072e-05,
3.989e-08,
1.208e-09,
-7.402e-06,
4.62e-06
],
"d1": [
2,
10,
12,
2,
1
],
"d2": [
1
],
"f": [
0
],
"g": [
1
],
"gamma": [
0,
-1,
-1,
-1,
-1
],
"h": [
0
],
"l": [
0,
1,
1,
2,
3
],
"p": [
1
],
"q": [
0
],
"rhomolar_reduce": 11183.9,
"rhomolar_reduce_units": "mol/m^3",
"t1": [
0.1,
0.25,
3.2,
0.9,
0.3
],
"t2": [
0
],
"type": "modified_Batschinski_Hildebrand"
}
}
}
}

94
dev/fluids/Oxygen.json
View File

@@ -358,5 +358,97 @@
"rhoVtriple_units": "mol/m^3"
}
],
"NAME": "Oxygen"
"NAME": "Oxygen",
"TRANSPORT": {
"BibTeX": "Lemmon-IJT-2004",
"epsilon_over_k": 118.5,
"epsilon_over_k_units": "K",
"sigma_eta": 3.428e-10,
"sigma_eta_units": "m",
"viscosity": {
"dilute": {
"C": 2.66958e-08,
"a": [
0.431,
-0.4623,
0.08406,
0.005341,
-0.00331
],
"molar_mass": 0.0319988,
"molar_mass_units": "kg/mol",
"t": [
0,
1,
2,
3,
4
],
"type": "collision_integral"
},
"higher_order": {
"T_reduce": 154.581,
"T_reduce_units": "K",
"a": [
1.767e-05,
4.042e-07,
1.077e-10,
3.51e-07,
-1.367e-05
],
"d1": [
1,
5,
12,
8,
1
],
"d2": [
1
],
"f": [
0
],
"g": [
1
],
"gamma": [
0,
0,
0,
-1,
-1
],
"h": [
0
],
"l": [
0,
0,
0,
1,
2
],
"p": [
1
],
"q": [
0
],
"rhomolar_reduce": 13630,
"rhomolar_reduce_units": "mol/m^3",
"t1": [
0.05,
0.0,
2.1,
0.0,
0.5
],
"t2": [
0
],
"type": "modified_Batschinski_Hildebrand"
}
}
}
}

5
doc/transport_table/table.tex
View File

@@ -13,7 +13,7 @@
\hline\hline
Fluid & Reference & $\eta^0$ & $\eta^r$ \\
\hline
Nitrogen, argon, oxygen air & (data) Lemmon and Jacobsen 2004 & $\eta^0 = \dfrac{0.0266958\sqrt{MT}}{\sigma^2\Omega(T^*)}$\newline$\Omega(T^*)=\exp\left(\sum_{i=0}^{4}b_i[\ln T^*]^i\right)$ & $\eta^r = \sum_{i=1}^NN_i\tau^{\tau_i}\delta^{d_i}\exp(-\gamma_i\delta^{l_i})$\\\hline
Ammonia & (data) Fenghour 1995 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\Game(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta^r = B_{BV}\rho\eta^0(T) + \Delta\eta$\newline$\Delta\eta = \sum_i b_i(T)\rho^i$\\\hline
Dimethyl Ether & Meng 2012 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\mathfrak{S}(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta^r = \Delta\eta$\newline$\Delta\eta = \sum_{i=0}^{1}n_i\tau^{t_i}\delta^{d_i} + \sum_{i=2}^{6}n_i\tau^{t_i}\delta^{d_i}\exp(-\delta^{p_i})$\\\hline
Ethanol & Kiselev 2005 & $\eta^0 = \sum_i a_i T^{n_i}$ & $\eta^r = B_{RF}\rho\eta^0(T)+\Delta \eta$ \newline $\Delta\eta = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+f_1\left(\frac{\delta}{\delta_0(\tau)-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=g_2+g_3\sqrt{\tau}$\\\hline
@@ -36,7 +36,7 @@ R23 & Shan 2000 & $\eta^0 = \frac{5}{16}\sqrt{\frac{MkT}{1000\pi N}}\frac{10^{24
Helium & Arp 1998 & NASTY & NASTY \\\hline
n-Hexane & (data) Michailidou 2013 &$\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2S(T^*)}$\newline$S(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$& $\eta^r = \eta^0(T)\rho B_{RF} + \Delta\eta$\newline$\Delta\eta = (\rho_r^{2/3}T_r^{1/2})\left\lbrace\dfrac{c_0}{T_r}+\dfrac{c_1}{c_2+T_r+c_3\rho_r^2}+\dfrac{c_4(1+\rho_r)}{c_5 + c_6T_r+c_7\rho_r+\rho_r^2+c_8\rho_rT_r} \right\rbrace$ \\\hline
Hydrogen & Muzny 2013 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2S^*(T^*)}$\newline$S^*(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta^r = B_{RF}\rho\eta^0(T) + \Delta\eta$\newline$\Delta\eta = c_1\rho_r^2\left[c_2T_r+c_3/T_r+\frac{c_4\rho_r^2}{c_5+T_r}+c_6\rho_r^6\right]$\\\hline
\hline\hline\hline\hline
\hline\hline \multicolumn{4}{c}{DONE below this line}\\ \hline\hline
Propane & (data) Vogel 1998 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\mathfrak{S}(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta_h = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+f_1\left(\frac{\delta}{\delta_0(\tau)-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=g_1(1+g_2\tau^{1/2})$\\\hline
n-Butane & Vogel 1999 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\mathfrak{S}(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta_h = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+f_1\left(\frac{\delta}{\delta_0(\tau)-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=g_1(1+\displaystyle\sum_{l=2}g_l\tau^{(l-1)/2})$ \\\hline
Isobutane & Vogel 2000 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\mathfrak{S}(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta_h = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+f_1\left(\frac{\delta}{\delta_0(\tau)-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=g_1(1+\displaystyle\sum_{l=2}g_l\tau^{(l-1)/2})$ \\\hline
@@ -44,6 +44,7 @@ Octane, nonane, decane & (data) Huber 2004 & $\eta^0 = \dfrac{0.021357\sqrt{MT}
n-Dodecane & (data) Huber 2004 & $\eta^0 = \dfrac{0.021357\sqrt{MT}}{\sigma^2\mathfrak{S}(T^*)}$\newline$\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\eta^r = \eta^0(T)\rho B_{RF} + \Delta\eta$\newline$\Delta\eta = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+c_1\left(\frac{\delta}{\delta_0-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=c_2 +c_3\sqrt{\tau}$\\\hline
R125 & (data) Huber 2006 & $\eta^0 = \dfrac{5}{16}\sqrt{\dfrac{MkT}{\pi N}}\dfrac{1}{\sigma^2\Omega^*(T^*)}$\newline $\Omega(T^*)$ from Neufeld & $\eta^r = \eta^0(T)\rho B_{RF} + \Delta\eta$\newline$\Delta\eta = \displaystyle\sum_{i=2}^n\displaystyle\sum_{j=0}^me_{ij}\frac{\delta^i}{\tau_j}+c_1\left(\frac{\delta}{\delta_0-\delta}-\frac{\delta}{\delta_0(\tau)}\right)$ \newline $\delta_0(\tau)=c_2 +c_3\sqrt{\tau}$\\\hline
Carbon Dioxide & (data) Fenghour 1998 & $\eta^0 = \dfrac{1.00697\sqrt{T}}{\sigma^2\mathfrak{S}(T^*)}$ \newline $\mathfrak{S}(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$ & $\Delta\eta = d_{11}\rho + d_{21}\rho^2+\frac{d_{64}\rho^6}{(T^*)^3}+d_{81}\rho^8+\frac{d_{82}\rho^8}{T^*}$ \newline $\Delta\eta = d_{11}\rho_c\delta + d_{21}\rho_c^2\delta^2+\frac{d_{64}\rho_c^6(\varepsilon/k)^3\delta^6\tau^3}{(T_c)^3}+\rho_c^8d_{81}\delta^8+\frac{\rho_c^8d_{82}\delta^8(\varepsilon/k)\tau}{T_c}$\\\hline
Nitrogen, argon, oxygen air & (data) Lemmon and Jacobsen 2004 & $\eta^0 = \dfrac{0.0266958\sqrt{MT}}{\sigma^2\Omega(T^*)}$\newline$\Omega(T^*)=\exp\left(\sum_{i=0}^{4}b_i[\ln T^*]^i\right)$ & $\eta^r = \sum_{i=1}^NN_i\tau^{\tau_i}\delta^{d_i}\exp(-\gamma_i\delta^{l_i})$\\\hline
\hline\hline
\end{tabular}

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

@@ -53,30 +53,30 @@ vel("R125", "T", 400, "Dmolar", 30.631, "V", 17.070e-6, 1e-3),
//vel("Ammonia", "T", 398, "Dmolar", 21066.7, "V", 43.95e-6, 1e-3),
//// From Lemmon and Jacobsen, JPCRD, 2004
//vel("Nitrogen", "T", 100, "Dmolar", 1e-14, "V", 6.90349e-6, 1e-3),
//vel("Nitrogen", "T", 300, "Dmolar", 1e-14, "V", 17.8771e-6, 1e-3),
//vel("Nitrogen", "T", 100, "Dmolar", 25000, "V", 79.7418e-6, 1e-3),
//vel("Nitrogen", "T", 200, "Dmolar", 10000, "V", 21.0810e-6, 1e-3),
//vel("Nitrogen", "T", 300, "Dmolar", 5000, "V", 20.7430e-6, 1e-3),
//vel("Nitrogen", "T", 126.195, "Dmolar", 11180, "V", 18.2978e-6, 1e-3),
//vel("Argon", "T", 100, "Dmolar", 1e-14, "V", 8.18940e-6, 1e-3),
//vel("Argon", "T", 300, "Dmolar", 1e-14, "V", 22.7241e-6, 1e-3),
//vel("Argon", "T", 100, "Dmolar", 33000, "V", 184.232e-6, 1e-3),
//vel("Argon", "T", 200, "Dmolar", 10000, "V", 25.5662e-6, 1e-3),
//vel("Argon", "T", 300, "Dmolar", 5000, "V", 26.3706e-6, 1e-3),
//vel("Argon", "T", 150.69, "Dmolar", 13400, "V", 27.6101e-6, 1e-3),
//vel("Oxygen", "T", 100, "Dmolar", 1e-14, "V", 7.70243e-6, 1e-3),
//vel("Oxygen", "T", 300, "Dmolar", 1e-14, "V", 20.6307e-6, 1e-3),
//vel("Oxygen", "T", 100, "Dmolar", 35000, "V", 172.136e-6, 1e-3),
//vel("Oxygen", "T", 200, "Dmolar", 10000, "V", 22.4445e-6, 1e-3),
//vel("Oxygen", "T", 300, "Dmolar", 5000, "V", 23.7577e-6, 1e-3),
//vel("Oxygen", "T", 150.69, "Dmolar", 13600, "V", 24.7898e-6, 1e-3),
//vel("Air", "T", 100, "Dmolar", 1e-14, "V", 7.09559e-6, 1e-3),
//vel("Air", "T", 300, "Dmolar", 1e-14, "V", 18.5230e-6, 1e-3),
//vel("Air", "T", 100, "Dmolar", 28000, "V", 107.923e-6, 1e-3),
//vel("Air", "T", 200, "Dmolar", 10000, "V", 21.1392e-6, 1e-3),
//vel("Air", "T", 300, "Dmolar", 5000, "V", 21.3241e-6, 1e-3),
//vel("Air", "T", 132.64, "Dmolar", 10400, "V", 17.7623e-6, 1e-3),
vel("Nitrogen", "T", 100, "Dmolar", 1e-14, "V", 6.90349e-6, 1e-3),
vel("Nitrogen", "T", 300, "Dmolar", 1e-14, "V", 17.8771e-6, 1e-3),
vel("Nitrogen", "T", 100, "Dmolar", 25000, "V", 79.7418e-6, 1e-3),
vel("Nitrogen", "T", 200, "Dmolar", 10000, "V", 21.0810e-6, 1e-3),
vel("Nitrogen", "T", 300, "Dmolar", 5000, "V", 20.7430e-6, 1e-3),
vel("Nitrogen", "T", 126.195, "Dmolar", 11180, "V", 18.2978e-6, 1e-3),
vel("Argon", "T", 100, "Dmolar", 1e-14, "V", 8.18940e-6, 1e-3),
vel("Argon", "T", 300, "Dmolar", 1e-14, "V", 22.7241e-6, 1e-3),
vel("Argon", "T", 100, "Dmolar", 33000, "V", 184.232e-6, 1e-3),
vel("Argon", "T", 200, "Dmolar", 10000, "V", 25.5662e-6, 1e-3),
vel("Argon", "T", 300, "Dmolar", 5000, "V", 26.3706e-6, 1e-3),
vel("Argon", "T", 150.69, "Dmolar", 13400, "V", 27.6101e-6, 1e-3),
vel("Oxygen", "T", 100, "Dmolar", 1e-14, "V", 7.70243e-6, 1e-3),
vel("Oxygen", "T", 300, "Dmolar", 1e-14, "V", 20.6307e-6, 1e-3),
vel("Oxygen", "T", 100, "Dmolar", 35000, "V", 172.136e-6, 1e-3),
vel("Oxygen", "T", 200, "Dmolar", 10000, "V", 22.4445e-6, 1e-3),
vel("Oxygen", "T", 300, "Dmolar", 5000, "V", 23.7577e-6, 1e-3),
vel("Oxygen", "T", 154.6, "Dmolar", 13600, "V", 24.7898e-6, 1e-3),
vel("Air", "T", 100, "Dmolar", 1e-14, "V", 7.09559e-6, 1e-3),
vel("Air", "T", 300, "Dmolar", 1e-14, "V", 18.5230e-6, 1e-3),
vel("Air", "T", 100, "Dmolar", 28000, "V", 107.923e-6, 1e-3),
vel("Air", "T", 200, "Dmolar", 10000, "V", 21.1392e-6, 1e-3),
vel("Air", "T", 300, "Dmolar", 5000, "V", 21.3241e-6, 1e-3),
vel("Air", "T", 132.64, "Dmolar", 10400, "V", 17.7623e-6, 1e-3),
//// From Michailidou, JPCRD, 2013
//vel("Hexane", "T", 250, "Dmass", 1e-14, "V", 5.2584e-6, 1e-3),
@@ -107,12 +107,6 @@ vel("CO2", "T", 800, "Dmass", 407.828, "V", 48.74e-6, 1e-3),
//vel("R152A", "T", 384, "Dmass", 504.51, "V", 43.29e-6, 1e-3),
//vel("R152A", "T", 384, "Dmass", 239.35, "V", 21.01e-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),
//vel("R152A", "T", 384, "Dmass", 504.51, "V", 43.29e-6, 1e-3),
//vel("R152A", "T", 384, "Dmass", 239.35, "V", 21.01e-6, 1e-3),
//
//// Huber, JPCRD, 2008 and IAPWS
//vel("Water", "T", 298.15, "Dmass", 998, "V", 889.735100e-6, 1e-3),
//vel("Water", "T", 298.15, "Dmass", 1200, "V", 1437.649467e-6, 1e-3),