天然氣工程課程報(bào)告-在線瀏覽
2024-12-11 14:09本頁面
【正文】 39。 Dim i As Integer For i = 0 To 11 n(i) = Val(歡迎界面 .摩爾分?jǐn)?shù) (i)) / 100 Next i Dim p As Double, T As Double39。 p = Val(絕對(duì)工作壓力 ) T = Val(絕對(duì)工作溫度 ) Dim Tc(11) As Double, Pc(11) As Double39。 Tc(0) = …… Pc(11) = Dim Tpc As Double, Ppc As Double39。 Tpc = 0 Ppc = 0 For i = 0 To 11 Tpc = Tpc + n(i) * Tc(i) Ppc = Ppc + n(i) * Pc(i) Next i Dim Tpr As Double, ppr As Double39。 Tpr = T / Tpc ppr = p / Ppc Dim ρr As Double39。 Dim x As Double, x1 As Double, F As Double, F1 As Double absolution = 1 x = 0 Do While absolution F = () * (ppr / Tpr) * Exp(() * ((1 (1 / Tpr)) ^ 2)) + (x + x ^ 2 + x ^ 3 x ^ 4) / ((1 x) ^ 3) (( / Tpr) ( / (Tpr ^ 2)) + ( / (Tpr ^ 3))) * (x ^ 2) + (( / Tpr) ( / (Tpr ^ 2)) + ( / (Tpr ^ 3))) * (x ^ ( + ( / Tpr))) F1 = (1 + 4 * x + 4 * (x ^ 2) 4 * (x ^ 3) + x ^ 4) / ((1 x) ^ 4) (( / Tpr) ( / (Tpr ^ 2)) + ( / (Tpr ^ 3))) * x + ( + ( / Tpr)) * (( / Tpr) ( / (Tpr ^ 2)) + ( / (Tpr ^ 3))) * (x ^ ( + ( / Tpr))) x1 = x F / F1 absolution = Abs(x1 x) x = x1 Loop ρ r = x Dim z As Double39。 z = (1 + ρr + ρr ^ 2 ρr ^ 3) / ((1 ρr) ^ 3) ( / Tpr / (Tpr ^ 2) + / (Tpr ^ 3)) * ρr + ( / Tpr / (Tpr ^ 2) + / (Tpr ^ 3)) * (ρr ^ ( + ( / Tpr))) HY 模型 Z = Format(z, )‘HY 模型求得未校正的 Z’ 校正法 WA 校正法引入?yún)?shù)ε ,主要考慮了一些常見的極性分子 ( H2 S、 CO2 )的影響 ,希望用此參數(shù)來彌補(bǔ)常用計(jì)算方法的缺陷。 N 為氣體混合物中 H2 S 的摩爾分?jǐn)?shù)。 pci 為 i 組分的臨界壓力 ,kPa 。 pci′為 i 組分的校正臨界壓力 ,kPa。WA 校正 39。 該法應(yīng)用范圍是 :1. 0 ≤ Tpr ≤ 3 ,0. 2 ≤ ppr ≤ 30 。 部分代碼如下： 39。 Dim A1 As Double, A2 As Double, A3 As Double, A4 As Double, A5 As Double, A6 As Double, A7 As Double, A8 As Double, A9 As Double, A10 As Double, A11 As Double A1 = A2 = A3 = A4 = A5 = A6 = A7 = A8 = A9 = A10 = A11 = Dim ρr As Double Dim f11 As Double Dim f12 As Double Dim x As Double, x1 As Double, F As Double, F1 As Double absolution = 1 x = 0 Do While absolution F = () * ppr / Tpr + x + (A1 + A2 / Tpr + A3 / (Tpr ^ 3) + A4 / (Tpr ^ 4) + A5 / (Tpr ^ 5)) * (x ^ 2) + (A6 + A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 3) A9 * (A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 6) + A10 * (1 + A11 * (x ^ 2)) * ((x ^ 3) / (Tpr ^ 3)) * Exp(A11 * (x ^ 2)) f11 = 1 + 2 * (A1 + A2 / Tpr + A3 / (Tpr ^ 3) + A4 / (Tpr ^ 4) + A5 / (Tpr ^ 5)) * x + 3 * (A6 + A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 2) 6 * A9 * (A7 / Tpr + A8 / (Tpr ^ 2)) * (x ^ 5) f12 = (A10 / (Tpr ^ 3)) * (3 * (x ^ 2) + A11 * (3 * (x ^ 4) 2 * A11 * (x ^ 6))) * Exp(A11 * (x ^ 2)) F1 = f11 + f12 x1 = x F / F1 absolution = Abs(x1 x) x = x1 Loop ρ r = x Dim z As Double39。 z = 1 + (A1 + A2 / Tpr + A3 / (Tpr ^
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1