| 查看: 726 | 回复: 7 | |||
20062993新虫 (小有名气)
|
[交流]
【求助】气液平衡程序 已有2人参与
|
|
用C语言编写一个气液平衡程序，液相为二元混合物，气相为非理想状态，可选任何一种状态方程，要循环实现。各位高手请指点一下。 [ Last edited by 20062993 on 2010-8-9 at 22:11 ] |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 甲缩醛——甲酸甲酯气液平衡
已经有5人回复
- 求 水和甲苯的 二相 气液平衡相图
已经有3人回复
- Aspen Plus 软件中气液平衡问题
已经有9人回复
- 正丙醇-异戊醇的气液平衡数据问题，求大神帮忙
已经有8人回复
- REFPROP怎么根据P和熵求焓？
已经有9人回复
- Detherm数据库气液平衡数据中的参考文献代表啥意思啊？
已经有1人回复
- 求三元物系不同温度下的气液液平衡数据
已经有7人回复
- 关于气液平衡的简单问题
已经有24人回复
- 工艺开发中蒸汽压和气液平衡数据的具体作用
已经有3人回复
- 甲苯——间二甲苯二元物系气液平衡关系
已经有3人回复
- 汽液平衡数据怎么算
已经有15人回复
- “异丁醇-水”体系气液平衡数据（国外文献）
已经有4人回复
- 2011注册化工工程师基础考试大纲 是不是和往年的有改变？？？
已经有5人回复
- 求助 计算泡点的问题
已经有8人回复
- aspen plus练习三：根据实验数据利用aspen plus回归二元交互参数实例
已经有19人回复
- 注册化工工程师
已经有10人回复
- 【求助】气液平衡的数据
已经有5人回复
- 【求助】 DMF-甲醇两相体系的气液平衡数据
已经有8人回复
2楼2010-08-10 11:18:24
20062993
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1563.8
- 散金: 538
- 红花: 3
- 帖子: 240
- 在线: 140.3小时
- 虫号: 1068987
- 注册: 2010-08-03
- 性别: GG
- 专业: 化学反应工程
3楼2010-08-10 11:33:16
20062993
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1563.8
- 散金: 538
- 红花: 3
- 帖子: 240
- 在线: 140.3小时
- 虫号: 1068987
- 注册: 2010-08-03
- 性别: GG
- 专业: 化学反应工程
4楼2010-08-10 11:35:05
5楼2010-08-10 14:58:03
20062993
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1563.8
- 散金: 538
- 红花: 3
- 帖子: 240
- 在线: 140.3小时
- 虫号: 1068987
- 注册: 2010-08-03
- 性别: GG
- 专业: 化学反应工程
nono2009:找个求根程序吧，很多程序集里都有。 2010-08-11 07:52:20
resonant:学一下matlab之类的东西就简单多了... 2010-08-11 09:25:24
resonant:学一下matlab之类的东西就简单多了... 2010-08-11 09:25:24
|
以下是我编的一部分程序，后面就进行不下去了，希望能给予指点，先谢了。 #include "stdio.h" #include "math.h" void main() { double pc1,pc2,tc1,tc2,tr1,tr2,w1,w2,k11,k12,k22,p,t,x1,x2,y1,y2,k1,k2,a1,a2,a11,a12,a22,amx,amy,b1,b2,bmx,bmy,ax,ay,bx,by,m1,m2,r=8.3145,i=1; scanf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f",&pc1,&pc2,&tc1,&tc2,&w1,&w2,&k11,&k12,&k22,&t,&x1,&x2); tr1=t/tc1; tr2=t/tc2; k1=0.37464+1.54226*w1-0.26992*w1*w1; k2=0.37464+1.54226*w2-0.26992*w2*w2; m1=1+k1*(1-sqrt(tr1)); m2=1+k2*(1-sqrt(tr2)); a1=0.45724*r*r*tc1*tc1*m1*m1/pc1; a2=0.45724*r*r*tc2*tc2*m2*m2/pc2; a11=a1*(1-k11); a12=sqrt(a1*a2)*(1-k12); a22=a2*(1-k22); b1=0.0778*r*tc1/pc1; b2=0.0778*r*tc2/pc2; (气相用PR状态方程，液相用活度系数法) { p=1; y1=0.5; y2=0.5; {amy=y1*y1*a11+2*y1*y2*a12+y2*y2*a22; bmy=b1*y1+b2*y2; ax=amx*p/(r*r*t*t); bx=bmy*p/(r*t); 之后就不知怎么编写了，碰到涉及关于Z的一元三次方程就不知道怎么求根了。 } } |
6楼2010-08-10 22:45:07
7楼2010-08-11 08:38:13
20062993
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1563.8
- 散金: 538
- 红花: 3
- 帖子: 240
- 在线: 140.3小时
- 虫号: 1068987
- 注册: 2010-08-03
- 性别: GG
- 专业: 化学反应工程
8楼2010-08-11 10:27:50