二氧化硫转化率最优化探讨.doc

二氧化硫转化率最优化 题目：SO2＋1/2O2=SO3，四段绝热反应器，级间间接换热。 1.基础数据： 混合物恒压热容Cp ＝0.2549[kcal/kg·K] －ΔH =23135[kcal/kmol] 床层空隙率ρy＝554[kg/m3] 进口SO2浓度8.0%，O2浓度9.0%，其余为氮气。 处理量131[kmolSO2/hr]，要求最终转化率98％。 2.动力学方程： 式中： 3.基本要求： (1)在T－X图上，做出平衡线，至少4条等速率线； (2)以一维拟均相平推流模型为基础，在催化剂用量最少的前提下，总的及各段的催化剂装量；进出口温度、转化率；并在T-X图上标出折线； 4.讨论： (1)要求的最终转化率从97％变化到99％对催化剂用量的影响； (2)如果有关系：YO2＋YSO2＝21％，SO2进口浓度在7－9％之间变化，对催化剂装量的影响。 计算过程： （1）由于反应速率的表达式是温度和转化率的函数，所以在入口浓度一定的条件下，当反应速率等于零的时候，不同的温度对应不同的转化率。根据这些数据就可以作出此反应的平衡线。改变反应速率的大小，分别可以得到不同反应速率下的等速率线。 平衡线和等速率线数据的计算程序如下： #includemath.h #includestdio.h #define R 1.987 #define h 0.0001 double 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); } void main() { double xso2=0.08,x,t,r0=0,y; int i; for(i=0;i5;i++){ printf(R%d=%e\n,i+1,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.001); printf(%5.2f %5.4f\n,t,x); t=t+5; } r0=r0+pow(10,-6); } } 运行上述程序后，获得的数据见下表： 温度T/k 反应速率＝r×106 0.0 1.0 2.0 3.0 4.0 693.15 0.9809 　 　 　 　 698.15 0.9797 　 　 　 　 703.15 0.9782 0.0117 　 　 　 708.15 0.9762 0.2719 　 　 　 713.15 0.9740 0.4750 　 　 　 718.15 0.9715 0.6205 0.1844 　 　 723.15 0.9686 0.7214 0.3863 0.0784 　 728.15 0.9655 0.7901 0.5373 0.2903 0.0638 733.15 0.9620 0.8365 0.6473 0.4536 0.2689 738.15 0.9582 0.8674 0.7265 0.5761 0.4285 743.15 0.9541 0.8873 0.7828 0.6669 0.5500 748.15 0.9496 0.8995 0.822