化工热力学R-K方程

1、化工创新班化工热力学小组讨论,应用RK方程计算气体流量,组员：褚可凡、杨洋、周昱希 指导老师：郭绪强,1、问题重述,根据每个物质的临界参数分别计算其 值,2、计算混合气体的RK方程,在RK方程中：,2、计算混合气体的RK方程,各个物质的 值，如下表1所示。,2、计算混合气体的RK方程,根据混合物状态方程的混合规则,对RK方程,结合上表1中数据，计算各物质数据yi*ai0.5，yi*bi如下表2所示。,2、计算混合气体的RK方程,2、计算混合气体的RK方程,最终算得,则混合物的RK方程为：,3.计算摩尔体积,#include #include void equa(double f4,double

2、 z2) double a=f0,b=f1,c=f2,d=f3; double A,B,C,L,k,t,rad; double gs1,gs2,y1,y2,p1,p2; double x3,s,r3; int i,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;,3.计算摩尔体积,if(A=0,if(y20) gs2=pow(y2,1.0/3); else gs2=-pow(fabs(y2),1.0/3); p1=(-2*b+gs1+gs2)/(6*a); p2=sqrt(3)*(gs1-gs2)/(6*a); x0=(-b-gs

3、1-gs2)/(3*a); ,3.计算摩尔体积,else t=(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(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;,void rkeos(do

4、uble tc, double pc, double w,double t,double p) ,double a,b,r=8.3146,e4,z2; a=4.08; b=3.28*pow(10,-5); e0=1.0;e1=-(r*t/p);e2=a/(pow(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); ,void main() double tc=305.4,pc=4880000.0,w=0.099,t=303.15,p=101000.0; rkeos(tc,pc,w,t,p); ,3.计算摩尔体积,3.计算摩尔体积,计算出来50度时的摩尔体积为0.590667，单位是L/mol,3.计算摩尔体积,计算出来30度时的摩尔体积为24.896001，单位是L/mol,3.计算摩尔体积,计算出来标