| 990 | return (1 - psi_w) * mu_a / ((1 - psi_w) + psi_w * Phi_av) + psi_w * mu_w / (psi_w + (1 - psi_w) * Phi_va); |
| 991 | } |
| 992 | double Conductivity(double T, double p, double psi_w) { |
| 993 | /* |
| 994 | Using the method of: |
| 995 | |
| 996 | P.T. Tsilingiris, 2009, Thermophysical and transport properties of humid air at temperature range between 0 and 100 oC, Energy Conversion and Management, 49, 1098-1010 |
| 997 | |
| 998 | but using the detailed measurements for pure fluid from IAPWS formulations |
| 999 | */ |
| 1000 | double mu_a = NAN, mu_w = NAN, k_a = NAN, k_w = NAN, Phi_av = NAN, Phi_va = NAN, Ma = NAN, Mw = NAN; |
| 1001 | Mw = MM_Water(); |
| 1002 | Ma = MM_Air(); |
| 1003 | |
| 1004 | // Viscosity of dry air at dry-bulb temp and total pressure |
| 1005 | Air->update(CoolProp::PT_INPUTS, p, T); |
| 1006 | mu_a = Air->keyed_output(CoolProp::iviscosity); |
| 1007 | k_a = Air->keyed_output(CoolProp::iconductivity); |
| 1008 | // Conductivity of saturated pure water at total pressure |
| 1009 | Water->update(CoolProp::PQ_INPUTS, p, 1); |
| 1010 | mu_w = Water->keyed_output(CoolProp::iviscosity); |
| 1011 | k_w = Water->keyed_output(CoolProp::iconductivity); |
| 1012 | Phi_av = sqrt(2.0) / 4.0 * pow(1 + Ma / Mw, -0.5) * pow(1 + sqrt(mu_a / mu_w) * pow(Mw / Ma, 0.25), 2); //[-] |
| 1013 | Phi_va = sqrt(2.0) / 4.0 * pow(1 + Mw / Ma, -0.5) * pow(1 + sqrt(mu_w / mu_a) * pow(Ma / Mw, 0.25), 2); //[-] |
| 1014 | return (1 - psi_w) * k_a / ((1 - psi_w) + psi_w * Phi_av) + psi_w * k_w / (psi_w + (1 - psi_w) * Phi_va); |
| 1015 | } |
| 1016 | /** |
| 1017 | @param T Temperature in K |
| 1018 | @param p Pressure in Pa |
no test coverage detected