化工热力学RK方程

1、化工创新班化工热力学小组讨论化工创新班化工热力学小组讨论应用应用RKRK方程计算气体流量方程计算气体流量组员：褚可凡、杨洋、周昱希指导老师：郭绪强1 1、问题重述、问题重述 根据每个物质的临界参数分别计算其 值2 2、计算混合气体的、计算混合气体的RKRK方程方程iiba,在RK方程中：cCipTRa5.2272748.0cCipTRb08664.02 2、计算混合气体的、计算混合气体的RKRK方程方程 各个物质的 值，如下表1所示。iiba,组成组成yiTc/KPc/PaaibiCH40.8315190.44.6*1063.21 2.98E-05C2H60.0906305.44.88*106

2、9.87 4.51E-05C3H80.0279369.84.25*10618.28 6.27E-05iC40.0035408.23.65*10627.25 8.06E-05nC40.0075425.23.80*10628.99 8.06E-05nC50.0033469.73.37*10641.92 1.00E-04CO20.001304.17.38*1066.46 2.97E-05N20.0347126.23.39*1061.56 2.68E-052 2、计算混合气体的、计算混合气体的RKRK方程方程 根据混合物状态方程的混合规则对RK方程215.0NiiiayaNiiibyb1结合上表1中数

3、据，计算各物质数据yi*ai0.5，yi*bi如下表2所示。2 2、计算混合气体的、计算混合气体的RKRK方程方程组成组成yiyi*ai0.5yi*biCH40.83151.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.034

4、70.0433332989.30503E-072 2、计算混合气体的、计算混合气体的RKRK方程方程 最终算得则混合物的RK方程为：08.4a51028.3b)1028.3(08.41028.352/15vvTvRTp3.计算摩尔体积#include #include void equa(double f4,double 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

5、*d;L=B*B-4*A*C;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-sqrt(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);3

6、.计算摩尔体积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(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(double tc, double pc, double w,

7、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