| 查看: 1199 | 回复: 3 | ||
[求助]
闪蒸计算 已有1人参与
|
| 谁有matlab编写关于闪蒸计算的程序或者例子 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有160人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
 |
3楼2013-12-31 16:53:21
4楼2014-01-05 19:32:43
30