泡点温度c语言编程,PR泡点温度计算小程序

#include

#include

main()

{

FILE *fp;

Double Pc[4]={3390000,5050000,7390000,22090000},Tc[4]={126,154,304,647

.3},Vc[4]={0.0000899,0.0000734,0.000094,0.000056},w[4]={0.040,0.021,0.225,0.344},P=100000,

M[4],qz,K1,K2,K3,K4,H,T,

x[4]={0.7,0.05,0.15,0.1}，y1=0,y2=0,y3=0,y4=0,y5,E,a,z,f[4],A,B，

q1,q2,q3,q4，aT[4],b[4],aT0=0,b0=0,aTc,aTr,Tr,k0,v;

int i,j;

double k[4][4];

for(T=100;T

{

for(i=0;i

{

for(j=i+1;j

{

y5=pow(Vc[i],1.0/3)+pow(Vc[j],1.0/3);

k[i][j]=1-8*sqrt(Vc[i]*Vc[j])/pow(y5,3);

k[j][i]=k[i][j];

}

k[i][i]=0;

}

for(i=0;i

{

k0=0.37464+1.54226*w[i]-0.26992*w[i]*w[i];

Tr=1.0*T/Tc[i];

aTr=(1+k0*(1-sqrt(Tr)))*(1+k0*(1-sqrt(Tr)));

aTc=0.45724*8.3145*8.3145*Tc[i]*Tc[i]/Pc[i];

aT[i]=aTc*aTr;

b[i]=0.07780*8.3145*Tc[i]/Pc[i];

}

b0=0;aT0=0;

for(i=0;i

{ b0=b0+x[i]*b[i];

for(j=0;j

{

aT0=aT0+x[i]*x[j]*sqrt(aT[i]*aT[j])*(1-k[i][j]);

}

}

v=8.3145*T/P;

E=1.0;

while

(E>=0.0001)

{a=v;

y1=P*v*v*v+(P*b0-8.3145*T)*v*v+(aT0-3*P*b0*b0-2*b0*8.3145*T)*v+P*b0*b0*b0+8.3145*T*b0*b0-aT0*b0;

y2=3*P*v*v+2*(P*b0-8.3145*T)*v+(aT0-3*P*b0*b0-2*b0*8.3145*T);

v=v-y1/y2;

E=fabs((v-a)/a);

}

z=P*v/(8.3145*T*1.0);

A=P*aT0/(8.3145*8.3145*T*T);

B=P*b0/(8.3145*T);

for(i=0;i

{ y4=0;

for(j=0;j

y4=y4+x[j]*sqrt(aT[i]*aT[j])*(1-k[i][j]);

y3=b[i]*(z-1)/b0-log(z-B)-A/(2*sqrt(2)*B)*(2*y4/aT0-b[i]/b0)*log((z+2.414*B)/(z-0.414*B));

M[i]=exp(y3);

}

q1=x[0]*M[0]; q2=x[1]*M[1];q3=x[2]*M[2];q4=x[3]*M[3];

qz=1-(q1+q2+q3+q4);

printf(“%fn”,T);

if(fabs(qz)

break;

}

printf(“泡点温度为%fn”,T);

}

文章知识点与官方知识档案匹配，可进一步学习相关知识C技能树首页概览115003 人正在系统学习中 相关资源：经典刻录软件尼禄Nero6.6.1.4中文安装版及序列 -其它工具类资源…