| 969 | } |
| 970 | } |
| 971 | double Viscosity(double T, double p, double psi_w) { |
| 972 | /* |
| 973 | Using the method of: |
| 974 | |
| 975 | 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 |
| 976 | |
| 977 | but using the detailed measurements for pure fluid from IAPWS formulations |
| 978 | */ |
| 979 | double mu_a = NAN, mu_w = NAN, Phi_av = NAN, Phi_va = NAN, Ma = NAN, Mw = NAN; |
| 980 | Mw = MM_Water(); |
| 981 | Ma = MM_Air(); |
| 982 | // Viscosity of dry air at dry-bulb temp and total pressure |
| 983 | Air->update(CoolProp::PT_INPUTS, p, T); |
| 984 | mu_a = Air->keyed_output(CoolProp::iviscosity); |
| 985 | // Saturated water vapor of pure water at total pressure |
| 986 | Water->update(CoolProp::PQ_INPUTS, p, 1); |
| 987 | mu_w = Water->keyed_output(CoolProp::iviscosity); |
| 988 | 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); //[-] |
| 989 | 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); //[-] |
| 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: |
no test coverage detected