化工热力学R-K方程

化工创新班化工热力学小组讨论应用RK方程计算气体流量组员：褚可凡、杨洋、周昱希指导老师：郭绪强1、问题重述根据每个物质的临界参数分别计算其值2、计算混合气体的RK方程iiba,在RK方程中：cCipTRa5.2272748.0cCipTRb08664.02、计算混合气体的RK方程各个物质的值，如下表1所示。iiba,组成yiTc/KPc/PaaibiCH40.8315190.44.6*10^63.212.98E-05C2H60.0906305.44.88*10^69.874.51E-05C3H80.0279369.84.25*10^618.286.27E-05iC40.0035408.23.65*10^627.258.06E-05nC40.0075425.23.80*10^628.998.06E-05nC50.0033469.73.37*10^641.921.00E-04CO20.001304.17.38*10^66.462.97E-05N20.0347126.23.39*10^61.562.68E-052、计算混合气体的RK方程根据混合物状态方程的混合规则对RK方程215.0NiiiayaNiiibyb1结合上表1中数据，计算各物质数据yi*ai^0.5，yi*bi如下表2所示。2、计算混合气体的RK方程组成yiyi*ai^0.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.03470.0433332989.30503E-072、计算混合气体的RK方程最终算得则混合物的RK方程为：08.4a51028.3b)1028.3(08.41028.352/15vvTvRTp3.计算摩尔体积#includestdio.h#includemath.hvoidequa(doublef[4],doublez[2]){doublea=f[0],b=f[1],c=f[2],d=f[3];doubleA,B,C,L,k,t,rad;doublegs1,gs2,y1,y2,p1,p2;doublex[3],s,r[3];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;3.计算摩尔体积if(A==0&&B==0)x[0]=x[1]=x[2]=-b/(3*a);else{if(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);x[0]=(-b-gs1-gs2)/(3*a);}3.计算摩尔体积else{t=(2*A*b-3*a*B)/(2*pow(A,1.5));rad=acos(t);x[0]=(-b-2*sqrt(A)*cos(rad/3))/(3*a);x[1]=(-b+sqrt(A)*(cos(rad/3)+sqrt(3)*sin(rad/3)))/(3*a);x[2]=(-b+sqrt(A)*(cos(rad/3)-sqrt(3)*sin(rad/3)))/(3*a);}}for(i=0;i3;i++)if(x[i]0)r[j++]=x[i];for(i=0;ij-1;i++)if(r[i]r[i+1]){s=r[i];r[i]=r[i+1];r[i+1]=s;}z[0]=r[j-1];z[1]=r[0];}voidrkeos(doubletc,doublepc,doublew,doublet,doublep){doublea,b,r=8.3146,e[4],z[2];a=4.08;b=3.28*pow(10,-5);e[0]=1.0;e[1]=-(r*t/p);e[2]=a/(pow(t,0.5)*p)-b*r*t/p-pow(b,2);e[3]=-a*b/(pow(t,0.5)*p);equa(e,z);printf(%lf,%lf\n,z[1]*1000,z[0]*1000);}voidmain(){doubletc=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/mol3.计算摩尔体积计算出来30度时的摩尔体积为24.896001，单位是L/mol3.计算摩尔体积计算出来标况时的摩尔体积为22.410519，单位是L/mol4.计算流量假设质量流量不变5030敬请各位老师、同学批评指正谢谢大家！