Simulation of An Isothermal Catalytic Membrane Reactor For The Dehydrogenation of ETHYLBENZENE

Chemical and Process Engineering Research
ISSN 2224-7467 (Paper) ISSN 2225-0913 (Online)

Vol 3, 2012
www.iiste.org
Simulation of an Isothermal Catalytic Membrane Reactor for

the Dehydrogenation of ETHYLBENZENE
Akpa, Jackson Gunorubon
Department of Chemical/Petrochemical Engineering
Rivers State University of Science and Technology,
Port-Harcourt, Rivers State, Nigeria
Tel: (+234) 8036696088
E-mail: Jacksonakpa@yahoo.com
Abstract
Mathematical models for predicting the fractional conversion of ethylbenzene and yields of products in a catalytic
membrane reactor for the dehydrogenation of ethylbenzene were developed. The mathematical models developed
consisted of nonlinear simultaneous differential equations which were solved numerically using the 4 th order
Runge-kuta algorithm. Prediction by the models of fractional conversion and yields of product compare favorably
with outputs of an industrial reactor with maximum deviations of 0.175. The models were subsequently used to
simulate the effects of feed inlet temperature, feed molar ratio of steam and ethylbenzene and inlet pressure on the
reactor performance.
Keywords: Dehydrogenation of Ethylbenzene, Catalytic membrane reactor, Reactor modeling, Simulation
1. Introduction
Styrene is one of the most important monomers in the petrochemical industry due to its protective, insulative, and
synthetic ability when polymerized. The worlds production at present is approximately twenty million tons per year
(Dennis and Castor, 1992). The styrene process was developed in the nineteen thirties and was used for the
production of many different polymeric materials, the most important being poly-styrene, styrene-acrylonitrile,
styrene-butadiene latex and acrylonitrile-butadiene styrene resins (ABS). Styrene can be produced by the
dehydrogenation of ethyl-benzene in the presence of steam over iron oxide based catalyst, as a by-product in the
epoxidation of propene with ethyl-benzene hydro-peroxide and molybdenum complex base catalyst or by the
oxidative dehydrogenation of ethyl-benzene (Yee, et al., 2003); however the dehydrogenation of ethyl-benzene
accounts for over ninety percent of the worlds styrene production (Abashar, 2004). The dehydrogenation process is
an endothermic reversible reaction and can be operated industrially either adiabatically or isothermally over a fixed
bed. The demand for higher conversion of ethylbenzene, high yield and selectivity of the desired reaction products
especially styrene had led to new ingenious configuration and design of reactors for the dehydrogenation process. In
this regards, multifunctional reactors where reactions combined with separation have received much attention
(Collins an Way, (1993); Dixon, (1999) and Devoldere and froment, (1999)). Membrane reactors are one of such type
of multi-functional reactors. Recent advances in materials used at high temperatures have allowed the consideration
of membranes for integration into reactors for catalytic reactions.
To accomplish this, a membrane that is permeable to a particular reaction product, but impermeable to all other
species is placed around the reacting mixture, If reaction is equilibrium limited, the decreased activity of the species
being removed permits further conversion to occur beyond that which would be possible if no species were removed.
The permeation of gaseous component through the membrane takes place from a higher partial pressure zone to a
14

Vol 3, 2012
www.iiste.org
lower partial pressure zone for a given component. This can be achieved by either a difference in the total pressure,
or by diluting the permeate side with enough inert to lower the mole fraction of the permeating species. Hence the
membrane reactor can be used to achieve conversions greater than the original equilibrium value (Ahari et al., 2004).
The objectives of this study are to develop the mathematical model for the ethylbenzene dehydrogenation and to
investigate the effect of operating conditions on an industrial membrane reactor using the developed models.
1.1 Process Description
Figure 1a shows a schematic of a typical membrane reactor. The membrane divides the reactor into two zones
namely; a reaction zone (tube side) packed with catalyst particles where reactant (ethylbenzene) is introduced and a
permeate zone (shell side) where a non-reactive purge gas (steam) is introduced co-currently to the feed in order to
send out the permeate gas (hydrogen). Ethylbenzene is first preheated in a heat-exchanger, mixed with superheated
steam and sent into the membrane reactor. Hydrogen produced within the reactor permeates out of the tube through
the membrane due to pressure difference between the reaction zone (tube) and the permeate zone (shell). The
removal of hydrogen reduces side reactions of ethylbenzene with hydrogen. The permeated hydrogen gas is swept
out of the reactor using steam. In the reaction zone of the reactor, styrene, other aromatic hydrocarbons such as
benzene and toluene, and lean amount of hydrogen found around the reactor tube are collected at the exit or product
stream of the reactor.
Steam +
Steam
Hydrogen
Pure
R1
Ethylbenze
ne
Steam
Figure 1a: Transverse section of the catalytic packed bed membrane reactor
Permeate
Feed
N1
R3
Products
Steam +
Hydrogen
Purge
R2
P1 '
Shell side
P1
N1+1
Tube side
R1
R2
Reject
dz
L
Figure 1b: Schematic of differential volume of membrane reactor used for material balance.
15

Vol 3, 2012
www.iiste.org
2. Kinetic Model
Numerous works on the development, type, composition and activity of various catalysts used in the
dehydrogenation of ethylbenzene to styrene have been performed: Lee, (1973), Hirano, (1986), Muhler et al., (1992)
and Wu et al., (1993). Various type of reaction mechanism has been postulated to describe the kinetics of
ethylbenzene dehydrogenation, a uni-molecular Langmuir-Hinshelwood mechanism, where the reaction rate depends
on the adsorption-desorption equilibrium of ethylbenzene and styrene was earlier proposed by Carra and Forni
(1965) and Sheel and Crowe, (1969); later supported by Sheppard and Maier, (1986) and Pradeep and Elnashaie,
(2004). Others include the works of Clough and Ramirez, (1976), Wu and Liu, (1992), Wu et al., (1993), Elnashaie et
al., (1993). Elnashaie et al., (2000) and Ganji, et al., (2004) postulated that the kinetic model for the dehydrogenation
of ethylbenzene to styrene can be represented by six reversible reactions viz: one main reaction (ethylbenzene to
styrene) with five side reactions. Sheppard and Maier, (1986) observed that all side reactions were much slower in
comparison to the main reaction and are far from equilibrium at finite time where the main reaction may get
completed; therefore all the side reactions could be considered irreversible and the reverse rates excluded from the
kinetic model.
The dehydrogenation reaction occurring in the reactor leading to the production of styrene as described by Elnashaie
et al., (2000) incorporating the conclusions from the observations of Sheppard and Maier (1986) were adapted as the
reactions occurring in the reactor. The reaction paths for the dehydrogenation process with their respective reaction
rate expressions as follows:
Main reaction:
K1
C6H5C2H5
C6H5C2H3 + H2
1 = 1 (
(1)
(2)
Side reactions:
K2
C6H5C2H5 + H2
C6H5CH3 + CH4
(3)
2 = 2 ( 2 )
(4)
K3
C6H5C2H5
C6H6 + C2H4
(5)
3 = 3 ( )
(6)
K4
H2O + C2H4
CO + 2H2
)
4 = 4 (2 0.5
2 4
(8)
K5
H2O + CH4
(7)
CO + 3H2
5 = 5 (2 4 )
(9)
(10)
H2O + CO
K6
6 =
6 ( 2 )(2 )
CO2 + H2
(11)
(12)
The reaction rate constants in these equations are obtained as:

(1 1 ) = 103 (
= (
(13)
(14)
H = a + bT + cT2
(15)
16

Vol 3, 2012
www.iiste.org
a(jmol1 )
b(Jmol-1K-1 )
c( Jmol-1K-2)
122725.16
-126.27
-2.19E-03
3. Reactor Model
3.1 Model Assumption
The following conditions are imposed on the reactor and process in developing the mathematical model of the
reactor.
The mass transfer in radial direction through the membrane is negligible. An axial one-dimensional steady state
model was developed; axial diffusion of mass is negligible because the ratio of reactor length to the particle size is
large, Radial temperature gradient across the membrane is neglected (adiabatic operation), The catalyst pellet
equations are discarded because concentration gradients in catalyst pellet are neglected due to small pellet size
(Zeynali, 2010), the catalyst deactivation is minimal, thus, effectiveness factor is equal to unity, The temperature of
the fluid phase and the solid phase are the same (Babu and Gujarathi, 2010), The system operates isothermally, the
effect of pressure drop on the performance of the reactor is negligible (Wu and Liu, 1992). Hence, isobaric
conditions are assumed to prevail in the tube and shell sides, therefore P tT and PsT are constant.
Based on these assumptions a pseudo-homogeneous one-dimensional model was developed for the reacting species
by taking a component mole balance of a differential element of the reactor as shown in Figure 1b thus:
TUBE SIDE: Transportation of reacting components and products through the tube:
The mole balance for a reacting component i through reaction path j flowing through the tube is:
= 21 12
(16)
In terms of fractional conversion:

,
(21 + 12 )
(17)
The basic reactants to be considered for the six reactions are ethylbenzene and steam. These two are not permeating
through the membrane, thus the permeate term is zero and equation (17) becomes,
,
12
(18)
Equation (18) was used to write the model equation for each reactant in the six reaction scheme.
For ethylbenzene: i = 1; j = 1, 2, 3
For steam: i = 2
j = 4, 5, 6
Equation (18) is the steady state one-dimensional model equation for the membrane reactor describing the fractional
conversion of reactant i along the reactor length.
1,1
1,2
1,3
=
=
=
1
1 1
1
1 1
1
1 1
12 (1 (
))
(19)
12 (2 ( 2 ))
(20)
12 (3( ))
(21)
17

Vol 3, 2012
2,4
2,5
2,6
=
=
=
1
2 2
1
2 2
1
2 2
www.iiste.org
))
12 (4 (2 0.5
2 4
(22)
12 (5 (2 4 ))
(23)
12 (6 ( 2 )(2 ))
(24)
SHELL SIDE: Permeation of hydrogen gas through the membrane and transportation of steam plus hydrogen gas
through the shell.
Similarly since no reaction occurs within the shell of the reactor, a non-reactive component mole balance was
developed for hydrogen gas as:
2
22 2
(25)
The permeation flux of hydrogen (2 ) through a palladium-palladium alloy composite membrane had been
determined by Moustafa and Elnashaie, (2000) as:
2 =
(
)

2
2
(26)
The permeability of hydrogen (2 ) as a function of the equilibrium solubility and its diffusivity in palladium is
given by the expression of Elnashaie et al., (2000) as:
2
(27)
2 2
1
Substituting equations (26) and (27) into equation (25), the component mole balance of hydrogen through membrane
becomes:
2
2 2
(
)

2
2
(28)
Where n is a constant; DH (m2s-1) Ficks diffusion coefficient of hydrogen dissolved in palladium and C O (mole m-3)
the solubility or standard concentration of dissolved hydrogen in palladium.
These constants have been determined by Elnashaie et al., (2000) as:
= 2.3 107 (
2610
(29)
= 3.03 105 1.0358
(30)
PO is the pressure at permeate side in atm.

4. Methodology
The component partial pressures in the model equations were converted to mole concentrations of the components
using the following relationships:
Tube side
= = ( )
(31)
Shell side
18

Vol 3, 2012
www.iiste.org
= = ( )
(32)
Where: Yi = mole fraction of component i, Ni = molar concentration of component i

The molar concentrations of each component ( ) in equations (31) and (32) of the reaction path for the
hydrogenation process were obtained using expressions that are functions of the fractional conversion for each
reactant as given in Table 2.
If the components are assumed to behave as ideal gases, then the mole fractions of the reactants and products ( ) in
equations (31) and (32) can be obtained by the expressions given in Table 3 using the molar concentrations obtained
from Table 2. The component partial pressures (equations 31, 32) with the component mole fractions given in Table 3
were substituted into the model equations (equations 19 24) to obtain model equations in terms of molar
concentrations. The fourth order Runge-Kuta algorithm was adopted to develop a visual basic program to solve the
final model equations using industrial plant data of Elnashaie and Elshishini, (1994) given in Table 4.
The following boundary conditions apply:
At z = 0:
, = 0
for i = 1, 2;
= ( );
j = 1, 2,..,6
= ( )
The results of the model equations gave the fractional conversions of the reacting components along the reactor
length. Substitution of these values into the expressions in Table 2 gave the amount (moles) of reactants and yields of
products along the reactor length.
Industrial plant data was used to test the suitability of the models in predicting the conversion of ethylbemzene and
yield of the products. The effects of the flowing process variables: Feed Temperature, Feed molar ratio of steam to
ethylbenzene (2 /) and Pressure on the models developed was then investigated.
The major products of this process are styrene, benzene and toluene, therefore only the results for these products are
presented.
5. Discussion of Results
5.1 Model Validation
The results from the model equations for the yield of the products (styrene, benzene and toluene) and the conversion
of ethylbezene as predicted by the model equations in comparison with industrial results of Elnashaie and Elshishini,
(1994) are presented in Table 5.
The results showed a reasonable agreement between the model predictions and the industrial plant data. The models
predicted the conversion of ethylbenzene and the yield of styrene very accurately. The fractional conversion of
ethybenzene and yield of styrene, benzene and toluene along the reactor length as predicted by the models are shown
in Figure 2. The fractional conversion of ethylbenzene increased continuously along the reactor length, the yield of
styrene rose rapidly initially and gradually reached equilibrium towards the reactor exit, benzene and toluene
formation rate are much slower compared to that of styrene, toluene formation is slower than benzene formation
because hydrogen produced from subsequent reactions which permeates to the shell side are required for its
formation. Selectivity can be expressed as the ratio of product formation to the rate of ethylbenzene consumption.
Figure 3 is a plot of the yield of products against conversion of ethylbenzene and shows the relative selectivity of
19

Vol 3, 2012
www.iiste.org
styrene formation over benzene and toluene formation.

5.2 Model Simulation
The effects of temperature of the feed mixture at the reactor inlet, pressure at the entrance to the reactor and feed
ratio of steam to ethylbenzene (H2O/EB) on the models developed were investigated. The effects of varying each of
these parameters on the reactor performance are presented.
5.2.1 Effect of inlet Feed Temperature
Figure 4 depicts the effect of an increase in feed temperature on conversion of ethylbenzene and yield of styrene,
benzene and toluene. Dehydrogenation of ethylbenzene is an endothermic reaction and high temperatures have been
reported by Mousavi et al., (2012) as necessary for high ethylbenzene conversion because of its thermodynamics. In
the temperature range examined, ethylbenzene conversion increased continuously and linearly with corresponding
increase in the yields of styrene, benzene and toluene.
Table 6 shows the percentage increase in the conversion of ethylbenzene and the yield of products as the feed
temperature was increased. There was a decrease in the percentage increase in the conversion of ethylbenzene and
yield of products as the feed temperature was increased in the interval of 850K to 950K. At 950K the conversion of
ethylbenzene and yield of styrene were only 6.49% and 0.53% higher than at 925K compared to a 19.02% and 16.11%
increase in conversion of ethylbenzene and yield of styrene obtained by increasing feed temperature from 850K to
875K. Therefore beyond 925K the increase in products yields compared to increase in feed temperature might not
be economical.
5.2.2 Effect of Feed Molar Ratio of Steam to Ethylbenzene (H 2o/Eb)
Figure 5 shows the effect of the feed molar ratio of steam to ethylbenzene on the conversion of ethylbenzene and the
yield of the products. An increase in the steam to ethylbenzene ratio means a decrease in the feed rate (F1) keeping
the stem rate (F2) constant, figure 5 shows that when the steam to ethylbenzene ratio was increased, there was a
minimal increase in ethylbenzene conversion and a selective increase in styrene yield. This is in agreement with Le
chateliars principle which predicts that the system will attempt to oppose the change affected to the original state of
equilibrium, since ethylbenzene dehydrogenation into styrene is a reversible reaction with increasing number of
moles, a reduction in ethylbenzene feed rate will shift the equilibrium in the direction of increased concentration, that
is the equilibrium shifts to the right. A decrease in the steam to ethylbenzene ratio means an increase in the feed rate
(F1) keeping the stem rate (F2) constant. The reverse trend is expected in this case. Simulation results also predicted
a reduction in benzene and toluene yields. (Toluene formation requires hydrogen which is selectively withdrawn by
the membrane to the shell side of the reactor). Similar trends were also reported in the works of Shuka and Anand,
(2011). Figure 5 also indicates a drop in conversion of ethylbemzene and yield of styrene above a steam to
ethylbenzene ratio of 7. Minimal amount of steam is recommended in terms of plant economics as the cost of
producing steam is reduced.
5.2.3 Effect of Operating Pressure
Figure 6 shows the effect of variation of inlet pressure on the conversion of ethylbenzene and the yield of styrene,
benzene and toluene. Figure 6 predicted an initial minimal increase then a continuous gradual decrease in the
conversion and yield of ethylbenzene and styrene respectively. This trend is correct as the equilibrium conversion of
20

Vol 3, 2012
www.iiste.org
styrene from the rate equation is inversely proportional to the pressure (p). (The reverse reaction as indicated by rate
equation is proportional to p2, while the forward reaction of styrene formation is proportional to p). Figure 5 also
predicts that the yields of benzene and toluene also increased minimally with pressure. This trend Le chateliars
principle also predicts; that an increase in system pressure due to decreasing volume will favor the reaction which
involves a reduction in pressure (the reaction will shift to the side with fewer moles of gas). That is to the left,
resulting in a decrease in conversion of ethylbenzene as predicted by the simulation results.
6. Conclusion
Mathematical models that can be used to predict the performance of industrial membrane reactors used for the
dehydrogenation of ethylbenzene were successfully developed. The accuracy of the developed models was tested
using industrial plant data as inputs to solve the model equations. The results of the models (conversion of
ethylbenzene and yield of styrene, benzene and toluene) compared favorable with output values from the industrial
reactor. Therefore the models could be used to simulate industrial reactors for the dehydrogenation of ethylbenzene.
Simulation of the effects of feed inlet temperature, feed molar ratio of steam to ethylbenzene and inlet pressure was
performed. The simulation studies have shown possible new operating conditions with improved performance
(conversion and yield values). However, final selection of best operating conditions are usually based on the overall
process economics.
Nomenclature
DH
Ficks diffusivity coefficient of hydrogen, m 2/sec
Pmi
Permeability coefficient of component i, mole.m/m2.sec.atm
Total pressure of gaseous mixture, atm
Pti
Partial pressure of component i in tube, atm
Psi
Partial pressure of component i in shell, atm
Pi
Partial pressure of component i, atm
Radial coordinate, m
R1
Internal diameter of tube, m
R2
Outer diameter of tube, m
R3
Internal diameter of shell, m
rj
Rate of reaction j, mol/sec.kg of catalyst
Xtij
Fractional conversion of reactant i by reaction j in tube
spatial coordinate, m
Ni
Permeation rate of component i through membrane, mole/sec
Ji
Molar flux of component i through membrane, mole/sec.m 2
Ftio
Molar flow rate of component i in feed to tube side, mole/sec
Ki
Rate constant for reaction (i), mole/kg catalyst.sec.atm
Keb
Equilibrium constant for reaction (1), atm
Co
Solubility or Standard concentration of dissolved H2, mole/m3
Aj
Dimensionless pre-exponential factor
Ej
Activation energy of reaction j, J/mole
21

Vol 3, 2012
www.iiste.org
Greek Symbols
Thickness of membranes, m
Density of catalyst, kg/m3
References
Abashar, M. E. E., (2004), Coupling of Ethylbenzene Dehydrogenation and Benzene Hydrogenation Reactions in
Fixed Bed Catalytic Reactors. Chem. Eng. and Proc., 43, 1195-1202.
Ahari, J. S., Kakavard, M. and Farshi, A., (2004), Modeling of Radial Flow Reactors of Oxidative Rehaet Process
for Production of Styrene Monomer, Chemical Engineering & Technology, 27, 2, 139-145.
Babu, B.V. and Gujarathi, A. M., (2010), Multi-objective optimiztion of industrial styrene reactor: adiabatic and
pseudo isothermal operation, Chemical Engineering Science, 65, 2009-2026.
Carra, S. and Forni, L., (1965), "Kinetics of Catalytic Dehydrogenation of Ethylbenzene to Styrene", Industrial and
Engineering Chemical Process Design Development, 4, 281.
Clough, D. E., and Ramirez, W. F., (1976), "Mathematical Modeling and Optimization of the
Dehydrogenation of Ethylbenzene to form Styrene", AIChE J., 22, 1097-1105.
Collins, J. P. and Way, J. D., (1993), "Preparation and Characterization of a Composite Palladium Ceramic Membrane",
Industrial and Engineering Chemistry Research, 32, 3006-3013.
Dennis, H. J. and Castor, W. M., (1992), "Styrene" In Ullmann's Encyclopedia of Industrial Chemistry, Ed. by Elvers,
B., Hawkin, S., and Russey, W., New York: John Wiley and Sons, A25, 325-335.
Devoldere, K. B. and Froment, G. F., (1999), Coke formation and Gasification in the Catalytic
Dehydrogenation of Ethylbenzene, Industrial and Engineering Chemistry Research, 38, 2626-2623.
Dixon, A. G., (1999), "Innovations in Catalytic Inorganic Membrane Reactors", the Royal Society of Chemistry;
Catalysis, 14, 40-92.
Elnashaie, S. S. E. H., Abdulla, B. K., and Hughes, R., (1993), "Simulation of the Industrial Fixed Bed Catalytic
Reactor for the Dehydrogenation of Ethylbenzene to Styrene: Heterogeneous Dusty Gas Model", Industrial
and Engineering Chemistry Research, 32, 2537-2541.
Elnashaie, S. S. E. H. and Elshishini, S. S., (1994), Modeling, Simulation and Optimization of Industrial Fixed Bed
Catalytic Reactors. Gordon and Breach science Publisher, London, 364-379
Elnashaie, S. S. E. H., Moustafa, T., Alsoudani, T. and Elshishini, S.S., (2000), "Modeling and basic characteristics of
novel Integrated Dehydrogenation-hydrogenation Membrane Catalytic Reactors", Computers and Chemical
Engineering, 24, 1293-1300.
Ganji, H., Ahari, J. S., Farshi, A. and Kakavand, M., (2004), Modelling and Simulation of Benzene Alkylation
Process Reactors for Production of Ethylbenzene, Petroleum & Coal 46 (1), 55-63.
Hirano, T., (1986), "Active Phase in Potassium Promoted Iron-oxide Catalyst for Dehydrogenation of Ethylbenzene."
Journal of Applied Catalysis, 26, 65-90.
Lee, H. E., (1973), "Iron-oxide Catalysts for Dehydrogenation of Ethylbenzene in presence of Steam", Catal. Rev.
-Science Engineering, 8, 285-305.
Mousavi, S. M., Panahi, P. N., Niaei, A., Farzi, A. and Salari, D., (2012), Modeling and simulation of Styrene
Monomer Reactor: Mathrmatical and Artificial Neutral Network Model, International Journal of Scientific
& Engineering Research, 3, 3, 1-7.
22

Vol 3, 2012
www.iiste.org
Moustafa, T. M. and Elnashaie, S. S. E. H., (2000), "Simultaneous Production of Styrene and Cyclohexane in an
Integrated Reactor", Journal of Membrane Science, 178, 171-184.
Muhler, M., Schlogl, R. and Erh, G., (1992), The Nature of the Iron Oxide-Based Catalyst for Dehydrogenation of
Ethylbenzene to Styrene, Journal of Catalysis, 138, 413-444.
Pradeep, P. and Elnashaie, S. S. H., (2004), Novel Circulating Fluidized-Bed Membrane Reformer Using Carbon
Dioxide Sequestration. Ind. Eng. Chem. Res., 43, 494-501
Sheel, J. G. P., and Crowe, C.M., (1969), "Simulation and Optimization of an existing Ethylbenzene Dehydrogenation
Reactor", Canadian Journal of Chemical Engineering, 47, 183-187.
Sheppard, C. M., and Maier, E. E., (1986), "Ethylbenzene Dehydrogenation Reactor Model", Industrial and
Engineering Chemical Process Design Development, 25, 207-210.
Shuka, S. and Anand A., (2011), Multi-objective Optimization of an Industrial styrene reactor using Harmony search
Algorithm, International Journal of Computer & Communication Technology, 2, VIII, 1-7.
Wu, J. C. S., Liu, D. S. and Ko, A.N., (1993), "Dehydrogenation of Ethylbenzene over TiO 2-Fe2O3 and ZrO2Fe2O3
mixed oxide catalyst." Catalysis Letter, 20, 191 [CAS].
Wu, J. C. S. and Liu, K. T. P., (1992), "Mathematical Analysis on Catalytic Dehydrogenation of Ethylbenzene using
Ceramic Membranes", Industrial and Engineering Chemistry Research, 31, 322-327.
Yee, A. K. Y.; Ray, A. K. and Rangaiah, G. P., (2003), Multiobjective Optimization of an Industrial Styrene
Reactor, Comput. Chem. Eng., 27, 111-130.
Zeynali, M. E., (2010), Effects of Catalyst pore size on Styrene production rate, Diffusion Fundamentals.org, 13, 2,
1-8.
REACTION
KINETIC PARAMETER
PATH
Ei(J/mol)
Ai
83357
0.75
100555
0.77
188630
12.19
103996
0.12
65723
-3.21
73628
21.24
Table 1: Kinetic Parameters in Rate expression for the six reactions (Elnashaie et al., 2001)
Component
Molar concentration (mole/sec)
Ethylbenzene (C6H5C2H5)
1 (1 1,1 1,2 1,3 )
Styrene (C6H5C2H3)
1 1,1
Toluene (C6H5CH3)
1 1,2
Benzene (C6H6)
1 1,3
Hydrogen (H2)
1 (1,1 1,2 ) + 2 (22,4 + 32,5 + 2,6 )
23

Vol 3, 2012
www.iiste.org
Steam (H2O)
2 (1 2,4 2,5 2,6 )
Ethylene (C2H4)
1 1,3 2 2,4
Methane (CH4)
1 1,2 2 2,5
Carbon monoxide (CO)
2 (2,4 + 2,5 2,6 )
Carbon dioxide (CO2)
2 2,6
Total moles ( )
2 (1 + 1,1 + 1,3 ) + 2 (1 + 1.52,4 + 22,5 )
1
2
Table 2: Table of Molar Concentrations for Reactants and Products

Component
Ethylbenzene (C6H5C2H5)
Styrene (C6H5C2H3)
Toluene (C6H5CH3)
Benzene (C6H6)
Hydrogen (H2)
Steam (H2O)
Mole Fraction ( )
1 (11,1 1,2 1,3 )
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1 1,1
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1 1,2
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1 1,3
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1 (1,1 1,2 )+2 (22,4 +32,5 +2,6 )
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
2 (12,4 2,5 2,6 )
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1
Ethylene (C2H4)
Methane (CH4)
Carbon monoxide (CO)
Carbon dioxide (CO2)
1 1,3 22 2,4
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
1 1,2 2 2,5
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
2 (2,4 +2,5 2,6 )
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
2 2,6
1 (1+1,1 +1,3 )+2 (1+1.52,4 +22,5 )
Table 3: Table of Mole Fraction for Reactants and Products.

S/N
QUANTITY
NUMERICAL VALUE
1.
Reactor diameter
1.95m
2.
Reactor Length
1.7m
3.
Catalyst bulk density
4.
Catalyst particle diameter
5.
Bed void fraction
2146 kg/m3
0.0047m
0.445
24

Vol 3, 2012
6.
Catalyst composition
7.
Inlet pressure
8.
Inlet temperature
9.
Ethyl benzene in the feed
www.iiste.org
62% Fe2O3, 36% K2CO3, 2% Cr2O3

2.4 bar
922.59 K
36.87 kmol/h
Table 4: Operating Conditions for the Industrial reactor (Elnashaie and Elshishini, 1994)
PARAMETER
MODEL
INDUSTRIAL
PREDICTION
DATA
% DEVIATION
Ethylbenzene Conversion (%)
80.23
82.13
2.37
Yield of Styrene (%)
70.62
71.67
1.54
Yield of Benzene (%)
8.85
9.8
10.73
Yield of Toluene (%)
0.81
0.66
-17.5
Table 5: Comparison of Model Predictions with Industrial Data.

Temperature (K)
850
875
900
925
950
1
Yield of styrene (%)
53.30
61.89
68.37
71.92
72.30
% Increase
16.11
10.48
5.19
0.53
2
Yield of Benzene (%)
2.34
4.23
7.07
10.96
15.9
% Increase
80.74
67.17
55.17
45.01
3
Yield of Toluene (%)
0.35
0.53
0.72
0.89
1.02
% Increase
50.05
35.99
23.97
14.26
4
Conversion of Ethylbenzene (%)
55.99
66.64
76.16
83.78
89.22
% Increase
19.02
14.28
10.01
6.49
Table 6: Percentage increase in Conversion of Ethylbenzene and Yield of products with Feed Temperature.
S/No.
Parameters
25

Vol 3, 2012
www.iiste.org
0.9
CONVERSION ETHYLBENZENE
0.8
YIELD STYRENE
YIELD BENZENE
YIELD TOLUENE
CONVERSION AND YIELD
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
0.05
0.1
0.15
REACTOR LENGTH (m)
0.2
0.25
Figure 2: Fractional Conversion of Ethylbenzene and Yield of Styrene, Benzene

and Toluene along Reactor Length.
0.8
YIELD STYRENE
YIELD BENZENE
YIELD TOLUENE
0.7
0.6
YIELD
0.5
0.4
0.3
0.2
0.1
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
FRACTIONAL CONVERSION OF ETHYLBENZENE
Figure 3: Yield of Styrene, Benzene and Toluene with Fractional Conversion of

Ethylbenzene
26
0.9

Vol 3, 2012
www.iiste.org
1
0.9
YIELD STYRENE
YIELD BENZENE
YIELD TOLUENE
0.8
Conversion/Yield
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
840
860
880
900
920
940
960
Temperature (K)
Figure 4: Effect of Feed Temperature on the conversion of Ethylbenzene and yield

of Styrene, Benzene and Toluene
0.9
YIELD STYRENE
YIELD BENZENE
YIELD TOLUENE
0.8
Staem to Ethylbenzene ratio
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
Mole
10
Figure 5: Effect of Steam to Ethylbenzene ratio on Conversion of Ethylbenzene

and Yield of Styrene, Benzene and Toluene
27
12

Vol 3, 2012
www.iiste.org
0.9
YIELD STYRENE
YIELD BENZENE
YIELD TOLUENE
0.8
Conversion/Yield
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0.6
0.8
Pressure (bar)
1.2
1.4
1.6
Figure 6: Effect of Operating Pressure on the conversion of Ethylbenzene and the

yield of Styrene, Benzene and toluene
28
1.8
This academic article was published by The International Institute for Science,
Technology and Education (IISTE). The IISTE is a pioneer in the Open Access
Publishing service based in the U.S. and Europe. The aim of the institute is
Accelerating Global Knowledge Sharing.
More information about the publisher can be found in the IISTEs homepage:
http://www.iiste.org
The IISTE is currently hosting more than 30 peer-reviewed academic journals and
collaborating with academic institutions around the world. Prospective authors of
IISTE journals can find the submission instruction on the following page:
http://www.iiste.org/Journals/
The IISTE editorial team promises to the review and publish all the qualified
submissions in a fast manner. All the journals articles are available online to the
readers all over the world without financial, legal, or technical barriers other than
those inseparable from gaining access to the internet itself. Printed version of the
journals is also available upon request of readers and authors.
IISTE Knowledge Sharing Partners
EBSCO, Index Copernicus, Ulrich's Periodicals Directory, JournalTOCS, PKP Open
Archives Harvester, Bielefeld Academic Search Engine, Elektronische
Zeitschriftenbibliothek EZB, Open J-Gate, OCLC WorldCat, Universe Digtial
Library , NewJour, Google Scholar

Simulation of An Isothermal Catalytic Membrane Reactor For The Dehydrogenation of ETHYLBENZENE

Uploaded by

Copyright:

Available Formats

Simulation of An Isothermal Catalytic Membrane Reactor For The Dehydrogenation of ETHYLBENZENE

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Simulation of An Isothermal Catalytic Membrane Reactor For The Dehydrogenation of ETHYLBENZENE

Uploaded by

Copyright:

Available Formats

Chemical and Process Engineering Research

ISSN 2224-7467 (Paper) ISSN 2225-0913 (Online)