化工热力学R-K方程优秀教案.pptVIP

下载本文档

2
0
约2.44千字
约 16页
2018-12-05 发布于江苏
举报
版权申诉

化工热力学R-K方程优秀教案.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

化工热力学R-K方程优秀教案

化工创新班化工热力学小组讨论应用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-05 2、计算混合气体的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 01.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-07 2、计算混合气体的RK方程最终算得则混合物的RK方程为： 3.计算摩尔体积 #include?stdio.h#include?math.hvoid?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==0B==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); 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;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];} void?rkeos(double?tc,?double?pc,?double?w,double?t,double?p){ double a,b,r=8.3146,e[4],z[2]; a=4.08; b=3.28*pow(10