| 查看: 2118 | 回复: 4 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 xiegangmai 的 20 个金币 | |||
[交流]
【求助】PR方程计算空气分离氮组分逸度系数问题
|
|||
|
小弟初学气液相平衡 根据已知的一组相平衡数据，学习用PR方程计算组分逸度，但计算结果总是不对! 现在有一个问题迷惑，在用PR方程计算压缩因子时，求得三个结果 不知哪个是气态的压缩因子，哪个是液态的压缩因子？ 先贴上此问题的m文件，文件中有详细说明 望各位前辈指教一二，不胜感激！ clear clc format long %在以下压力、温度以及组分条件下，氮氧两组份满足相平衡，相平衡常数KN=1.14，这段程序是为了学习并验证PR方程 p=5.23e5; %压力 T=96.6; %温度 yN=0.904; %气相氮组分 yO=1-yN; %气态氧组分 xN=0.793; %液态氮组分 xO=1-xN; %液态氧组分 %物性基本常数 %N2 TcN=126.2; %临界温度 pcN=34e5; %临界压力 wN=0.03726; %偏心因子 RN=296.8; %气体常数 mN=0.37464+1.54226*wN-0.26992*wN^2; %用来计算PR方程中的a(T) %O2 TcO=154.58; %同氮气 pcO=50.43e5; wO=0.022; RO=259.84; mO=0.37464+1.54226*wO-0.26992*wO^2; %PR方程参数求取 %N2 TrN=T/TcN; %对比温度 alphaN=(1+mN*(1-TrN^0.5))^2; %用来求a(T) acN=0.45724*RN^2*TcN^2/pcN; %用来求a(T) aN=acN*alphaN; %a(T)计算式 bN=0.0778*RN*TcN/pcN; %b计算式 AN=aN*p/(RN^2*T^2); %A计算式 BN=bN*p/RN/T; %B计算式 aaN=BN-1; bbN=AN-2*BN-3*BN^2; ccN=BN^3+BN^2-AN*BN; %以上三个参数用做求体积压缩率的多项式系数 %O2 注释同氮气 TrO=T/TcO; alphaO=(1+mO*(1-TrO^0.5))^2; acO=0.45724*RO^2*TcO^2/pcO; aO=acO*alphaO; bO=0.0778*RO*TcO/pcO; AO=aO*p/(RO^2*T^2); BO=bO*p/RO/T; aaO=BO-1; bbO=AO-2*BO-3*BO^2; ccO=BO^3+BO^2-AO*BO; %混合规则 %气相 kNO=-0.023644+0.15427e-3; %两元交互因子 am=yN*yN*aN+yO*yO*aO+2*yN*yO*sqrt(aN*aO)*(1-kNO); %混合a % am=(yN*sqrt(aN)+yO*sqrt(aO))^2; bm=yN*bN+yO*bO; %混合物的b Am=am*p/(RN^2*T^2); Bm=bm*p/RN/T; aam=Bm-1; bbm=Am-2*Bm-3*Bm^2; ccm=Bm^3+Bm^2-Am*Bm; %以上三个参数用做求体积压缩率的多项式系数 %液相 am1=xN*xN*aN+xO*xO*aO+2*xN*xO*sqrt(aN*aO)*(1-kNO); bm1=xN*bN+(1-xN)*bO; Am1=am1*p/(RN^2*T^2); Bm1=bm1*p/RN/T; aam1=Bm1-1; bbm1=Am1-2*Bm1-3*Bm1^2; ccm1=Bm1^3+Bm1^2-Am1*Bm1; %求体积压缩因子 n=[1,aam,bbm,ccm]; z=roots(n) n1=[1,aam1,bbm1,ccm1]; z1=roots(n1) %求两相N2的逸度 phi_yN=exp(bN*(z-1)/bm-log(z-Bm)-Am/(2*sqrt(2))*((2*(yO*sqrt(aN*aO)+yN*sqrt(aN*aN)))/am-bN/bm)*(log((z+2.414*Bm)./(z+0.414*Bm)))) phi_xN=exp(bN*(z1-1)/bm1-log(z1-Bm1)-Am1/(2*sqrt(2))*((2*((1-xN)*sqrt(aN*aO)+xN*sqrt(aN*aN)))/am1-bN/bm1)*(log((z1+2.414*Bm1)./(z1+0.414*Bm1)))) % 无论气相还是液相都取得了三个体积压缩因子，根据相平衡的原则，KN=1.14应该等于phi_xN/phi_yN，分别按照三个压缩因子求得逸度，但无 % 论怎么组合，都不能满足以上式子。 %问题：可能是因为压缩因子求取有问题?另问：当求得三个压缩因子时，到底哪个用作液相的压缩因子，哪个用作气相压缩因子? |
回复此楼» 本帖已获得的红花（最新10朵）
moirakaola» 猜你喜欢
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有12人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有8人回复
-
我不理解！
已经有24人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
国自然面上复盘~欢迎讨论
已经有19人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求数学高手 下面的微分方程的数值解 急急急
已经有3人回复
- 解一个类似的泊松方程
已经有11人回复
- 求大侠帮忙解个方程
已经有4人回复
- 热压缩试验的流变应力模型和本构方程
已经有7人回复
- 用quantum espresso 软件来计算染料分子吸附在cluster 表面后的光谱性能
已经有11人回复
- 用均方位移求扩散系数的疑问？
已经有9人回复
- Fluent计算时的迭代问题
已经有7人回复
- XRD图用布拉格方程计算晶面间距的θ
已经有4人回复
- 二维波动方程正演模拟，平面波解析解？
已经有5人回复
- 微分方程解的个数与方程的阶数关系
已经有5人回复
- 请问如何根据XRD计算多晶粉末的ordering parameter
已经有5人回复
- 有关castep能量计算，请高手赐教！
已经有12人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- NEB计算过渡态和反应路径
已经有9人回复
- TPR图能否说明助剂的添加改善了催化剂的还原性能
已经有24人回复
- 柴、蜡油组分分离
已经有11人回复
- 如何用软件画出这个方程的图形
已经有13人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
2楼2010-09-21 09:46:11
3楼2016-06-17 16:13:36
4楼2016-06-17 16:16:17
5楼2016-06-23 09:41:59
5