化工熱力學(xué)方程

化工熱力學(xué)方程第一頁，共十五頁，2022年，8月28日1、問題重述第二頁，共十五頁，2022年，8月28日根據(jù)每個物質(zhì)的臨界參數(shù)分別計算其值2、計算混合氣體的RK方程在RK方程中：第三頁，共十五頁，2022年，8月28日2、計算混合氣體的RK方程各個物質(zhì)的值，如下表1所示。組成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-05第四頁，共十五頁，2022年，8月28日2、計算混合氣體的RK方程根據(jù)混合物狀態(tài)方程的混合規(guī)則對RK方程結(jié)合上表1中數(shù)據(jù)，計算各物質(zhì)數(shù)據(jù)yi*ai^0.5，yi*bi如下表2所示。第五頁，共十五頁，2022年，8月28日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-07第六頁，共十五頁，2022年，8月28日2、計算混合氣體的RK方程最終算得則混合物的RK方程為：第七頁，共十五頁，2022年，8月28日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;

第八頁，共十五頁，2022年，8月28日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);

else

gs1=-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); }第九頁，共十五頁，2022年，8月28日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)

{

第十頁，共十五頁，2022年，8月28日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.計算摩爾體積第十一頁，共十五頁，2022年，8月28日3.計算摩爾體積計算出來50度時的摩爾體積為0.590667，單位是L/mol

第十二頁，共十五頁，2022年，8月28日3.計算摩爾體積計算出來30度時的摩爾體積為24.896001，單位是L/mol