热力学作业.doc_第1页
热力学作业.doc_第2页
热力学作业.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

Estimate the solubility of naphthalene in Sc-CO2 at 308.5K and pressure up to 300bar.（=bar，=125，混合性质参数为k12=0.068；l12=-0.218）解答:1.CO2-1组分，naphthalene(萘)-2组分，=bar，=125很低，接近于理想气体，所以=1所以上式可以简化为： 组分2在流体相中的摩尔分数 纯组分2在温度T下的饱和蒸汽压, Pa 纯组分2的摩尔体积, 组分2在流体相中的逸度系数 纯组分2在其饱和蒸汽压下的逸度系数。2.利用A、B、Z值求解萘在超临界CO2相中的逸度其中， ；物质临界压力临界温度压缩因子超临界CO2Naphthalene(萘)在给定的T=308.15K和给定的某一个压力P (1300 bar)下分别求出两种物质的对比温度和对比压力； ； ； 由上面两式可分别求得两种物质的Ai和Bi值。将混合性质参数k12=0.068；l12=-0.218代入混合规则表达式（呈一阶线性关系，因此A与a的混合规则相同，B与b的混合规则相同） 和 推导可知：；代入 ，所求y是未知数，所以用迭代方法。以y2=0为初值来对下面的所有公式进行迭代并得到一定精度的y2值，由y2=0代入得A、B即可得Z,因此得到，最后求得y2，与初值对比，否：重新迭代，是：相同就为所求值。3.迭代：在command 窗口输入压力p的取值 “p =1 20.5 68.7 100.4 152.3 190.5 230.6 267.3 300”并用命令“y=PR(p)”来调用下面的程序求解对应的y值。function y=PR(p)n=length(p);pc=73.82 40.5;tc=304.2 748.4;w=0.225 0.302;t=308.15;for j=1:n for i=1:2 tr(i)=t/tc(i); pr(i)=p(j)/pc(i); alpha(i)=(1+(0.3746+1.54226*w(i)-0.26992*w(i)2)*(1-sqrt(tr(i)2; A(i)=0.45724*alpha(i)*pr(i)/(tr(i)2; B(i)=0.07779*pr(i)/tr(i); end A(3)=sqrt(A(1)*A(2)*(1-0.068); B(3)=0.5*(B(1)+B(2)*(1+0.218); y2=0;y(j)=1;M=1; while (abs(y(j)-y2)=1e-6 & M=1e-8 & N p =1 20.5 68.7 100.4 152.3 190.5 230.6 267.3 300; y=PR(p); y2=0,0,0.007
人人文库> 全部分类> 专业文献 > 生活休闲

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论