天然气压缩因子计算
按物性值
GB17747.3
Function B(T, HCHu, x1, x2, x3, x4, x5) As Double
Dim Bij() As Double
ReDim Preserve Bij(1 To 5, 1 To 5)
On Error Resume Next
Bij(1, 1) = -0.425468 + 0.002865 * T - 0.00000462073 * T * T
Bij(1, 1) = Bij(1, 1) + (0.000877118 - 0.00000556281 * T + 0.0000000088151 * T * T) *
HCHu
Bij(1, 1) = Bij(1, 1) + (-0.000000824747 + 0.00000000431436 * T - 6.08319E-12 * T * T)
* HCHu * HCHu
Bij(2, 2) = -0.1446 + 0.00074091 * T - 0.00000091195 * T * T
Bij(3, 3) = -0.86834 + 0.0040376 * T - 0.0000051657 * T * T
Bij(4, 4) = -0.00110596 + 0.0000813385 * T - 0.000000098722 * T * T
Bij(5, 5) = -0.13082 + 0.00060254 * T - 0.0000006443 * T * T
Bij(1, 2) = (0.72 + 0.00001875 * (320 - T) * (320 - T)) * (Bij(1, 1) + Bij(2, 2)) / 2
Bij(1, 3) = -0.865 * (Bij(1, 1) * Bij(3, 3)) ^ (1 / 2)
Bij(1, 4) = -0.052128 + 0.00027157 * T - 0.00000025 * T * T
Bij(1, 5) = -0.068729 - 0.00000239381 * T + 0.000000518195 * T * T
Bij(2, 3) = -0.339693 + 0.00161176 * T - 0.00000204429 * T * T
Bij(2, 4) = 0.012
B = x1 * x1 * Bij(1, 1) + 2 * x1 * x2 * Bij(1, 2) + 2 * x1 * x3 * Bij(1, 3)
B = B + 2 * x1 * x4 * Bij(1, 4) + 2 * x1 * x5 * Bij(1, 5) + x2 * x2 * Bij(2, 2)
B = B + 2 * x2 * x3 * Bij(2, 3) + 2 * x2 * x4 * Bij(2, 4)
B = B + x3 * x3 * Bij(3, 3) + x4 * x4 * Bij(4, 4) + x5 * x5 * Bij(5, 5)
End Function
Function C(T, HCH, x1, x2, x3, x4, x5) As Double
Dim Cijk() As Double
ReDim Preserve Cijk(1 To 5, 1 To 5, 1 To 5)
评论0
最新资源