xúc tác mcm22 trong lọc hóa dầu

MODELING OF ETHYLBENZENE PRODUCTION PROCESS WITH INTERNAL RECYCLE Peter Burian **, Marcell Horváth *, Peter Mizsey * and Alois Mészáros ** ** P.. 3-9, H-111 Budapest, Hungary E-mail: m

Slovak University of Technology in Bratislava

PROCEEDINGS

33 rd International Conference of Slovak Society of Chemical Engineering

Hotel Hutník Tatranské Matliare, Slovakia May 22 – 26, 2006

Editors: J Markoš and V Štefuca

ISBN 80-227-2409-2

MODELING OF ETHYLBENZENE PRODUCTION PROCESS WITH

INTERNAL RECYCLE

Peter Burian **, Marcell Horváth *, Peter Mizsey * and Alois Mészáros **

(**) P Burian, A Mészáros Slovak University of Technology, Faculty of Chemical and Food Technology, Radlinského 9, 81237 Bratislava, Slovak Republic

E-Mail: peter.burian@stuba.sk

(*) M Horváth, P Mizsey, Budapest University of Technology and Economics,

Faculty of Chemical Engineering, Műegyetem rkp 3-9, H-111 Budapest, Hungary

E-mail: mhorvath@mail.bme.hu

Keywords: recycle system, dynamic analysis, separation process

Introduction:

Ethylbenzene is a very important industrial substance: it is mainly used as an intermediate in the production of styrene Ethylbenzene is also a major component of mixed xylenes used as solvents in agricultural and home insecticide sprays, rubber and chemical manufacturing, and household degreasers, paints, adhesives, and rust preventives [1]

The overall dynamics of chemical processing plants with material recycle can be very different from the dynamics of the individual processing units Material recycle may dramatically alter the overall gain and time constants of the plant, and may give rise to oscillatory or instable behaviour, even when the individual processing units are stable by themselves In most cases, recycle leads to positive feedback effects For example, increasing the concentration of a chemical species in a process stream will normally increase the amount of this species in the recycle stream, and, thus, lead to a reinforcement of the original increase It refers to a self-reinforcing mechanism associated with the recycle [3,5]

The aim of this paper is to give a comprehensive picture on dynamic behaviour of processes with internal recycle The analysis is carried out applying theoretical and simulation experimental tools

on the model of the ethylbenzene production process

1 Ethylbenzene Production

Ethylbenzene is a clear, colourless liquid with a characteristic sweet, gasoline-like, aromatic odour It is a flammable and combustible liquid Its combustion may produce irritants and toxic gases It is exclusively used as an intermediate for the manufacture of styrene Currently, the primary source of ethylbenzene is the alkylation of benzene with ethylene [4]

Process Description

The ethylbenzene process consists of two subsystems (Fig 1):

• a reactor, where ethylene is reacted with benzene

• a distillation section, where the rest of benzene and polyethylbenzenes are separated from the reactor effluent to produce ethylbenzene of high purity

Fig 1 Ehtylbenzene production system

Benzene is alkylated with ethylene to form ethylbenzene in the reactor in present of AlCl3 catalyst Diethylbenzenes, triethylbenzenes and other heavier polyethylbenzenes are also formed The benzene feed is heated to reaction temperature

To maximize ethylbenzene formation, the diethylbenzene is usually transalkylated with benzene to form more ethylbenzene

The reactor effluent, a mixture of benzene, ethylbenzene and polyethylbenzenes, flows to a distillation block Untreated benzene is recovered in the first distillation column, which is called benzene column The benzene from the top of the column is recycled to the reactor The bottoms containing ethylbenzene, diethylbenzene and triethylbenzene is fed to the second column called ethylbenzene column It separates ethylbenzene from the polyethylbenzenes as a final product The ethylbenzene in a high purity (99,9%) is led away from the top of the second column The mixture

of diethylbenzene and triethylbenzene takes out from the bottom of the second column and it is fed into the third column Diethylbenzene from the third, the so-called polyethylbenzene column, is recycled back to the reactor [1,2]

2 Dynamic model

The system of the reactor and the first distillation column is studied in the first step

a Reactor

Benzene, with flow rate of 327,612 kmol/h and ethylene with flow rate of 121,802 kmol/h is

fed to the reactor (at the temperature of 180 ºC) The temperature in the reactor is 180 ºC and the

pressure is 10 bars Three exothermic reactions occur in liquid phase, in presence of AlCl3 catalyst

+ CH2=CH2

k1

k-1

∆ r H 1

+ CH2=CH2

k2

k-2

∆ r H 2

+ CH2=CH2

k3

k-3

∆ r H 3

where k i stands for reaction rate

RT E i

i k e

The conversion of the benzene in the reactor is approximately 28% The reactor effluent contains

benzene, ethylbenzene, diethylbenzene and triethylbenzene This mixture is led out from the bottom

of the reactor

The model of the process can be derived from the material and energy balance laws

Material balance:

Benzene:

[ ]

dt

Vx d x F V V

Ethylene:

[ ]

dt

Vx d x F V V

V V

V V

F2 −νEξV1 +νEBξV−1 −νEξV2 +νDEBξV−2 −νEξV3 +νTEBξV−3 = 3 E + E (2)

Ethylbenzene:

dt

Vx d x F V V

V

EB V

DEB V

E V

EB V

Eξ 1 −ν ξ −1 −ν ξ 2 +ν ξ −2 = 3 +

Diethylbenzene:

dt

Vx d x

F V V

V

V

Eξ 2 −ν ξ −2 −ν ξ 3 +ν ξ −3 = 3 +

Triethylbenzene:

+

=

−ν ξ

ξ

ν

Energy balance:

Reactor effluent:

Ak T c F T c

F

T

c

F1 p1 1+ 2 p2 2= 3 p3 + − c +ξV1 ∆r 1+ξV2 ∆r 2 +ξV3 ∆r 3+ p3 (6)

Coolant:

dt

T c dV T c F T T Ak

T

c

b First column

The first column (see Fig 2) is used to separate benzene from the other product components

The column consists of 12 trays, a reboiler and a condenser The top-product of the column is

benzene with purity higher than 99% The other three components leave the first column at the

bottom, and are fed into the second column. The feed is saturated liquid Equimolal overflow and

constant relative volatility are assumed

The model of the column is derived from the material balance law

Condenser:

dt

x H d x n x F y

n v 1= D 0+ L 0 + 0 0

1 st tray:

dt

x H d x n y n y n x

n L 0 + V 2 = V 1+ L 1+ 1 1

( )1

1 f x

y =

i th tray:

dt

x H d x n y n y n x

n L i=1 + V i+1 = V i + L i + i i

( )i

i f x

y =

7th (feed) tray:

dt

x H d x n y n x F y n x

7 7

8

( )7

7 f x

y =

k th tray:

dt

x H d x n y n y n x

n

( )k

k f x

Reboiler:

dt

x H d x F y n x

n

W W n V n

L

=

n f x

Fig 2 Distillation column

List of symbols:

F 1

F 2

F 3

ν i

V

∆ r H j

c pi

T

T C0

T C

A

– flow rate of benzene

– flow rate of ethylene

– flow rate of reactor mixture

– stoichiometric coefficient

– volume of the reactor

– reaction enthalpy

– heat capacity

– temperature of the reactor effluent

– input temperature of the coolant

– output temperature of the coolant

– heat transfer surface

n L

n V

F F

F W

F D

x i

y i

– flow rate of liquid phase in column – flow rate of vapour phase in column – fresh feed to column

(flow rate of reactor mixture) – flow rate of bottoms

– flow rate of distillate – composition in liquid phase – composition in vapour phase

A simulation schema, created in MATLAB SIMULINK environment can be seen in Fig 3

Fig 3 Simulation model in Matlab - Simulink environment

3 Simulation results

Simulation experiments, applying input variables step changes to the process, were carried out to demonstrate the dynamic behaviour of the system First, process transient behaviour with recycle, as responses to step-wise changes in inlet benzene flow rate (change of + 7%), ethylene flow rate (–10%), bottoms flow rate (+10%) and coolant temperature (+10%), are shown in Fig 4 a, b, c and d, respectively.

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

x

w

0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09

x w

0 5 10 15 20 25 30

1.4

1.6

1.8

2

2.2

2.4

2.6

2.8

3x 10

-4

time (h)

x

w

Tcoolant+10%

0 0.02 0.04 0.06 0.08 0.1 0.12

time (h)

x w

F

bottoms +10%

c) d)

Fig 4 Dinamic behaviour of the system with recycle

In order to demonstrate the effect of recycle on process dynamics, simulation experiments in form of

transient responses were carried out and compared, see Fig 5 It is transparent that the recycle modifies the

time constant as well as the static gain of the overall system

0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045

time (h)

x w

with recycle without recycle

Fig 5 Comparison of the response of the system with and without recycle

4 Process identification in linearized form

Process identification, on basis of transient dynamic responses, was carried out in order to obtain the

linearized models of the system with and without recycle, in form of linear transfer functions

The transfer functions of the process with recycle, GS, and without recycle, GM (G1, G2), take then

the following forms:

s

1 ,0122s 3 ,3511s 3 1,1806s

+ +

+

s

e

1 s 4552 , 8 23,8301s 22,3876s

+ +

+

e

1

0,912 −

Fig 6 Transfer function block diagraml of the ethylbenzeneproduction plant

The transfer function of the recycle, (23), was derived from the model configuration in Fig 6

2

1

G

G G

G G

S

S

(22)

0,0447 0,3776s

1,0648s s

0444 , 0 0,4002s ,1522s

1 1,0911s

2 3

2 3

+ +

+

+ +

+

=

R

Fig 7 shows the comparison of the transient responses of the nonlinear system and the linearized

model of system with recycle The modeling mismatch problem leads to an oscillatory behaviour of

the linearized model

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07

x w

nonlinear system linear system

5 Conclusion

Modeling problem of a realistic industrial scaled coupled reactor/column process to produce ethylbenzene is studied in this paper The model was created using Matlab – Simulink environment

Simulation experiments were carried out for various step-wise changes to demonstrate the dynamic behaviour of the system It is transparent that the recycle modifies the time constant as well as the static gain

of the overall system Process identification was carried out in order to obtain the linearized model of the ethylbenzene producing system

In further study, the linear transfer function models, obtained above, will be utilized to eliminate the recycle effect by introduction of a recycle compensator To control the overall process with robust recycle elimination, an adaptive controller with gain-scheduling performance can extend the compensator

Acknowledgements:

The authors highly acknowledge the financial support of the Scientific Grant Agency of the Slovak Republic under grants No 1/1046/04 and 1/3081/06

References

[1] Ullmann’s Encyclopedia of Industrial Chemistry, Vol A 10, 1987 VCH Verlagsgellschaft

mbH, D-6940 Weinheim, Germany

[2] Encyclopedia of Chemical Processing and Design 371.698/20, 1984 Marcel Dekker Inc.,

270 Madison Avenue, New York, 10016

[3] Alois Mészáros, Peter Mizsey, Marcell Horváth, Zsolt Fonyo, Dynamic analysis and control

of recycle processes 2004, SSCHE

[4] Toxicology and Carcinogenesis Studies of Ethylbenzene (CAS No 100-41-4) in F344/N

Rats and B6C3F1 Mice (Inhalation Studies),

http://ntp-server.niehs.nih.gov/ntpweb/printFriendly.cfm?objectid=070A7F1A-CD44-B3EE-05596CD253F268B6

[5] MORUD, J.; SKOGESTAD S.: Effects of recycle on dynamics and control of chemical processing

plants.1994, Computers Chem Engng, 18, Suppl., S529-S534.