diff --git a/dev/fluids/R404A.json b/dev/fluids/R404A.json index c3ffd700..ce9057b6 100644 --- a/dev/fluids/R404A.json +++ b/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" + } + } } \ No newline at end of file diff --git a/dev/fluids/R407C.json b/dev/fluids/R407C.json index 02c648e5..4a0a7599 100644 --- a/dev/fluids/R407C.json +++ b/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" + } + } } \ No newline at end of file diff --git a/dev/fluids/R410A.json b/dev/fluids/R410A.json index b57c38f7..e3f96ec5 100644 --- a/dev/fluids/R410A.json +++ b/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" + } + } } \ No newline at end of file diff --git a/dev/fluids/R507A.json b/dev/fluids/R507A.json index 9c6c302f..95b83bf9 100644 --- a/dev/fluids/R507A.json +++ b/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" + } + } } \ No newline at end of file diff --git a/doc/transport_table/table.tex b/doc/transport_table/table.tex index 79a6c994..e058d70a 100644 --- a/doc/transport_table/table.tex +++ b/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} diff --git a/src/Tests/CoolProp-Tests.cpp b/src/Tests/CoolProp-Tests.cpp index 30307e4e..76863757 100644 --- a/src/Tests/CoolProp-Tests.cpp +++ b/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