化工热力学C语言程序.docVIP

利用C++程序解决中压下两组分的汽液平衡问题： （活度系数+状态方程法） 已知压力P和液相组成xi求温度T和气相组成yi #include #include void main 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 x[] ; 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