利用C++程序解决中压下两组分的汽液平衡问题： （活度系数+状态方程法）
已知压力P和液相组成xi求温度T和气相组成yi
#include
{int i,j,q;
double ps,pc[3],pc0[2],m[2],y1[100],y2[100],y[2],p[2],A[2],C[2],a[2],b[2],c[2],v[2]; double k,s,na1,na2,B0[3],B1[3],tc[3],Vc[2],vc[3],w0[2],w[3],Bm[2],Bi0,Bi1,B[3],z[3]; double x[2],T,Tc[2],dt,l[2],h[2],Z[2],ha;
cout<<\输入压力ps: \ //输入已知压力 cin>>ps;
for(i=0;i<2;i++) {cout<<\
cin>>x[i]; //输入已知组成 m[i]=1; //设初始汽相逸度系数为1 cout<<\输入组分的A,B,C: \
cin>>A[i]; cin>>Bm[i];
cin>>C[i]; //输入安托因因子 cout<<\输入组分i的Tc,pc0,Vc,Z,w0: \cin>>Tc[i]; cin>>pc0[i];
cin>>Vc[i]; cin>>Z[i];
cin>>w0[i]; //输入组分的临界参数 cout<<\认为V与温度成V=a+bT+cT^2的关系，输入组分1的a,b,c: \cin>>a[i];
cin>>b[i];
cin>>c[i]; //体积与温度的关系 }
cout<<\输入Wilson方程能量参数k,s: \cin>>k; cin>>s;
ha=25.0;y1[0]=0.0;
for(q=0;fabs(y1[q]-1)>=1e-5;q++)
{
cout<<\假设T=\ //未满足精度则重新输入温度 cin>>T;
for(i=0;i<2;i++)
{
p[i]=exp(A[i]-Bm[i]/(T+C[i])); //求取饱和蒸汽压 v[i]=a[i]+b[i]*T+c[i]*pow(T,2); //求液体体积 Bi0=0.083-0.422/(pow((T/Tc[i]),1.6));
Bi1=0.139-0.172/(pow((T/Tc[i]),4.2));
h[i]=exp(p[i]*Tc[i]*(Bi0+w0[i]*Bi1)/pc0[i]/T); //用两项的维里方程求液体逸度系数 }
na1=v[1]/v[0]*exp(-k/8.314/T);
na2=v[0]/v[1]*exp(-s/8.314/T);
l[0]=exp(-log(x[0]+na1*x[1])+x[1]*(na1/(x[0]+na1*x[1])-na2/(x[1]+na2*x[0]))); l[1]=exp(-log(x[1]+na2*x[0])-x[0]*(na1/(x[0]+na1*x[1])-na2/(x[1]+na2*x[0])));
//利用威尔森方程求活度系数 y2[0]=0.0; for(j=0;;j++)
{
for(i=0;i<2;i++)
y[i]=x[i]*l[i]*p[i]*h[i]/ps/m[i]; //得到各组分的y的值
y2[j+1]=y[1]+y[0];
if(j==0) {
tc[0]=Tc[0];tc[1]=Tc[1];tc[2]=sqrt(Tc[0]*Tc[1]);
vc[0]=Vc[0];vc[1]=Vc[1];vc[2]=pow(((pow(Vc[0],(1/3))+pow(Vc[1],(1/3)))/2.0),3); w[0]=w0[0];w[1]=w0[1];w[2]=(w0[1]+w0[0])/2.0; z[0]=Z[0];z[1]=Z[1];z[2]=(Z[1]+Z[0])/2.0; for(i=0;i<3;i++)
{ pc[i]=z[i]*8.314*tc[i]/vc[i]; B0[i]=0.083-0.422/(pow((T/tc[i]),1.6));
B1[i]=0.139-0.172/(pow((T/tc[i]),4.2)); B[i]=8.314*tc[i]*(B0[i]+B1[i]*w[i])/pc[i];
}
dt=2*B[2]-B[0]-B[1];
}
else if(fabs(y2[j+1]-y2[j])>=1e-5) {y[0]=y[0]/(y[1]+y[0]); y[1]=y[1]/(y[1]+y[0]); }
else break;
m[0]=exp(ps/8.314/T*(B[0]+pow(y[1],2.0)*dt));
m[1]=exp(ps/8.314/T*(B[1]+pow(y[0],2.0)*dt)); //求气体逸度系数 }y1[q+1]=y2[j+1]; //y1[]指各组分y之和 cout<<\输出T,y值:T=\
}
cout<<\输出最后的T,y值:T=\}
运行结果：
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库化工热力学C语言程序在线全文阅读。
相关推荐：