 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1779 | 回復(fù): 4 | |||
| 【有獎交流】積極回復(fù)本帖子，參與交流，就有機會分得作者 xiegangmai 的 20 個金幣 | |||
[交流]
【求助】PR方程計算空氣分離氮組分逸度系數(shù)問題
|
|||
|
小弟初學(xué)氣液相平衡 根據(jù)已知的一組相平衡數(shù)據(jù)，學(xué)習(xí)用PR方程計算組分逸度，但計算結(jié)果總是不對! 現(xiàn)在有一個問題迷惑，在用PR方程計算壓縮因子時，求得三個結(jié)果 不知哪個是氣態(tài)的壓縮因子，哪個是液態(tài)的壓縮因子？ 先貼上此問題的m文件，文件中有詳細(xì)說明 望各位前輩指教一二，不勝感激！ clear clc format long %在以下壓力、溫度以及組分條件下，氮氧兩組份滿足相平衡，相平衡常數(shù)KN=1.14，這段程序是為了學(xué)習(xí)并驗證PR方程 p=5.23e5; %壓力 T=96.6; %溫度 yN=0.904; %氣相氮組分 yO=1-yN; %氣態(tài)氧組分 xN=0.793; %液態(tài)氮組分 xO=1-xN; %液態(tài)氧組分 %物性基本常數(shù) %N2 TcN=126.2; %臨界溫度 pcN=34e5; %臨界壓力 wN=0.03726; %偏心因子 RN=296.8; %氣體常數(shù) mN=0.37464+1.54226*wN-0.26992*wN^2; %用來計算PR方程中的a(T) %O2 TcO=154.58; %同氮氣 pcO=50.43e5; wO=0.022; RO=259.84; mO=0.37464+1.54226*wO-0.26992*wO^2; %PR方程參數(shù)求取 %N2 TrN=T/TcN; %對比溫度 alphaN=(1+mN*(1-TrN^0.5))^2; %用來求a(T) acN=0.45724*RN^2*TcN^2/pcN; %用來求a(T) aN=acN*alphaN; %a(T)計算式 bN=0.0778*RN*TcN/pcN; %b計算式 AN=aN*p/(RN^2*T^2); %A計算式 BN=bN*p/RN/T; %B計算式 aaN=BN-1; bbN=AN-2*BN-3*BN^2; ccN=BN^3+BN^2-AN*BN; %以上三個參數(shù)用做求體積壓縮率的多項式系數(shù) %O2 注釋同氮氣 TrO=T/TcO; alphaO=(1+mO*(1-TrO^0.5))^2; acO=0.45724*RO^2*TcO^2/pcO; aO=acO*alphaO; bO=0.0778*RO*TcO/pcO; AO=aO*p/(RO^2*T^2); BO=bO*p/RO/T; aaO=BO-1; bbO=AO-2*BO-3*BO^2; ccO=BO^3+BO^2-AO*BO; %混合規(guī)則 %氣相 kNO=-0.023644+0.15427e-3; %兩元交互因子 am=yN*yN*aN+yO*yO*aO+2*yN*yO*sqrt(aN*aO)*(1-kNO); %混合a % am=(yN*sqrt(aN)+yO*sqrt(aO))^2; bm=yN*bN+yO*bO; %混合物的b Am=am*p/(RN^2*T^2); Bm=bm*p/RN/T; aam=Bm-1; bbm=Am-2*Bm-3*Bm^2; ccm=Bm^3+Bm^2-Am*Bm; %以上三個參數(shù)用做求體積壓縮率的多項式系數(shù) %液相 am1=xN*xN*aN+xO*xO*aO+2*xN*xO*sqrt(aN*aO)*(1-kNO); bm1=xN*bN+(1-xN)*bO; Am1=am1*p/(RN^2*T^2); Bm1=bm1*p/RN/T; aam1=Bm1-1; bbm1=Am1-2*Bm1-3*Bm1^2; ccm1=Bm1^3+Bm1^2-Am1*Bm1; %求體積壓縮因子 n=[1,aam,bbm,ccm]; z=roots(n) n1=[1,aam1,bbm1,ccm1]; z1=roots(n1) %求兩相N2的逸度 phi_yN=exp(bN*(z-1)/bm-log(z-Bm)-Am/(2*sqrt(2))*((2*(yO*sqrt(aN*aO)+yN*sqrt(aN*aN)))/am-bN/bm)*(log((z+2.414*Bm)./(z+0.414*Bm)))) phi_xN=exp(bN*(z1-1)/bm1-log(z1-Bm1)-Am1/(2*sqrt(2))*((2*((1-xN)*sqrt(aN*aO)+xN*sqrt(aN*aN)))/am1-bN/bm1)*(log((z1+2.414*Bm1)./(z1+0.414*Bm1)))) % 無論氣相還是液相都取得了三個體積壓縮因子，根據(jù)相平衡的原則，KN=1.14應(yīng)該等于phi_xN/phi_yN，分別按照三個壓縮因子求得逸度，但無 % 論怎么組合，都不能滿足以上式子。 %問題：可能是因為壓縮因子求取有問題?另問：當(dāng)求得三個壓縮因子時，到底哪個用作液相的壓縮因子，哪個用作氣相壓縮因子? |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
moirakaola» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求數(shù)學(xué)高手 下面的微分方程的數(shù)值解 急急急
已經(jīng)有3人回復(fù)
- 解一個類似的泊松方程
已經(jīng)有11人回復(fù)
- 求大俠幫忙解個方程
已經(jīng)有4人回復(fù)
- 熱壓縮試驗的流變應(yīng)力模型和本構(gòu)方程
已經(jīng)有7人回復(fù)
- 用quantum espresso 軟件來計算染料分子吸附在cluster 表面后的光譜性能
已經(jīng)有11人回復(fù)
- 用均方位移求擴散系數(shù)的疑問？
已經(jīng)有9人回復(fù)
- Fluent計算時的迭代問題
已經(jīng)有7人回復(fù)
- XRD圖用布拉格方程計算晶面間距的θ
已經(jīng)有4人回復(fù)
- 二維波動方程正演模擬，平面波解析解？
已經(jīng)有5人回復(fù)
- 微分方程解的個數(shù)與方程的階數(shù)關(guān)系
已經(jīng)有5人回復(fù)
- 請問如何根據(jù)XRD計算多晶粉末的ordering parameter
已經(jīng)有5人回復(fù)
- 有關(guān)castep能量計算，請高手賜教！
已經(jīng)有12人回復(fù)
- 一階微分方程求解，并優(yōu)化參數(shù)值
已經(jīng)有14人回復(fù)
- NEB計算過渡態(tài)和反應(yīng)路徑
已經(jīng)有9人回復(fù)
- TPR圖能否說明助劑的添加改善了催化劑的還原性能
已經(jīng)有24人回復(fù)
- 柴、蠟油組分分離
已經(jīng)有11人回復(fù)
- 如何用軟件畫出這個方程的圖形
已經(jīng)有13人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較
已經(jīng)有4人回復(fù)
2樓2010-09-21 09:46:11
moirakaola
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 90.5
- 帖子: 205
- 在線: 108.4小時
- 蟲號: 1660588
- 注冊: 2012-03-03
- 專業(yè): 理論環(huán)境化學(xué)
3樓2016-06-17 16:13:36
moirakaola
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 90.5
- 帖子: 205
- 在線: 108.4小時
- 蟲號: 1660588
- 注冊: 2012-03-03
- 專業(yè): 理論環(huán)境化學(xué)
4樓2016-06-17 16:16:17
moirakaola
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 90.5
- 帖子: 205
- 在線: 108.4小時
- 蟲號: 1660588
- 注冊: 2012-03-03
- 專業(yè): 理論環(huán)境化學(xué)
5樓2016-06-23 09:41:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiegangmai 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +3 | 煙汐憶海 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 293求調(diào)劑 +7 | zjl的號 2026-03-16 | 12/600 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[基金申請]
有必要更換申報口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|