二氧化碳粘度计算-1998版

最新推荐文章于 2025-05-14 15:59:36 发布

CODE-boy1

最新推荐文章于 2025-05-14 15:59:36 发布

阅读量515

点赞数 9

文章标签： excel 算法 microsoft visual studio code

本文链接：https://blog.csdn.net/weixin_43441028/article/details/136774419

版权

Function yita0(T) As Double


Dim ke, Tx, Ge, a0, a1, a2, a3, a4, e As Double


a0 = 0.235156
a1 = -0.491266
a2 = 0.05211155
a3 = 0.05347906
a4 = -0.01537102
ke = 1 / 251.196
Tx = ke * T
e = a0 + a1 * Log(Tx) + a2 * Log(Tx) ^ 2 + a3 * Log(Tx) ^ 3 + a4 * Log(Tx) ^ 4
Ge = 2.718282 ^ e

yita0 = 1.00697 * T ^ (1 / 2) / Ge


End Function

Function dyita(rou, T) As Double
Dim ke, Tx, d11, d21, d64, d81, d82 As Double
ke = 1 / 251.196
Tx = ke * T
d11 = 0.004071119
d21 = 0.00007198037
d64 = 2.411697E-17
d81 = 2.971072E-23
d82 = -1.627888E-23
 
dyita = d11 * rou + d21 * rou ^ 2 + d64 * (rou ^ 6) / (Tx * Tx * Tx) + d81 * (rou ^ 8) + d82 * (rou ^ 8) / Tx


End Function

Function yita(rou, T) As Double
yita = yita0(T) + dyita(rou, T)


End Function



Sub viscosity()

Dim R, T, vis As Double

R = 407.824
T = 800

vis = yita(R, T)

Debug.Print "vis="; vis; "Pa.s"


End Sub