54 {
55 const double rho_c = 314;
56 const double CVF = 14.058;
57
58 const double sigma = 0.36502496e-09;
59 const double k = 1.38062e-23;
60 const double eps = 138.08483e-23;
61 const double c1 = 0.3125;
62 const double c2 = 2.0442e-49;
63
64 const double T_star =
T * k / eps;
66
67 double Omega = loop1_term<0>(T_star);
68 Omega += loop1_term<1>(T_star);
69 Omega += loop1_term<2>(T_star);
70 Omega += loop1_term<3>(T_star);
71 Omega += loop1_term<4>(T_star);
72
73 Omega = std::exp(Omega);
74
75
76 const double eta_0 = c1 * std::sqrt(c2 * T) / (sigma * sigma * Omega);
77
78 double sum = loop2_term<2>(rho);
79 sum += loop2_term<3>(rho);
80 sum += loop2_term<4>(rho);
81
82
83 const double eta_r =
84 CVF * 1e-6 * (
C[0] / (
rho -
C[1]) +
C[0] /
C[1] + sum);
85
86 return eta_0 + eta_r;
87 }
@ rho
density. For some models, rho substitutes p