化工创新班化工热力学小组讨论应用RK方程计算气体流量组员：褚可凡、杨洋、周昱希指导老师：郭绪强1、问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
重述 根据每个物质的临界参数分别计算其值2、计算混合气体的RK方程在RK方程中：2、计算混合气体的RK方程 各个物质的值，如下表1所示。 组成 yi Tc/K Pc/Pa ai bi CH4 0.8315 190.4 4.6*10^6 3.21 2.98E-05 C2H6 0.0906 305.4 4.88*10^6 9.87 4.51E-05 C3H8 0.0279 369.8 4.25*10^6 18.28 6.27E-05 iC4 0.0035 408.2 3.65*10^6 27.25 8.06E-05 nC4 0.0075 425.2 3.80*10^6 28.99 8.06E-05 nC5 0.0033 469.7 3.37*10^6 41.92 1.00E-04 CO2 0.001 304.1 7.38*10^6 6.46 2.97E-05 N2 0.0347 126.2 3.39*10^6 1.56 2.68E-052、计算混合气体的RK方程 根据混合物状态方程的混合规则对RK方程结合上表1中数据，计算各物质数据yi*ai^0.5，yi*bi如下表2所示。2、计算混合气体的RK方程 组成 yi yi*ai^0.5 yi*bi CH4 0.8315 1.490510152 2.47913E-05 C2H6 0.0906 0.284624642 4.08419E-06 C3H8 0.0279 0.11929886 1.74868E-06 iC4 0.0035 0.018271777 2.81953E-07 nC4 0.0075 0.040381175 6.04504E-07 nC5 0.0033 0.021366066 3.31309E-07 CO2 0.001 0.002541033 2.96817E-08 N2 0.0347 0.043333298 9.30503E-072、计算混合气体的RK方程 最终算得则混合物的RK方程为：3.计算摩尔体积#include <stdio.h>#include <math.h>void equa(double f[4],double z[2]){double a=f[0],b=f[1],c=f[2],d=f[3];double A,B,C,L,k,t,rad;double gs1,gs2,y1,y2,p1,p2;double x[3],s,r[3];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&&B==0)x[0]=x[1]=x[2]=-b/(3*a);else{if(L>0){y1=A*b+(3*a*(-B+sqrt(L)))/2;y2=A*b+(3*a*(-B-sqrt(L)))/2;if(y1>0)gs1=pow(y1,1.0/3);elsegs1=-pow(fabs(y1),1.0/3);if(y2>0) 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); 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;i<3;i++) if(x[i]>0) r[j++]=x[i]; for(i=0;i<j-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];}void rkeos(double tc, double pc, double w,double t,double p){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敬请各位老师、同学批评指正 谢谢大家！此课件下载可自行编辑修改，此课件供参考！部分
内容
财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容
来源于网络，如有侵权请与我联系删除！感谢你的观看！
浙公网安备 33021202002483