化學(xué)反應(yīng)器大作業(yè)

1、化學(xué)反應(yīng)器大作業(yè)化學(xué)反應(yīng)器理論大作業(yè)二氧化硫轉(zhuǎn)換器最優(yōu)化班 級(jí): 化研1612 學(xué) 號(hào): 62 指導(dǎo)教師: 文利雄 姓 名: 閆曉宇 1818 / 19二氧化硫轉(zhuǎn)換器最優(yōu)化反應(yīng)方程式：so21/2o2=so3 （放熱反應(yīng)）四段絕熱反應(yīng)器，級(jí)間間接換熱，常壓下反應(yīng)。1. 基礎(chǔ)數(shù)據(jù) 混合物恒壓熱容cp kcal/kgk h =23135 kcal/kmol 催化劑堆密度 b554 kg/m3 進(jìn)口so2濃度 mol%，o2濃度 mol%，其余為氮?dú)?處理量131 kmol so2/hr，要求最終轉(zhuǎn)化率98 2.動(dòng)力學(xué)方程式中： 3. 基本要求 在tx 圖上，做出平衡線；至少4條等速率線。 以一維擬

2、均相平推流模型為基礎(chǔ)，在催化劑用量最少的前提下，求總的及各段的催化劑裝量，進(jìn)出口溫度、轉(zhuǎn)化率并在tx 圖上標(biāo)出折線。 程序用c、fortran、basic語言之一編制。4.討論 要求的最終轉(zhuǎn)化率從98變化到99對(duì)催化劑用量的影響； yo2yso221，so2進(jìn)口濃度在79之間變化，對(duì)催化劑裝量的影響。一t-x圖繪制平衡線與等反應(yīng)速率線 本次大作業(yè)計(jì)算程序，使用matlab編程實(shí)現(xiàn)。 表1.平衡線所需數(shù)據(jù)溫度/k xe 使用matlab導(dǎo)出的數(shù)據(jù)作平衡線圖，如圖所示。圖1 t-x平衡線圖圖2. 完整范圍內(nèi)的tx圖圖中rso2的適宜反應(yīng)范圍是420600，但在更高的溫度范圍內(nèi)也是會(huì)有反應(yīng)的，即使反

3、應(yīng)曲線在適宜溫度范圍以外精確度低，或者反應(yīng)體系發(fā)生變化，此圖權(quán)當(dāng)?shù)玫礁用烙^、完整的圖，以期反映出整個(gè)t-x圖的趨勢(shì)。計(jì)算數(shù)據(jù)如下表：表2 等反應(yīng)速率線數(shù)據(jù)r=r=r=r=r=r=r=r=t/ xso2t/ xso2t/ xso2t/ xso2t/ xso2t/ xso2t/ xso2t/ xso2620620620620620 620 620 620610610610610610 610 610 610600600600600600 600 600 600590590590590590 590 590 590580580580580580 580 580 580570570570570570

4、 570 570 570560560560560560 560 560 560550550550550550 550 550 550540540540540540 540 540 540530530530530530 530 530 530520520520520520 520 520 520510510510510510 510 510 510500500500500500 500 500 500490490490490490 490 490 490480480480480480 480 480 480470470470470470 470 470 470460460460460460460

5、460 460450450450450440二.四段反應(yīng)器數(shù)據(jù)的計(jì)算及優(yōu)化在工業(yè)實(shí)踐中，對(duì)于任何化學(xué)反應(yīng)，要保證反應(yīng)在盡量高的反應(yīng)速率下進(jìn)行，即意味著減小反應(yīng)器的體積，減少設(shè)備投資，減少催化劑的填量，即意味著減少操作費(fèi)用，更快的反應(yīng)速率對(duì)應(yīng)著更小的反應(yīng)器體積，也意味著更大的收益。對(duì)于放熱反應(yīng)，反應(yīng)所放熱量抑制反應(yīng)向正方向進(jìn)行。隨著溫度的升高，雖然正反應(yīng)速率增大，但逆反應(yīng)速率增大更快，整個(gè)反應(yīng)便在更低的轉(zhuǎn)化率下達(dá)到平衡，放熱反應(yīng)高溫對(duì)應(yīng)著低轉(zhuǎn)化率。而且隨著反應(yīng)溫度升高，會(huì)出現(xiàn)反應(yīng)體系高溫下出現(xiàn)的問題、能耗、設(shè)備要求等其他方面，一般不在過高的溫度下進(jìn)行生產(chǎn)。從多方面考慮，但反應(yīng)放熱到一定溫度，會(huì)使

6、物料抽出與冷源換熱降低物料溫度，以達(dá)到更快的反應(yīng)速率、更大的轉(zhuǎn)化率。多級(jí)反應(yīng)器的級(jí)間換熱即是出于這個(gè)目的。對(duì)于四級(jí)反應(yīng)器，若使反應(yīng)所需催化劑用量最少。由把wcat分別對(duì)各段求x和t的微分，使其等于0；即即下一段入口溫度點(diǎn)的反應(yīng)速率和上一段出口的反應(yīng)速率相等。假設(shè)第一段入口溫度tin(1),進(jìn)行操作線計(jì)算。1. 操作線線斜率：已知入口溫度、組成，出口組成，求出口溫度：so2反應(yīng)所放出的熱量nso2*xso2 * (-h) =m*cp*txso2/ t= m*cp/ nso2* (-h) 操作線斜率即可求得。 由function t1 = txxtot (t0,x0,x1 )實(shí)現(xiàn) xso2/ t=

7、1/2.求反應(yīng)器出口轉(zhuǎn)化率xout(i)、出口溫度tout(i)和所用催化劑的量wcat(i)：已知tout（i-1）,xout（i-1）,求tin（i）當(dāng)滿足條件 （1）時(shí)，反應(yīng)器催化劑用量在此條件下取得極值，求微分、積分可以用matlab中的求偏導(dǎo)函數(shù)diff、積分函數(shù)int，亦可用其他高精度的微積分函數(shù)，但我所使用對(duì)(1/r)的ti偏導(dǎo)數(shù)計(jì)算時(shí)間長(zhǎng)，計(jì)算一萬次可花費(fèi)數(shù)分鐘（本人所使用的計(jì)算機(jī)），對(duì)偏微分的積分耗時(shí)更長(zhǎng)，一次可達(dá)數(shù)分鐘。因此選用梯形法求微分、積分也可以用一x的微小偏差，獲得的y的增量，用(（y+y）-y）/x即為此處導(dǎo)數(shù)，求積分則用梯形法，選取足夠小的步長(zhǎng)，算出每一步長(zhǎng)對(duì)應(yīng)

8、的面積，作為積分值。用梯形法計(jì)算對(duì)偏微分的積分，十萬次只需數(shù)秒，具有實(shí)踐性。求出口t見函數(shù)function daera.求催化劑用量見函數(shù) function dwcat.應(yīng)當(dāng)注意的是，xso2滿足條件（1）時(shí)，不可大于600對(duì)應(yīng)的轉(zhuǎn)化率，如果大于600時(shí)的轉(zhuǎn)化率，則此段出口轉(zhuǎn)化率為600對(duì)應(yīng)的轉(zhuǎn)化率，出口溫度即為600。已知入口溫度，出口轉(zhuǎn)化率，已知操作線斜率，即可算出出口溫度、出口時(shí)的反應(yīng)速率r.見函數(shù) function txxtot。 3求反應(yīng)器入口溫度已知上一段出口轉(zhuǎn)化率xout（i-1）、出口溫度tout(i-1)、反應(yīng)速率rout(i-1)。下一段入口轉(zhuǎn)化率xin（i）=xout（i

9、-1）,由條件可知r(i)=r(i-1). 已知xso2、r，求t。已知第二段入口反應(yīng)速率r(tin(2),xin(2) 利用function t = rxtot2( r,x)（變步長(zhǎng)搜索） 或者 function t = rxtot( r,x)（割線法）對(duì)于求第二段入口溫度時(shí)，如果使用割線法/牛頓法，因?yàn)樵?20475時(shí)，導(dǎo)數(shù)值很大，而且有重根的情況（因?yàn)?r(tin(2),xin(1)=r(tout(1)，xin(1)），有時(shí)候所得不是想要的解，甚至得到不到解。有時(shí)需手動(dòng)改變初值。比較麻煩在已知溫度在420到tout(i-1)之間情況下，考慮采用變步長(zhǎng)依次搜索法，計(jì)算50次，即可達(dá)到t的7

10、位小數(shù)的精確度，且屢試不爽。見函數(shù) function rxtot2.function dw = dwcat (xd,tin,xin) 已知入口溫度，出口轉(zhuǎn)化率，已知操作線斜率，即可算出出口溫度、出口時(shí)的反應(yīng)速率r.見函數(shù) function txxtot。給定一個(gè)第一入口溫度，求得四級(jí)反應(yīng)器的數(shù)據(jù)，xso2步長(zhǎng)如果使用matalb， （從0，計(jì)算9800次）積分計(jì)算一次需要數(shù)分鐘，如果使用梯形法計(jì)算微分積分，步長(zhǎng)選擇（計(jì)算近10萬次），計(jì)算一次四級(jí)反應(yīng)器數(shù)據(jù)只需不到兩秒。梯形法雖然簡(jiǎn)單，但計(jì)算速率高了2、3個(gè)數(shù)量級(jí)。給定第一段進(jìn)口溫度，運(yùn)行一次即可在excel中得到四段反應(yīng)器數(shù)據(jù)，例子中當(dāng)?shù)谝欢?/p>

11、進(jìn)口溫度設(shè)為時(shí)，計(jì)算四段反應(yīng)器數(shù)據(jù)如下：表3初設(shè)第一段進(jìn)口溫度設(shè)為，反應(yīng)器數(shù)據(jù)xintin/xouttout/wcat/kg總10234.對(duì)第一段入口溫度的優(yōu)化對(duì)第一段入口溫度tin（1）的優(yōu)化,現(xiàn)在溫度范圍內(nèi)得到tin（1）對(duì)催化劑總量和最終轉(zhuǎn)化率的影響趨勢(shì)，再在小范圍內(nèi)搜索最優(yōu)值。根據(jù)以上程序，使第一段入口溫度tin(1)從變化，得到第一段入口溫度tin(1)對(duì)催化劑總量和最終轉(zhuǎn)化率的影響，得到數(shù)據(jù)列表如下表4不同第一段進(jìn)口溫度反應(yīng)器數(shù)據(jù)xintin/xouttout/wcat/kg1234123412341234123412341234matlab源程序見附錄。根據(jù)上表數(shù)據(jù)作圖如下，最優(yōu)

12、值區(qū)間圖4最終轉(zhuǎn)化率和催化劑總量隨著tin(1)變化趨勢(shì)圖由上圖可以看出，在420450范圍內(nèi)，存在第一段入口溫度使得最后轉(zhuǎn)化率大于，并且在440445之間。下一步從第一段入口溫度445向440依次搜索，步長(zhǎng)為。所得結(jié)果如下表5轉(zhuǎn)化率達(dá)標(biāo)、催化劑最少用量時(shí)，反應(yīng)器數(shù)據(jù)xintin/xouttout/wcat/kg總10. 600. 2 3 4 matlab程序見附錄由上表可知滿足最終轉(zhuǎn)化率時(shí)，所有催化劑量為.根據(jù)上表數(shù)據(jù)用origin作操作折線圖如下：圖4轉(zhuǎn)化率達(dá)標(biāo)、催化劑最少用量時(shí)，反應(yīng)器操作線三.討論： 1.要求的最終轉(zhuǎn)化率從98變化到99對(duì)催化劑用量的影響；根據(jù)平衡線圖，平衡轉(zhuǎn)化率xe=

13、時(shí)，對(duì)應(yīng)的平衡溫度te=，而平衡轉(zhuǎn)化率隨著溫度升高而減小，所以在催化劑適用范圍420600內(nèi)，最終轉(zhuǎn)化率達(dá)不到。求在420600攝氏度溫度范圍內(nèi)，最大轉(zhuǎn)化率時(shí)，四段反應(yīng)器數(shù)據(jù)結(jié)果如下表所示：表6最大轉(zhuǎn)化率時(shí)，反應(yīng)器數(shù)據(jù)xintin/xouttout/wcat/kg總1 2 3 4 計(jì)算matlab源程序見附錄. yso221，so2進(jìn)口濃度在79之間變化，對(duì)催化劑裝量的影響。計(jì)算so2進(jìn)料量在保持131kmol,求滿足轉(zhuǎn)化率條件下催化劑用量最少表7催化劑用量隨著進(jìn)口so2濃度變化數(shù)據(jù)yso2, mol%tin(1)/450wcat/kg可以看出，隨著反應(yīng)器進(jìn)口yso2濃度升高。所需催化劑填裝量

14、增大。以yso2為變量，t=450，求r表7 t=450，反應(yīng)速率隨著so2濃度變化數(shù)據(jù)xso2以yso2為變量，t=500，求r表7 t=500，反應(yīng)速率隨著so2濃度變化數(shù)據(jù)xso2由上表可以看出，隨著so2入口濃度的增大，反應(yīng)速率的確降低了。符合以上結(jié)論。附錄1 t-x圖繪制平衡線與等反應(yīng)速率線matlab代碼clc,clear,yso2=;yo2=;yn2=1-yso2-yo2;epsilong=-yso2*;pso20=*yso2; xe(26)=;for i=1:26 if i<=20 xe(i)=+i*;elseif i<26 xe(i)=+(i-21)*;else

15、xe(i)=; end pso2=pso20*(1-xe(i)/(1+epsilong*xe(i);po2=(yo2-yso2*xe(i)/yso2/(1-xe(i)*pso2;pso3=pso2*xe(i)/(1-xe(i); kp=pso3/(pso2*po2;te(i)=log(kp/;endte=;texe=te;xe;texe=texe',xlswrite('',texe,'sheet2'); xso2(1)=0;for k=1:8j=2;t(1)=530+; xso2(k)=k*; rso2=txtorso2(t(1),xso2(k); tx

16、result(1,1)=rso2;t=1300+;for i=1:90 t=t-10; a=trxtox(t,rso2,; if a<1 txresult(j,1)=;txresult(j,2)=a; fprintf('rso2(%d)=%5g ,t= % ,xso2= % ,n',k,rso2,a); j=j+1; else continue; endendfprintf('n');if k=1 xlswrite('',txresult,'sheet2','c1:d90');elseif k=2 xlswr

17、ite('',txresult,'sheet2','e1:f90');elseif k=3 xlswrite('',txresult,'sheet2','g1:h90');elseif k=4 xlswrite('',txresult,'sheet2','i1:j90');elseif k=5 xlswrite('',txresult,'sheet2','k1:l90');elseif k=6 xlswr

18、ite('',txresult,'sheet2','m1:n90');elseif k=7 xlswrite('',txresult,'sheet2','o1:p90');elseif k=8 xlswrite('',txresult,'sheet2','q1:r90');end ;clear txresult;end附錄 計(jì)算四段反應(yīng)器各段數(shù)據(jù)及優(yōu)化for k=1:8tin=415+5*k;xin=0;dx=;fprintf(2,'t0=%,

19、x0=%',xin);xso2=*dx;wcat=0;for i=1:4sum=0;sumw=0;x600=(600+(i)/ +xin(i);while sum<=0 da=daera(xso2,tin(i),xin(i); dw=dwcat(xso2,tin(i),xin(i); sum=sum+da;sumw=sumw+dw; xso2=xso2+dx; if xso2>=x600 xso2=x600; break; endend xso2(i)=xso2;wcat(i)=sumw;wcat=wcat+wcat(i); tout(i)=txxtot(tin(i),xi

20、n(i), xso2(i); if tout(i)> tout(i)=; end r(i)=txtorso2(tout(i), xso2(i); fprintf(2,'xout(%d)= % ,tout(%d)= %¡æ wcat(%d)=%;n',i,xso2(i),i,tout(i),i,wcat(i); tin(i+1)=rxtot2(r(i),xso2(i); if tin(i+1)< tin(i+1)=; end xin(i+1)=xso2; fprintf(2,'tin(%d)=%¡æn',i+1,

21、tin(i+1);q2result(i,:)=i,xin(i),tin(i),xso2(i),tout(i),wcat(i);endif k=1 xlswrite('',q2result,'sheet2','s1:x4');elseif k=2 xlswrite('',q2result,'sheet2','s6:x9');elseif k=3 xlswrite('',q2result,'sheet2','s11:x14');elseif k=4 xls

22、write('',q2result,'sheet2','s16:x19');elseif k=5 xlswrite('',q2result,'sheet2','s21:x24');elseif k=6 xlswrite('',q2result,'sheet2','s26:x29');elseif k=7 xlswrite('',q2result,'sheet2','s31:x34')endend附錄 計(jì)算

23、四段反應(yīng)器各段數(shù)據(jù)及優(yōu)化for k=1:50tin=445+'t0=%,x0=%',xin);xso2=*dx;wcat=0;for i=1:4sum=0;sumw=0;x600=(600+(i)/ +xin(i);while sum<=0 da=daera(xso2,tin(i),xin(i); dw=dwcat(xso2,tin(i),xin(i); sum=sum+da;sumw=sumw+dw; xso2=xso2+dx; if xso2>=x600 xso2=x600; break; endend xso2(i)=xso2;wcat(i)=sumw;wca

24、t=wcat+wcat(i); tout(i)=txxtot(tin(i),xin(i), xso2(i); if tout(i)> tout(i)=; end r(i)=txtorso2(tout(i), xso2(i); fprintf(2,'k=%d xout(%d)= % ,tout(%d)= %¡æ wcat(%d)=%;n',k,i,xso2(i),i,tout(i),i,wcat(i); tin(i+1)=rxtot2(r(i),xso2(i); if tin(i+1)< tin(i+1)=; end xin(i+1)=xso2;

25、fprintf(2,'tin(%d)=%¡æn',i+1,tin(i+1);q2result(i,:)=i,xin(i),tin(i),xso2(i),tout(i),wcat(i);endif i=4&&xso2> xlswrite('',q2result,'sheet2','s36:x39'); break;endend附錄3 求在420600范圍內(nèi)，催化劑用量最少情況下所達(dá)到的最大轉(zhuǎn)化率maxso2=0for k=1:50tin=+ i=1:4sum=0;sumw=0;x600=(6

26、00+(i)/+xin(i);while sum<=0 da=daera(xso2,tin(i),xin(i); dw=dwcat(xso2,tin(i),xin(i); sum=sum+da;sumw=sumw+dw; xso2=xso2+dx; if xso2>=x600 xso2=x600; break; endend xso2(i)=xso2;wcat(i)=sumw;wcat=wcat+wcat(i); tout(i)=txxtot(tin(i),xin(i), xso2(i); if tout(i)> tout(i)=; end r(i)=txtorso2(tou

27、t(i), xso2(i); tin(i+1)=rxtot2(r(i),xso2(i); if tin(i+1)< tin(i+1)=; end xin(i+1)=xso2; q2result(i,:)=i,xin(i),tin(i),xso2(i),tout(i),wcat(i);end if i=4 maxso2(k+1)=xso2,wcat if maxso2(k+1)<maxso2(k) xlswrite('',q2result,'sheet2','s36:x39'); break; endendend附錄4 計(jì)算過程中所用到

28、的函數(shù)function rso2 = txtorso2( tx,xso2 )%untitled3 summary of this function goes hereyso2=;yo2=;yn2=1-yo2-yso2;epsilong=-yso2*;pso20=*yso2; r=;pso2=pso20*(1-xso2)/(1+epsilong*xso2);po2=(yo2-yso2*xso2)/yso2/(1-xso2)*pso2;pso3=pso2*xso2/(1-xso2);kp=pso3/(pso2*po2;te=log(kp/;if <470 keff=*exp(-76062/r

29、/tx);else keff=*exp(-35992/r/tx);endb=48148*exptx);k=*exp(27200/r/tx);kp=*exptx);xi=pso3/(kp*pso2*sqrt(po2);rso2=-keff*po2*k*pso2/pso3*(1-xi2)/(sqrt(b+(b-1)*pso2/pso3)+sqrt(k*pso2/pso3)2;end function rxt = rxtot( r,x)t=415+;t(1)=t;t(2)=t+5;f(1)=-txtorso2( t(1),x )+r;f(2)=-txtorso2( t(2),x )+r;j=2;wh

30、ile abs(t(j)-t(j-1)>&&j<500 t(j+1)=t(j)-f(j)*(t(j)-t(j-1)/(f(j)-f(j-1); j=j+1; f(j)=-txtorso2(t(j),x)+r;endif j>=499 fprintf('error in function rxtot,n') rxt=nan;else rxt=t(j);endend function xr = trxtox(t3, r3,x3 )%trtox summary of this function goes here x3so3(1)=x3; x3so3(2)=x3+; if isnan(trxtof(t3,r3,x3so3(1)|isnan(trxtof(t3,r3,x3so3(2) xr=nan;return; else fx(1)=trxtof(t3,r3,x3so3(1);fx(2)=trxtof(t3,r3,x3so3(2); endj=2;format long;while j<500&abs(x3so3(j)-