HFC pseudo-pure mixtures viscosity

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>

dev/fluids/R404A.json

@@ -287,5 +287,95 @@
       "rhoVtriple_units": "mol/m^3"
     }
   ], 
-  "NAME": "R404A"
+  "NAME": "R404A", 
+  "TRANSPORT": {
+    "viscosity": {
+      "BibTeX": "Geller-PURDUE-2000", 
+      "dilute": {
+        "a": [
+          9.766e-07, 
+          3.676e-08, 
+          2.938e-12
+        ], 
+        "t": [
+          0, 
+          1, 
+          2
+        ], 
+        "type": "powers_of_T"
+      }, 
+      "epsilon_over_k": 1e+99, 
+      "epsilon_over_k_units": "K", 
+      "higher_order": {
+        "T_reduce": 345.27, 
+        "T_reduce_units": "K", 
+        "a": [
+          1.08968786472e-06, 
+          4.152109565230651e-05, 
+          -4.710175334443676e-05, 
+          4.588082670553865e-05, 
+          -2.2593706031401724e-05, 
+          4.4806744468118286e-06
+        ], 
+        "d1": [
+          1, 
+          2, 
+          3, 
+          4, 
+          5, 
+          6
+        ], 
+        "d2": [
+          0
+        ], 
+        "f": [
+          0.0
+        ], 
+        "g": [
+          1
+        ], 
+        "gamma": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "h": [
+          0
+        ], 
+        "l": [
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1
+        ], 
+        "p": [
+          1
+        ], 
+        "q": [
+          0
+        ], 
+        "rhomolar_reduce": 4939.999999999999, 
+        "rhomolar_reduce_units": "mol/m^3", 
+        "t1": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "t2": [
+          0
+        ], 
+        "type": "modified_Batschinski_Hildebrand"
+      }, 
+      "sigma_eta": 1e+99, 
+      "sigma_eta_units": "m"
+    }
+  }
 }

dev/fluids/R407C.json

@@ -283,5 +283,95 @@
       "rhoVtriple_units": "mol/m^3"
     }
   ], 
-  "NAME": "R407C"
+  "NAME": "R407C", 
+  "TRANSPORT": {
+    "viscosity": {
+      "BibTeX": "Geller-PURDUE-2000", 
+      "dilute": {
+        "a": [
+          -1.507e-06, 
+          4.894e-08, 
+          -9.305e-12
+        ], 
+        "t": [
+          0, 
+          1, 
+          2
+        ], 
+        "type": "powers_of_T"
+      }, 
+      "epsilon_over_k": 1e+99, 
+      "epsilon_over_k_units": "K", 
+      "higher_order": {
+        "T_reduce": 359.345, 
+        "T_reduce_units": "K", 
+        "a": [
+          -1.3775231957200001e-06, 
+          6.0179007998135324e-05, 
+          -8.911399521418381e-05, 
+          7.350962141560484e-05, 
+          -2.788808852966893e-05, 
+          4.133412533885438e-06
+        ], 
+        "d1": [
+          1, 
+          2, 
+          3, 
+          4, 
+          5, 
+          6
+        ], 
+        "d2": [
+          0
+        ], 
+        "f": [
+          0.0
+        ], 
+        "g": [
+          1
+        ], 
+        "gamma": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "h": [
+          0
+        ], 
+        "l": [
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1
+        ], 
+        "p": [
+          1
+        ], 
+        "q": [
+          0
+        ], 
+        "rhomolar_reduce": 5260.000046401775, 
+        "rhomolar_reduce_units": "mol/m^3", 
+        "t1": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "t2": [
+          0
+        ], 
+        "type": "modified_Batschinski_Hildebrand"
+      }, 
+      "sigma_eta": 1e+99, 
+      "sigma_eta_units": "m"
+    }
+  }
 }

dev/fluids/R410A.json

@@ -283,5 +283,95 @@
       "rhoVtriple_units": "mol/m^3"
     }
   ], 
-  "NAME": "R410A"
+  "NAME": "R410A", 
+  "TRANSPORT": {
+    "viscosity": {
+      "BibTeX": "Geller-PURDUE-2000", 
+      "dilute": {
+        "a": [
+          -2.695e-06, 
+          5.85e-08, 
+          -2.129e-11
+        ], 
+        "t": [
+          0, 
+          1, 
+          2
+        ], 
+        "type": "powers_of_T"
+      }, 
+      "epsilon_over_k": 1e+99, 
+      "epsilon_over_k_units": "K", 
+      "higher_order": {
+        "T_reduce": 344.4943434, 
+        "T_reduce_units": "K", 
+        "a": [
+          4.1528450396712e-06, 
+          1.2187385701457372e-05, 
+          1.2660855545258757e-05, 
+          -1.0753223728015742e-05, 
+          1.9206178288211474e-06, 
+          3.6793471063956824e-07
+        ], 
+        "d1": [
+          1, 
+          2, 
+          3, 
+          4, 
+          5, 
+          6
+        ], 
+        "d2": [
+          0
+        ], 
+        "f": [
+          0.0
+        ], 
+        "g": [
+          1
+        ], 
+        "gamma": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "h": [
+          0
+        ], 
+        "l": [
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1
+        ], 
+        "p": [
+          1
+        ], 
+        "q": [
+          0
+        ], 
+        "rhomolar_reduce": 6323.999999999999, 
+        "rhomolar_reduce_units": "mol/m^3", 
+        "t1": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "t2": [
+          0
+        ], 
+        "type": "modified_Batschinski_Hildebrand"
+      }, 
+      "sigma_eta": 1e+99, 
+      "sigma_eta_units": "m"
+    }
+  }
 }

dev/fluids/R507A.json

@@ -287,5 +287,95 @@
       "rhoVtriple_units": "mol/m^3"
     }
   ], 
-  "NAME": "R507A"
+  "NAME": "R507A", 
+  "TRANSPORT": {
+    "viscosity": {
+      "BibTeX": "Geller-PURDUE-2000", 
+      "dilute": {
+        "a": [
+          -2.53e-06, 
+          5.626e-08, 
+          -2.323e-11
+        ], 
+        "t": [
+          0, 
+          1, 
+          2
+        ], 
+        "type": "powers_of_T"
+      }, 
+      "epsilon_over_k": 1e+99, 
+      "epsilon_over_k_units": "K", 
+      "higher_order": {
+        "T_reduce": 343.765, 
+        "T_reduce_units": "K", 
+        "a": [
+          2.60484792e-07, 
+          5.380047201384e-05, 
+          -7.967886221772443e-05, 
+          8.183382443771656e-05, 
+          -3.9504517246163185e-05, 
+          7.366291094257731e-06
+        ], 
+        "d1": [
+          1, 
+          2, 
+          3, 
+          4, 
+          5, 
+          6
+        ], 
+        "d2": [
+          0
+        ], 
+        "f": [
+          0.0
+        ], 
+        "g": [
+          1
+        ], 
+        "gamma": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "h": [
+          0
+        ], 
+        "l": [
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1
+        ], 
+        "p": [
+          1
+        ], 
+        "q": [
+          0
+        ], 
+        "rhomolar_reduce": 4964.029650250053, 
+        "rhomolar_reduce_units": "mol/m^3", 
+        "t1": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "t2": [
+          0
+        ], 
+        "type": "modified_Batschinski_Hildebrand"
+      }, 
+      "sigma_eta": 1e+99, 
+      "sigma_eta_units": "m"
+    }
+  }
 }

doc/transport_table/table.tex

@@ -14,16 +14,12 @@
 Fluid & Reference & $\eta^0$ & $\eta^r$ \\
 \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
-R404A, R410A, R507, R407 & Geller 2000 & $\eta^0 = \sum_i A_iT^i$&$\eta^r = \sum_j b_j\rho^j$ \\\hline
-\hline\hline
 
-SF6 & Quinones-Cisneros 2012 & $\eta^0 = \sum_i d_i T_r^{n_i}$ & FRICTION THEORY\\\hline
+\hline\hline
 H2S & Quinones-Cisneros 2012 & $\eta^0 = 8.7721\dfrac{\sqrt{T}}{S^*(T^*)}$ \newline $S^*(T^*) = \sum_i \frac{\alpha_i}{T^{*i}}$ & FRICTION THEORY\\\hline
 R23 & Shan 2000 & $\eta^0 = \frac{5}{16}\sqrt{\frac{MkT}{1000\pi N}}\frac{10^{24}}{\sigma^2\Omega^*(T^*)}$\newline $\Omega(T^*)=\exp\left(\sum_{i=0}^{4}a_i[\ln T^*]^i\right)$  & \\\hline
 Helium & Arp 1998 & NASTY & NASTY \\\hline
 Ethane & Friend 1991 & $\eta^0 = \dfrac{12.0085\sqrt{t}}{\Omega^{(2,2)*}(t)}$ \newline $\Omega^{(2,2)*}(t) = \left[\sum_i C_it^{(i-1)/3-1} \right]^{-1}$& $\Delta\eta = 15.977\left[\displaystyle\sum_i g_i\delta^{r_i}\tau^{s_i}\right]\left[1+\displaystyle\sum_{i=10}^{11}g_i\delta^{r_i}\tau^{s_i}\right]^{-1}$\\\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
@@ -41,6 +37,8 @@ Ethanol & Kiselev 2005 & $\eta^0 = \sum_i a_i T^{n_i}$ & $\eta^r = B_{RF}\rho\et
 Water & (data) Huber 2009 & & \\\hline
 Hydrogen & (from v5) 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\exp\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
 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
+SF6 & (data) Quinones-Cisneros 2012 & $\eta^0 = \sum_i d_i T_r^{n_i}$ & FRICTION THEORY\\\hline
+R404A, R410A, R507, R407 & (data) Geller 2000 & $\eta^0 = \sum_i A_iT^i$&$\eta^r = \sum_j b_j\rho^j$ \\\hline
 \hline\hline
 \end{tabular}
 

src/Tests/CoolProp-Tests.cpp

@@ -140,14 +140,39 @@ vel("Water", "T", 1173.15, "Dmass", 1, "V", 44.217245e-6, 1e-3),
 vel("Water", "T", 1173.15, "Dmass", 100, "V", 47.640433e-6, 1e-3),
 vel("Water", "T", 1173.15, "Dmass", 400, "V", 64.154608e-6, 1e-3),
 
-// 
-vel("SF6", "T", 300, "Dmass", 1e-14, "V", 15.2887e-6, 1e-3),
-vel("SF6", "T", 300, "Dmass", 5.92, "V", 15.3043e-6, 1e-3),
-vel("SF6", "T", 300, "Dmass", 1345.1, "V", 117.417e-6, 1e-3),
-vel("SF6", "T", 400, "Dmass", 1e-14, "V", 19.6796e-6, 1e-3),
-vel("SF6", "T", 400, "Dmass", 278.47, "V", 24.4272e-6, 1e-3),
-vel("SF6", "T", 400, "Dmass", 1123.8, "V", 84.7835e-6, 1e-3),
+// Quinones-Cisneros, JPCRD, 2012
+vel("SF6", "T", 300, "Dmass", 1e-14, "V", 15.2887e-6, 1e-4),
+vel("SF6", "T", 300, "Dmass", 5.92, "V", 15.3043e-6, 1e-4),
+vel("SF6", "T", 300, "Dmass", 1345.1, "V", 117.417e-6, 1e-4),
+vel("SF6", "T", 400, "Dmass", 1e-14, "V", 19.6796e-6, 1e-4),
+vel("SF6", "T", 400, "Dmass", 278.47, "V", 24.4272e-6, 1e-4),
+vel("SF6", "T", 400, "Dmass", 1123.8, "V", 84.7835e-6, 1e-4),
 
+// Vogel, JCED, 2011
+vel("H2S", "T", 405.76, "Dmass", 0.740, "V", 15.56e-6, 1e-2),
+// Quinones-Cisneros, JCED, 2012
+vel("H2S", "T", 373.15, "P", 1e6, "V", 15.9e-6, 1e-3),
+vel("H2S", "T", 373.15, "P", 10e6, "V", 48.3e-6, 1e-3),
+vel("H2S", "T", 423.15, "P", 1e6, "V", 16.5e-6, 1e-3),
+vel("H2S", "T", 423.15, "P", 10e6, "V", 21.2e-6, 1e-3),
+
+// Geller, Purdue Conference, 2000
+vel("R410A", "T", 243.15, "Q", 0, "V", 238.61e-6, 5e-2),
+vel("R410A", "T", 243.15, "Q", 1, "V", 10.37e-6, 5e-2),
+vel("R410A", "T", 333.15, "Q", 0, "V", 70.71e-6, 5e-2),
+vel("R410A", "T", 333.15, "Q", 1, "V", 19.19e-6, 5e-2),
+vel("R407C", "T", 243.15, "Q", 0, "V", 304.18e-6, 1e-2),
+vel("R407C", "T", 243.15, "Q", 1, "V", 9.83e-6, 1e-2),
+vel("R407C", "T", 333.15, "Q", 0, "V", 95.96e-6, 1e-2),
+vel("R407C", "T", 333.15, "Q", 1, "V", 16.38e-6, 1e-2),
+vel("R404A", "T", 243.15, "Q", 0, "V", 264.67e-6, 1e-2),
+vel("R404A", "T", 243.15, "Q", 1, "V", 10.13e-6, 1e-2),
+vel("R404A", "T", 333.15, "Q", 0, "V", 73.92e-6, 1e-2),
+vel("R404A", "T", 333.15, "Q", 1, "V", 18.56e-6, 1e-2),
+vel("R507A", "T", 243.15, "Q", 0, "V", 284.59e-6, 3e-2),
+vel("R507A", "T", 243.15, "Q", 1, "V", 9.83e-6, 1e-2),
+vel("R507A", "T", 333.15, "Q", 0, "V", 74.37e-6, 1e-2),
+vel("R507A", "T", 333.15, "Q", 1, "V", 19.35e-6, 1e-2),
 };
 
 class ViscosityValidationFixture