| 5 | 1/1 | 返回列表 |
| 查看: 1134 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
zhuzhuge銅蟲 (小有名氣)
|
[求助]
閃蒸計(jì)算 已有1人參與
|
|
| 誰(shuí)有matlab編寫關(guān)于閃蒸計(jì)算的程序或者例子 |

銅蟲 (小有名氣)

|
% Authors: Nasri Zakia and Housam Binous % Flash distillation using the SRK equation of state function f=zakia2(x) z=[0.025 0.35 0.60 0.025]; Pc=[709.8 617.4 550.4 489.5]; Tc=[550.0 665.9 765.3 845.6]; w=[0.1064 0.1538 0.1954 0.2387]; T=197.33+459.67; P=20*14.696; Tre=T./Tc; Pre=P./Pc; m=0.480+1.574.*w-0.176.*w.^2; a=(1+m.*(1-Tre.^0.5)).^2; Ap=0.42747.*a.*Pre./Tre.^2; Bp=0.08664.*Pre./Tre; for i=1:4 for j=1:4 Ab(i,j)=(Ap(i)*Ap(j))^0.5; end end Av=0; for i=1:4 for j=1:4 Av=Av+x(i+4)*x(j+4)*Ab(i,j); end end Bv=0; for i=1:4 Bv=Bv+x(i+4)*Bp(i); end Bl=0; for i=1:4 Bl=Bl+x(i)*Bp(i); end Al=0; for i=1:4 for j=1:4 Al=Al+x(i)*x(j)*Ab(i,j); end end Zv=max(roots([1 -1 Av-Bv-Bv^2 -Av*Bv])); Zl=min(roots([1 -1 Al-Bl-Bl^2 -Al*Bl])); phiv=exp((Zv-1).*Bp/Bv-log(Zv-Bv)... -Av/Bv*log((Zv+Bv)/Zv).*(2.*Ap.^0.5./Av^0.5-Bp./Bv)); phil=exp((Zl-1).*Bp/Bl-log(Zl-Bl)... -Al/Bl*log((Zl+Bl)/Zl).*(2.*Ap.^0.5./Al^0.5-Bp./Bl)); K=phil./phiv; for i=1:4 f(i)=x(i+4)-K(i)*x(i); end for i=1:4 f(i+4)=x(i)-z(i)/(1+x(9)*(K(i)-1)); end f(9)=0; for i=1:4 f(9)=f(9)+z(i)*(K(i)-1)/(1+x(9)*(K(i)-1)); end SRK方程計(jì)算乙烷、丙烷、丁烷、戊烷閃蒸程序,地址http://www.mathworks.com/matlabc ... te/content/zakia2.m |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +13 | Mr. Z 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +6 | 101次希望 2026-03-28 | 6/300 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|