化工热力学方程优秀课件.ppt

《化工热力学方程优秀课件.ppt》由会员分享，可在线阅读，更多相关《化工热力学方程优秀课件.ppt（15页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、化工热力学方程第1页，本讲稿共15页1 1、问题重述、问题重述第2页，本讲稿共15页根据每个物质的临界参数分别计算其 值2 2、计算混合气体的、计算混合气体的RKRK方程方程在RK方程中：第3页，本讲稿共15页2 2、计算混合气体的、计算混合气体的RKRK方程方程各个物质的 值，如下表1所示。组成组成yiTc/KPc/PaaibiCH40.8315190.44.6*1063.212.98E-05C2H60.0906305.44.88*1069.874.51E-05C3H80.0279369.84.25*10618.286.27E-05iC40.0035408.23.65*10627.258.0

2、6E-05nC40.0075425.23.80*10628.998.06E-05nC50.0033469.73.37*10641.921.00E-04CO20.001304.17.38*1066.462.97E-05N20.0347126.23.39*1061.562.68E-05第4页，本讲稿共15页2 2、计算混合气体的、计算混合气体的RKRK方程方程根据混合物状态方程的混合规则对RK方程结合上表1中数据，计算各物质数据yi*ai0.5，yi*bi如下表2所示。第5页，本讲稿共15页2 2、计算混合气体的、计算混合气体的RKRK方程方程组成组成yiyi*ai0.5yi*biCH40.831

3、51.4905101522.47913E-05C2H60.09060.2846246424.08419E-06C3H80.02790.119298861.74868E-06iC40.00350.0182717772.81953E-07nC40.00750.0403811756.04504E-07nC50.00330.0213660663.31309E-07CO20.0010.0025410332.96817E-08N20.03470.0433332989.30503E-07第6页，本讲稿共15页2 2、计算混合气体的、计算混合气体的RKRK方程方程最终算得则混合物的RK方程为：第7页，本讲稿共

4、15页3.计算摩尔体积#include#includevoidequa(doublef4,doublez2)doublea=f0,b=f1,c=f2,d=f3;doubleA,B,C,L,k,t,rad;doublegs1,gs2,y1,y2,p1,p2;doublex3,s,r3;inti,j=0;A=b*b-3*a*c;B=b*c-9*a*d;C=c*c-3*b*d;L=B*B-4*A*C;第8页，本讲稿共15页3.计算摩尔体积if(A=0&B=0)x0=x1=x2=-b/(3*a);elseif(L0)y1=A*b+(3*a*(-B+sqrt(L)/2;y2=A*b+(3*a*(-B-s

5、qrt(L)/2;if(y10)gs1=pow(y1,1.0/3);elsegs1=-pow(fabs(y1),1.0/3);if(y20)gs2=pow(y2,1.0/3);elsegs2=-pow(fabs(y2),1.0/3);p1=(-2*b+gs1+gs2)/(6*a);p2=sqrt(3)*(gs1-gs2)/(6*a);x0=(-b-gs1-gs2)/(3*a);第9页，本讲稿共15页3.计算摩尔体积elset=(2*A*b-3*a*B)/(2*pow(A,1.5);rad=acos(t);x0=(-b-2*sqrt(A)*cos(rad/3)/(3*a);x1=(-b+sqrt

6、(A)*(cos(rad/3)+sqrt(3)*sin(rad/3)/(3*a);x2=(-b+sqrt(A)*(cos(rad/3)-sqrt(3)*sin(rad/3)/(3*a);for(i=0;i0)rj+=xi;for(i=0;iri+1)s=ri;ri=ri+1;ri+1=s;z0=rj-1;z1=r0;voidrkeos(doubletc,doublepc,doublew,doublet,doublep)第10页，本讲稿共15页doublea,b,r=8.3146,e4,z2;a=4.08;b=3.28*pow(10,-5);e0=1.0;e1=-(r*t/p);e2=a/(po

7、w(t,0.5)*p)-b*r*t/p-pow(b,2);e3=-a*b/(pow(t,0.5)*p);equa(e,z);printf(%lf,%lfn,z1*1000,z0*1000);voidmain()doubletc=305.4,pc=4880000.0,w=0.099,t=303.15,p=101000.0;rkeos(tc,pc,w,t,p);3.计算摩尔体积第11页，本讲稿共15页3.计算摩尔体积计算出来计算出来50度时的摩尔体积为度时的摩尔体积为0.590667，单位是，单位是L/mol第12页，本讲稿共15页3.计算摩尔体积计算出来计算出来30度时的摩尔体积为度时的摩尔体积为24.896001，单位是，单位是L/mol第13页，本讲稿共15页3.计算摩尔体积计算出来标况时的摩尔体积为计算出来标况时的摩尔体积为22.410519，单位是，单位是L/mol第14页，本讲稿共15页4.计算流量假设质量流量不变假设质量流量不变5030第15页，本讲稿共15页