【文章內容簡介】 059926011 ????? ?② 第一次迭代 ③ ???x ???x④ 由 wilson方程求得 4 2 7 ?? ??⑤ k P a.......PyPsiii73096384745581405992427716011 ???????? ?⑥ 010722901677396 ....PP ??????第二次迭代 P= 43 900599242 771 7396601111 .....PPyxs ??????③ 5 6 1 5 5 2222 ?????sPPyx?⑥ 01000407309672696 ....PP ??????④ 由 wilson方程求得 416811 .?? ??⑤ k P a.......P 7269638474645140599241681601 ??????5580384746451739640442059924168173966021.....x.....x????????⑦ P= 3 等壓泡點計算 已知 P 與 { xi }， 求 T與 { yi }。 siP② 任選一個 k組分， 由 飽和蒸汽壓方程求出各 由 活度系數(shù)方程求各 γi skisiiisk P/PxPP???③ ??isiiTxT 0①取溫度初值 kkk CPlnABT ?????????2,600 返回轉TTNoY e sTT ? ⑤ skP④由 和 k組分的 Antoine方程 計算出改進后的 溫度 T ? ?N,iPPxy siiii ?21?????iii yyy輸出 T和各 yi 將 yi 值歸一化 siP⑥由 飽和蒸汽壓方程求出各 由 活度系數(shù)方程求各 γi 計算甲醇 （ 1） 水 （ 2） 體系在 下的汽液平衡 。 已知 wilson方程能量參數(shù) g12g11=, g21g22=甲醇 、 水的 Antoine方程及液相摩爾體積與溫度的關系 V1= 102T+ 104T2 V2= 102T+ 104T2 單位 Pis, bar。 Vi , cm3/mol。 T, K。 29345536 2696 73111 .T..Pln s???13464438 1668 34112 .T..Pln s???例 74 ???????????????211221122112212211 ??????xxxxx)xxl n (ln???????????????211221122112121122 ??????xxxxx)xxl n (ln?????? ???RTgge x pVVll11121212??????? ???RTgge x pVVll22212121?? ?21 ,iPxPy siiii ?? ?解 : 計算 P=、 x1= T、 y y2 ① 取溫度初值 2934553 62 69 67 3110 1311 .T...lns ???1346443 81 66 83 4110 1312 .T...lns ???663371 .T s ?153732 .T s ? ss 9535822110 ???選 k=2 第一次迭代 ② 由 T0= b a ,b a ss 5 9 6 502 1 8 92 21 ??3 5 5 512 9 6 9064186543 211221 .,.,.V,.V ???? ??1 4 611 5 51 21 .,. ?? ??421 00146160596 50218 92155140013122..../....P/PxPPsisiiis???????? ?③ 第二次迭代 T0= ② 由 T0= ? ?b a ,b a ss 4 2 1 006 3 2 01 21 ??3 2 3 812 9 7 2054189742 211221 .,.,.V,.V ???? ??1 5 211 6 51 21 .,. ?? ??69895358263500 ...TT ????⑤ 263501346 44381668341142100 . ...ln ????④ 第三次迭代 T0= ② 由 T0= ? ?, ss 4056057931 21 ??90026350363490 ...TT ????⑤ 363491346 44381668341140560 . ...ln ????④ 4 0 5601521604 2 1006321165140 01312 ..../... .P s ??????③ 320712972053189142 211221 .,.,.V,.V ???? ??1525116611 21 .,. ?? ?? 4 4 90 ????? TT⑤ 273491346 44381668341140400 . ...ln ????④ 40400152516040560579311661140 01312 ..../... .P s ??????③ 7 2 4 800131 5 7 4 111 6 6 11401111 .. ...P Pxys????? ?b a rPb a rP ss 4 0 4 ,5 7 4 21 ??3 2 0 ,2 9 7 , , 211221 ?????? VV1526116611 21 .,. ?? ??275800131 4040015261602222 .. ...P Pxys????? ?⑥ 由 T= ? ??? 0 0 0 612 7 5 807 2 4 80 ...y i72440000617248011 ...yyyi??? ? 12 ??? yy4 等壓露點計算 已知 P 與 { yi }， 求 T與 { xi }。 siiii PPyx????iii xxx③ siP② 任選一個 k組分， 由 飽和蒸汽壓方程求出各 令各 γi =1 ??isii TyT 0① 取溫度初值 kskkk CPABT ???lnskP⑥ 由 和 k組分的 Antoine方程計算出改進后的溫度 T ??isksiiisk P/PyPP?⑤ ④ 由 活度系數(shù)方程求各 γi siiii PPyx????iii xxx輸出 T和 { xi } siP⑧由 飽和蒸汽壓方程求出各 ??????返回2,TTN轉8YεTT00⑦ 利用簡化 K進行汽液平衡計算 iii xyK ?Ki 應是 T、 P、 yi、 xi的函數(shù)，但對輕烴類體系，組成對 Ki值影響較小，可將 Ki看作 T、 P的函數(shù)，利用 PTK列線圖查 Ki值。 1 泡點計算 ?? iiiii xKxKy已知 P,xi 取 T初值 (已知 T,xi 取 P初值 ) 由 PTK 圖查 K ??? ? ii xK1是 否 調整 T（ 或 P) 1?? ii xK1?? ii xK所設溫度偏高（或壓力偏低） 所設溫度偏低（或壓力偏高） 2 露點計算 已知 P,yi 取 T初值 (已知 T,yi 取 P初值 ) 由 PTK 圖查 K ??? ? iiKy1是 ??iiiiiKyKyx否 調整 T（ 或 P) 1??iiKy1??iiKy所設溫度偏高 （ 或壓力偏低 ） 所設溫度偏低 （ 或壓力偏高 ） 給定 T、 P條件下 ， 若總組成為 zi的混合物進行減壓 ， 則可能分為相互成平衡的汽 、 液兩相的過程 閃蒸計算的目的 是確定汽 、 液相組成 (yi， xi)及汽化分率 (e＝ V/F)。 進料 F, zi P,T 汽相 V, yi 液相 L, xi 3 等溫閃蒸計算 相平衡方程 ， yi＝ Kixi (i＝ l， 2， …… N) 物料衡算方程 ， Fzi＝ Vyi+Lxi (i＝ 1， 2， …… N) Fzi＝ VKixi+(FV)xi ? ? 11 ??? eKzxiii ? ? 11 ??? eKKzyiiii? ? iii xKFVz ?????? ??? 11求解等溫閃蒸問題中 ,收斂性質比較好的一種函數(shù)形式是 0?? ?? ii xy? ? ? ?? ? 011 1 ??? ?? ? eK Kzefiii 在進行等溫閃蒸計算之前，可先用上式考察在規(guī)定的溫度、壓力下體系是否存在汽、液兩相。當 f(0) 和 f(1) 異號時 , 則體系存在汽、液兩相。 割線法解方程 f(x)=0 迭代公式 ? ? ? ?? ? ? ?? ?? ? ? ?? ?? ?? ?kkkkkkk xfxfxfxxxx111??????? 從迭代公式可以看出，割線法在各輪迭代中只需計算函數(shù)值，在這一點上它比牛頓法簡單，但是在做每一輪計算時，它卻需要前兩輪的信息，需要兩個初始點才能使算法起步。 例題 等溫閃蒸計算 進閃蒸器的物流組成 (摩爾分率 )為 20％ 正丁烷 (1)，50％ 正戊烷 (2)和 30％ 正己烷 (3)， 閃蒸壓力 l MPa， 閃蒸溫度 132℃ 。 計算汽液兩相分率和組成 。 已知在上述條件下 ， 各組分的汽液平衡常數(shù)值為 Kl＝ ， K2＝ ， K3＝ 。 f(0)＝ ＞ 0， f(1)＝－ ＜ 0， 所以在規(guī)定條件下 ， 該物流將存在于汽 、 液兩相 。 ? ? 0 .2 1 .1 3 0 .5 0 .1 0 .3 0 .4 11 1 .1 3 1 0 .1 1 0 .4 1fe e e e? ? ? ?? ? ?? ? ?? ? ? ?? ? 111iiizKfeKe?????取 e(0)＝ 0, e(1)＝ 1 ? ? ? ?? ? ? ?? ?? ? ? ?? ?? ?? ?kkkkkkk efefefeeee111???????e= ? ? ? ? ? ? ? ?? ?? ? ? ?? ?? ?? ? ? ? 72890056901530056900111010112 ....efefefeeee ???????????? ?? ? 00490