Determination of adsorption and kinetic parameters for methyl tert-butyl ether synthesis from tert-butyl alcohol and methanol

Article abstract of DOI:10.1006/jcat.2001.3180

Synthesis of MTBE by direct reaction between methanol and tert-butyl alcohol (TBA) on acid ion-exchange resin, Amberlyst 15, packed in a reactor was studied at 318-328 K and at different flow rates. A mathematical model based on a quasi-homogeneous kinetics was developed, assuming the reaction in the polymer phase to be homogeneous. H₂O traveled more slowly than MTBE. The reaction rate increased with increasing reaction temperature, and the conversion of the limiting reactant, TBA, was favored at high temperatures and at low flow rates of TBA. The adsorption constant of H₂O was almost 10 times greater than that of MTBE and as a result, the desired product MTBE always elutes faster than H₂O. The effects of temperature on forward reaction constant, and reaction equilibrium constant were also explored. The activation energy for the reaction was 130.1 kJ/mole. The fitted parameters obtained for enthalpy and entropy of adsorption from Arrhenius plots were consistent with thermodynamics when the standard states were changed to liquid at 1 mole to pure gas at 1 atm. The accuracy of the proposed mathematical model was further verified when it was observed that the model could predict experimental results at different feed concentrations and flow rates quite well.

Full text of DOI:10.1006/jcat.2001.3180

Journal of Catalysis 200, 209–221 (2001)
doi:10.1006/jcat.2001.3180, available online at http://www.idealibrary.com on
Determination of Adsorption and Kinetic Parameters for Methyl tert-Butyl
Ether Synthesis from tert-Butyl Alcohol and Methanol
Zhang Ziyang, K. Hidajat, and Ajay K. Ray1
Department of Chemical and Environmental Engineering, The National University of Singapore, 10 Kent Ridge Crescent, Singapore 119260
Received July 18, 2000; revised January 22, 2001; accepted January 23, 2001; published online April 24, 2001
MTBE content in U.S. gasoline) to unleaded gasoline in the
U.S.A. in 1979, rapid growth took place almost immediately
since MTBE could be manufactured and blended at the re-
finery complex. The production and import of MTBE in
the U.S.A. reached the highest levels in 1998 (3) although
presently there is a ban on the use of MTBE in California.
MTBE has highly favorable performance qualities from the
refiner’sperspective, includinghigh octane number, lowsul-
fur content, acceptable blending vapor pressure, high mis-
cibility in gasoline, moderate boiling point, and stability
in storage. Additionally, MTBE’s cost, compared to those
of other high-octane components, made it the economical
choice in the refinery marketplace. Commercially, MTBE
is synthesized today by reversible etherification reaction of
isobutene (IB, a feedstock derived from natural gas or as a
by-product ofpetroleum refining) with methanolat an oper-
ating temperature between 313 and 353 K under a pressure
high enough to maintain the reaction system in liquid phase.
Sulfonated acidic ion-exchange resin (such as Lewatit SPC
118, Amberlyst 15 or Purolite CT-115) is typically used as
catalyst. The reversible reaction is represented by
Methyl tert-butyl ether (MTBE), a high-octane blending agent
for motor gasoline, is produced by reacting directly tert-butyl alco-
hol (TBA) with methanol using Amberlyst 15 ion-exchange resin,
which acts as both catalyst and adsorbent. Experiments were car-
ried out in a fixed-bed reactor in the temperature range 318–328 K
using rectangular pulse input and measuring the elution profiles of
TBA, MTBE, and H2O. A mathematical model based on a quasi-
homogeneous kinetics was developed, which assumes the reaction
in the polymer phase to be homogeneous. The kinetic parameters,
as well as the adsorption equilibrium constants of water, TBA, and
MTBE in methanol together with their dependence on temperature,
were determined by tuning the simulation results to fit the experi-
mental data using a state-of-the-art optimization technique,the ge-
netic algorithm. The model was further validated using the tuned
adsorption and rate parameters to predict other experimental re-
sults. The kinetics reported in the present study were obtained un-
der conditions free of both external and internal mass transfer re-
sistance. The enthalpy and entropy of adsorption obtained from
Arrhenius plots were found to be consistent with thermodynamics.
c
2001 Academic Press
Key Words: MTBE; kinetic parameters; adsorption parameters;
direct synthesis; TBA; methanol.
*
)
(CH ) CCH + CH OH
(CH ) COCH
3 3 3
3 2
2
3
[1]
INTRODUCTION
1H = 37.3 kJ/mol.
Pollution from motor vehicles is responsible in industri-
alized nations for ozone-forming smog, hazardous carbon
monoxide pollution, and other toxic air pollutants. Methyl
tert-butyl ether (MTBE) is widely used as an oxygenate for
gasoline (1) not only to enhance the octane number but
also to make motor vehicle fuel burn more cleanly (replac-
ing toxic additives like lead), thereby significantly reducing
toxic tailpipe pollution (2).
MTBE was first manufactured commercially in the early
1970s in Europe by Chemische Werke Huls in West
Germany and by ANIC in Italy. After the Environmental
Protection Agency granted approval for MTBE blending
(although recently the U.S. EPA is in favor of reduction of
Numerous publications (4–15) report thermodynamics, ki-
netics, and reaction mechanism studiesofthisreaction using
ion-exchange resins as catalyst under different experimen-
tal conditions. Several attempts have also been made to
synthesize MTBE using a variety of zeolite catalysts such
as ZSM-5 (16–19) and ZCIC-10 (20). However, a problem
concerning MTBE production from isobutene is that the
source of isobutene (IB) is limited to catalytic cracking and
steam cracking fractions of petroleum refining. Another
possible source of isobutene is isobutane dehydrogenation.
Approximately 30% of MTBE production worldwide at
present is made from this feedstock.
Commercially, MTBE is also produced from tert-butyl
alcohol (TBA), a by-product of propylene oxide synthe-
sis. There are two ways to produce MTBE from TBA,
1 To whom correspondence should be addressed. Fax: +65 779 1936.
E-mail: cheakr@nus.edu.sg.
209
0021-9517/01 $35.00
c
Copyright
2001 by Academic Press
All rights of reproduction in any form reserved.

210
ZIYANG, HIDAJAT, AND RAY
an indirect and a direct method. In the indirect method device for carrying out chemical reaction and separation si-
(ARCO process), TBA is dehydrated to IB in the first re- multaneously in a fixed bed (24, 25). These are novel reac-
actor according to Eq. [2],
tors in which separation takes place at the site of chemical
reaction to improve product purities and conversions be-
yond those prescribed by thermodynamic equilibrium. The
integration of reaction and separation not only reduces the
capital and operating costs, but also results in better yields
and selectivity, reduces adsorbent requirement, and mini-
mizes problems associated with product inhibition in addi-
tion to the added advantages of performing reactions at low
pressure and temperature. The primary objective of our re-
search in studying the direct synthesis of MTBA from TBA
and methanol is to determine to what extent the moving-
bed reactor advantages of high purity and favorable equi-
librium shifts are retained in SCMCR operations. In this pa-
per, we report an experimental investigation to determine
the adsorption and kinetic parameters of this reaction.
*
)
(CH ) COH
(CH ) CCH + H O,
[2]
3 3
3 2
2
2
followed byreaction ofIB with methanolto produce MTBE
in the second reactor (10) according to Eq. [1]. In the direct
method, MTBE can be produced by reacting TBA directly
with methanol in one reactor in the presence of an acid
catalyst in which water is also formed as a by-product:
*
)
(CH ) COH + CH OH
(CH ) COCH + H O. [3]
3 3
3
3 3
3
2
The choices of solid acid catalyst for the above reaction
(Eq. [3]) are (a) a variety of heteropoly inorganic solid
acid catalysts (21, 22), and (b) hydrogen ion-exchange acid
resins, such as Amberlyst 15 (23). To the best of our knowl-
edge, there are no reported works on the direct synthesis of
MTBE from TBA using acid resin catalyst other than the
work of Matouq and Shigeo (23). Recently, another brand
of Amberlyst resin (Amberlyst 35) has been found to be
REACTION KINETICS AND ADSORPTION ISOTHERM
Numerous publications (4–15) exist on the reaction
suitable as catalyst for many reactions. But, in this study we mechanism and kinetic and thermodynamic studies of the
have used Amberlyst 15 only. reaction given by Eq. [1] using ion-exchange resins as cata-
In contrast to kinetic studies of MTBE synthesis from lyst. Ancillotti et al. (4) concluded from experimental stud-
isobutene and methanol, on which numerous papers have ies on the MTBE synthesis that the reaction rate is negative
been published during the last 20 years, direct reaction of order dependent on methanol concentration below 4 mol/l,
TBA and methanol to produce MTBE drew little atten- while it is zero-order dependent above 4 mol/l due to a shift
tion. In the open literature, Matouq and Shigeo’s work in the reaction mechanism. At low methanol concentra-
(23) seems to be the only one reported. They studied the tion, the undissociated sulfonic acid group is active, while
direct synthesis of MTBE from methanol and TBA in a the solvated protons become the catalytic agents at increas-
stirred batch reactor at different temperatures under at- ing methanol concentration as the network of hydrogen-
mospheric pressure using Amberlyst 15 ion-exchange resin bonded SO3H groups is broken down by the excess alcohol
as catalyst. Their experimental results showed that cata- present in the pores of the resin owing to selective swelling.
lyst sizes and rotation speed had no significant effects on In addition, they concluded that the olefin protonation was
reaction rates. They found that the liquid film diffusion re- the rate-limiting step and the rate is first-order dependent
sistance and particle internal diffusion resistance are negli- on the concentration of IB (9). However, no kinetic expres-
gible. They proposed a reaction mechanism in which three sion was reported.
reactions take place and determined rate expressions based
As for the kinetic rate models for the reaction (Eq. [1])
on simple kinetics. However, in developing the kinetics, the in the presence of resin catalysts, several papers have been
polymer resin catalyst was treated as a common dissolved published recently based on either the homogeneous (8, 26)
electrolyte, and the reaction was regarded as homogeneous or the heterogeneous (6, 10–12, 14, 15) model, which were
in the liquid phase. However, it is well-known that resin derived mainly from either the Langmuir–Hinshelwood
catalyst can offer a degree of adsorption selectivity that is (LH) or the Eley–Rideal (ER) mechanism. The rate equa-
unattainable in homogeneous acid or base catalysis. There- tions were described as a function of either concentrations
fore, the advantages of solid resin catalysts over the tradi- or activities. It can be generalized from the previous works
tional liquid acid or base catalysts were ignored.
that the reaction proceeds quasi-homogeneously when
The economics of many industrial chemical processes are methanol is present in a large excess, and the rate-limiting
unfavorably influenced by the equilibrium limitations of the step is the protonation of isobutene. The resin is com-
reactions involved, which results in low yield and/or selec- pletely swollen by the polar protophilic species methanol,
tivity. This can be seen in the high additional costs required and the polymer-fixed acid is dissociated completely. How-
for the separation of nonconverted reactants from the re- ever, when methanol concentration is decreased, the cata-
actor outlet product and their recycling to the reactor inlet. lytic agent is the sulfonic acid group itself and a quasi-
This can be circumvented by the use of a simulated counter- heterogeneous phenomenon begins to conquest. In most
current moving-bed chromatographic reactor (SCMCR), a heterogeneous kinetic models, the rate-limiting step is the

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
211
surface reaction either between methanol and isobutene methanol consumed and MTBE produced. Based on the
both adsorbed on the resin (LH mechanism), or between above, the overall reaction can be described by the follow-
the methanol adsorbed on the resin and the isobutene ing equation,
present in the solution. In a recent article, Tejero et al. (15)
suggested that a transition from the ER to the faster LH
mechanism takes place, depending on the methanol con-
centration in the liquid phase.
n(CH3)3COH + CH3OH
*
)
(CH ) COCH + nH O + (n 1)(CH3)2CCH2, [7]
3 3
3
2
where n (usually greater than 1) is an unknown parameter,
which indicates the ratio of reaction rates between the first
(Eq. [4]) and the second (Eq. [5]) reactions and the amount
of isobutene produced. However, if isobutene produced is
neglected, n is equal to 1, and Eq. [7] reduces to Eq. [6].
It should also be noted that although methanol is one of
the reactants, it also acts as carrier solvent and is usually
present in excess. Its concentration varies very little in the
entire reaction process and, therefore, can be regarded as
constant. The concentration of isobutene in the solid phase
can also be neglected for its low boiling point, being present
in a trace amount, and most importantly, it has very low
affinity toward the ion-exchange resin. To the best of our
knowledge, there is no reported work in the literature on
the kinetics of this reaction. Determination of adsorption
and kinetic parameters of Eq. [7] is the main objective of
this work.
Ancillotti et al. (4, 5) proposed a reaction mechanism for
the reaction (depending on the concentration of methanol
and its proportion to other reactants) between methanol
and IB (Eq. [1]) to produce MTBE catalyzed by the same
ion-exchange resin (Amberlyst 15). They concluded that
the reaction is catalyzed mainly by solvated proton when
methanol is present in large excess compared to IB. The
network of hydrogen-bonded SO3H groups can be broken
by methanol, which is usually present at high concentra-
tion in the pore liquid. The dissociated proton is distributed
evenly in the pore liquid. Likewise, the reaction given by
Eq. [7] can also be visualized as a quasi-homogeneous reac-
tion, or homogeneous in the polymer phase. The adsorbed
TBA reacts with methanol to produce MTBE, H2O, and IB
catalyzed by the dissociated proton present in the pore of
the ion-exchange resin particle, which is initially saturated
with methanol. The concentrations of MTBE and H2O in
the polymer phase depend on their adsorption equilibrium,
while the nonadsorbed isobutene desorbs into the liquid
phase as soon as it is produced, and therefore, it has no
impact on the total reaction rate. Based on the above, we
propose the following kinetic expression,
Mazzotti et al. (27) reported an interesting kinetic study
for a similar reaction system, esterification of ethanol and
acetic acid with ethyl acetate and water as products using
Amberlyst 15 ion-exchange resin as catalyst. They (27) put
forward a quasi-homogeneous model based on the compo-
sition in the polymer phase. The polymer phase concentra-
tions were related to those in the liquid phase by equating
the activities of each component in the two phases. The ac-
tivities were estimated using UNIFAC for the liquid phase,
and the extended Flory–Huggins model for the polymer
phase. The parameters were fitted to adsorption equilib-
rium experimental results of four binary systems where no
reactions were involved. Coupled with their kinetic model,
their phase equilibrium model was successful in predict-
ing the breakthrough curves of individual components in
a fixed-bed chromatographic reactor (28). However, their
phase equilibrium model is impractical for most adsorption
systems, as nonreactive binary mixtures are scarce, besides
its complexity and inconvenience in computation. Hence,
the method is not suitable for predicting the phase equilib-
rium of the reacting system and is not used in the present
study.
In order to develop an overall kinetic expression for the
reaction system studied in this work, analysis of the three
reactions proposed by Matouq and Shiego (23) is indispens-
able. For the TBA–methanol–acid ion-exchange system, the
following three reactions may take place:
*
)
TBA
IB + MeOH
TBA + MeOH
IB + H O
[4]
[5]
[6]
2
*
)
MTBE
*
)
MTBE + H O.
2
If IB produced in the first reaction (Eq. [4]) is consumed in-
stantaneously and completely according to the second reac-
tion (Eq. [5]), then the overall reaction can be described by
Eq. [6]. However, in our experimental studies, gas bubbles
(isobutene) were detected in the effluent from the packed-
bed reactor as well as when experiments conducted in a
well-stirred batch reactor using the same catalyst. Hence,
the intermediate product, isobutene, is not consumed com-
pletely, possibly because the rate of formation of IB in the
first reaction (Eq. [4]) is faster than that of the consumption
n
H2O
q
MTBE q
R = kf qTnBA
,
[8]
Ke
by the second reaction (Eq. [5]). In this case, the third re- where R is the reaction rate, qi is the concentration of com-
action (Eq. [6]) takes place in addition to reactions 1 and 2 ponent i in the polymer phase, kf is the forward reaction rate
(Eqs. [4] and [5]).However, it can be noticed from the above constant, and Ke is the reaction equilibrium constant. The
three reaction expressions that the overall TBA consumed adsorbed phase concentrations, qi, in Eq. [8] are assumed
equals the overall rate of H2O produced, and similarly for to follow a linear adsorption relation (Henry’s law) in

212
ZIYANG, HIDAJAT, AND RAY
equilibrium with those in liquid phase:
HETP for the corresponding compounds by (28)
qi = Ki Ci .
[9]
Hi u
Lu
Di =
=
.
[15]
2
2Nap,i
The linearity assumption is generally valid at sufficiently
low concentration of the reactants and products, as adopted
in our experiments. However, at higher concentrations, the
multicomponent Langmuir model (29) or the more compli-
cated Flory–Huggins and UNIFAC models (27, 30) must be
used in predicting the phase equilibrium on ion-exchange
resin.
In linear chromatography, HETP is related to the axial dis-
persion, the adsorption equilibrium, and the coefficients
of resistance to mass transfer as described by Van Deemter
etal. (31). However, determination ofthe value ofthe height
equivalent to theoretical plate, Hi, or the apparent plate
number, Nap,i, is a tedious, lengthy process. Hence, in this
work, D values are obtained by fitting the experimental
elution profiles for each component to the solution of the
above model equations (Eqs. [8]–[14]). We have made an
additional assumption that the apparent dispersion coeffi-
cient for TBA is the same as that for MTBE due to the
experimental limitation as there is no way to prevent TBA
from undergoing dehydration reaction (Eq. [4]).
The PDE in Eq. [10] together with the initial and bound-
ary conditions (Eqs. [11]–[14]), kinetic model equation
(Eq. [8]), and adsorption equilibrium (Eq. [9]) was solved
using the method of lines (32). In this technique the PDE is
first discretized in space using the finite difference method
(FDM) to convert it into a set of several coupled ODE-
IVPs. The numerical method of lines combines a numerical
method for the initial value problems (IVPs) of ordinary
differential equations (ODEs) and a numerical method for
the boundary value problems (BVPs). In this work, the
resultant stiff ODEs of the initial value kind were solved
using the subroutine, DIVPAG (which is based on Gear’s
method), in the IMSL library. The simulation results ob-
tained from the model describing the concentration profiles
of reactants and products are discussed below.
MATHEMATICAL MODEL
As discussed earlier, the direct etherification reaction
between methanol and TBA was visualized as a quasi-
homogeneous reaction in the polymer phase considering
the large excess of methanol used in the reaction mixture.
A new kinetic expression, which combines the three main
reactions into one (Eq. [7]), is put forward for convenience
in our subsequent studies being conducted at present with
a simulated countercurrent moving-bed chromatographic
reactor.
A mathematical model based on quasi-homogeneous ki-
netics was developed, which assumes the reaction in the
polymer phase to be homogeneous. The behavior of reac-
tants and products in the fixed-bed reactor was described
by a kinetic model, which assumes that the mobile and the
stationary phases are always in equilibrium, and the contri-
butions of all the nonequilibrium effects are lumped into an
apparent axial dispersion coefficient, D (28). Mass balance
equations for each component i (reactants and products)
can be written as follows:
∂Ci
∂t
l
ε
∂qi u ∂Ci
1
ε
∂2Ci
∂z2
+
+
i R = Di
.
[10]
EXPERIMENTAL DETAILS
ε
∂t
ε ∂z
ε
The initial and boundary conditions are given by
HPLC grade methanol, MTBE, and TBA have been
used in all the experimental runs. The catalyst used is a
macroporous sulfonic ion-exchange acid resin, Amberlyst
15 from Rohm & Haas Co. These are cross-linked three-
dimensional structures of polymeric material obtained by
sulfonation of a copolymer of polystyrene and divinyl
benzene. These resin are heat sensitive and lose activity
above 393 K. Macroporous resins are better catalysts com-
pared to microporous resins, particularly in nonaqueous
Ci [t = 0] = Ci0
Ci [0 < t tp]z=0 = Cf,i
Ci [t > tp]z=0 = 0
[11]
[12]
[13]
∂Ci (t)
= 0,
[14]
∂z
z=0
where i = TBA, MTBE, H2O, and u is the superficial fluid media where the latter resins do not swell appreciably. In
phase velocity, which was assumed to be constant. Equa- the MTBE process, the catalyst is poisoned by basic and
tion [10] is the overall mass balance equation of each com- cationic compounds, which neutralize the active acid groups
ponent i in a single packed column and similar to an equa- of the resin. The main properties of the ion-exchange resin
tion described elsewhere (24). The terms in Eq. [10] are (Amberlyst 15) are listed in Table 1. New fresh resin cata-
the unsteady state term (the first two terms denote the un- lysts before use were kept at 368 K in a vacuum oven
steady state term in the fluid and solid phases, respectively), overnight to get rid of any moisture.
the convective term (the third term), the reaction term
A HPLC column 25 cm long, with an overall volume of
(the fourth term), and the diffusion term (the last term). 17.35 cm3, was used in the experimental study. It was packed
The apparent dispersion coefficient, Di, is related to the with 9.66 g of dry resin. A water bath with a temperature

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
213
TABLE 1
nents in methanol (solvent) was injected into the column
packed with the Amberlyst 15 resin. No reaction occurs be-
tween H2O and methanol, but MTBE may decompose to
isobutene (IB) and methanol on the resin catalyst accord-
ing to Eq. [6]. However, from the reaction rate expression
proposed by Ali and Bhatia (10), the rate of MTBE decom-
position is negligible at a low concentration of MTBE and
when temperature of the system is low. This assumption was
further verified in this work by injecting a rectangular pulse
input of MTBE in methanol to the column for 10 min. When
a mass balance was performed between the total effluent
(area in the breakthrough profile of MTBE) and feed, it
was observed that practically no decomposition of MTBE
occurred as the discrepancy was only 1.38% . Moreover, no
gas bubble (IB) was detected in the effluent. Therefore, ad-
Typical Properties of Amberlyst 15 Dry Ion-Exchange Resin
Appearance
Hard, dry, spherical particle
Typical particle size distribution
16 mesh
% retained on US standard screens
2–5
16–20 mesh
20–30mesh
30–40 mesh
40–50 mesh
20–30
45–55
15–25
5–10
>50 mesh
1.0
Bulk density, lb/ft3
Moisture, by weight
Hydrogen ion concentration,
mequiv/g dry
38 (608 g/l)
<1%
4.7
Surface area, m2/g
Porosity, ml pore/ml bead
Average pore diameter, A
50
0.36
240
˚
sorption constants, KMTBE and KH O , can be determined by
2
fitting model equations described by Eqs. [9]–[14] to the
breakthrough curves of MTBE and H2O.
In order to obtain best-fit (tuned) values of the parame-
ters in the model, an error function, F, is defined. This com-
prised of a sum of square errors between the experimental
and model predicted values,
controller was used to maintain a constant temperature. A
binary, Series 200 LC pump from Perkin-Elmer was con-
nected to the packed bed to provide a rectangular pulse
input of width tp. Effluent from the exit of the column was
collected manually at fixed time (2 min) intervals.
Analysis for methanol, TBA, and MTBE was carried out
in a HP6890 gas chromatograph equipped with a 30-m-long
OV-1 fused silica capillary column from Ohio Valley. A vol-
umetric Karl Fischer titrator with Model 100 titration con-
n
m
X X
F(p) =
[Cij,exp Cij,m]2,
[16]
i=1 j=1
troller from Denver Instrument was used to measure the where Cij is the concentration of the ith component for
concentration of water. jth data point, p is a vector of parameters tuned, and the
Experiments were conducted at three different tempera- subscripts, exp and m, represent experimental values and
tures (318, 323, and 328 K) and flow rates. The upper limit of values predicted by the model, respectively. The state-of-
temperature (328 K) was used to maintain MTBE in the liq- the-art optimization method, the genetic algorithm (GA)
uid phase. Before feeding, the packed bed was first washed (33–35), was used to obtain the values of the model pa-
with solvent, pure methanol, for about 30 min to drive out rameters, which minimize the error function, F. GA is a
the impurities and moisture adsorbed on the resin particle. search technique developed by Holland (33) that mimics
Then LC pump connected to the feed reservoir was then the process of natural selection and natural genetics. GA
switched on to continuously provide a rectangular input is noted for its robustness and the algorithm is superior to
pulse of width 5–10 min. Afterward, the chemicals adsorbed traditional optimization algorithms in many aspects (34),
on the resin catalyst were washed off the column by solvent. and has become quite popular in recent years (35). A short
It is to be noted that for the reaction studied in this work, note on GA is added in Appendix A.
methanol, used in the concentration range of 20–24 mol/l,
acts as both solvent and reactant. Two sets of experiments
were performed. In the first set of experiments, binary mix-
tures of MTBE and H2O in methanol were used as feed
to study their adsorption and desorption equilibrium in the
column. The second set of experiments was conducted to in-
vestigate the reaction kinetics. The elution (breakthrough)
profiles of the various compounds from the exit of the col-
umn were monitored continuously.
DETERMINATION OF ADSORPTION CONSTANTS
MTBE AND KH O
K
2
In the experiments to determine KMTBE and KH O , a rect-
2
angular pulse input (of width 5 min) of binary mixtures of
MTBE and H2O at low concentration was adopted. Separa-
tion of different components takes place due to the differ-
ence in their affinitytoward the adsorbent, and each compo-
nent elutes (breakthrough) from the column at a different
time. Figure 1 shows experimental as well as model pre-
dicted breakthrough curves of MTBE and H2O at different
RESULTS AND DISCUSSIONS
The first set of experiments conducted was to measure temperatures. Experimental results show that H2O travels
breakthrough curves of MTBE and H2O when a rect- more slowly than MTBE (less strongly adsorbed) and there
angular pulse of a binary mixture of these two compo- is some band broadening. The reasons for the broadening

214
ZIYANG, HIDAJAT, AND RAY
TABLE 2
Adsorption Equilibrium Constants and Apparent Dispersion
Coefficients for MTBE and H2O
T, K
KMTBE
KH
106 DMTBE
106 DH
F
O
O
2
2
318
323
328
0.375
0.330
0.300
2.846
2.800
2.750
1.948
2.333
2.350
7.092
7.708
8.167
0.014
0.030
0.014
tabulated in Table 2 and are compared to experimental re-
sults in Fig. 1 for three different temperatures. Figure 1a
shows that model prediction is quite good when D = 0
compared to when D is set equal to zero. The minimum
values of the error function F in Eq. [16] obtained are 0.547
and 0.014 for D = 0 and D = 0 respectively. It is obvious
that model prediction is worse when the dispersion term is
neglected.The presence of a dispersion term in Eq. [17] is,
therefore, indispensable in predicting breakthrough curves
accurately. The adsorption equilibrium constant of water
was found out to be about 10 times greater than the corre-
sponding value for MTBE (Table 2). This is expected due
to the stronger polarity of H2O. Table 2 also shows that the
adsorption equilibrium constants for both MTBE and H2O
decrease with temperature while the dispersion coefficients
increase.
At low concentrations, the model results fit the exper-
imental breakthrough curve very well. However, the lin-
ear adsorption isotherm is less accurate in predicting the
elution profile when the solute is present at high con-
centration. When a mixture of [MTBE]f = 0.54 mol/l and
[H2O]f = 1.328 mol/l was injected with a pulse width of
20 min, which is comparatively larger than the exper-
imental conditions for Fig. 1a, it was observed that a
FIG. 1. Effect of temperature on breakthrough curve of the MTBE–
H2O system. Symbols, experiment; lines, model prediction (for all the fig-
ures, if not specified). Experimental conditions: flow rate Q = 1 ml/min, model using the K and D values of Table 2 cannot pre-
tp = 5 min. (a) T = 318 K, Cf, MTBE = 0.5488 mol/l, Cf, H O = 1.34 mol/l; (b)
2
dict accurately the experimentally observed breakthrough
curves, particularly for H2O. This is shown in Fig. 2. Al-
though model prediction for MTBE is quite good, it shows
T = 323 K, Cf, MTBE = 0.6114 mol/l, Cf, H O = 0.79 mol/l; (c) T = 328 K,
2
Cf, MTBE = 0.4104 mol/l, Cf, H O = 0.799 mol/l.
2
are manifold: mass transfer resistance, adsorption, and ax-
ial dispersion are among the most important factors. In this
work, for simplicity, these factors were lumped into one pa-
rameter, namely, the apparent axial dispersion coefficient,
D (see Eq. [10]). For the nonreactive case (MTBE–H2O
system), Eq. [10] simplifies to
∂Ci
∂t
1
ε
∂qi
∂t
u ∂Ci
ε ∂z
∂2Ci
∂z2
+
+
= Di
,
[17]
ε
where i = MTBE, H2O. The objective function F (Eq. [16])
was minimized by tuning four parameters, p, namely,
KMTBE , KH O , DMTBE , and DH O , to match the model pre-
2
2
diction to the experimental breakthrough curves. The tuned
adsorption and dispersion parameters for MTBE and H2O
obtained by fitting the model to experimental results are
FIG. 2. Effect of nonlinearity on the breakthrough curve of the
MTBE–H2O system. Experimental conditions: Q = 1 ml/min; tp = 20 min;
T = 318 K, Cf, MTBE = 0.54 mol/l, Cf, H O = 1.328 mol/l.
2

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
215
that due to the tailing effect prediction for H2O is poor,
and a nonlinear adsorption isotherm must be considered.
Figure 2 also reveals that at higher concentrations, the mi-
gration velocity obtained experimentally is larger than that
of the model predicted value (using linear isotherm con-
stants) as can be seen from the difference in the retention
time of the concentration peak between model and experi-
mental profiles. These effects indicate that adsorption of
H2O on the Amberlyst 15 resin follows a convex upward
isotherm, where dq/dC decreases with increasing concen-
tration.
DETERMINATION OF ADSORPTION CONSTANT, KTBA
,
AND REACTION PARAMETERS, kf, Ke, AND n
Experiments were conducted at three different temper-
atures by feeding a rectangular pulse (10 min width) in-
put of TBA of a certain concentration in methanol. Pure
methanol (solvent) was continuously pumped into the col-
umn to wash away any unconverted TBA and products
MTBE and H2O out of the packed-bed reactor. Adsorp-
tion constant KTBA and reaction kinetic parameters kf, Ke,
and n were determined from these reactive experiments, by
minimizing Eq. [16]. The model equations used for tuning
four parameters (KTBA , kf, Ke, and n) are Eqs. [10]–[14] to-
gether with the kinetic expression (Eq. [8]) and the linear
adsorption isotherm (Eq. [9]). In solving the model equa-
tions, the values of KMTBE , KH O , DMTBE and DH O obtained
2
2
earlier (given in Table 2) were used unchanged. In addition,
DTBA was assumed to be equal to DMTBE . The four param-
eters KTBA , kf, Ke and n) were tuned at each of the three
temperatures using the genetic algorithm that minimizes
the error function defined in Eq. [16]. The tuned values
of the four parameters at different temperatures are listed
in Table 3. The equilibrium conversion of TBA and the
yield, selectivity, and purity of MTBE for three different
temperatures are also reported in Table 3. The experimen-
tal and model predicted results (using the values of Tables 2
and 3) are shown in Fig. 3. Once again, the proposed model
can predict the experimental breakthrough curves reason-
ably well when D = 0.
FIG. 3. Effect of temperature on the breakthrough curves of the
TBA–H2O–MTBE system. Experimental conditions: Q
=
1 ml/min,
tp = 10 min, Cf,TBA = 2.189 mol/l. (a) T = 318 K; (b) T = 323 K; (c)
T = 328 K.
the genetic algorithm. A gene pool of 50 chromosomes was
considered and GA operations were carried out for 50 gen-
erations, after which it was observed that all 50 chromo-
somes converged to a single global optimum. The CPU time
The single objective function optimization problem in-
volving minimization of the error function was solved using
TABLE 3
Adsorption Equilibrium Constant, KTBA, and Kinetic Parameters, kf, Ke, and n
a
a
a
a
KTBA
(–)
kf,
l(n 1)/min
Ke,
mol/l
F,
mol/l
XTBA
(–)
,
YMTBE
(–)
,
SMTBE
(–)
,
PMTBE,
(–)
n)
T, K
mol(l
n, (–)
318
323
328
0.460
0.440
0.460
0.025
0.060
0.112
24.682
20.943
18.202
1.018
1.092
1.120
0.056
0.026
0.048
0.8839
0.8586
0.8524
0.8683
0.7863
0.7611
0.4955
0.4780
0.4717
0.4648
0.4402
0.4322
a Calculation is based on [TBA]o = 1.5 mol/l. XTBA = 1
[TBA]out/[TBA]o; YMTBE = [MTBE]out/[TBA]o; SMTBE = [MTBE]out/([MTBE]out +
[H2O]out); PMTBE = [MTBE]out/([MTBE]out + [H2O]out + [TBA]out).

216
ZIYANG, HIDAJAT, AND RAY
required for 50 generations was 600 min on the CRAY J916
supercomputer.
In the heterogeneous reaction sequence, mass transfer of
reactants first takes place from the bulk fluid to the exter-
nal surface of the pellet. The reactants then diffuse from
the external surface into and through the pores within the
pellet. In order to determine intrinsic kinetic parameters,
the effect of bulk diffusion resistance and pore diffusion
resistance must be estimated first.
Estimation of Bulk (External) Diffusion Resistance
Mears’s criterion (36, 41), which uses the measured rate
of reaction, helps to determine if external diffusion is limit-
ing the reaction. Mears’s criterion states that external mass
transfer can be neglected if
( rA0 b)Rn
< 0.15,
[18]
kcCA
where ( rA0 b) is the measured rate of reaction, in kilo-
moles per cubic meter peer second, R is the average radius
of catalyst particles (3.75 10 4 m), n is the order of reac-
tion, CA is the bulk concentration of the reactant (TBA), in
kilomoles per cubic meter, and kc is the mass transfer coef-
ficient, in meters per second. The measured rate of reaction
FIG. 5. Effect of feed concentration on the breakthrough curves of
the TBA–H2O–MTBE system. Experimental conditions: T = 323 K, Q =
1 ml/min, tp = 10 min. (a) Cf,TBA = 2.766 mol/l; (b) Cf,TBA = 1.398 mol/l.
can be determined from Eq. [8] as 0.381 kmol/m3-catalyst
bed/s, and the mass transfer coefficient, kc, can be estimated
from the Dwidevi–Upadhyay mass transfer correlation (42)
as 3.49 10 5 m/s, which results in Mears’s criterion param-
eter value as 2.98 10 3, which is less than 0.15. Therefore,
bulk diffusion can be neglected. It should also be noted that
estimated kinetic parameters as reported in Table 3 can pre-
dict the breakthrough curves very well when experiments
were performed at other flow rates, confirming further that
external mass transfer resistance is negligible. Detailed cal-
culation is shown in Appendix B.
Estimation of Pore Diffusion Resistance
The Weisz–Prater criterion is used to determine whether
internal mass transfer is limiting the reaction.The Weisz–
Prater criterion (36) states that internal pore diffusion is
negligible if
[ rA0
]
b obs L2
< 1,
[19]
DeCAs
where CAs is the concentration of TBA on the resin surface,
De is the effective diffusivity of TBA in methanol, and L for
a spherical pellet is given by R/3, where R is the average
radius of the resin particles. The Weisz–Prater parameter
FIG. 4. Effect of feed concentration on the breakthrough curves of
the TBA–H2O–MTBE system. Experimental conditions: T = 318 K, Q =
1 ml/min, tp = 10 min, (a) Cf,TBA = 1.7117 mol/l; (b) Cf,TBA = 1.203 mol/l.

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
217
FIG. 7. Temperature dependence of adsorption constants.
tration changes inside the packed-bed reactor. As a result,
the local voidage will change, which subsequently changes
the interfacial flow velocity. However, considering the low
concentration range adopted in this work, the concentra-
tion change is relatively small, and the error caused is in-
significant.
The dependence of adsorption constants, KMTBE, KH O
,
2
and KTBA , on temperature can be determined from the
equation,
1Hi
FIG. 6. Effect of flow rate on the breakthrough curves of the
Ki = Ki,o exp
,
i = MTBE, H2O, TBA, [20]
TBA–H2O–MTBE system. Experimental conditions: T = 323 K, Cf, TBA
=
RT
2.189 mol/l, tp = 10 min. (a) Q = 2 ml/min; (b) Q = 0.5 ml/min.
where 1H is the heat of adsorption. The heat of adsorption
is exothermic (1Hi < 0), and therefore, Ki decreases with
increase of temperature. The values of Ki,o and 1Hi were
obtained for each component (MTBE, H2O, and TBA) by
least-squares fitting as shown in Fig. 7. It was observed that
the effect of temperature is not significant in the temper-
ature range under study. The values of Ki,o and 1Hi for
TBA, MTBE, and H2O are given in Table 4. In order to
verify the thermodynamic consistency, enthalpy, and en-
tropy of adsorption of the components, the standard states
were changed to 1 atm in the gas phase. The calculations
made were identical to those reported by Boudart (37) and
Singh and Vannice (38). The standard states for 1H and 1S
were changed from liquid at 1 mol to pure gas at 1 atm by
following the path shown by Singh and Vannice (38). The
enthalpyand entropyofadsorption ofMTBE and methanol
were found to be negative as required by thermodynamics
calculated (see Appendix B) for the given system yields a
value of 0.012, signifying that internal pore diffusion is also
negligible.
The tuned adsorption and kinetic parameters of the three
components at three different temperatures were used next
to verify the validity of the model by checking whether it
can correctly predict experimental breakthrough curves of
the three components when experiments were performed
at different flow rates and feed concentrations. Figure 4
shows results obtained when feed concentrations were var-
ied keeping the temperature at 318 K and flow rate at
1 ml/min, the same values as in Fig. 3a. Figure 5 shows
a similar effect of feed concentrations at T = 323 K while
Fig. 6 compares model predicted values to experimental re-
sults at different flow rates. These figures show that when
adsorption and kinetic parameter values of Tables 2 and 3
are used the model can predict quite well the breakthrough
curves of all three components, particularly for MTBE and
TBA. From the figures it is also apparent that when reaction
occurs broadening of the elution peaks is less significant as
the overall rate is controlled by kinetics (rate-determining
step) rather than by axial dispersion. The main reason for
the small error is believed to be that we have neglected the
varied degree to which the resin gets swollen when concen-
TABLE 4
Pre-exponential Factors and Heat of Adsorption
Ki,o
1Hi, J/mol
TBA
0.446
MTBE
H2O
TBA
39.9
MTBE
H2O
4
2.461 10
0.924
19363.3
2974.7

218
ZIYANG, HIDAJAT, AND RAY
TABLE 5
TABLE 6
Comparison of Thermodynamic Parameters for Liquid
and Vapor Phase Studies at Standard State of 1 atm
Reaction Kinetic Parameters
kf,o
,
Ef,
1SRo ,
1HRo ,
kJ/mol
mol(1
l(n 1)/min
kJ/mol
J/mol/K
n)
Liquid phase
Vapor phase
6.028 1019
130.147
56.444
26.418
1Ho,
1So,
1Ho,
1So,
kJ/mol
kJ/mol/K
kJ/mol
kJ/mol/K
KMTBE
KTBA
KH
18.6
0.03
3
59
16
38
20.3
0.04
1.7
48
8
29
The dependence of the forward reaction rate constant
Kf on temperature was determined from the Arrhenius
equation,
O
2
Ef
kf = kf,o exp
.
[21]
RT
and are given in Table 5. It is interesting to note the common
trends between the liquid and vapor phase values reported
in Table 5 and the fitted values obtained as reported in
Table 4. Some of the values do differ slightly between the
liquid and vapor phases but they all satisfy thermodynamic
constraints. It should also be noted that this is the only study
where synthesis of MTBE is done from TBA and methanol.
All other MTBE synthesis studies reported in the literature
are from isobutene (IB) and methanol and Refs. (39, 40)
did not report enthalpy and entropy of adsorption in their
reported work. In addition, in these works Lewatit SPC
118 BG was used as adsorbent and catalyst instead of the
Amberlyst 15 used in this work.
The values of activation energy, Ef, and pre-exponential
factor, kf,o were obtained from the slope of the trend line
(shown in Fig. 8) and the values are given in Table 6. The
reaction equilibrium constant, Ke, is related by
kf
kf,o
(Ef Eb)
1GoR
ke =
=
exp
= exp
kb
kb,o
RT
RT
1SRo
1HRo
= exp
exp
.
[22]
R
RT
The values of 1SRo [ R ln( kf,o )] and 1HRo [ (Ef Eb)] were
kb,o
obtained by least-squares fitting and the values are given
in Table 6. The effect of temperature on kinetic parameter
n, which represents the ratio of the production rate of H2O
to that of MTBE, is shown in Table 3. n increases with in-
creasing temperature, indicating that more H2O and IB is
produced than MTBE at higher reaction temperature. The
1HRo and 1SRo values obtained in this work (see Table 6) dif-
fer from the value reported by Izqulerdo et al. (39, 40). They
reported 1Ho = 37.1 kJ/mol and 1So = 79.3 J/mol/K
compared to our valuesof26.42kJ/moland 56.44J/mol/K,
respectively. However, it should be noted that they con-
sidered MTBE synthesis reaction from isobutene (IB) and
methanol, whereas in this work we considered synthesis of
MTBE by a different route, namely, direct reaction of TBA
with methanol. Moreover, they reported large discrepan-
cies of values of 1HRo at different temperatures.
CONCLUSIONS
Synthesis of MTBE by direct reaction between methanol
and TBA on acid ion-exchange resin, Amberlyst 15, packed
in a reactor is studied in the temperature range of 318–
328 K and at different flow rates. Methanol also acts as sol-
vent and is used in excess. The reaction was visualized as a
quasi-homogeneous reaction in the polymer phase, and is
assumed to be in (adsorption) equilibrium with that of the
liquid phase. The adsorption equilibrium was described by
a linear isotherm as only the low concentration range was
FIG. 8. Effects of temperature on (a) kf and (b) Ke.

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
219
considered. A new kinetic expression was established. The
Three common operators are used in GA to obtain an
breakthrough curves of the reactants and products were improved (next) generation of chromosomes. These are re-
experimentally measured at different temperatures, feed ferred to as reproduction, crossover, and mutation. Repro-
concentrations, and flow rates. Experimental results show duction is the generation of the mating pool, where the
that H2O travels more slowly than MTBE (less strongly chromosomes are copied probabilistically on the basis of
adsorbed), the reaction rate increases with an increase in their fitness values. However, no new strings are formed
reaction temperature, and the conversion of the limiting in the reproduction phase. New strings are created using
reactant, TBA, is favored at high temperatures and at low the crossover operator by exchanging information among
flow rates of TBA. The adsorption equilibrium constants pairs of strings in the mating pool. A pair of daughter chro-
and reaction kinetic parameters together with their depen- mosomes are produced by selecting a crossover site (cho-
dence on temperature were determined by least-squares sen randomly) and exchanging the two parts of the pair
fitting of the proposed model to the experimental results of parent chromosomes (selected randomly from the mat-
using a state-of-the-art optimization technique, the genetic ing pool). The effect of crossover may be detrimental or
algorithm. Pure kinetic parameters were obtained as it was beneficial. It is hoped that the daughter strings are supe-
found that under the experimental conditions used, both rior. If they are worse than the parent chromosomes, they
external and internal mass resistances are negligible. The will slowly die a natural death over the next few genera-
adsorption constant of H2O was found to be almost 10 tions (the Darwinian principle at work). In order to pre-
times greater than that of MTBE and as a result, the de- serve some of the good strings that are already present
sired product MTBE always elutes faster than H2O. The in the mating pool, not all strings in the pool are used in
effects of temperature on forward reaction constant, kf, and crossover. A crossover probability, Pcross, is used, where
reaction equilibrium constant, Ke, were also explored. The only 100Pcross percent of the strings in the mating pool are
activation energy for this reaction was determined to be involved in crossover while the rest continue unchanged to
130.1 kJ/mol. The fitted parameters obtained for enthalpy the next generation. After a crossover is performed, muta-
and entropy of adsorption from Arrhenius plots were found tion takes place. The mutation operator changes a binary
to be consistent with thermodynamics when the standard number at any location in a chromosome from a 1 to a 0
states were changed to liquid at 1 mol to pure gas at 1 atm. and vice versa, with a small probability, Pmute. Mutation is
The accuracy of the proposed mathematical model was fur- needed to create a point in the neighborhood of the cur-
ther verified when it was observed that the model could rent point, thereby achieving a local search around the cur-
predict experimental results at different feed concentra- rent solution and maintaining diversity in the population.
tions and flow rates quite well.
The entire process is repeated until some termination cri-
terion is met (the specified maximum number of genera-
tions is attained, or the improvements in the values of the
objective functions become lower than a specified toler-
ance). A more elaborate description of GA is available in
the works Holland (33), Goldberg (34), and Bhaskar et al.
(35).
APPENDIX A: A NOTE ON GENETIC ALGORITHM
GA is a search technique developed by Holland (33) that
mimics the process of natural selection and natural genet-
ics. In this algorithm, a set of decision variables are first
coded in the form of a set of randomly generated binary
numbers (0 and 1), called strings or chromosomes, thereby
creating a “population (gene pool)” of such binary strings.
Each chromosome is then mapped into a set of real values
of the decision variables, using the upper and lower bounds
of each of these. A model of the process is then used to pro-
vide values of the objective function for each chromosome.
The value of the objective function of any chromosome re-
flects its “fitness.” The Darwinian principle of “survival of
the fittest” is used to generate a new and improved gene
pool (new generation). This is done by preparing a “mating
pool,” comprising copies of chromosomes, the number of
copies of any chromosome being proportional to its fitness
(Darwin’s principle). Pairs of chromosomes are then se-
lected randomly, and pairs of daughter chromosomes gen-
erated using operations similar to those in genetic repro-
duction. The gene pool evolves, with the fitness improving
over the generations.
APPENDIX B. ESTIMATION OF EXTERNAL AND
INTERNAL MASS TRANSFER RESISTANCE
Estimation of External Mass Transfer
Mears’s criterion (41) is given by
( rA0 b)Rn
< 0.15,
[B1]
kcCA
where ( rA0 b) is the measured rate of the reaction based
on the bulk volume of packed bed, in kilomoles per cubic
meter of catalyst bed per second. The initial reaction rate
is given by
( rA000)obs = kfqTnB = (0.06)[(0.44)(1.5)]1.092[103/60]
= 0.635 mol/m3 catalyst/s
= 0.381 mol/(m3 catalyst bed)/s.

220
ZIYANG, HIDAJAT, AND RAY
The mass transfer coefficient was calculated from the L = R/3, for a spherical pellet where R is the average
Dwidevi and Upadhyay correlation (42) which is valid for radius of the resin particles. Substituting in Eq. [B3], we
liquids (Re > 0.01) in fixed beds and is given by
have Weisz–Prater parameter = 0.012, signifying that inter-
nal pore diffusion is negligible.
εb Sh
0.765
0.365
εb JD
=
=
+
,
[B2]
Sc1/3 Re
Re0.82
Re0.386
NOTATION
where
C
D
E
F
1G
H
1H
IB
k
liquid phase concentration, mol/l
apparent axial dispersion coefficient, m2/s
activation energy, J/mol
error (objective) function for optimization
change in Gibbs free energy of reaction, J/mol
height equivalent to a theoretical plate, m
heat of adsorption, J/mol
U
MeOH dp
(0.00024 m/s)(790 kg/m3)(0.00075 m)
(0.0006 kg/m s)
Re =
=
MeOH
= 0.237 > 0.01,
where U = superficial velocity of fluid through bed =
0.00024 m/s, and dp = average catalyst particle diameter =
0.00075 m.
isobutene
reaction rate constant
K
L
equilibrium constant
length of the packed bed reactor, m
Schmidt number
(0.0006 kg/m s)
MeOH methanol
MTBE methyl tert-butyl ether
Sc =
=
= 393.52.
(790 kg/m3)(1.93 10 9 m2/s)
DAB
n
N
P
q
R
1S
t
moles of TBA reacted per mole of methanol
equivalent plate number
Substituting into Eq. [B2], we have Sh = kcdp/DAB
=
13.578. DAB can be obtained from the Wilke–Chang equa-
vector of parameters tuned
concentration in the polymer phase, mol/l
reaction rate, mol/min/l, universal gas constant
entropy change of reaction, J/mol/K
time, s
tion (36) as
(117.3 10 18)(ϕMB)0.5T
DAB
=
m2/s,
0.6
A
T
TBA
u
z
temperature, K
tert-butyl alcohol
superficial fluid phase flow rate, m/s
axial coordinate, m
where ϕ = 1.9 (methanol), MB = 32 kg/kmol, T = 323 K,
= 0.0006 kg/ms (methanol), A = molar volume at normal
boiling point of TBA = 0.285(Vc)1.048 = 0.10263 m3/kmol.
9
Substituting, we have DAB = 1.93 10 m2/s. Therefore,
5
the mass transfer coefficient, kc = 3.49 10 m/s. Substi-
Greek Letters
tuting in Mears’s criterion parameter, we have
ε
void fraction
(rA0 b)Rn (0.381)(0.000375)(1.092)
stoichiometric coefficient of component
=
= 2.98 10 3 < 0.15.
kcCA
(3.49 10 5)(1500)
Subscripts/Superscripts
Therefore, external mass transfer (bulk diffusion) can be
neglected.
o
ap
b
initial, pre-exponential, standard
apparent
backward
equilibrium
Estimation of Internal Mass Transfer
e
exp experiment
The Wiesz–Prater criterion is used to determine whether
internal mass transfer is limiting the reaction.The Wiesz–
Prater criterion (36) states that internal pore diffusion is
negligible if
f
feed, forward
i
component i
j
data point
m
n
p
model, number of data points
exponent, number of components
width of rectangular pulse
[rA0
]
b obs L2
< 1,
[B3]
DeCAs
where CAs is the concentration of TBA on the resin sur-
face. Since bulk diffusion is negligible, CAs can be taken as
CA = 1500 mol/m3. De is the effective diffusivity of TBA
in methanol, and is given by [ε/ ] DAB, where ε is particle
porosity = 0.36 (Table 1), t is tortuosity factor = 1.3 (43).
REFERENCES
1. Peeples, J. E., Fuel Reform. 1(1), 27 (1991).
2. Hadder, G. R., Energy 21, 118 (1992).
3. Swain, E. J., Oil Gas J. 97(24), 99 (1999).

ADSORPTION AND KINETIC PARAMETERS FOR MTBE SYNTHESIS
221
4. Ancillotti, F., Massi Mauri, M., and Pescarollo, E., J. Catal. 46, 49
(1977).
5. Ancillotti, F., Massi Mauri, M., Pescarollo, E., and Romagnoni, L.,
J. Mol. Catal. 4, 37 (1978).
6. Gicquel, A., and Torck, B., J. Catal. 83, 9 (1983).
7. Voloch, M., Ladisch, M. R., and Tsao, G. T., React. Polym. 4, 91 (1986).
8. Subramanium, C., and Bhatia, S., Can. J. Chem. Eng. 65, 613 (1987).
9. Al-Jarallah, A. M., Siddiqui, M. A. B., and Lee, A. K. K., Can. J. Chem.
Eng. 66, 802 (1988).
10. Ali, A., and Bhatia, S., Chem. Eng. J. 44, 97 (1990).
11. Rehfinger, A., and Hoffmann, U., Chem. Eng. Sci. 45, 1605 (1990).
12. Parra, D., Tejero, J., Cunill, F., Iborra, M., and Izquierdo, J. F., Chem.
Eng. Sci. 49, 4563 (1994).
13. Izquierdo, J. F., Cunill, F., Vila, M., Iborra, M., and Tejero, J., Ind. Eng.
Chem. Res. 33, 2830 (1994).
14. Caetano, N. S., Loureiro, J. M., and Rodrigues, A. E., Chem. Eng. Sci.
49, 4589 (1994).
15. Tejero, J., Cunill, F., Izquierdo, J. F., Iborra, M., Fite´, C., and Parra, D.,
Appl. Catal. A-Gen. 134, 21 (1996).
16. Chu, P., and Kuhl, G. H., Ind. Eng. Chem. Res. 26, 365 (1987).
17. Kogelbauer, A., Nikolopoulos, A. A., Goodwin, Jr., J. G., and
Marcelin, G., J. Catal. 152, 122 (1995).
25. Ray, A. K., and Carr, R. W., Chem. Eng. Sci. 50(14), 2195 (1995).
26. Colombo, F., Corl, L., Dalloro, L., and Delogu, P., Ind. Eng. Chem.
Fundam. 22(2), 219 (1983).
27. Mazzotti, M., Kruglov, A., Neri, B., Gelosa, D., and Morbidelli, M.,
Chem. Eng. Sci. 51, 1827 (1996).
28. Guiochon, G., Shirazi, S. G., and Katti, A. M., “Fundamentals of
preparative and nonlinear chromatography.” Academic Press, Boston,
1994.
29. Migliorini, C., Fillinger, M., Mazzotti, M., and Morbidelli, M., Chem.
Eng. Sci. 54, 2475 (1999).
30. Mazzotti, M., Neri, B., Gelosa, D., Kruglov, A., and Morbidelli, M.,
Ind. Eng. Chem. Res. 36, 3 (1997).
31. Van Deemter, J. J., Zuiderweg, F. J., and Klinkenberg, A., Chem. Eng.
Sci. 5, 271 (1956).
32. Schiesser, W. E., “The Numerical Method of Lines.” Academic Press,
New York, 1991.
33. Holland, J. H., “Adaptation in natural and artificial systems.” Univ.
Michigan Press, Ann Arbor, MI, 1975.
34. Goldberg, D. E., “Genetic algorithms in search, optimization and ma-
chine learning.” Addison-Wesley, Reading, MA, 1989.
35. Bhaskar, V., Gupta, S. K., and Ray, A. K., Rev. Chem. Eng. 16(1), 1
(2000).
18. Collignon, F., Loenders, R., Martens, J. A., Jacobs, P. A., and Poncelet,
G., J. Catal. 182, 302 (1999).
36. Fogler, H. C., “Elements of Chemical Reaction Engineering.”
Prentice-Hall, Englewood Cliffs, NJ, 1986.
19. Horvath, T., Seiler, M., and Hunger, M., Appl. Catal. A-Gen. 193, 227
(2000).
20. Ali, M. A., Brisdon, B. J., and Thomas, W. J., Appl. Catal. A-Gen. 197,
303 (2000).
37. Boudart, M., AIChE J. 18(3), 465 (1972).
38. Singh, U. K., and Vannice, M. A., AIChE J. 45(5), 1059 (1999).
39. Izqulerdo, J. F., Cunill, F., Vlla, M., Tejero, J., and Iborra, M., J. Chem.
Eng. Data 37, 339 (1992).
21. Sugiyama, K., Kato, K., Miura, H., and Matsuda, T., J. Jpn. Petrol. Inst. 40. Izqulerdo, J. F., Cunill, F., Vlla, M., Iborra, M., and Tejero, J., Ind. Eng.
26, 243 (1983).
Chem. Res. 33, 2830 (1994).
22. Knifton, J. F., and Edwards, J. C., Appl. Catal. A-Gen. 183, 1 (1999).
23. Matouq, M. H., and Shigeo, G., Int. J. Chem. Kinet. 25, 825 (1993).
24. Ray, A. K., Carr, R. W., and Aris, R., Chem. Eng. Sci. 49(4), 469
(1994).
41. Mears, D. F., Ind. Eng. Chem. Process Des. Dev. 10, 541 (1971).
42. Dwidevi, P. N., and Upadhyay, S. N., Ind. Eng. Chem. Process Des.
Dev. 16, 157 (1977).
43. Rehfinger, A., and Hoffmann, U., Chem. Eng. Sci. 45, 1619 (1990).