| 24小時(shí)熱門(mén)版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 1125 | 回復(fù): 3 | ||
zhuzhuge銅蟲(chóng) (小有名氣)
|
[求助]
閃蒸計(jì)算 已有1人參與
|
| 誰(shuí)有matlab編寫(xiě)關(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 |
|
銅蟲(chóng) (小有名氣)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南航材料專(zhuān)317分求調(diào)劑 +4 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研]
|
pk3725069 2026-03-19 | 16/800 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
[考研]
|
孅華 2026-03-22 | 4/200 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 材料工程專(zhuān)碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|