化反大作业--二氧化硫转化器最优化.docx

二氧化硫转化器最优化题目要求及数据题目:SO2＋1/2O2=SO3，四段绝热反应器，级间间接换热。基础数据：混合物恒压热容Cp ＝0.2549 [kcal/kgK]－ΔH =23135 [kcal/kmol]催化剂堆密度ρb＝554 [kg/m3]进口SO2浓度8.0%，O2浓度9.0%，其余为氮气处理量131 [kmol SO2/hr]，要求最终转化率98％ 动力学方程式中（1）已知反应速率是温度和转化率的函数，在入口浓度及反应速率一定的条件下，通过改变反应温度可以得到相对应的转化率。因此，设置5个不同的反应速率，即r=0、1×10-6、2×10-6、3×10-6、4×10-6，通过编程进行计算机模拟，根据运行结果数据，可作出作T-X线。现编写计算机c++程序，如下：#includemath.h#includestdio.h#define R 1.987#define h 0.00001double r(double x, double t, double xso2){double r,r1,r2,r3, B,keff,K, Kp,Pso2,Pso3,Po2; if(t=693.15 t748.15)keff = 7.6915*pow(10,18)*exp(-76062 / (R*t)); if(t=748.15 t=873.15) keff = 1.5128*pow(10,7)*exp(-35992 / (R*t)); K = 2.3*pow(10,-8)*exp(27200/(R*t)); Kp=2.26203*pow(10,-5)*exp(11295.3 / t); Pso2=(xso2-xso2*x) / (1-xso2*x/2); Pso3=xso2*x / (1-xso2*x/2); Po2=(0.17-xso2-xso2*x/2) / (1-xso2*x/2); r1 =Po2*Pso2/Pso3; r2 =Pso3/(Pso2*sqrt(Po2)*Kp); B =48148*exp(-7355.5/t); r3 =sqrt(B+(B-1)*(1-x)/x)+sqrt(K*(1-x)/x); r =keff*K*r1*(1-r2*r2)/(r3*r3); return(r);}int main(){double xso2=0.08,x,t,r0=0,y; int i; FILE *fp;fp=fopen(data.txt,w); for(i=0;i5;i++){fprintf(fp,R=%e\n,r0);for(t=693.15;t=873.15;){x=0.01;do{x=x+h;if(x=1) break;y=fabs(pow(10,5)*r(x,t,xso2)-pow(10,5)*r0);}while(y0.0001);fprintf(fp,%5.2f %5.4f\n,t,x);t=t+5;}r0=r0+pow(10,-6);}fclose(fp);return(0);}运行程序，所得数据如表所示温度T/k反应速率r=0r=1*10-6r=2*10-6r=3*10-6r=4*10-6693.150.9838698.150.9820703.150.97990.0197708.150.97760.2789713.150.97510.4804718.150.97240.62460.1880723.150.96940.72430.38920.0809728.150.96610.79220.53950.29240.0656733.150.96250.83800.64900.45530.2705738.150.95860.86840.72780.57740.4298743.150.95440.88800.78380.66790.5510748.150.94990.90000.82310.73430.6424753.150.94500.89690.82440.74050.6536758.150.93970.89320.82450.74500.6624763.150.93400.88900.82370.74810.6692768.150.92800.88410.82190.74980.6743773.150.92160.87880.81930.75030.6778778.150.91470.87290.81590.74980.6800783.150.90740.86650.81180.74820.6810788.150.89970.85970.80700.74580.6809793.150.89150.85230.80140.74260.6798798.150.88290.84440.79530.73850