【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-
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为基础,依照该网格温度变化
μg(IJK)=1.7×10−4×(Tg273)1.5×(383Tg+110) \mu_g(IJK) = 1.7\times10^{-4} \times (\frac{T_g}{273})^{1.5}\times(\frac{383}{T_g+110}) μg(IJK)=1.7×10−4×(273Tg)1.5×(Tg+110383)
其中Tg是当地网格的气体温度
如果温度为1043K则计算出来为0.00042(4.2E-4)
可见温度的影响是十分明显的
不同气体的sutherland常数是不同的,MFiX中用的 是110
其他气体参考http://www.360doc.com/content/12/0216/14/492439_187081813.shtml

DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。
更多推荐
所有评论(0)