【MFiX源代码】sutherland公式计算粘度

MFiX采用sutherland公式计算粘度只能适用于空气，因为源代码里只给了空气

源代码位置\mfix-19.3.1\model\calc_mu_g.f

其中F2O3为2/3

最主要的语句是136行

MU_G(IJK) = to_SI*1.7D-4 * &
   (T_G(IJK)/273.0D0)**1.5D0 * (383.D0/(T_G(IJK)+110.D0))

动力粘度以273K 时 1.7E-4为基础，依照该网格温度变化

$${\mu }_g(IJK)=1.7\times 10^{-4}\times (\frac{T_g}{273}{)}^{1.5}\times (\frac{383}{T_g+110})$$
其中Tg是当地网格的气体温度

如果温度为1043K则计算出来为0.00042（4.2E-4）
可见温度的影响是十分明显的

不同气体的sutherland常数是不同的，MFiX中用的 是110
其他气体参考http://www.360doc.com/content/12/0216/14/492439_187081813.shtml