| 查看: 1578 | 回复: 1 | |||
[交流]
【求助】PR方程计算空气分离氮组分逸度系数问题 已有1人参与
|
|
在仿真模拟版发了这个求助帖（https://muchong.com/bbs/viewthread.php?tid=2343000&fpage=1），好几天了都没高手相助。再到这发一次，请达人指点，谢谢！ 小弟初学气液相平衡 根据已知的一组相平衡数据，学习用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，分别按照三个压缩因子求得逸度，但无 % 论怎么组合，都不能满足以上式子。 %问题：可能是因为压缩因子求取有问题?另问：当求得三个压缩因子时，到底哪个用作液相的压缩因子，哪个用作气相压缩因子? |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程数值解
已经有6人回复
- 纳氏试剂法和水杨酸法测氨氮
已经有39人回复
- 求助：关于肋片效率问题以及换热系数测量误差分析
已经有14人回复
- 南京工业大学ASPEN
已经有96人回复
- 纳米级载体负载活性组分，再压片磨碎成20-40目，与直接的20-40目的载体负载 有区别吗
已经有16人回复
- 最新ACS有关锂空气电池催化剂的文章
已经有112人回复
- 升温速率与扩散系数的关系
已经有5人回复
- 有关磁流体方程计算的询问
已经有8人回复
- [原创]一种基于能带计算结果的载流子有效质量计算方法
已经有174人回复
- 拟合非线性方程系数的精度问题
已经有5人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 2012年理论与高性能计算化学国际会议
已经有60人回复
- PR系列杂志投稿的copy right transfer。
已经有7人回复
- 如何处理自旋计算产生的能带
已经有3人回复
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求助这个非线性方程改怎么解啊？？
已经有4人回复
2楼2016-06-17 11:19:29
