| 4 | 1/1 | 返回列表 |
| 查看: 1126 | 回復(fù): 3 | |||
zhuzhuge銅蟲 (小有名氣)
|
[求助]
閃蒸計(jì)算 已有1人參與
|
| 誰有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 |
|
銅蟲 (小有名氣)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|