化工熱力學(xué)RK方程

化工新班化工力學(xué)小用RK方程算氣體流量：褚可凡、洋、周昱希指老：郭1、重述2、算混合氣體的RK方程§

根據(jù)每個(gè)物的界參數(shù)分算其在RK方程中：2、算混合氣體的RK方程，如下表1所示?！?/p>

各個(gè)物的CH4C2H6C3H8iC40.83150.09060.02790.00350.00750.00330.001190.4305.4369.8408.2425.2469.7304.1126.24.6*10^63.212.98E-054.51E-056.27E-058.06E-058.06E-051.00E-042.97E-052.68E-054.88*10^64.25*10^63.65*10^63.80*10^63.37*10^67.38*10^63.39*10^69.8718.2827.2528.9941.926.46nC4nC5CO2N20.03471.562、算混合氣體的RK方程§

根據(jù)混合物狀方程的混合RK方程合上表1中數(shù)據(jù)，算各物數(shù)據(jù)yi*ai^0.5，yi*bi如下表2所示。2、算混合氣體的RK方程CH4C2H6C3H8iC40.83150.09060.02790.00350.00750.00330.0011.4905101520.2846246420.119298860.0182717770.0403811750.0213660660.0025410330.0433332982.47913E-054.08419E-061.74868E-062.81953E-076.04504E-073.31309E-072.96817E-089.30503E-07nC4nC5CO2N20.03472、算混合氣體的RK方程§

最算得混合物的RK方程：3.算摩體#include

<stdio.h>#include

<math.h>void

equa(double

f[4],doublez[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)if(y2>0)x[0]=x[1]=x[2]=-b/(3*a);else{gs2=pow(y2,1.0/3);elseif(L>0)gs2={-pow(fabs(y2),1.0/3);p1=(-2*b+gs1+gs2)/(6*a);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);p2=sqrt(3)*(gs1-gs2)/(6*a);-gs2)/(3*a);x[0]=(-b-gs1}3.算摩體else{for(i=0;i<3;i++)if(x[i]>0)r[j++]=x[i];t=(2*A*b-3*a*B)/(2*pow(A,1.5));rad=acos(t);x[0]=(-b-for(i=0;i<j-1;i++)if(r[i]>r[i+1])2*sqrt(A)*cos(rad/3))/(3*a);x[1]=(-{s=r[i];r[i]=r[i+1];r[i+1]=s;}z[0]=r[j-1];b+sqrt(A)*(cos(rad/3)+sqrt(3)*sin(rad/3)))/(3*a);z[1]=r[0];}x[2]=(-b+sqrt(A)*(cos(rad/3)-sqrt(3)*sin(rad/3)))/(3*a);}}void

rkeos(double

tc,

double

pc,

double

w,double

t,double

p){3.算摩體double

a,b,r=8.3146,e[4],z[2];a=4.08;void

main(){doubletc=305.4,pc=4880000.0,w=0.099,t=303.15,p=101000.0;rkeos(tc,pc,w,t,p);}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-p