| 查看: 1082 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
闪蒸计算已有1人参与
|
||
| 谁有matlab编写关于闪蒸计算的程序或者例子 |
回复此楼» 猜你喜欢
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有154人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
4楼2014-01-05 19:32:43
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhuzhuge: 金币+20, ★★★很有帮助 2014-01-05 19:27:42
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2014-01-06 08:20:38
zhuzhuge: 金币+20, ★★★很有帮助 2014-01-05 19:27:42
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2014-01-06 08:20:38
|
% 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方程计算乙烷、丙烷、丁烷、戊烷闪蒸程序，地址http://www.mathworks.com/matlabc ... te/content/zakia2.m |
2楼2013-12-30 16:59:04