2024年反应工程大作业.doc

[

[反应工程學期作业]

[有关管式换热器的计算]

化工0802王風景

J

/5/20

一、绝热管式反应器

題设：

在内径為1.22m的绝热管式反应器中進行乙苯催化脱氢反应

進料温度為898K，反应压力為0.12MPa。進料為乙苯与水蒸气的混合物，其摩尔比為1:20。

反应的速率方程為

式中，p為分压，下標A、S及H分别代表乙苯、苯乙烯及氢。

反应速率常数与温度的关系為

不一样温度下的化學平衡常数值可根据下列近似式估算

反应混合物的平均比热容為，反应热等于。催化剂床层的堆密度為。试计算：

催化剂用量；

反应器的轴向温度及反应率分布。

解析：

求解基础：将管式反应器按活塞流反应器处理；将反应按拟均相模型处理；忽视副反应影响。

求解思绪：求解公式為

将式中的用的函数表述，而实际上等于單位時间單位反应器截面积上流過的乙苯摩尔数，可以求得详细数值。因此由求解公式可以得到的函数，积分函数式即可求转化率的轴向分布，從而求出反应体积。

绝热反应過程的温度与转化率成线性关系，因此可由转化率的轴向分布求得轴向温度分布。

求解過程：

反应為变摩尔反应，则

将各反应组分的分压变為转化率的函数：

（1）

（2）

（3）

（4）

将式（1）、（2）、（3）、（4）代入反应速率方程，成果為

（5）

本題為绝热操作，属变温状况，需将題給的k及K与温度的关系代入式子（5）

（6）

進料中乙苯与水蒸气的摩尔比1/20，因此

绝热温升為

由值可得反应過程的温度与转化率的关系為

（7）

而

将此值及式子（6），（7）代入求解公式可得轴向转化率分布方程

（8）

（8）

式子（8）的初值条件為Z=0，XA=0。积分式（8）即可求转化率的轴向分布，即不一样反应器高度所能到达的转化率，從而求出反应体积。這裏用龙格-库塔-吉尔法，用计算机求解。

求解的程序為：

functionzODE

clearall;clc,formatlong

XA0=0;

[z,XA]=ode45(@f,0:0.1:3.5,XA0)

plot(z,XA,k-),xlabel(z),ylabel(XA)

%------------------------------------------------------

functiondXAdz=f(z,XA)

dXAdz=4.104e6/(21+XA)*exp(-10983/(898-137*XA))*(1-XA-3.03e-7*XA^2/(21+XA)*exp(14520/(898-137*XA)))

由于绝热反应過程的温度与转化率成线性关系，由式子（7）便可算出給定转化率下的温度。

求解成果列表如下：

乙苯脱氢反应器轴向转化率及温度分布

Z/m

XA

T/K

Z/m

XA

T/K

Z/m

XA

T/K

0

0

898

1.2

0.460269

834.9431

2.4

0.576101

819.0741

0.1

0.083959

886.4977

1.3

0.475537

832.8514

2.5

0.581332

818.3575

0.2

0.150586

877.3697

1.4

0.489379

830.9551

2.6

0.586147

817.6979

0.3

0.205259

869.8796

1.5

0.501958

829.2318

2.7

0.59058

817.0905

0.4

0.251122

863.5963

1.6

0.513414

827.6623

2.8

0.594664

816.531

0.5

0.290359

858.2208

1.7

0.523865

826.2305

2.9

0.598428

816.0153

0.6

0.324376

853.5605

1.8

0.533414

824.9223

3.0

0.601899

815.5398

0.7

0.354143

849.4824

1.9

0.542151

823.7253

3.1

0.605100

815.1013

0.8

0.380419

845.8827

2.0

0.550154

822.6289

3.2

0.608053

814.6968

0.9

0.403776

842.6827

2.1

0.557492

8